Source author record

R. Rebolo

R. Rebolo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

215works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

215 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

We present a detailed analysis of the old, extended stellar component of the Local Group dwarf galaxy IC 10 using deep resolved-star photometry in the VIS and NISP bands of the Euclid Early Release Observations. Leveraging Euclid's unique combination of wide field of view and high spatial resolution, we trace red giant branch (RGB) stars out to $\sim$8 kpc from the galaxy centre, reaching azimuthally-averaged surface brightness levels as faint as $μ_{HE}\sim$29 mag arcsec$^{-2}$. Our analysis reveals that IC 10's stellar distribution is significantly more extended than previously thought. After correcting for foreground extinction and subtracting contamination from Milky Way stars and background galaxies, we derive a radial stellar density profile from RGB star counts. The profile shows a marked flattening beyond $\sim$5 kpc, and is best fit by a two-component (Sersic + exponential) model, yielding a total stellar mass in old (age $>$1 Gyr) stars of $M_{\star}=(6.7-8.1)\times10^8 M_{\odot}$. The origin of the outer stellar component is unclear. It may be accreted, even possibly associated with the counter-rotating HI gas in the outer regions of IC 10, or it may represent an ancient in-situ stellar halo. We tentatively detect two symmetric stellar overdensities at the edge of our imagery. These roughly align with the direction of IC 10's orbit around M31, suggesting that they may be signatures of tidal stripping. As part of our analysis, we derive a new distance to IC 10 based on the RGB tip, finding $D=762\pm 20$ kpc and the distance modulus is $(m-M)_0=24.41\pm 0.05$.

preprint2025arXiv

Unveiling the nature of HE 0107-5240

The vast majority of the most iron-poor stars in the Galaxy exhibit a strong carbon enhancement, with C/H ratios only about two orders of magnitude below solar. This unusual chemical composition likely reflects the properties of the gas cloud from which these stars formed, having been enriched by one, or at most a few, supernovae. A remarkable member of this stellar class, HE 0107-5240 with [Fe/H]=-5.56, has been identified as part of a binary system. To constrain its orbital parameters, radial velocity monitoring has been carried out using the ESPRESSO spectrograph. Radial velocities were derived using cross-correlation with a template, taking advantage of the strong G-band feature. Combining all observations yielded a high signal-to-noise spectrum, which has been used to refine our understanding of the stellar chemical composition. Additionally, a co-added UVES spectrum in the blue was used to complement the wavelength coverage of ESPRESSO. Observations of HE 0107-5240 over a span of more than four years have yielded a revised orbital period of about 29 years. Updated elemental abundances have been determined for Sc, Cr, Co, and, tentatively, Al, along with a new upper limit for Be. The iron abundance has been derived from ionised Fe lines. Significant upper limits have been established for Li, Si, and Sr. The star is confirmed to be a long-period binary. Iron abundances derived from neutral and ionised lines are consistent with local thermodynamical equilibrium (LTE) assumption, casting doubt on published deviation from LTE corrections for Fe for this star. The heavy elements Sr and Ba remain undetected, confirming the classification of HE 0107-5240 as a carbon enhanced metal-poor and non enhanced in heavy elements (CEMP-no) star and supporting the absence of an n-capture element plateau at the lowest metallicities.

preprint2023arXiv

QUIJOTE scientific results -- IV. A northern sky survey in intensity and polarization at 10-20GHz with the Multi-Frequency Instrument

We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17 and 19GHz, and covering approximately 29000 deg$^2$, including most of the Northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE instrument (MFI), and have angular resolutions of around $1^\circ$, and sensitivities in polarization within the range 35-40 $μ$K per 1-degree beam, being a factor $\sim 2$-$4$ worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterise the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5%. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A and 3C274) which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11GHz and WMAP 23GHz map peaks at $β=-3.09$ with a standard deviation of 0.14. The measured BB/EE ratio at scales of $\ell=80$ is $0.26\pm 0.07$ for a Galactic cut $|b|>10^\circ$. We find a positive TE correlation for 11GHz at large angular scales ($\ell \lesssim 50$), while the EB and TB signals are consistent with zero in the multipole range $30 \lesssim \ell \lesssim 150$. The maps discussed in this paper are publicly available.

preprint2023arXiv

QUIJOTE scientific results -- IX. Radio sources in the QUIJOTE-MFI wide survey maps

We present the catalogue of Q-U-I JOint TEnerife (QUIJOTE) Wide Survey radio sources extracted from the maps of the Multi-Frequency Instrument compiled between 2012 and 2018. The catalogue contains 786 sources observed in intensity and polarization, and is divided into two separate sub-catalogues: one containing 47 bright sources previously studied by the \emph{Planck} collaboration and an extended catalogue of 739 sources either selected from the \emph{Planck} Second Catalogue of Compact Sources or found through a blind search carried out with a Mexican Hat 2 wavelet. A significant fraction of the sources in our catalogue (38.7 per cent) are within the $|b| \leq 20^\circ$ region of the Galactic plane. We determine statistical properties for those sources that are likely to be extragalactic. We find that these statistical properties are compatible with currently available models, with a $\sim$1.8 Jy completeness limit at 11 GHz. We provide the polarimetric properties of (38, 33, 31, 23) sources with P detected above the $99.99\%$ significance level at (11, 13, 17, 19) GHz, respectively. Median polarization fractions are in the $2.8$-$4.7$\% range in the 11-19 GHz frequency interval. We do not distinguish between Galactic and extragalactic sources here. The results presented here are consistent with those reported in the literature for flat- and steep-spectrum radio sources.

preprint2023arXiv

QUIJOTE scientific results -- V. The microwave intensity and polarisation spectra of the Galactic regions W49, W51 and IC443

We present new intensity and polarisation maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim$ 1 deg angular resolution, with a sensitivity in the range 35-79 $μ$K/beam for total intensity and 13-23 $μ$K/beam for polarisation. For each region, we combine QUIJOTE maps with ancillary data at frequencies ranging from 0.4 to 3000 GHz, reconstruct the spectral energy distribution and model it with a combination of known foregrounds. We detect anomalous microwave emission (AME) in total intensity towards W49 at 4.7$σ$ and W51 at 4.0$σ$ with peak frequencies $ν_{AME}$ = (20.0 $\pm$ 1.4) GHz and $ν_{AME}$ = (17.7 $\pm$ 3.6) GHz respectively; this is the first detection of AME towards W51. The contamination from ultra-compact HII regions to the residual AME flux density is estimated at 10% in W49 and 5% in W51, and does not rule out the AME detection. The polarised SEDs reveal a synchrotron contribution with spectral indices $α_s$ = -0.67 $\pm$ 0.10 in W49 and $α_s$ = -0.51 $\pm$ 0.07 in W51, ascribed to the diffuse Galactic emission and to the local supernova remnant respectively. Towards IC443 in total intensity we measure a broken power-law synchrotron spectrum with cut-off frequency $ν_{0,s}$ = (114 $\pm$ 73) GHz, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarised AME emission is detected in this study.

preprint2023arXiv

QUIJOTE scientific results -- VI. The Haze as seen by QUIJOTE

The Haze is an excess of microwave intensity emission surrounding the Galactic centre. It is spatially correlated with the $γ$-ray Fermi bubbles, and with the S-PASS radio polarization plumes, suggesting a possible common provenance. The models proposed to explain the origin of the Haze, including energetic events at the Galactic centre and dark matter decay in the Galactic halo, do not yet provide a clear physical interpretation. In this paper we present a re-analysis of the Haze including new observations from the Multi-Frequency Instrument (MFI) of the Q-U-I JOint TEnerife (QUIJOTE) experiment, at 11 and 13 GHz. We analyze the Haze in intensity and polarization, characterizing its spectrum. We detect an excess of diffuse intensity signal ascribed to the Haze. The spectrum at frequencies 11$\,\leqν\leq\,$70 GHz is a power-law with spectral index $β^{\rm H}=-2.79\pm0.08$, which is flatter than the Galactic synchrotron in the same region ($β^{\rm S}=-2.98\pm0.04$), but steeper than that obtained from previous works ($β^{\rm H}\sim-2.5$ at 23$\,\leq\,ν\leq\,$70 GHz). We also observe an excess of polarized signal in the QUIJOTE-MFI maps in the Haze area. This is a first hint detection of polarized Haze, or a consequence of curvature of the synchrotron spectrum in that area. Finally, we show that the spectrum of polarized structures associated with Galactic centre activity is steep at low frequencies ($β\sim -3.2$ at 2.3 $\leqν\leq$ 23 GHz), and becomes flatter above 11 GHz.

preprint2023arXiv

QUIJOTE Scientific Results -- VII. Galactic AME sources in the QUIJOTE-MFI Northern Hemisphere Wide-Survey

The QUIJOTE-MFI Northern Hemisphere Wide-Survey has provided maps of the sky above declinations $-30^\circ$ at 11, 13, 17 and 19$\,$GHz. These data are combined with ancillary data to produce Spectral Energy Distributions in intensity in the frequency range 0.4--3\,000$\,$GHz on a sample of 52 candidate compact sources harbouring anomalous microwave emission (AME). We apply a component separation analysis at 1$^\circ$ scale on the full sample from which we identify 44 sources with high AME significance. We explore correlations between different fitted parameters on this last sample. QUIJOTE-MFI data contribute to notably improve the characterisation of the AME spectrum, and its separation from the other components. In particular, ignoring the 10--20\,GHz data produces on average an underestimation of the AME amplitude, and an overestimation of the free-free component. We find an average AME peak frequency of 23.6 $\pm$ 3.6$\,$GHz, about 4$\,$GHz lower than the value reported in previous studies. The strongest correlation is found between the peak flux density of the thermal dust and of the AME component. A mild correlation is found between the AME emissivity ($A_{\rm AME}/τ_{250}$) and the interstellar radiation field. On the other hand no correlation is found between the AME emissivity and the free-free radiation Emission Measure. Our statistical results suggest that the interstellar radiation field could still be the main driver of the intensity of the AME as regards spinning dust excitation mechanisms. On the other hand, it is not clear whether spinning dust would be most likely associated with cold phases of the interstellar medium rather than with hot phases dominated by free-free radiation.

preprint2023arXiv

QUIJOTE scientific results -- VIII. Diffuse polarized foregrounds from component separation with QUIJOTE-MFI

We derive linearly polarized astrophysical component maps in the Northern Sky from the QUIJOTE-MFI data at 11 and 13 GHz in combination with the WMAP K and Ka bands (23 and 33 GHz) and all Planck polarized channels (30-353 GHz), using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the parameter estimation of the low-frequency foregrounds, especially the estimation of the synchrotron spectral index, $β_s$. We present the first detailed $β_s$ map of the Northern Celestial Hemisphere at a smoothing scale of $2^{\circ}$. We find statistically significant spatial variability across the sky. We obtain an average value of $-3.08$ and a dispersion of $0.13$, considering only pixels with reliable goodness-of-fit. The power law model of the synchrotron emission provides a good fit to the data outside the Galactic plane but fails to track the complexity within this region. Moreover, when we assume a synchrotron model with uniform curvature, $c_s$, we find a value of $c_s = -0.0797 \pm 0.0012$. However, there is insufficient statistical significance to determine which model is favoured, either the power law or the power law with uniform curvature. Furthermore, we estimate the thermal dust spectral parameters in polarization. Our CMB, synchrotron, and thermal dust maps are highly correlated with the corresponding products of the PR4 Planck release, although some large-scale differences are observed in the synchrotron emission. Finally, we find that the $β_s$ estimation in the high signal-to-noise synchrotron emission areas is prior-independent while, outside these regions, the prior governs the $β_s$ estimation.

preprint2023arXiv

The pristine nature of SMSS 1605$-$1443 revealed by ESPRESSO

SMSS J160540.18$-$144323.1 is the carbon-enhanced metal-poor (CEMP) star with the lowest iron abundance ever measured, [Fe/H]=-6.2, which was first reported with the SkyMapper telescope. The carbon abundance is A(C)~6.1 in the low-C band, as the majority of the stars in this metallicity range. Yet, constraining the isotopic ratio of key species, such as carbon, sheds light on the properties and origin of these elusive stars. We performed high-resolution observations of SMSS1605$-$1443 with the ESPRESSO spectrograph to look for variations in the radial velocity ($v_{rad}$) with time. These data have been combined with older MIKE and UVES archival observations to enlarge the temporal baseline. The $^{12}$C/$^{13}$C isotopic ratio is also studied to explore the possibility of mass transfer from a binary companion. A cross-correlation function against a natural template was applied to detect $v_{rad}$ variability and a spectral synthesis technique was used to derive $^{12}$C/$^{13}$C in the stellar atmosphere. We confirm previous indications of binarity in SMSS1605$-$1443 and measured a lower limit $^{12}$C/$^{13}$C$>60$ at more than a 3$σ$ confidence level, proving that this system is chemically unmixed and that no mass transfer from the unseen companion has happened so far. Thus, we confirm the CEMP-no nature of SMSS1605$-$1443 and show that the pristine chemical composition of the cloud from which it formed is currently imprinted in its stellar atmosphere free of contamination.

preprint2022arXiv

A candidate short-period sub-Earth orbiting Proxima Centauri

Proxima Centauri is the closest star to the Sun. This small, low-mass, mid M dwarf is known to host an Earth-mass exoplanet with an orbital period of 11.2 days within the habitable zone, as well as a long-period planet candidate with an orbital period of close to 5 years. We report on the analysis of a large set of observations taken with the ESPRESSO spectrograph at the VLT aimed at a thorough evaluation of the presence of a third low-mass planetary companion, which started emerging during a previous campaign. Radial velocities (RVs) were calculated using both a cross-correlation function (CCF) and a template matching approach. The RV analysis includes a component to model Proxima's activity using a Gaussian process (GP). We use the CCF's full width at half maximum to help constrain the GP, and we study other simultaneous observables as activity indicators in order to assess the nature of any potential RV signals. We detect a signal at 5.12 $\pm$ 0.04 days with a semi-amplitude of 39 $\pm$ 7 cm/s. The analysis of subsets of the ESPRESSO data, the activity indicators, and chromatic RVs suggest that this signal is not caused by stellar variability but instead by a planetary companion with a minimum mass of 0.26 $\pm$ 0.05 $M_\oplus$ (about twice the mass of Mars) orbiting at 0.029 au from the star. The orbital eccentricity is well constrained and compatible with a circular orbit.

preprint2022arXiv

A novel framework for semi-Bayesian radial velocities through template matching

The detection and characterization of an increasing variety of exoplanets has been in part possible thanks to the continuous development of high-resolution, stable spectrographs, and using the Doppler radial-velocity (RV) method. The Cross Correlation Function (CCF) method is one of the traditional approaches for RV extraction. More recently, template matching was introduced as an advantageous alternative for M-dwarf stars. In this paper, we describe a new implementation of template matching within a semi-Bayesian framework, providing a more statistically principled characterization of the RV measurements. In this context, a common RV shift is used to describe the difference between each spectral order of a given stellar spectrum and a template built from the available observations. Posterior probability distributions are obtained for the relative RV associated with each spectrum, after marginalizing with respect to the continuum. This methodology was named S-BART: Semi-Bayesian Approach for RVs with Template-matching, and it can be applied to HARPS and ESPRESSO. The application of our method to HARPS archival observations of Barnard's star allowed us to validate our implementation against HARPS-TERRA and SERVAL. Then, we applied it to 33 ESPRESSO targets, evaluating its performance and comparing it with the CCF method. We found a decrease in the median RV scatter of \sim 10\% and \sim 4\% for M- and K-type stars, respectively. S-BART yields more precise RV estimates than the CCF method, particularly in the case of M-type stars where a median uncertainty of \sim 15 cm/s is achieved over 309 observations. Further, we estimated the nightly zero point (NZP) of ESPRESSO, finding a weighted NZP scatter below \sim 0.7 m/s. As this includes stellar variability, photon noise, and potential planetary signals, it should be taken as an upper limit of the RV precision attainable with ESPRESSO data.

preprint2022arXiv

CaRM: Exploring the chromatic Rossiter-McLaughlin effect. The cases of HD 189733b and WASP-127b

In this paper we introduce CaRM, a semi-automatic code for the retrieval of broadband transmission spectra of transiting planets through the chromatic Rossiter-McLaughlin method. We applied it to HARPS and ESPRESSO observations of two exoplanets to retrieve the transmission spectrum and we analyze its fitting transmission models. We used the strong radius dependence of the Rossiter-McLaughlin (RM) effect amplitude, caused by planetary companions, to measure the apparent radius change caused by the exoplanet atmosphere. In order to retrieve the transmission spectrum, the radial velocities, which were computed over wavelength bins that encompass several spectral orders, were used to simultaneously fit the Keplerian motion and the RM effect. From this, the radius ratio was computed as a function of the wavelength, which allows one to retrieve the low-resolution broadband transmission spectrum of a given exoplanet. CaRM offers the possibility to use two Rossiter-McLaughlin models taken from ARoME and PyAstronomy, associated with a Keplerian function to fit radial velocities during transit observations automatically. Furthermore it offers the possibility to use some methods that could, in theory, mitigate the effect of perturbation in the radial velocities during transits. The CaRM code allows one to retrieve the transmission spectrum of a given exoplanet using minimal user interaction. We demonstrate that it allows one to compute the low-resolution broadband transmission spectra of exoplanets observed using high-resolution spectrographs such as HARPS and ESPRESSO.

preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

preprint2022arXiv

Euclid: Constraining ensemble photometric redshift distributions with stacked spectroscopy

The ESA Euclid mission will produce photometric galaxy samples over 15000 square degrees of the sky that will be rich for clustering and weak lensing statistics. The accuracy of the cosmological constraints derived from these measurements will depend on the knowledge of the underlying redshift distributions based on photometric redshift calibrations. A new approach is proposed to use the stacked spectra from Euclid slitless spectroscopy to augment broad-band photometric information to constrain the redshift distribution with spectral energy distribution fitting. The high spectral resolution available in the stacked spectra complements the photometry and helps to break the colour-redshift degeneracy and constrain the redshift distribution of galaxy samples. We modelled the stacked spectra as a linear mixture of spectral templates. The mixture may be inverted to infer the underlying redshift distribution using constrained regression algorithms. We demonstrate the method on simulated Vera C. Rubin Observatory and Euclid mock survey data sets based on the Euclid Flagship mock galaxy catalogue. We assess the accuracy of the reconstruction by considering the inference of the baryon acoustic scale from angular two-point correlation function measurements. We selected mock photometric galaxy samples at redshift z>1 using the self-organising map algorithm. Considering the idealised case without dust attenuation, we find that the redshift distributions of these samples can be recovered with 0.5% accuracy on the baryon acoustic scale. The estimates are not significantly degraded by the spectroscopic measurement noise due to the large sample size. However, the error degrades to 2% when the dust attenuation model is left free. We find that the colour degeneracies introduced by attenuation limit the accuracy considering the wavelength coverage of Euclid near-infrared spectroscopy.

preprint2022arXiv

HADES RV Programme with HARPS-N at TNG XV. Planetary occurrence rates around early-M dwarfs

We present the complete Bayesian statistical analysis of the HArps-n red Dwarf Exoplanet Survey (HADES), which monitored the radial velocities of a large sample of M dwarfs with HARPS-N at TNG, over the last 6 years. The targets were selected in a narrow range of spectral types from M0 to M3, $0.3$ M$_\odot < M_\star < 0.71$ M$_\odot$, in order to study the planetary population around a well-defined class of host stars. We take advantage of Bayesian statistics to derive an accurate estimate of the detectability function of the survey. Our analysis also includes the application of Gaussian Process approach to take into account stellar activity induced radial velocity variations, and improve the detection limits, around the most-observed and most-active targets. The Markov chain Monte Carlo and Gaussian process technique we apply in this analysis has proven very effective in the study of M-dwarf planetary systems, helping the detection of most of the HADES planets. From the detectability function we can calculate the occurrence rate of small mass planets around early-M dwarfs, either taking into account only the 11 already published HADES planets or adding also the 5 new planetary candidates discovered in this analysis, and compare them with the previous estimates of planet occurrence around M-dwarf or Solar-type stars: considering only the confirmed planets, we find the highest frequency for low-mass planets ($1$ M$_\oplus < m_p \sin i < 10$ M$_\oplus$) with periods $10$ d$ < P < 100$ d, $f_\text{occ} = 85^{+5}_{-19}\%$, while for short-period planets ($1$ d$ < P < 10$ d) we find a frequency of $f_\text{occ} = 10.3^{+8.4}_{-3.3}\%$, significantly lower than for later-M dwarfs. These results, and their comparison with other surveys focused on different stellar types, confirms the central role that stellar mass plays in the formation and evolution of planetary systems.

preprint2022arXiv

HORuS transmission spectroscopy and revised planetary parameters of KELT-7b

We report on the high-resolution spectroscopic observations of two planetary transits of the hot Jupiter KELT-7b ($M_{\rm p}$ $=$ 1.28 $\pm$ 0.17 M$_{\rm Jup}$, $T_{\rm eq}$ $=$ 2028 K) observed with the High Optical Resolution Spectrograph (HORuS) mounted on the 10.4-m Gran Telescopio Canarias (GTC). A new set of stellar parameters are obtained for the rapidly rotating parent star from the analysis of the spectra. Using the newly derived stellar mass and radius, and the planetary transit data of the Transiting Exoplanet Survey Satellite (TESS) together with the HORuS velocities and the photometric and spectroscopic data available in the literature, we update and improve the ephemeris of KELT-7b. Our results indicate that KELT-7 has an angle $λ$ = $-$10.55 $\pm$ 0.27 deg between the sky projections of the star's spin axis and the planet's orbital axis. By combining this angle and our newly derived stellar rotation period of 1.38 $\pm$ 0.05 d, we obtained a 3D obliquity $ψ$ = 12.4 $\pm$ 11.7 deg (or 167.6 deg), thus reinforcing that KELT-7 is a well-aligned planetary system. We search for the presence of H$α$, Li I, Na I, Mg I, and Ca II features in the transmission spectrum of KELT-7b but we are only able to determine upper limits of 0.08-1.4 % on their presence after accounting for the contribution of the stellar variability to the extracted planetary spectrum. We also discuss the impact of stellar variability in the planetary data. Our results reinforce the importance of monitoring the parent star when performing high-resolution transmission spectroscopy of the planetary atmosphere in the presence of stellar activity.

preprint2022arXiv

Radio emission in a nearby ultracool dwarf binary: a multi-frequency study

The substellar triple system VHS J125601.92$-$125723.9 is composed by an equal-mass M7.5 brown dwarf binary and a L7 low-mass substellar object. In this work, we aim to identify the origin of the radio emission occurring in the central binary of VHS 1256$-$1257 while discussing the expected mechanisms involved in the radio emission of ultracool dwarfs (UCDs). We observed this system with the Karl G. Jansky Very Large Array, the European very-long-baseline interferometry (VLBI) Network, the enhanced Multi Element Remotely Linked Interferometer Network, the NOrthern Extended Millimeter Array, and the Atacama Large Millimetre Array at frequencies ranging from 5 GHz up to 345 GHz in several epochs during 2017, 2018, and 2019. We have found radio emission at 6 GHz and 33 GHz coincident with the expected position of the central binary of VHS~1256$-$1257. The Stokes I density flux detected were 73 $\pm$ 4 $μ$Jy and 83 $\pm$ 13 $μ$Jy, respectively, with no detectable circular polarisation or pulses. No emission is detected at higher frequencies (230 GHz and 345 GHz) nor at 5 GHz with VLBI arrays. The emission appears to be stable over almost 3 years at 6 GHz. To explain the constraints obtained both from the detections and non-detections we considered multiple scenarios including thermal and non-thermal emission, and different contributions from each component of the binary. Our results can be well explained by non-thermal gyrosynchrotron emission originating at radiation belts with a low plasma density (n$_e$ = 300$-$700 cm$^{-3}$), a moderate magnetic field strength (B $\approx$ 140 G), and an energy distribution of electrons following a power-law ($d N / d E \propto E^{- δ}$) with $δ$ fixed at 1.36. These radiation belts would need to be present in both components and also be viewed equatorially.

preprint2022arXiv

The early afterglow of GRB 190829A

GRB 190829A at z=0.0785 is the fourth closest long GRB ever detected by the Neil Gehrels Swift observatory, and the third confirmed case with a very high energy component. We present our multi-wavelength analysis of this rare event, focusing on its early stages of evolution, and including data from Swift, the MASTER global network of optical telescopes, ALMA, and ATCA. We report sensitive limits on the linear polarization of the optical emission, disfavouring models of off-axis jets to explain the delayed afterglow peak. The study of the multi-wavelength light curves and broadband spectra supports a model with at least two emission components: a bright reverse shock emission, visible at early times in the optical and X-rays and, later, in the radio band; and a forward shock component dominating at later times and lower radio frequencies. A combined study of the prompt and afterglow properties shows many similarities with cosmological long GRBs, suggesting that GRB 190829A is an example of classical GRBs in the nearby universe.

preprint2022arXiv

Transmission spectroscopy of MASCARA-1b with ESPRESSO: Challenges of overlapping orbital and Doppler tracks

Atmospheric studies at high spectral resolution have shown the presence of molecules, neutral and ionised metals, and hydrogen in the transmission spectrum of ultra-hot Jupiters, and have started to probe the dynamics of their atmospheres. We analyse the transmission spectrum of MASCARA-1b, one of the densest ultra-hot Jupiters orbiting a bright (V=8.3) star. We focus on the CaII H&K, NaI, LiI, H$α$, and KI D1 spectral lines and on the cross-correlated FeI, FeII, CaI, YI, VI, VII, CaH, and TiO lines. For those species that are not present in the stellar spectrum, no detections are reported, but we measure upper limits with an excellent precision ($\sim10$ ppm for particular species). For those species that are present in the stellar spectrum and whose planet-occulted spectral lines induce spurious features in the planetary transmission spectrum, an accurate modelling of the Rossiter-McLaughlin effect (RM) and centre-to-limb variations (CLV) is necessary to recover possible atmospheric signals. In the case of MASCARA-1b, this is difficult due to the overlap between the radial velocities of the stellar surface regions occulted by MASCARA-1b and the orbital track along which the planet atmospheric signal is expected to be found. To try to disentangle a possible planetary signal, we compare our results with models of the RM and CLV effects, and estimate the uncertainties of our models depending on the different system parameters. Unfortunately, more precise measurements of the spin-orbit angle are necessary to better constrain the planet-occulted track and correct for the transit effects in the transmission spectrum with enough precision to be able to detect or discard possible planetary absorptions. Finally, we discuss the possibility that non-detections are related to the low absorption expected for a high surface gravity planet such as MASCARA-1b.

preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

preprint2021arXiv

Detection of Spectral Variations of Anomalous Microwave Emission with QUIJOTE and C-BASS

Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range $10$-$60\,$GHz and a new window into the properties of sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the $\approx10^{\circ}$ diameter $λ$ Orionis ring by combining intensity data from the QUIJOTE experiment at $11$, $13$, $17$ and $19\,$GHz and the C-Band All Sky Survey (C-BASS) at $4.76\,$GHz, together with 19 ancillary datasets between $1.42$ and $3000\,$GHz. Maps of physical parameters at $1^{\circ}$ resolution are produced through Markov Chain Monte Carlo (MCMC) fits of spectral energy distributions (SEDs), approximating the AME component with a log-normal distribution. AME is detected in excess of $20\,σ$ at degree-scales around the entirety of the ring along photodissociation regions (PDRs), with three primary bright regions containing dark clouds. A radial decrease is observed in the AME peak frequency from $\approx35\,$GHz near the free-free region to $\approx21\,$GHz in the outer regions of the ring, which is the first detection of AME spectral variations across a single region. A strong correlation between AME peak frequency, emission measure and dust temperature is an indication for the dependence of the AME peak frequency on the local radiation field. The AME amplitude normalised by the optical depth is also strongly correlated with the radiation field, giving an overall picture consistent with spinning dust where the local radiation field plays a key role.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond.

preprint2021arXiv

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

preprint2021arXiv

Euclid: Forecasts from redshift-space distortions and the Alcock-Paczynski test with cosmic voids

Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection of cosmic voids therein and the unique benefit they provide for cosmology. In particular, we study the imprints of dynamic and geometric distortions of average void shapes and their constraining power on the growth of structure and cosmological distance ratios. To this end, we made use of the Flagship mock catalog, a state-of-the-art simulation of the data expected to be observed with Euclid. We arranged the data into four adjacent redshift bins, each of which contains about 11000 voids and estimated the stacked void-galaxy cross-correlation function in every bin. Fitting a linear-theory model to the data, we obtained constraints on $f/b$ and $D_M H$, where $f$ is the linear growth rate of density fluctuations, $b$ the galaxy bias, $D_M$ the comoving angular diameter distance, and $H$ the Hubble rate. In addition, we marginalized over two nuisance parameters included in our model to account for unknown systematic effects. With this approach, Euclid will be able to reach a relative precision of about 4% on measurements of $f/b$ and 0.5% on $D_M H$ in each redshift bin. Better modeling or calibration of the nuisance parameters may further increase this precision to 1% and 0.4%, respectively. Our results show that the exploitation of cosmic voids in Euclid will provide competitive constraints on cosmology even as a stand-alone probe. For example, the equation-of-state parameter $w$ for dark energy will be measured with a precision of about 10%, consistent with previous more approximate forecasts.

preprint2021arXiv

Observations with KIDs Interferometer Spectrum Survey(KISS)

We describe the preliminary on-sky results of the KIDs Interferometer Spectrum Survey (KISS), a spectral imager with a 1 deg field of view (FoV). The instrument operates in the range 120-180 GHz from the 2.25 m Q-U-I JOint TEnerife telescope in Teide Observatory (Tenerife, Canary Islands), at 2 395 m altitude above sea level. Spectra at low resolution, up to 1.45 GHz, are obtained using a fast (3.72 Hz mechanical frequency) Fourier transform spectrometer, coupled to a continuous dilution cryostat with a stabilized temperature of 170 mK that hosts two 316-pixel arrays of lumped-element kinetic inductance detectors. KISS generates more than 3 000 spectra per second during observations and represents a pathfinder to demonstrate the potential for spectral mapping with large FoV. We give an overall description of the spectral mapping paradigm and we present recent results from observations, in this paper.

preprint2021arXiv

Retrieving the transmission spectrum of HD 209458b using CHOCOLATE: A new chromatic Doppler tomography technique

Multiband photometric transit observations or low-resolution spectroscopy (spectro-photometry) are normally used to retrieve the broadband transmission spectra of transiting exoplanets in order to assess the chemical composition of their atmospheres. In this paper, we present an alternative approach for recovering the broadband transmission spectra using chromatic Doppler tomography based on physical modeling through the SOAP tool: CHOCOLATE (CHrOmatiC line prOfiLe tomogrAphy TEchnique). To validate the method and examine its performance, we use observational data recently obtained with the ESPRESSO instrument to retrieve the transmission spectra of the archetypal hot Jupiter HD 209458b. Our findings indicate that the recovered transmission spectrum is in good agreement with the results presented in previous studies, which used different methodologies to extract the spectrum, achieving similar precision. We explored several atmospheric models and inferred from spectral retrieval that a model containing H2O and NH3 is the preferred scenario. The CHOCOLATE methodology is particularly interesting for future studies of exoplanets around young and active stars or moderate to fast rotating stars, considering SOAP's ability to model stellar active regions and the fact that the rotational broadening of spectral lines favors its application. Furthermore, CHOCOLATE will allow the broad transmission spectrum of a planet to be retrieved using high signal-to-noise ratio, high-resolution spectroscopy with the next generation of Extremely Large Telescopes (ELTs), where low-resolution spectroscopy will not always be accessible.

preprint2021arXiv

Six transiting planets and a chain of Laplace resonances in TOI-178

Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth's density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes.

preprint2021arXiv

Study of the thermal and nonthermal emission components in M31: the Sardinia Radio Telescope view at 6.6 GHz

The Andromeda galaxy is the best-known large galaxy besides our own Milky Way. Several images and studies exist at all wavelengths from radio to hard X-ray. Nevertheless, only a few observations are available in the microwave range where its average radio emission reaches the minimum. In this paper, we want to study the radio morphology of the galaxy, decouple thermal from nonthermal emission, and extract the star formation rate. We also aim to derive a complete catalog of radio sources for the mapped patch of sky. We observed the Andromeda galaxy with the Sardinia Radio Telescope at 6.6 GHz with very high sensitivity and angular resolution, and an unprecedented sky coverage. Using new 6.6 GHz data and Effelsberg radio telescope ancillary data, we confirm that, globally, the spectral index is $\sim 0.7-0.8$, while in the star forming regions it decreases to $\sim 0.5$. By disentangling (gas) thermal and nonthermal emission, we find that at 6.6 GHz, thermal emission follows the distribution of HII regions around the ring. Nonthermal emission within the ring appears smoother and more uniform than thermal emission because of diffusion of the cosmic ray electrons away from their birthplaces. This causes the magnetic fields to appear almost constant in intensity. Furthermore, we calculated a map of the star formation rate based on the map of thermal emission. Integrating within a radius of $R_{max}=15$ kpc, we obtained a total star formation rate of $0.19 \pm 0.01$ $M_{\odot}$/yr in agreement with previous results in the literature. Finally, we correlated our radio data with infrared images of the Andromeda galaxy. We find an unexpectedly high correlation between nonthermal and mid-infrared data in the central region, with a correlation parameter $r=0.93$.

preprint2021arXiv

The atmosphere of HD 209458b seen with ESPRESSO. No detectable planetary absorptions at high resolution

We observed two transits of the iconic gas giant HD 209458b between 380 and 780 nm, using the high-resolution ESPRESSO spectrograph. The derived planetary transmission spectrum exhibits features at all wavelengths where the parent star shows strong absorption lines, for example, NaI, MgI, FeI, FeII, CaI, VI, H$α$, and KI. We interpreted these features as the signature of the deformation of the stellar line profiles due to the Rossiter-McLaughlin effect, combined with the centre-to-limb effects on the stellar surface, which is in agreement with similar reports recently presented in the literature. We also searched for species that might be present in the planetary atmosphere but not in the stellar spectra, such as TiO and VO, and obtained a negative result. Thus, we find no evidence of any planetary absorption, including previously reported NaI, in the atmosphere of HD 209458b. The high signal-to-noise ratio in the transmission spectrum allows us to compare the modelled deformation of the stellar lines in assuming different one-dimensional stellar atmospheric models. We conclude that the differences among various models and observations remain within the precision limits of the data. However, the transmission light curves are better explained when the centre-to-limb variation is not included in the computation and only the Rossiter-McLaughlin deformation is considered. This demonstrates that ESPRESSO is currently the best facility for spatially resolving the stellar surface spectrum in the optical range using transit observations and carrying out empirical validations of stellar models.

preprint2020arXiv

Atmospheric Rossiter-McLaughlin effect and transmission spectroscopy of WASP-121b with ESPRESSO

WASP-121b is one of the most studied Ultra-hot Jupiters: many recent analyses of its atmosphere report interesting features at different wavelength ranges. In this paper we analyze one transit of WASP-121b acquired with the high-resolution spectrograph ESPRESSO at VLT in 1-telescope mode, and one partial transit taken during the commissioning of the instrument in 4-telescope mode. We investigate the anomalous in-transit radial velocity curve and study the transmission spectrum of the planet. By analysing the in-transit radial velocities we were able to infer the presence of the atmospheric Rossiter-McLaughlin effect. We measured the height of the planetary atmospheric layer that correlates with the stellar mask (mainly Fe) to be 1.052$\pm$0.015 Rp and we also confirmed the blueshift of the planetary atmosphere. By examining the planetary absorption signal on the stellar cross-correlation functions we confirmed the presence of a temporal variation of its blueshift during transit, which could be investigated spectrum-by-spectrum. We detected significant absorption in the transmission spectrum for Na, H, K, Li, CaII, Mg, and we certified their planetary nature by using the 2D tomographic technique. Particularly remarkable is the detection of Li, with a line contrast of $\sim$0.2% detected at the 6$σ$ level. With the cross-correlation technique we confirmed the presence of FeI, FeII, CrI and VI. H$α$ and CaII are present up to very high altitudes in the atmosphere ($\sim$1.44 Rp and $\sim$2 Rp, respectively), and also extend beyond the transit-equivalent Roche lobe radius of the planet. These layers of the atmosphere have a large line broadening that is not compatible with being caused by the tidally-locked rotation of the planet alone, and could arise from vertical winds or high-altitude jets in the evaporating atmosphere.

preprint2020arXiv

Broadband transmission spectroscopy of HD209458b with ESPRESSO: Evidence for Na, TiO, or both

The detection and characterization of exoplanet atmospheres is currently one of the main drivers pushing the development of new observing facilities. In this context, high-resolution spectrographs are {proving} their potential and showing that high-resolution spectroscopy will be paramount in this field. We aim to make use of ESPRESSO high-resolution spectra, which cover two transits of HD209458b, to probe the broadband transmission optical spectrum of the planet. We applied the chromatic Rossiter-McLaughin method to derive the transmission spectrum of HD209458b. We compared the results with previous HST observations and with synthetic spectra. We recover a transmission spectrum of HD209458b similar to the one obtained with HST data. The models suggest that the observed signal can be explained by only Na, only TiO, or both Na and TiO, even though none is fully capable of explaining our observed transmission spectrum. Extra absorbers may be needed to explain the full dataset, though modeling approximations and observational errors can also be responsible for the observed mismatch. Using the chromatic Rossiter-McLaughlin technique, ESPRESSO is able to provide broadband transmission spectra of exoplanets from the ground, in conjunction with space-based facilities, opening good perspectives for similar studies of other planets.

preprint2020arXiv

ESPRESSO high resolution transmission spectroscopy of WASP-76b

Aims. We report on ESPRESSO high-resolution transmission spectroscopic observations of two primary transits of the highly-irradiated, ultra-hot Jupiter-size planet WASP-76b. We investigate the presence of several key atomic and molecular features of interest that may reveal the atmospheric properties of the planet. Methods. We extracted two transmission spectra of WASP-76b with R approx 140,000 using a procedure that allowed us to process the full ESPRESSO wavelength range (3800-7880 A) simultaneously. We observed that at a high signal-to-noise ratio, the continuum of ESPRESSO spectra shows wiggles that are likely caused by an interference pattern outside the spectrograph. To search for the planetary features, we visually analysed the extracted transmission spectra and cross-correlated the observations against theoretical spectra of different atomic and molecular species. Results. The following atomic features are detected: Li I, Na I, Mg I, Ca II, Mn I, K I, and Fe I. All are detected with a confidence level between 9.2 sigma (Na I) and 2.8 sigma (Mg I). We did not detect the following species: Ti I, Cr I, Ni I, TiO, VO, and ZrO. We impose the following 1 sigma upper limits on their detectability: 60, 77, 122, 6, 8, and 8 ppm, respectively. Conclusions. We report the detection of Li I on WASP-76b for the first time. In addition, we found the presence of Na I and Fe I as previously reported in the literature. We show that the procedure employed in this work can detect features down to the level of ~ 0.1 % in the transmission spectrum and ~ 10 ppm by means of a cross-correlation method. We discuss the presence of neutral and singly ionised features in the atmosphere of WASP-76b.

preprint2020arXiv

ESPRESSO@VLT -- On-sky performance and first results

ESPRESSO is the new high-resolution spectrograph of ESO's Very-Large Telescope (VLT). It was designed for ultra-high radial-velocity precision and extreme spectral fidelity with the aim of performing exoplanet research and fundamental astrophysical experiments with unprecedented precision and accuracy. It is able to observe with any of the four Unit Telescopes (UT) of the VLT at a spectral resolving power of 140,000 or 190,000 over the 378.2 to 788.7 nm wavelength range, or with all UTs together, turning the VLT into a 16-m diameter equivalent telescope in terms of collecting area, while still providing a resolving power of 70,000. We provide a general description of the ESPRESSO instrument, report on the actual on-sky performance, and present our Guaranteed-Time Observation (GTO) program with its first results. ESPRESSO was installed on the Paranal Observatory in fall 2017. Commissioning (on-sky testing) was conducted between December 2017 and September 2018. The instrument saw its official start of operations on October 1st, 2018, but improvements to the instrument and re-commissioning runs were conducted until July 2019. The measured overall optical throughput of ESPRESSO at 550 nm and a seeing of 0.65 arcsec exceeds the 10% mark under nominal astro-climatic conditions. We demonstrate a radial-velocity precision of better than 25 cm/s during one night and 50 cm/s over several months. These values being limited by photon noise and stellar jitter show that the performanceis compatible with an instrumental precision of 10 cm/s. No difference has been measured across the UTs neither in throughput nor RV precision. The combination of the large collecting telescope area with the efficiency and the exquisite spectral fidelity of ESPRESSO opens a new parameter space in RV measurements, the study of planetary atmospheres, fundamental constants, stellar characterisation and many other fields.

preprint2020arXiv

HADES RV Programme with HARPS-N at TNG XII. The abundance signature of M dwarf stars with planets

Most of our current knowledge on planet formation is still based on the analysis of main-sequence, solar-type stars. Conversely, detailed chemical studies of large samples of M-dwarf planet hosts are still missing. We develop for the first time a methodology to determine stellar abundances of elements others than iron for M dwarf stars from high-resolution, optical spectra. Our methodology is based on the use of principal component analysis and sparse Bayesian's methods. We made use of a set of M dwarfs orbiting around an FGK primary with known abundances to train our methods. We applied our methods to derive stellar metalliticies and abundances of a large sample of M dwarfs observed within the framework of current radial velocity surveys. We then used a sample of nearby FGK stars to cross-validate our technique by comparing the derived abundance trends in the M dwarf sample with those found on the FGK stars. The metallicity distribution of the different subsamples shows that M dwarfs hosting giant planets show a planet-metallicity correlation as well as a correlation with the stellar mass. M dwarfs hosting low-mass planets do not seem to follow the planet-metallicity correlation. We also found that the frequency of low-mass planets does not depend on the mass of the stellar host. These results seem in agreement with previous works. However, we note that for giant planet hosts our metallicities predict a weaker planet metallicity correlation but a stronger mass-dependency than photometric values. We show, for the first time, that there seems to be no differences in the abundance distribution of elements different from iron between M dwarfs with and without known planets. Our data shows that low-mass stars with planets follow the same metallicity, mass, and abundance trends than their FGK counterparts.

preprint2020arXiv

KISS: a spectrometric imager for millimetre cosmology

Clusters of galaxies are used to map the large-scale structures in the universe and as probe of universe evolution. They can be observed through the Sunyaev-Zel'dovich (SZ) effect. At this respect the spectro-imaging at low resolution frequency is an important tool, today, for the study of cluster of galaxies. We have developed KISS (KIDs-Interferometer-Spectrum-Survey), a spectrometric imager dedicated to the secondary anisotropies of the Cosmic Microwave Background (CMB). The multi-frequency approach permits to improve the component separation with respect to predecessor experiments. In this paper, firstly, we provide a description of the scientific context and the state of the art of SZ observations. Secondly, we describe the KISS instrument. Finally, we show preliminary results of the ongoing commissioning campaign.

preprint2020arXiv

Lowly polarized light from a highly magnetized jet of GRB 190114C

We report multi-color optical imaging and polarimetry observations of the afterglow of the first TeV- detected gamma-ray burst, GRB 190114C, using RINGO3 and MASTER II polarimeters. Observations begin 31 s after the onset of the GRB and continue until $\sim 7000\,$s post-burst. The light curves reveal a chromatic break at $\sim 400- 500\,$s, with initial temporal decay $α= 1.669 \pm 0.013$ flattening to $α\sim 1$ post-break, which we model as a combination of reverse and forward-shock components, with magnetization parameter $R_{\rm B} \sim 70$. The observed polarization degree decreases from $7.7 \pm 1.1\%$ to $2-4\%$ during $52-109\,$s post-burst and remains steady at this level for the subsequent $\sim 2000$-s, at constant position angle. Broadband spectral energy distribution modeling of the afterglow confirms GRB 190114C is highly obscured (A$_{\rm v, HG} = 1.49 \pm 0.12 \,$mag; N$_{\rm H, HG}= (9.0 \pm 0.3) \times 10^{22}\,$cm$^{-2}$). We interpret the measured afterglow polarization as intrinsically low and dominated by dust, in contrast to ${\rm P} >10\%$ measured previously for other GRB reverse shocks, with a small contribution from polarized prompt photons in the first minute. We test whether 1st and higher-order inverse Compton scattering in a magnetized reverse shock can explain the low optical polarization and the sub-TeV emission but conclude neither is explained in the reverse shock Inverse Compton model. Instead, the unexpectedly low intrinsic polarization degree in GRB 190114C can be explained if large-scale jet magnetic fields are distorted on timescales prior to reverse shock emission.

preprint2020arXiv

Measuring and characterizing the line profile of HARPS with a laser frequency comb

Aims. We study the 2D spectral line profile of HARPS (High Accuracy Radial Velocity Planet Searcher), measuring its variation with position across the detector and with changing line intensity. The characterization of the line profile and its variations are important for achieving the precision of the wavelength scales of 10^{-10} or 3.0 cm/s necessary to detect Earth-twins in the habitable zone around solar-like stars. Methods. We used a laser frequency comb (LFC) with unresolved and unblended lines to probe the instrument line profile. We injected the LFC light (attenuated by various neutral density filters) into both the object and the reference fibres of HARPS, and we studied the variations of the line profiles with the line intensities. We applied moment analysis to measure the line positions, widths, and skewness as well as to characterize the line profile distortions induced by the spectrograph and detectors. Based on this, we established a model to correct for point spread function distortions by tracking the beam profiles in both fibres. Results. We demonstrate that the line profile varies with the position on the detector and as a function of line intensities. This is consistent with a charge transfer inefficiency (CTI) effect on the HARPS detector. The estimate of the line position depends critically on the line profile, and therefore a change in the line amplitude effectively changes the measured position of the lines, affecting the stability of the wavelength scale of the instrument. We deduce and apply the correcting functions to re-calibrate and mitigate this effect, reducing it to a level consistent with photon noise.

preprint2020arXiv

Optical Observations Reveal Strong Evidence for High Energy Neutrino Progenitor

We present the earliest astronomical observation of a high energy neutrino error box in which its variability was discovered after high-energy neutrinos detection. The one robotic telescope of the MASTER global international network (Lipunov et al. 2010) automatically imaged the error box of the very high-energy neutrino event IceCube-170922A. Observations were carried out in minute after the IceCube-170922A neutrino event was detected by the IceCube observatory at the South Pole. MASTER found the blazar TXS 0506+056 to be in the off-state after one minute and then switched to the on-state no later than two hours after the event. The effect is observed at a 50-sigma significance level. Also we present own unique 16-years light curve of blazar TXS 0506+056 (518 data set).

preprint2020arXiv

Revisiting Proxima with ESPRESSO

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days.

preprint2020arXiv

WASP-127b: A misaligned planet with a partly cloudy atmosphere and tenuous sodium signature seen by ESPRESSO

The study of exoplanet atmospheres is essential to understand the formation, evolution and composition of exoplanets. The transmission spectroscopy technique is playing a significant role in this domain. In particular, the combination of state-of-the-art spectrographs at low- and high-spectral resolution is key to our understanding of atmospheric structure and composition. Two transits of the close-in sub Saturn-mass planet,WASP-127b, have been observed with ESPRESSO in the frame of the Guaranteed Time Observations Consortium. Transit observations allow us to study simultaneously the system architecture and the exoplanet atmosphere. We found that this planet is orbiting its slowly rotating host star (veq sin(i)=0.53+/-0.07 km/s) on a retrograde misaligned orbit (lambda=-128.41+/-5.60 deg). We detected the sodium line core at the 9-sigma confidence level with an excess absorption of 0.3+/-0.04%, a blueshift of 2.7+/-0.79 km/s and a FWHM of 15.18+/-1.75 km/s. However, we did not detect the presence of other atomic species but set upper-limits of only few scale heights. Finally, we put a 3-sigma upper limit, to the average depth of the 1600 strongest water lines at equilibrium temperature in the visible band, of 38 ppm. This constrains the cloud-deck pressure between 0.3 and 0.5 mbar by combining our data with low-resolution data in the near-infrared and models computed for this planet. To conclude, WASP-127b, with an age of about 10 Gyr, is an unexpected exoplanet by its orbital architecture but also by the small extension of its sodium atmosphere (~7 scale heights). ESPRESSO allows us to take a step forward in the detection of weak signals, thus bringing strong constraints on the presence of clouds in exoplanet atmospheres. The framework proposed in this work can be applied to search for molecular species and study cloud-decks in other exoplanets.

preprint2019arXiv

A transiting super-Earth close to the inner edge of the habitable zone of an M0 dwarf star

We present a super-Earth orbiting close to the inner edge of the habitable zone of the cool dwarf star K2-286 (EPIC 249889081), detected with data from the K2 mission in its $15^{th}$ campaign. The planet has radius of $2.1\pm0.2$ R$_{\oplus}$, near the 1.5 - 2.0 R$_{\oplus}$ gap in the radii distribution. The equilibrium temperature is $347^{+21}_{-11}$ K, cooler than most of the small planets with well measured masses, and the orbital period is $27.359\pm0.005$ days. K2-286, located at a distance of $76.3\pm0.3$ pc, is an M0V star with estimated effective temperature of $3926\pm100$ K, less active than other M dwarf stars hosting exoplanets. The expected radial velocity semi-amplitude induced by the planet on the star is $1.9^{+1.3}_{-1.2}$ m$\cdot$s$^{-1}$, and the amplitude of signals in transit transmission spectroscopy is estimated at $5.0\pm3.0$ ppm. Follow-up observations for mass measurements and transit spectroscopy should be desirable for this relatively bright target ($m_V=12.76, m_{Ks}=9.32$) hosting a transiting super-Earth within the inner edge of the habitable zone.

preprint2019arXiv

The KISS experiment

Mapping millimetre continuum emission has become a key issue in modern multi-wavelength astrophysics. In particular, spectrum-imaging at low frequency resolution is an asset for characterizing the clusters of galaxies via the Sunyaev Zeldovich (SZ) effect. In this context, we have built a ground-based spectrum-imager named KIDs Interferometer Spectrum Survey (KISS). This instrument is based on two 316-pixel arrays of Kinetic Inductance Detectors (KID) cooled to 150 mK by a custom dilution refrigerator-based cryostat. By using Ti-Al and Al absorbers, we can cover a wide frequency range between 80 and 300 GHz. In order to preserve a large instantaneous Field of View (FoV) 1 degree the spectrometer is based on a Fourier Transform interferometer. This represents a technological challenge due to the fast scanning speed that is needed to overcome the effects of background atmospheric fluctuations. KISS is installed at the QUIJOTE 2.25 m telescope in Tenerife since February 2019 and is currently in its commissioning phase. In this proceeding we present an overview of the instrument and the latest results.

preprint2016arXiv

Extremely fast orbital decay of the black hole X-ray binary Nova Muscae 1991

We present new medium-resolution spectroscopic observations of the black hole X-ray binary Nova Muscae 1991 taken with X-Shooter spectrograph installed at the 8.2m-VLT telescope. These observations allow us to measure the time of inferior conjunction of the secondary star with the black hole in this system that, together with previous measurements, yield an orbital period decay of $\dot P=-20.7\pm12.7$ ms yr$^{-1}$ ($-24.5\pm15.1$ $μ$s per orbital cycle). This is significantly faster than those previously measured in the other black hole X-ray binaries A0620-00 and XTE J1118+480. No standard black hole X-ray binary evolutionary model is able to explain this extremely fast orbital decay. At this rate, the secondary star would reach the event horizon (as given by the Schwarzschild radius of about 32 km) in roughly 2.7 Myr. This result has dramatic implications on the evolution and lifetime of black hole X-ray binaries.

preprint2016arXiv

Follow-up observations of extremely metal-poor stars identified from SDSS

The most metal-poor stars in the Milky Way witnessed the early phases of formation of the Galaxy, and have chemical compositions that are close to the pristine mixture from Big Bang nucleosynthesis, polluted by one or few supernovae. Only two dozen stars with ([Fe/H]< -4) are known, and they show a wide range of abundance patterns. It is therefore important to enlarge this sample. We present the first results of an effort to identify new extremely metal-poor stars in the Milky Way halo. Our targets have been selected from low-resolution spectra obtained as part of the Sloan Digital Sky Survey, and followed-up with medium resolution spectroscopy on the 4.2 m William Herschel Telescope and, in a few cases, at high resolution on the the 9.2 m Hobby-Eberly Telescope. Stellar parameters and the abundances of magnesium, calcium, iron, and strontium have been inferred from the spectra using classical model atmospheres. We have also derived carbon abundances from the G band. We find consistency between the metallicities estimated from SDSS and those from new data at the level of 0.3 dex. The analysis of medium resolution data obtained with ISIS on the WHT allow us to refine the metallicities and in some cases measure other elemental abundances. Our sample contains 11 new metal-poor stars with [Fe/H] < -3.0, one of them with an estimated metallicity of [Fe/H] ~ -4.0. We also discuss metallicity discrepancies of some stars in common with previous works in the literature. Only one of these stars is found to be C-enhanced at about [C/Fe]~+1, whereas the other metal-poor stars show C abundances at the level of [C/Fe]~+0.45.

preprint2016arXiv

High spatial resolution optical imaging of the multiple T Tauri system LkHα 262/LkHα 263

We report high spatial resolution i' band imaging of the multiple T Tauri system LkH$α$ 262/LkH$α$ 263 obtained during the first commissioning period of the Adaptive Optics Lucky Imager (AOLI) at the 4.2 m William Herschel Telescope, using its Lucky Imaging mode. AOLI images have provided photometry for each of the two components LkH$α$ 263 A and B (0.41 arcsec separation) and marginal evidence for an unresolved binary or a disc in LkH$α$ 262. The AOLI data combined with previously available and newly obtained optical and infrared imaging show that the three components of LkH$α$ 263 are co-moving, that there is orbital motion in the AB pair, and, remarkably, that LkH$α$ 262-263 is a common proper motion system with less than 1 mas/yr relative motion. We argue that this is a likely five-component gravitationally bounded system. According to BT-settl models the mass of each of the five components is close to 0.4 M$_{\odot}$ and the age is in the range 1-2 Myr. The presence of discs in some of the components offers an interesting opportunity to investigate the formation and evolution of discs in the early stages of multiple very low-mass systems. In particular, we provide tentative evidence that the disc in 263C could be coplanar with the orbit of 263AB.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Magnetic Cycles and Rotation Periods of Late Type Stars from photometric time series

We investigate the photometric modulation induced by magnetic activity cycles and study the relationship between rotation period and activity cycle(s) in late-type (FGKM) stars. We analyse light-curves spanning up to 9 years of 125 nearby stars provided by the ASAS survey. The sample is mainly conformed by low-activity main sequence late A to mid M-type stars. A search is performed for short (days) and long-term (years) periodic variations in the photometry. We modelled with combinations of sinusoids the light-curves to measure the properties of these periodic signals. To provide a better statistical interpretation of our results we complement them with the results from previous similar works. We have been able to measure long-term photometric cycles of 47 stars. Rotational modulation was also detected and rotational periods measured in 36 stars. For 28 stars we have simultaneous measurements of both, activity cycles and rotational periods, being 17 of them M-type stars. From sinusoidal fits we measured both photometric amplitudes and periods. The measured cycle periods range from 2 up to 14 yr with photometric amplitudes in the range of 5-20 mmag. We have found that the distribution of cycle lengths for the different spectral types is similar. On the other hand the distribution of rotation periods is completely different, trending to longer periods for later type stars. The amplitudes induced by magnetic cycles and rotation show a clear correlation. A trend of photometric amplitudes with rotation period is also outlined in the data. For a given activity index the amplitudes of the photometric variability induced by activity cycles of main sequence GK stars are lower than those of early and mid-M dwarfs. Using spectroscopic data we also provide an update in the empirical relationship between the level of chromospheric activity as given by log(Rhk) and the rotation periods.

preprint2016arXiv

Near-infrared colors of minor planets recovered from VISTA - VHS survey (MOVIS)

The Sloan Digital Sky Survey (SDSS) and Wide-field Infrared Survey Explorer (WISE) provide information about the surface composition of about 100,000 minor planets. The resulting visible colors and albedos enabled us to group them in several major classes, which are a simplified view of the diversity shown by the few existing spectra. We performed a serendipitous search in VISTA-VHS observations using a pipeline developed to retrieve and process the data that corresponds to solar system objects (SSo). The colors and the magnitudes of the minor planets observed by the VISTA survey are compiled into three catalogs that are available online: the detections catalog (MOVIS-D), the magnitudes catalog (MOVIS-M), and the colors catalog (MOVIS-C). They were built using the third data release of the survey (VISTA VHS-DR3). A total of 39,947 objects were detected, including 52 NEAs, 325 Mars Crossers, 515 Hungaria asteroids, 38,428 main-belt asteroids, 146 Cybele asteroids, 147 Hilda asteroids, 270 Trojans, 13 comets, 12 Kuiper Belt objects and Neptune with its four satellites. The colors found for asteroids with known spectral properties reveal well-defined patterns corresponding to different mineralogies. The distributions of MOVIS-C data in color-color plots shows clusters identified with different taxonomic types. All the diagrams that use (Y-J) color separate the spectral classes more effectively than the (J-H) and (H-Ks) plots used until now: even for large color errors (<0.1), the plots (Y-J) vs (Y-Ks) and (Y-J) vs (J-Ks) provide the separation between S-complex and C-complex. The end members A, D, R, and V-types occupy well-defined regions.

preprint2016arXiv

Planck 2015 results. II. Low Frequency Instrument data processing

We present an updated description of the Planck Low Frequency Instrument (LFI) data processing pipeline, associated with the 2015 data release. We point out the places where our results and methods have remained unchanged since the 2013 paper and we highlight the changes made for the 2015 release, describing the products (especially timelines) and the ways in which they were obtained. We demonstrate that the pipeline is self-consistent (principally based on simulations) and report all null tests. For the first time, we present LFI maps in Stokes Q and U polarization. We refer to other related papers where more detailed descriptions of the LFI data processing pipeline may be found if needed.

preprint2016arXiv

Planck 2015 results. XIII. Cosmological parameters

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

preprint2016arXiv

Planck 2015 results. XIV. Dark energy and modified gravity

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

preprint2016arXiv

Planck 2015 results. XIX. Constraints on primordial magnetic fields

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfvén waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

preprint2016arXiv

Planck 2015 results. XV. Gravitational lensing

We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 sigma), using temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400 and an associated likelihood for cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in the best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level measurement of the parameter combination $σ_8 Ω_m^{0.25} = 0.591\pm 0.021$. We combine our determination of the lensing potential with the E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 sigma, confirming Planck's sensitivity to this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.

preprint2016arXiv

Planck 2015 results. XVI. Isotropy and statistics of the CMB

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

preprint2016arXiv

Planck 2015 results. XVII. Constraints on primordial non-Gaussianity

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

preprint2016arXiv

Planck 2015 results. XVIII. Background geometry & topology

Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance $χ_{rec}$). We specialize to flat spaces with toroidal and slab topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $Δ\ln{L}>-5$ relative to a simply-connected flat Planck best-fit model) are $R_i>0.97χ_{rec}$ for the cubic torus and $R_i>0.56χ_{rec}$ for the slab. The limit for the cubic torus from the matched-circles search is numerically equivalent, $R_i>0.97χ_{rec}$ (99% CL) from polarisation data alone. We also perform a Bayesian search for a Bianchi VII$_h$ geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.1\pm0.04$ compared to +1 if the model were to be correct. In the physical setting where the Bianchi parameters are fit simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models to $(ω/H)_0<7.6\times10^{-10}$ (95% CL). [Abridged]

preprint2016arXiv

Planck 2015 results. XXI. The integrated Sachs-Wolfe effect

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,σ$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,σ$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $Ω_Λ$ at more than $3\,σ$. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, polarization data is used to study the anomalously large ISW signal previously reported through the stacking of CMB features at the locations of known superstructures. We find that the current Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the LSS traced by the 2MPZ survey by directly inverting the density field into the gravitational potential field.

preprint2016arXiv

Planck 2015 results. XXV. Diffuse low-frequency Galactic foregrounds

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

preprint2016arXiv

Planck 2015 results. XXVI. The Second Planck Catalogue of Compact Sources

The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.

preprint2016arXiv

Planck intermediate results. XLIII. The spectral energy distribution of dust in clusters of galaxies

Although infrared (IR) overall dust emission from clusters of galaxies has been statistically detected using data from the Infrared Astronomical Satellite (IRAS), it has not been possible to sample the spectral energy distribution (SED) of this emission over its peak, and thus to break the degeneracy between dust temperature and mass. By complementing the IRAS spectral coverage with Planck satellite data from 100 to 857 GHz, we provide new constraints on the IR spectrum of thermal dust emission in clusters of galaxies. We achieve this by using a stacking approach for a sample of several hundred objects from the Planck cluster sample; this procedure averages out fluctuations from the IR sky, allowing us to reach a significant detection of the faint cluster contribution. We also use the large frequency range probed by Planck, together with component-separation techniques, to remove the contamination from both cosmic microwave background anisotropies and the thermal Sunyaev-Zeldovich effect (tSZ) signal, which dominate below 353 GHz. By excluding dominant spurious signals or systematic effects, averaged detections are reported at frequencies between 353 and 5000 GHz. We confirm the presence of dust in clusters of galaxies at low and intermediate redshifts, yielding an SED with a shape similar to that of the Milky Way. Planck's beam does not allow us to investigate the detailed spatial distribution of this emission (e.g., whether it comes from intergalactic dust or simply the dust content of the cluster galaxies), but the radial distribution of the emission appears to follow that of the stacked SZ signal, and thus the extent of the clusters. The recovered SED allows us to constrain the dust mass responsible for the signal, as well as its temperature. We additionally explore the evolution of the IR emission as a function of cluster mass and redshift.

preprint2016arXiv

Planck intermediate results. XLIV. The structure of the Galactic magnetic field from dust polarization maps of the southern Galactic cap

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($ψ$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $ψ$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $ψ$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

preprint2016arXiv

Planck intermediate results. XXIX. All-sky dust modelling with Planck, IRAS, and WISE observations

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

preprint2016arXiv

Planck Intermediate Results. XXXVI. Optical identification and redshifts of Planck SZ sources with telescopes in the Canary Islands Observatories

We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with telescopes at the Canary Islands observatories as part of the general optical follow-up programme undertaken by the Planck collaboration. In total, 78 SZ sources are discussed. Deep-imaging observations were obtained for most of these sources; spectroscopic observations in either in long-slit or multi-object modes were obtained for many. We effectively used 37.5 clear nights. We found optical counterparts for 73 of the 78 candidates. This sample includes 53 spectroscopic redshift determinations, 20 of them obtained with a multi-object spectroscopic mode. The sample contains new redshifts for 27 Planck clusters that were not included in the first Planck SZ source catalogue (PSZ1).

preprint2016arXiv

QUIJOTE Scientific Results. II. Polarisation Measurements of the Microwave Emission in the Galactic molecular complexes W43 and W47 and supernova remnant W44

We present Q-U-I JOint TEnerife (QUIJOTE) intensity and polarisation maps at 10-20 GHz covering a region along the Galactic plane 24<l<45 deg, |b|<8 deg. These maps result from 210 h of data, have a sensitivity in polarisation of ~40 muK/beam and an angular resolution of ~1 deg. Our intensity data are crucial to confirm the presence of anomalous microwave emission (AME) towards the two molecular complexes W43 (22 sigma) and W47 (8 sigma). We also detect at high significance (6 sigma) AME associated with W44, the first clear detection of this emission towards a SNR. The new QUIJOTE polarisation data, in combination with WMAP, are essential to: i) Determine the spectral index of the synchrotron emission in W44, beta_sync =-0.62 +/-0.03, in good agreement with the value inferred from the intensity spectrum once a free-free component is included in the fit. ii) Trace the change in the polarisation angle associated with Faraday rotation in the direction of W44 with rotation measure -404 +/- 49 rad/m2. And iii) set upper limits on the polarisation of W43 of Pi_AME <0.39 per cent (95 per cent C.L.) from QUIJOTE 17~GHz, and <0.22 per cent from WMAP 41 GHz data, which are the most stringent constraints ever obtained on the polarisation fraction of the AME. For typical physical conditions (grain temperature and magnetic field strengths), and in the case of perfect alignment between the grains and the magnetic field, the models of electric or magnetic dipole emissions predict higher polarisation fractions.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2015arXiv

A Joint Analysis of BICEP2/Keck Array and Planck Data

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5°$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $μ$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$Λ$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,σ$ significance.

preprint2015arXiv

A search for lithium in metal-poor L dwarfs

The aim of the project is to search for lithium in absorption at 6707.8 Angstroms to constrain the nature and the mass of the brightest low-metallicity L-type dwarfs (refered to as L subdwarfs) identified in large-scale surveys. We obtained low- to intermediate-resolution (R~2500-9000) optical (~560-770 nm) spectra of two mid-L subdwarfs, SDSSJ125637.13-022452.4 (SDSS1256; sdL3.5) and 2MASSJ162620.14+392519.5 (2MASS1626; sdL4) with spectrographs on the European Southern Observatory Very Large Telescope and the Gran Telescopio de Canarias. We report the presence of a feature at the nominal position of the lithium absorption doublet at 6707.8 Angstroms in the spectrum of SDSS1256, with an equivalent width of 66+/-27 Angstroms at 2.4 sigma, which we identify as arising from a CaH molecular transition based on atmosphere models. We do not see any feature at the position of the lithium feature in the spectrum of 2MASS1626. The existence of overlapping molecular absorption sets a confusion detection limit of [Li/H]=-3 for equivalently-typed L subdwarfs. We provided improved radial velocity measurements of -126+/-10 km/s and -239+/-12 km/s for SDSS1256 and 2MASS1626, respectively, as well as revised Galactic orbits. We implemented adjusting factors for the CaH molecule in combination with the NextGen atmosphere models to fit the optical spectrum of SDSS1256 in the 6200-7300 Angstroms range. We also estimate the expected Li abundance from interstellar accretion ([Li/H]=-5), place limits on circumstellar accretion (10^9 g/yr), and discuss the prospects of Li searches in cooler L and T subdwarfs.

preprint2015arXiv

Adaptive Optics and Lucky Imager (AOLI): presentation and first light

In this paper we present the Adaptive Optics Lucky Imager (AOLI), a state-of-the-art instrument which makes use of two well proved techniques for extremely high spatial resolution with ground-based telescopes: Lucky Imaging (LI) and Adaptive Optics (AO). AOLI comprises an AO system, including a low order non-linear curvature wavefront sensor together with a 241 actuators deformable mirror, a science array of four 1024x1024 EMCCDs, allowing a 120x120 down to 36x36 arcseconds field of view, a calibration subsystem and a powerful LI software. Thanks to the revolutionary WFS, AOLI shall have the capability of using faint reference stars ({\it I\/} $\sim$ 16.5-17.5), enabling it to be used over a much wider part of the sky than with common Shack-Hartmann AO systems. This instrument saw first light in September 2013 at William Herschel Telescope. Although the instrument was not complete, these commissioning demonstrated its feasibility, obtaining a FWHM for the best PSF of 0.151$\pm$0.005 arcsec and a plate scale of 55.0$\pm$0.3 mas/pixel. Those observations served us to prove some characteristics of the interesting multiple T Tauri system LkH$α$ 262-263, finding it to be gravitationally bounded. This interesting multiple system mixes the presence of proto-planetary discs, one proved to be double, and the first-time optically resolved pair LkH$α$ 263AB (0.42 arcsec separation).

preprint2015arXiv

An eclipsing double-line spectroscopic binary at the stellar/substellar boundary in the Upper Scorpius OB association

We aim at constraining evolutionary models at low mass and young ages by identifying interesting transiting system members of the nearest OB association to the Sun, Upper Scorpius, targeted by the Kepler mission. We produced light curves for M dwarf members of the USco region surveyed during the second campaign of the Kepler K2 mission. We identified 'by eye' a transiting system, UScoJ161630.68-251220.1 (=EPIC203710387) with a combined spectral type of M5.25 whose photometric, astrometric, and spectroscopic properties makes it a member of USco. We conducted an extensive photometric and spectroscopic follow-up of this transiting system with a suite of telescopes and instruments to characterise the properties of each component of the system. We calculated a transit duration of about 2.42 hours occuring every 2.88 days with a slight difference in transit depth and phase between the two components. We estimated a mass ratio of 0.922+/-0.015 from the semi-amplitudes of the radial velocity curves for each component. We derived masses of 0.091+/-0.005 Msun and 0.084+/-0.004 Msun,radii of 0.388+/-0.008 Rsun and 0.380+/-0.008 Rsun, luminosities of log(L/Lsun)=-2.020 (-0.121+0.099) dex and -2.032 (-0.121+0.099) dex, and effective temperatures of 2901 (-172+199) K and 2908 (-172+199) K for the primary and secondary, respectively. We present a complete photometric and radial velocity characterisation of the least massive double-line eclipsing binary system in the young USco association with two components close to the stellar/substellar limit. This system fills in a gap between the least massive eclipsing binaries in the low-mass and substellar regimes at young ages and represents an important addition to constrain evolutionary models at young ages.

preprint2015arXiv

An equatorial ultra iron-poor star identified in BOSS

We report the discovery of SDSS J131326.89-001941.4, an ultra iron-poor red giant star ([Fe/H] ~ -4.3) with a very high carbon abundance ([C/Fe]~ +2.5). This object is the fifth star in this rare class, and the combination of a fairly low effective temperature (Teff ~ 5300 K), which enhances line absorption, with its brightness (g=16.9), makes it possible to measure the abundances of calcium, carbon and iron using a low-resolution spectrum from the Sloan Digital Sky Survey. We examine the carbon and iron abundance ratios in this star and other similar objects in the light of predicted yields from metal-free massive stars, and conclude that they are consistent. By way of comparison, stars with similarly low iron abundances but lower carbon-to-iron ratios deviate from the theoretical predictions.

preprint2015arXiv

Constraints on the substellar companions in wide orbits around the Barnard's Star from CanariCam mid-infrared imaging

We have performed mid-infrared imaging of Barnard's Star, one of the nearest stars to the Sun, using CanariCam on the 10.4 m Gran Telescopio Canarias. We aim to investigate an area within 1-10 arcsec separations, which for the 1.83 pc distance of the star translates to projected orbital separations of 1.8-18 AU (P > 12 yr), which have not been explored yet with astrometry or radial velocity programs. It is therefore an opportunity to enter the domain of distances where most giant planets are expected to form. We performed deep imaging in the N-band window (Si-2 filter, 8.7 μm) reaching a 3σ detection limit of 0.85+/-0.18 mJy and angular resolution of 0.24 arcsec, close to the diffraction limit of the telescope at this wavelength. A total of 80 min on-source integration time data were collected and combined for the deepest image. We achieved a dynamical range of 8.0+/-0.1 mag in the 8.7 μm band, at angular separations from ~2 to 10 arcsec and of ~6-8 mag at 1-2 arcsec. No additional sources were found. Our detectability limits provide further constraints to the presence of substellar companions of the Barnard's Star. According to solar metallicity evolutionary models, we can exclude companions of masses larger than 15 MJup (Teff > 400 K), ages of a few Gyr, and located in ~3.6-18 AU orbits with a 3σ confidence level. This minimum mass is approximately 5 MJup smaller than any previous imaging survey that explored the surroundings of Barnard's Star could restrict.

preprint2015arXiv

Liverpool Telescope 2: a new robotic facility for rapid transient follow-up

The Liverpool Telescope is one of the world's premier facilities for time domain astronomy. The time domain landscape is set to radically change in the coming decade, with surveys such as LSST providing huge numbers of transient detections on a nightly basis; transient detections across the electromagnetic spectrum from other facilities such as SVOM, SKA and CTA; and the era of `multi-messenger astronomy', wherein events are detected via non-electromagnetic means, such as gravitational wave emission. We describe here our plans for Liverpool Telescope 2: a new robotic telescope designed to capitalise on this new era of time domain astronomy. LT2 will be a 4-metre class facility co-located with the LT at the Observatorio del Roque de Los Muchachos on the Canary island of La Palma. The telescope will be designed for extremely rapid response: the aim is that the telescope will take data within 30 seconds of the receipt of a trigger from another facility. The motivation for this is twofold: firstly it will make it a world-leading facility for the study of fast fading transients and explosive phenomena discovered at early times. Secondly, it will enable large-scale programmes of low-to-intermediate resolution spectral classification of transients to be performed with great efficiency. In the target-rich environment of the LSST era, minimising acquisition overheads will be key to maximising the science gains from any follow-up programme. The telescope will have a diverse instrument suite which is simultaneously mounted for automatic changes, but it is envisaged that the primary instrument will be an intermediate resolution, optical/infrared spectrograph for scientific exploitation of transients discovered with the next generation of synoptic survey facilities. In this paper we outline the core science drivers for the telescope, and the requirements for the optical and mechanical design.

preprint2015arXiv

Measurements of the Intensity and Polarization of the Anomalous Microwave Emission in the Perseus molecular complex with QUIJOTE

Anomalous microwave emission (AME) has been observed in numerous sky regions, in the frequency range ~10-60 GHz. One of the most scrutinized regions is G159.6-18.5, located within the Perseus molecular complex. In this paper we present further observations of this region (194 hours in total over ~250 deg^2), both in intensity and in polarization. They span four frequency channels between 10 and 20 GHz, and were gathered with QUIJOTE, a new CMB experiment with the goal of measuring the polarization of the CMB and Galactic foregrounds. When combined with other publicly-available intensity data, we achieve the most precise spectrum of the AME measured to date, with 13 independent data points being dominated by this emission. The four QUIJOTE data points provide the first independent confirmation of the downturn of the AME spectrum at low frequencies, initially unveiled by the COSMOSOMAS experiment in this region. We accomplish an accurate fit of these data using models based on electric dipole emission from spinning dust grains, and also fit some of the parameters on which these models depend. We also present polarization maps with an angular resolution of ~1 deg and a sensitivity of ~25 muK/beam. From these maps, which are consistent with zero polarization, we obtain upper limits of Pi<6.3% and <2.8% (95% C.L.) respectively at 12 and 18 GHz, a frequency range where no AME polarization observations have been reported to date. These constraints are compatible with theoretical predictions of the polarization fraction from electric dipole emission originating from spinning dust grains. At the same time, they rule out several models based on magnetic dipole emission from dust grains ordered in a single magnetic domain, which predict higher polarization levels. Future QUIJOTE data in this region may allow more stringent constraints on the polarization level of the AME.

preprint2015arXiv

Planck 2013 results. XXII. Constraints on inflation

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $σ$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $χ^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $χ^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

preprint2015arXiv

Planck 2013 results. XXIX. The Planck catalogue of Sunyaev-Zeldovich sources: Addendum

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

preprint2015arXiv

Planck 2015 results. I. Overview of products and scientific results

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

preprint2015arXiv

Planck 2015 results. IV. Low Frequency Instrument beams and window functions

This paper presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). The structure of the paper is similar to that presented in the 2013 Planck release; the main differences concern the beam normalization and the delivery of the window functions to be used for polarization analysis. The in-flight assessment of the LFI main beams relies on measurements performed during observations of Jupiter. By stacking data from seven Jupiter transits, the main beam profiles are measured down to -25 dB at 30 and 44 GHz, and down to -30 dB at 70 GHz. It has been confirmed that the agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band (within the 20 dB contour from the peak, the rms values are 0.1% at 30 and 70 GHz; 0.2% at 44 GHz). Simulated polarized beams are used for the computation of the effective beam window functions. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer band shapes. The total uncertainties in the effective beam window functions are 0.7% and 1% at 30 and 44 GHz, respectively (at $\ell \approx 600$); and 0.5% at 70 GHz (at $\ell \approx 1000$).

preprint2015arXiv

Planck 2015 results. IX. Diffuse component separation: CMB maps

We present foreground-reduced CMB maps derived from the full Planck data set in both temperature and polarization. Compared to the corresponding Planck 2013 temperature sky maps, the total data volume is larger by a factor of 3.2 for frequencies between 30 and 70 GHz, and by 1.9 for frequencies between 100 and 857 GHz. In addition, systematic errors in the forms of temperature-to-polarization leakage, analogue-to-digital conversion uncertainties, and very long time constant errors have been dramatically reduced, to the extent that the cosmological polarization signal may now be robustly recovered on angular scales $\ell\gtrsim40$. On the very largest scales, instrumental systematic residuals are still non-negligible compared to the expected cosmological signal, and modes with $\ell < 20$ are accordingly suppressed in the current polarization maps by high-pass filtering. As in 2013, four different CMB component separation algorithms are applied to these observations, providing a measure of stability with respect to algorithmic and modelling choices. The resulting polarization maps have rms instrumental noise ranging between 0.21 and 0.27$\,μ\textrm{K}$ averaged over 55 arcmin pixels, and between 4.5 and 6.1$\,μ\textrm{K}$ averaged over 3.4 arcmin pixels. The cosmological parameters derived from the analysis of temperature power spectra are in agreement at the $1σ$ level with the Planck 2015 likelihood. Unresolved mismatches between the noise properties of the data and simulations prevent a satisfactory description of the higher-order statistical properties of the polarization maps. Thus, the primary applications of these polarization maps are those that do not require massive simulations for accurate estimation of uncertainties, for instance estimation of cross-spectra and cross-correlations, or stacking analyses.

preprint2015arXiv

Planck 2015 results. V. LFI calibration

We present a description of the pipeline used to calibrate the Planck Low Frequency Instrument (LFI) timelines into thermodynamic temperatures for the Planck 2015 data release, covering four years of uninterrupted operations. As in the 2013 data release, our calibrator is provided by the spin-synchronous modulation of the cosmic microwave background dipole, but we now use the orbital component, rather than adopting the Wilkinson Microwave Anisotropy Probe (WMAP) solar dipole. This allows our 2015 LFI analysis to provide an independent Solar dipole estimate, which is in excellent agreement with that of HFI and within $1σ$ (0.3% in amplitude) of the WMAP value. This 0.3% shift in the peak-to-peak dipole temperature from WMAP and a global overhaul of the iterative calibration code increases the overall level of the LFI maps by 0.45% (30 GHz), 0.64% (44 GHz), and 0.82% (70 GHz) in temperature with respect to the 2013 Planck data release, thus reducing the discrepancy with the power spectrum measured by WMAP. We estimate that the LFI calibration uncertainty is now at the level of 0.20% for the 70 GHz map, 0.26% for the 44 GHz map, and 0.35% for the 30 GHz map. We provide a detailed description of the impact of all the changes implemented in the calibration since the previous data release.

preprint2015arXiv

Planck 2015 results. VI. LFI mapmaking

This paper describes the mapmaking procedure applied to Planck LFI (Low Frequency Instrument) data. The mapmaking step takes as input the calibrated timelines and pointing information. The main products are sky maps of $I,Q$, and $U$ Stokes components. For the first time, we present polarization maps at LFI frequencies. The mapmaking algorithm is based on a destriping technique, enhanced with a noise prior. The Galactic region is masked to reduce errors arising from bandpass mismatch and high signal gradients. We apply horn-uniform radiometer weights to reduce effects of beam shape mismatch. The algorithm is the same as used for the 2013 release, apart from small changes in parameter settings. We validate the procedure through simulations. Special emphasis is put on the control of systematics, which is particularly important for accurate polarization analysis. We also produce low-resolution versions of the maps, and corresponding noise covariance matrices. These serve as input in later analysis steps and parameter estimation. The noise covariance matrices are validated through noise Monte Carlo simulations. The residual noise in the map products is characterized through analysis of half-ring maps, noise covariance matrices, and simulations.

preprint2015arXiv

Planck 2015 results. X. Diffuse component separation: Foreground maps

Planck has mapped the microwave sky in nine frequency bands between 30 and 857 GHz in temperature and seven bands between 30 and 353 GHz in polarization. In this paper we consider the problem of diffuse astrophysical component separation, and process these maps within a Bayesian framework to derive a consistent set of full-sky astrophysical component maps. For the temperature analysis, we combine the Planck observations with the 9-year WMAP sky maps and the Haslam et al. 408 MHz map to derive a joint model of CMB, synchrotron, free-free, spinning dust, CO, line emission in the 94 and 100 GHz channels, and thermal dust emission. Full-sky maps are provided with angular resolutions varying between 7.5 arcmin and 1 deg. Global parameters (monopoles, dipoles, relative calibration, and bandpass errors) are fitted jointly with the sky model, and best-fit values are tabulated. For polarization, the model includes CMB, synchrotron, and thermal dust emission. These models provide excellent fits to the observed data, with rms temperature residuals smaller than 4 uK over 93% of the sky for all Planck frequencies up to 353 GHz, and fractional errors smaller than 1% in the remaining 7% of the sky. The main limitations of the temperature model at the lower frequencies are degeneracies among the spinning dust, free-free, and synchrotron components; additional observations from external low-frequency experiments will be essential to break these. The main limitations of the temperature model at the higher frequencies are uncertainties in the 545 and 857 GHz calibration and zero-points. For polarization, the main outstanding issues are instrumental systematics in the 100-353 GHz bands on large angular scales in the form of temperature-to-polarization leakage, uncertainties in the analog-to-digital conversion, and very long time constant corrections, all of which are expected to improve in the near future.

preprint2015arXiv

Planck 2015 results. XII. Full Focal Plane simulations

We present the 8th Full Focal Plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising $10^4$ mission realizations reduced to about $10^6$ maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects; remaining subdominant effects will be included in future updates. Generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used for the validation and verification of analysis algorithms, as well as their implementations, and for removing biases from and quantifying uncertainties in the results of analyses of the real data.

preprint2015arXiv

Planck 2015 results. XXVII. The Second Planck Catalogue of Sunyaev-Zeldovich Sources

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

preprint2015arXiv

Planck 2015 results. XXVIII. The Planck Catalogue of Galactic Cold Clumps

We present the Planck Catalogue of Galactic Cold Clumps (PGCC), an all-sky catalogue of Galactic cold clump candidates detected by Planck. This catalogue is the full version of the Early Cold Core (ECC) catalogue, which was made available in 2011 with the Early Release Compact Source Catalogue (ERCSC) and contained 915 high S/N sources. It is based on the Planck 48 months mission data that are currently being released to the astronomical community. The PGCC catalogue is an observational catalogue consisting exclusively of Galactic cold sources. The three highest Planck bands (857, 545, 353 GHz) have been combined with IRAS data at 3 THz to perform a multi-frequency detection of sources colder than their local environment. After rejection of possible extragalactic contaminants, the PGCC catalogue contains 13188 Galactic sources spread across the whole sky, i.e., from the Galactic plane to high latitudes, following the spatial distribution of the main molecular cloud complexes. The median temperature of PGCC sources lies between 13 and 14.5 K, depending on the quality of the flux density measurements, with a temperature ranging from 5.8 to 20 K after removing sources with the 1% largest temperature estimates. Using seven independent methods, reliable distance estimates have been obtained for 5574 sources, which allows us to derive their physical properties such as their mass, physical size, mean density and luminosity. The PGCC sources are located mainly in the solar neighbourhood, up to a distance of 10.5 kpc towards the Galactic centre, and range from low-mass cores to large molecular clouds. Because of this diversity and because the PGCC catalogue contains sources in very different environments, the catalogue is useful to investigate the evolution from molecular clouds to cores. Finally, the catalogue also includes 54 additional sources located in the SMC and LMC.

preprint2015arXiv

Planck intermediate results. XXI. Comparison of polarized thermal emission from Galactic dust at 353 GHz with interstellar polarization in the visible

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, $p_V$ , and the optical depth in the V band to the star, $τ_V$. Toward these stars we measure the submillimetre polarized intensity, $P_S$, and total intensity, $I_S$, in the Planck 353 GHz channel. For those lines of sight through the diffuse interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, $R_{S/V} = (P_S/I_S)/(p_V/τ_V)$ and $R_{P/p} = P_S/p_V$ , the latter focusing directly on the polarization properties of the aligned grain population alone. We find $R_{S/V}$ = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and $R_{P/p}$ = 5.4 MJy sr$^{-1}$, with uncertainties 0.2 and 0.3 MJy sr$^{-1}$, respectively. Our estimate of $R_{S/V}$ is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. However, our estimate of $R_{P/p}$ is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, $R_{P/p}$, indicates that changes to the optical properties in the models of the aligned grain population are required.

preprint2015arXiv

Planck intermediate results. XXIII. Galactic plane emission components derived from Planck with ancillary data

Planck data when combined with ancillary data provide a unique opportunity to separate the diffuse emission components of the inner Galaxy. The purpose of the paper is to elucidate the morphology of the various emission components in the strong star-formation region lying inside the solar radius and to clarify the relationship between the various components. The region of the Galactic plane covered is l=300-0-60deg where star-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1deg and 2deg, which correspond to FWHM z-widths of 100-200pc at a typical distance of 6kpc. The four emission components studied here are synchrotron, free-free, anomalous microwave emission (AME), and thermal (vibrational) dust emission. These components are identified by constructing spectral energy distributions (SEDs) at positions along the Galactic plane using the wide frequency coverage of Planck (28.4-857GHz) in combination with low-frequency radio data at 0.408-2.3GHz plus WMAP data at 23-94GHz, along with far-infrared (FIR) data from DIRBE and IRAS. The free-free component is determined from radio recombination line (RRL) data. AME is found to be comparable in brightness to the free-free emission on the Galactic plane in the frequency range 20-40GHz with a width in latitude similar to that of the thermal dust; it comprises 45+/-1% of the total 28.4GHz emission in the longitude range l=300-0-60deg. The free-free component is the narrowest, reflecting the fact that it is produced by current star-formation as traced by the narrow distribution of OB stars. It is the dominant emission on the plane between 60 and 100GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM ~100pc...(abridged)

preprint2015arXiv

Planck intermediate results. XXVIII. Interstellar gas and dust in the Chamaeleon clouds as seen by Fermi LAT and Planck

Shortened abstract: Observations of the nearby Chamaeleon clouds in gamma rays with the Fermi Large Area Telescope and in thermal dust emission with Planck and IRAS have been used with the HI and CO radio data to (i) map the gas column densities in the different phases and at the dark neutral medium (DNM) transition between the HI-bright and CO-bright media; (ii) constrain the CO-to-$H_2$ conversion factor, $X_{CO}$; (iii) probe the dust properties per gas nucleon in each gas phase and spatially across the clouds. We have separated clouds in velocity in HI and CO emission and modelled the 0.4-100 GeV intensity, the dust optical depth at 353 GHz, the thermal radiance of the large grains, and an estimate of the dust extinction empirically corrected for the starlight intensity, $A_{VQ}$. The gamma-ray emissivity spectra confirm that the GeV-TeV cosmic rays uniformly permeate all gas phases up to the CO cores. The dust and cosmic rays reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the HI-DNM-CO transitions and CO-dark $H_2$ fractions for separate clouds. The corrected extinction provides the best fit to the total gas traced by the gamma rays, but we find evidence for a rise in $A_{VQ}/N_H$ and a steep rise in opacity, with increasing $N_H$ and $H_2$ fraction, and with decreasing dust temperature. We observe less variations for the specific power of the grains, except for a decline by half in the CO cores. This combined information suggests grain evolution. The gamma rays and dust radiance yield consistent $X_{CO}$ estimates near $0.7\times10^{20}$ cm$^{-2}$ (K km/s)$^{-1}$. The other dust tracers yield biased values because of the grain opacity rise in the CO clouds. These results also confirm a factor of 2 difference between $X_{CO}$ estimates at pc and kpc scales.

preprint2015arXiv

Planck intermediate results. XXXIII. Signature of the magnetic field geometry of interstellar filaments in dust polarization maps

Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial information to separate Stokes I, Q, and U of the filaments from those of their backgrounds, an essential step to measure the intrinsic polarization fraction (p) and angle (psi) of each emission component. We find that the polarization angles in the three filaments (psi_fil) are coherent along their lengths and not the same as in their backgrounds (psi_bg). The differences between psi_fil and psi_bg are 12deg, 6deg, and 54deg for Musca, B211, and L1506, respectively. These differences for Musca and L1506 are larger than the dispersions of psi, both along the filaments and in their backgrounds. The observed changes of psi are direct evidence for variations of the orientation of the plane of the sky (POS) projection of the B-field. As in previous studies, we find a decrease of several percent of p with NH. We show that the drop in p cannot be explained by random fluctuations of the orientation of B within the filaments because they are too small (sigma_psi<10deg). We recognize the degeneracy between dust alignment efficiency and the structure of B in causing variations in p, but we argue that the decrease of p from the backgrounds to the filaments results in part from depolarization associated with the 3D structure of B: both its orientation in the POS and with respect to the POS. We do not resolve the inner structure of the filaments, but at the smallest scales accessible with Planck (~0.2pc), the observed changes of psi and p hold information on the B-field structure within filaments.

preprint2015arXiv

Planck intermediate results. XXXIV. The magnetic field structure in the Rosette Nebula

Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding HII region on the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line-of-sight component of the magnetic field of about 3microG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 12cm-3. The dust shell that surrounds the Rosette HII region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky orientation of the magnetic field in the Rosette's parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for a given orientation of the field. This comparison leads to an upper limit of about 45degr on the angle between the line of sight and the magnetic field in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4%. This field direction can reproduce the RM values detected in the ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5--9microG.

preprint2015arXiv

Planck intermediate results. XXXVIII. E- and B-modes of dust polarization from the magnetized filamentary structure of the interstellar medium

The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and B-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to 2° (corresponding to 3.5 pc in length for a typical distance of 100 pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes I, Q, U, E, and B, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be 11 %. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the E and B asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck 353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.

preprint2015arXiv

Rotation periods of late-type dwarf stars from time-series high-resolution spectroscopy of chromospheric indicators

We determine rotation periods of a sample of 48 late F-type to mid-M dwarf stars using time-series high-resolution spectroscopy of the Ca II H&K and H-alpha chromospheric activity indicators. We find good agreement between the rotation periods obtained from each of these two indicators. An empirical relationship between the level of chromospheric emission measured by log (R'HK) and the spectroscopic rotation periods is reported. This relation is largely independent of the spectral type and the metallicity of the stars and can be used to make a reliable prediction of rotation periods for late K to mid-M dwarfs with low levels of activity. For some stars in the sample, the measured spectroscopic rotation periods coincide, or are very close, to the orbital periods of postulated planets. In such cases, further studies are needed to clarify whether the associated periodic radial velocity signals reveal the existence of planets or are due to magnetic activity.

preprint2015arXiv

Stellar parameters of early M dwarfs from ratios of spectral features at optical wavelengths

(Abridged) Low-mass stars have been recognised as promising targets in the search for rocky, small planets with the potential of supporting life. Doppler search programmes using high-resolution spectrographs like HARPS or HARPS-N are providing huge quantities of optical spectra of M dwarfs. We aim to calibrate empirical relationships to determine stellar parameters for early M dwarfs (spectral types M0-M4.5) using the same spectra that are used for the radial velocity determinations. Our methodology consists in the use of ratios of pseudo equivalent widths of spectral features as a temperature diagnostic. Stars with effective temperatures obtained from interferometric estimates of their radii are used as calibrators. Empirical calibrations for the spectral type are also provided. Combinations of features and ratios of features are used to derive calibrations for the stellar metallicity. Our methods are then applied to a large sample of M dwarfs that are being observed in the framework of the HARPS search for extrasolar planets.The derived temperatures and metallicities are used together with photometric estimates of mass, radius, and surface gravity to calibrate empirical relationships for these parameters. A total of 112 temperature sensitive ratios have been calibrated over the range 3100-3950 K, providing Teff values with typical uncertainties of the order of 70 K. Eighty-two ratios of pseudo equivalent widths of features were calibrated to derive spectral types. Regarding stellar metallicity, 696 combinations of pseudo equivalent widths of individual features and temperature-sensitive ratios have been calibrated, over the metallicity range from -0.54 to +0.24 dex, with estimated uncertainties in the range of 0.07-0.10 dex. We provide our own empirical calibrations for stellar mass, radius, and surface gravity.

preprint2015arXiv

The QUIJOTE experiment: project overview and first results

QUIJOTE (Q-U-I JOint TEnerife) is a new polarimeter aimed to characterize the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range 10-40 GHz. The multi-frequency (10-20~GHz) instrument, mounted on the first QUIJOTE telescope, saw first light on November 2012 from the Teide Observatory (2400~m a.s.l). During 2014 the second telescope has been installed at this observatory. A second instrument at 30~GHz will be ready for commissioning at this telescope during summer 2015, and a third additional instrument at 40~GHz is now being developed. These instruments will have nominal sensitivities to detect the B-mode polarization due to the primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

preprint2015arXiv

VLT X-shooter spectroscopy of the nearest brown dwarf binary

The aim of the project is to characterise both components of the nearest brown dwarf sytem to the Sun, WISE J104915.57-531906.1 (=Luhman16AB) at optical and near-infrared wavelengths. We obtained high signal-to-noise intermediate-resolution (R~6000-11000) optical (600-1000 nm) and near-infrared (1000-2480nm) spectra of each component of Luhman16AB, the closest brown dwarf binary to the Sun, with the X-Shooter instrument on the Very Large Telescope. We classify the primary and secondary of the Luhman16 system as L6-L7.5 and T0+/-1, respectively, in agreement with previous measurements published in the literature. We present measurements of the lithium pseudo-equivalent widths, which appears of similar strength on both components (8.2+/-1.0 Angstroms and 8.4+/-1.5 Angstroms for the L and T components, respectively). The presence of lithium (Lithium 7) in both components imply masses below 0.06 Msun while comparison with models suggests lower limits of 0.04 Msun. The detection of lithium in the T component is the first of its kind. Similarly, we assess the strength of other alkali lines (e.g. pseudo-equivalent widths of 6-7 Angstroms for RbI and 4-7 Angstroms for CsI) present in the optical and near-infrared regions and compare with estimates for L and T dwarfs. We also derive effective temperatures and luminosities of each component of the binary: -4.66+/-0.08 dex and 1305(+180)(-135) for the L dwarf and -4.68+/-0.13 dex and 1320(+185)(-135) for the T dwarf, respectively. Using our radial velocity determinations, the binary does not appear to belong to any of the well-known moving group. Our preliminary theoretical analysis of the optical and J-band spectra indicates that the L- and T-type spectra can be reproduced with a single temperature and gravity but different relative chemical abundances which impact strongly the spectral energy distribution of L/T transition objects.

preprint2014arXiv

Binary frequency of planet-host stars at wide separations: A new brown dwarf companion to a planet-host star

The aim of the project is to improve our knowledge on the multiplicity of planet-host stars at wide physical separations. We cross-matched approximately 6200 square degree area of the Southern sky imaged by the Visible Infrared Survey Telescope for Astronomy (VISTA) Hemisphere Survey (VHS) with the Two Micron All Sky Survey (2MASS) to look for wide common proper motion companions to known planet-host stars. We complemented our astrometric search with photometric criteria. We confirmed spectroscopically the co-moving nature of seven sources out of 16 companion candidates and discarded eight, while the remaining one stays as a candidate. Among these new wide companions to planet-host stars, we discovered a T4.5 dwarf companion at 6.3 arcmin (~9000 au) from HIP70849, a K7V star which hosts a 9 Jupiter mass planet with an eccentric orbit. We also report two new stellar M dwarf companions to one G and one metal-rich K star. We infer stellar and substellar binary frequencies for our complete sample of 37 targets of 5.4+/-3.8% and 2.7+/-2.7% (1 sigma confidence level), respectively, for projected physical separations larger than ~60-160 au assuming the range of distances of planet-host stars (24-75 pc). These values are comparable to the frequencies of non planet-host stars. We find that the period-eccentricity trend holds with a lack of multiple systems with planets at large eccentricities (e > 0.2) for periods less than 40 days. However, the lack of planets more massive than 2.5 Jupiter masses and short periods (<40 days) orbiting single stars is not so obvious due to recent discoveries by ground-based transit surveys and space missions.

preprint2014arXiv

Chemical abundances of the secondary star in the Neutron star X-ray binary Cygnus X-2

We present UES@WHT high-resolution spectra of the low-mass X-ray binary (LMXB) Cygnus X-2. We have derived the stellar parameters of the secondary star using $χ^{2}$ minimisation procedure, and taking into account any possible veiling from the accretion disk. We determine a metallicity higher than solar ([Fe/H]$=0.27\pm0.19$), as seen also in the neutron star X-ray binary Centaurus X-4. The high-quality of the secondary's spectrum allow us to determine the chemical abundances of O, Mg, Si, Ca, S, Ti, Fe and Ni. We found that some $α$-elements (Mg, Si, S, Ti) are enhanced, consistent with a scenario of contamination of the secondary star during the supernova event. Surprisingly oxygen appears to be under-abundant, whereas enhanced abundances of Fe and Ni are measured. Assuming that these abundances come from matter that has been processed in the SN and then captured by the secondary star, we explore different SN explosion scenarios with diverse geometries. A non-spherically symmetric SN explosion, with a low mass cut, seems to reproduce better the observed abundance pattern of the secondary star compared to the spherical case.

preprint2014arXiv

ESPRESSO: The next European exoplanet hunter

The acronym ESPRESSO stems for Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations; this instrument will be the next VLT high resolution spectrograph. The spectrograph will be installed at the Combined-Coudé Laboratory of the VLT and linked to the four 8.2 m Unit Telescopes (UT) through four optical Coudé trains. ESPRESSO will combine efficiency and extreme spectroscopic precision. ESPRESSO is foreseen to achieve a gain of two magnitudes with respect to its predecessor HARPS, and to improve the instrumental radial-velocity precision to reach the 10 cm/s level. It can be operated either with a single UT or with up to four UTs, enabling an additional gain in the latter mode. The incoherent combination of four telescopes and the extreme precision requirements called for many innovative design solutions while ensuring the technical heritage of the successful HARPS experience. ESPRESSO will allow to explore new frontiers in most domains of astrophysics that require precision and sensitivity. The main scientific drivers are the search and characterization of rocky exoplanets in the habitable zone of quiet, nearby G to M-dwarfs and the analysis of the variability of fundamental physical constants. The project passed the final design review in May 2013 and entered the manufacturing phase. ESPRESSO will be installed at the Paranal Observatory in 2016 and its operation is planned to start by the end of the same year.

preprint2014arXiv

Exoplanet Science with the European Extremely Large Telescope. The Case for Visible and Near-IR Spectroscopy at High Resolution

Exoplanet science is booming. In 20 years our knowledge has expanded considerably, from the first discovery of a Hot Jupiter, to the detection of a large population of Neptunes and super-Earths, to the first steps toward the characterization of exoplanet atmospheres. Between today and 2025, the field will evolve at an even faster pace with the advent of several space-based transit search missions, ground-based spectrographs, high-contrast imaging facilities, and the James Webb Space Telescope. Especially the ESA M-class PLATO mission will be a game changer in the field. From 2024 onwards, PLATO will find transiting terrestrial planets orbiting within the habitable zones of nearby, bright stars. These objects will require the power of Extremely Large Telescopes (ELTs) to be characterized further. The technique of ground-based high-resolution spectroscopy is establishing itself as a crucial pathway to measure chemical composition, atmospheric structure and atmospheric circulation in transiting exoplanets. A high-resolution spectrograph covering the visible and near-IR domains, mounted on the European ELT, will be able to detect molecules such as water vapour, carbon dioxide and oxygen in the atmospheres of habitable planets under favourable circumstances. E-ELT HiRES is the perfect ground-based match to the PLATO space mission and represents a unique opportunity for Europe to lead the world into the era of exploration of exoplanets with habitable conditions. HiRES will also be extremely complementary to other E-ELT planned instruments specialising in different kinds of planets, such as METIS and EPICS.

preprint2014arXiv

Planck 2013 results. I. Overview of products and scientific results

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

preprint2014arXiv

Planck 2013 results. II. Low Frequency Instrument data processing

We describe the data processing pipeline of the Planck Low Frequency Instrument (LFI) data processing centre (DPC) to create and characterize full-sky maps based on the first 15.5 months of operations at 30, 44 and 70 GHz. In particular, we discuss the various steps involved in reducing the data, starting from telemetry packets through to the production of cleaned, calibrated timelines and calibrated frequency maps. Data are continuously calibrated using the modulation induced on the mean temperature of the cosmic microwave background radiation by the proper motion of the spacecraft. Sky signals other than the dipole are removed by an iterative procedure based on simultaneous fitting of calibration parameters and sky maps. Noise properties are estimated from time-ordered data after the sky signal has been removed, using a generalized least square map-making algorithm. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices, required to compute statistical uncertainties on LFI and Planck products, are also produced. Main beams are estimated down to the -20 dB level using Jupiter transits, which are also used for the geometrical calibration of the focal plane.

preprint2014arXiv

Planck 2013 results. III. LFI systematic uncertainties

We present the current estimate of instrumental and systematic effect uncertainties for the Planck-Low Frequency Instrument relevant to the first release of the Planck cosmological results. We give an overview of the main effects and of the tools and methods applied to assess residuals in maps and power spectra. We also present an overall budget of known systematic effect uncertainties, which are dominated sidelobe straylight pick-up and imperfect calibration. However, even these two effects are at least two orders of magnitude weaker than the cosmic microwave background (CMB) fluctuations as measured in terms of the angular temperature power spectrum. A residual signal above the noise level is present in the multipole range $\ell<20$, most notably at 30 GHz, and is likely caused by residual Galactic straylight contamination. Current analysis aims to further reduce the level of spurious signals in the data and to improve the systematic effects modelling, in particular with respect to straylight and calibration uncertainties.

preprint2014arXiv

Planck 2013 results. IV. Low Frequency Instrument beams and window functions

This paper presents the characterization of the in-flight beams, the beam window functions and the associated uncertainties for the Planck Low Frequency Instrument (LFI). Knowledge of the beam profiles is necessary for determining the transfer function to go from the observed to the actual sky anisotropy power spectrum. The main beam distortions affect the beam window function, complicating the reconstruction of the anisotropy power spectrum at high multipoles, whereas the sidelobes affect the low and intermediate multipoles. The in-flight assessment of the LFI main beams relies on the measurements performed during Jupiter observations. By stacking the data from multiple Jupiter transits, the main beam profiles are measured down to -20 dB at 30 and 44 GHz, and down to -25 dB at 70 GHz. The main beam solid angles are determined to better than 0.2% at each LFI frequency band. The Planck pre-launch optical model is conveniently tuned to characterize the main beams independently of any noise effects. This approach provides an optical model whose beams fully reproduce the measurements in the main beam region, but also allows a description of the beams at power levels lower than can be achieved by the Jupiter measurements themselves. The agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band. The simulated beams are used for the computation of the window functions for the effective beams. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer bandshapes. The total uncertainties in the effective beam window functions are: 2% and 1.2% at 30 and 44 GHz, respectively (at $\ell \approx 600$), and 0.7% at 70 GHz (at $\ell \approx 1000$).

preprint2014arXiv

Planck 2013 results. V. LFI calibration

We discuss the methods employed to photometrically calibrate the data acquired by the Low Frequency Instrument on Planck. Our calibration is based on a combination of the Orbital Dipole plus the Solar Dipole, caused respectively by the motion of the Planck spacecraft with respect to the Sun and by motion of the Solar System with respect to the CMB rest frame. The latter provides a signal of a few mK with the same spectrum as the CMB anisotropies and is visible throughout the mission. In this data release we rely on the characterization of the Solar Dipole as measured by WMAP. We also present preliminary results (at 44GHz only) on the study of the Orbital Dipole, which agree with the WMAP value of the Solar System speed within our uncertainties. We compute the calibration constant for each radiometer roughly once per hour, in order to keep track of changes in the detectors' gain. Since non-idealities in the optical response of the beams proved to be important, we implemented a fast convolution algorithm which considers the full beam response in estimating the signal generated by the dipole. Moreover, in order to further reduce the impact of residual systematics due to sidelobes, we estimated time variations in the calibration constant of the 30GHz radiometers (the ones with the largest sidelobes) using the signal of a reference load. We have estimated the calibration accuracy in two ways: we have run a set of simulations to assess the impact of statistical errors and systematic effects in the instrument and in the calibration procedure, and we have performed a number of consistency checks on the data and on the brightness temperature of Jupiter. Calibration errors for this data release are expected to be about 0.6% at 44 and 70 GHz, and 0.8% at 30 GHz. (Abriged.)

preprint2014arXiv

Planck 2013 results. VI. High Frequency Instrument data processing

We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12 and 39 microKelvin in HFI four lowest frequency channel (100--353 GHz) and 13 and 14 kJy/sr for the 545 and 857 GHz channels. Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5% and the 353 GHz relative calibration is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50 < l <2500), are intercalibrated at better than 0.2 %.

preprint2014arXiv

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

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

preprint2014arXiv

Planck 2013 results. XIV. Zodiacal emission

The Planck satellite provides a set of all-sky maps at nine frequencies from 30 GHz to 857 GHz. Planets, minor bodies, and diffuse interplanetary dust emission (IPD) are all observed. The IPD can be separated from Galactic and other emissions because Planck views a given point on the celestial sphere multiple times, through different columns of IPD. We use the Planck data to investigate the behaviour of zodiacal emission over the whole sky at sub-millimetre and millimetre wavelengths. We fit the Planck data to find the emissivities of the various components of the COBE zodiacal model -- a diffuse cloud, three asteroidal dust bands, a circumsolar ring, and an Earth-trailing feature. The emissivity of the diffuse cloud decreases with increasing wavelength, as expected from earlier analyses. The emissivities of the dust bands, however, decrease less rapidly, indicating that the properties of the grains in the bands are different from those in the diffuse cloud. We fit the small amount of Galactic emission seen through the telescope's far sidelobes, and place limits on possible contamination of the CMB results from both zodiacal and far-sidelobe emission. When necessary, the results are used in the Planck pipeline to make maps with zodiacal emission and far sidelobes removed. We show that the zodiacal correction to the CMB maps is small compared to the Planck CMB temperature power spectrum and give a list of flux densities for small Solar System bodies.

preprint2014arXiv

Planck 2013 results. XVI. Cosmological parameters

We present the first results based on Planck measurements of the CMB temperature and lensing-potential power spectra. The Planck spectra at high multipoles are extremely well described by the standard spatially-flat six-parameter LCDM cosmology. In this model Planck data determine the cosmological parameters to high precision. We find a low value of the Hubble constant, H0=67.3+/-1.2 km/s/Mpc and a high value of the matter density parameter, Omega_m=0.315+/-0.017 (+/-1 sigma errors) in excellent agreement with constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent-level precision using Planck CMB data alone. We present results from an analysis of extensions to the standard cosmology, using astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured significantly over standard LCDM. The deviation of the scalar spectral index from unity is insensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find a 95% upper limit of r<0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles. Using BAO and CMB data, we find N_eff=3.30+/-0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the summed neutrino mass. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of N_eff=3.046. We find no evidence for dynamical dark energy. Despite the success of the standard LCDM model, this cosmology does not provide a good fit to the CMB power spectrum at low multipoles, as noted previously by the WMAP team. While not of decisive significance, this is an anomaly in an otherwise self-consistent analysis of the Planck temperature data.

preprint2014arXiv

Planck 2013 results. XVII. Gravitational lensing by large-scale structure

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

preprint2014arXiv

Planck 2013 results. XX. Cosmology from Sunyaev-Zeldovich cluster counts

We present constraints on cosmological parameters using number counts as a function of redshift for a sub-sample of 189 galaxy clusters from the Planck SZ (PSZ) catalogue. The PSZ is selected through the signature of the Sunyaev--Zeldovich (SZ) effect, and the sub-sample used here has a signal-to-noise threshold of seven, with each object confirmed as a cluster and all but one with a redshift estimate. We discuss the completeness of the sample and our construction of a likelihood analysis. Using a relation between mass $M$ and SZ signal $Y$ calibrated to X-ray measurements, we derive constraints on the power spectrum amplitude $σ_8$ and matter density parameter $Ω_{\mathrm{m}}$ in a flat $Λ$CDM model. We test the robustness of our estimates and find that possible biases in the $Y$--$M$ relation and the halo mass function are larger than the statistical uncertainties from the cluster sample. Assuming the X-ray determined mass to be biased low relative to the true mass by between zero and 30%, motivated by comparison of the observed mass scaling relations to those from a set of numerical simulations, we find that $σ_8=0.75\pm 0.03$, $Ω_{\mathrm{m}}=0.29\pm 0.02$, and $σ_8(Ω_{\mathrm{m}}/0.27)^{0.3} = 0.764 \pm 0.025$. The value of $σ_8$ is degenerate with the mass bias; if the latter is fixed to a value of 20% we find $σ_8(Ω_{\mathrm{m}}/0.27)^{0.3}=0.78\pm 0.01$ and a tighter one-dimensional range $σ_8=0.77\pm 0.02$. We find that the larger values of $σ_8$ and $Ω_{\mathrm{m}}$ preferred by Planck's measurements of the primary CMB anisotropies can be accommodated by a mass bias of about 40%. Alternatively, consistency with the primary CMB constraints can be achieved by inclusion of processes that suppress power on small scales relative to the $Λ$CDM model, such as a component of massive neutrinos (abridged).

preprint2014arXiv

Planck 2013 results. XXI. All-sky Compton parameter power spectrum and high-order statistics

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim θ\lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $σ_8(Ω_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $σ_{8}$.

preprint2014arXiv

Planck 2013 results. XXIII. Isotropy and statistics of the CMB

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

preprint2014arXiv

Planck 2013 results. XXIX. Planck catalogue of Sunyaev-Zeldovich sources

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

preprint2014arXiv

Planck intermediate results. XIII. Constraints on peculiar velocities

Using \Planck\ data combined with the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we address the study of peculiar motions by searching for evidence of the kinetic Sunyaev-Zeldovich effect (kSZ). By implementing various filters designed to extract the kSZ generated at the positions of the clusters, we obtain consistent constraints on the radial peculiar velocity average, root mean square (rms), and local bulk flow amplitude at different depths. For the whole cluster sample of average redshift 0.18, the measured average radial peculiar velocity with respect to the cosmic microwave background (CMB) radiation at that redshift, i.e., the kSZ monopole, amounts to $72 \pm 60$ km s$^{-1}$. This constitutes less than 1% of the relative Hubble velocity of the cluster sample with respect to our local CMB frame. While the linear $Λ$CDM prediction for the typical cluster radial velocity rms at $z=0.15$ is close to 230km s$^{-1}$, the upper limit imposed by \Planck\ data on the cluster subsample corresponds to 800 km s$^{-1}$ at 95% confidence level, i.e., about three times higher. \Planck\ data also set strong constraints on the local bulk flow in volumes centred on the Local Group. There is no detection of bulk flow as measured in any comoving sphere extending to the maximum redshift covered by the cluster sample. A blind search for bulk flows in this sample has an upper limit of 254 km s$^{-1}$ (95% confidence level) dominated by CMB confusion and instrumental noise, indicating that the Universe is largely homogeneous on Gpc scales. In this context, in conjunction with supernova observations, \Planck\ is able to rule out a large class of inhomogeneous void models as alternatives to dark energy or modified gravity. The \Planck\ constraints on peculiar velocities and bulk flows are thus consistent with the $Λ$CDM scenario.

preprint2014arXiv

Planck intermediate results. XIV. Dust emission at millimetre wavelengths in the Galactic plane

We use Planck HFI data combined with ancillary radio data to study the emissivity index of the interstellar dust emission in the frequency range 10 - 353 GHz, or 3 - 0.8 mm, in the Galactic plane. We analyse the region l=20 degr - 44 degr and |b| \leq 4 degr where the free-free emission can be estimated from radio recombination line data. We fit the spectra at each sky pixel with a modified blackbody model and two spectral indices, beta_mm and beta_FIR, below and above 353 GHz respectively. We find that beta_mm is smaller than beta_FIR and we detect a correlation between this low frequency power-law index and the dust optical depth at 353 GHz, tau_353. The opacity spectral index beta_mm increases from about 1.54 in the more diffuse regions of the Galactic disk, |b| = 3 degr - 4 degr and tau_353 ~ 5 x 10^{-5}, to about 1.66 in the densest regions with an optical depth of more than one order of magnitude higher. We associate this correlation with an evolution of the dust emissivity related to the fraction of molecular gas along the line of sight. This translates into beta_mm ~ 1.54 for a medium that is mostly atomic and beta_mm ~ 1.66 when the medium is dominated by molecular gas. We find that both the Two-Level System model and the emission by ferromagnetic particles can explain the results. The results improve our understanding of the physics of interstellar dust and lead towards a complete model of the dust spectrum of the Milky Way from far-infrared to millimetre wavelengths.

preprint2014arXiv

Planck intermediate results. XIX. An overview of the polarized thermal emission from Galactic dust

This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%), in particular in some of the intermediate dust column density (AV < 1mag) regions. There is a systematic decrease in the dust polarization fraction with increasing dust column density, and we interpret the features of this correlation in light of both radiative grain alignment predictions and fluctuations in the magnetic field orientation. We also characterize the spatial structure of the polarization angle using the angle dispersion function and find that, in nearby fields at intermediate latitudes, the polarization angle is ordered over extended areas that are separated by filamentary structures, which appear as interfaces where the magnetic field sky projection rotates abruptly without apparent variations in the dust column density. The polarization fraction is found to be anti-correlated with the dispersion of the polarization angle, implying that the variations are likely due to fluctuations in the 3D magnetic field orientation along the line of sight sampling the diffuse interstellar medium.We also compare the dust emission with the polarized synchrotron emission measured with the Planck LFI, with low-frequency radio data, and with Faraday rotation measurements of extragalactic sources. The two polarized components are globally similar in structure along the plane and notably in the Fan and North Polar Spur regions. A detailed comparison of these three tracers shows, however, that dust and cosmic rays generally sample different parts of the line of sight and confirms that much of the variation observed in the Planck data is due to the 3D structure of the magnetic field.

preprint2014arXiv

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

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

preprint2014arXiv

Planck intermediate results. XXIV. Constraints on variation of fundamental constants

Any variation of the fundamental physical constants, and more particularly of the fine structure constant, $α$, or of the mass of the electron, $m_e$, would affect the recombination history of the Universe and cause an imprint on the cosmic microwave background angular power spectra. We show that the Planck data allow one to improve the constraint on the time variation of the fine structure constant at redshift $z\sim 10^3$ by about a factor of 5 compared to WMAP data, as well as to break the degeneracy with the Hubble constant, $H_0$. In addition to $α$, we can set a constraint on the variation of the mass of the electron, $m_{\rm e}$, and on the simultaneous variation of the two constants. We examine in detail the degeneracies between fundamental constants and the cosmological parameters, in order to compare the limits obtained from Planck and WMAP and to determine the constraining power gained by including other cosmological probes. We conclude that independent time variations of the fine structure constant and of the mass of the electron are constrained by Planck to ${Δα}/α= (3.6\pm 3.7)\times10^{-3}$ and ${Δm_{\rm e}}/{m_{\rm e}}= (4 \pm 11)\times10^{-3}$ at the 68% confidence level. We also investigate the possibility of a spatial variation of the fine structure constant. The relative amplitude of a dipolar spatial variation of $α$ (corresponding to a gradient across our Hubble volume) is constrained to be $δα/α=(-2.4\pm 3.7)\times 10^{-2}$.

preprint2014arXiv

Planck intermediate results. XXX. The angular power spectrum of polarized dust emission at intermediate and high Galactic latitudes

The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $α^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $β_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2π)$ of $1.32\times10^{-2}μ$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.

preprint2014arXiv

Spectroscopic follow-up of L- and T-type proper-motion member candidates in the Pleiades

(Abridged) We report on the near-infrared low-resolution spectroscopy and red optical (Z-band) photometry of seven proper-motion, very low-mass substellar member candidates of the Pleiades cluster with magnitudes in the interval J=17.5-20.8 and K=16.1-18.5 mag. Spectra were acquired for six objects with the LIRIS and NIRSPEC instruments mounted on the 4.2-m WHT and the 10-m Keck II telescopes. Z-band images of two of the faintest candidates were collected with ACAM/WHT. The new data confirm the low temperatures of all seven Pleiades candidates. From the imaging observations, we find extremely red Z-J and Z-K colors that suggest that the faintest target, Calar Pleiades 25, has a Galactic rather than extragalactic nature. We tentatively classify the spectroscopic targets from early-L to ~T0 and suggest that the L/T transition, which accounts for the onset of methane absorption at 2.1 micron, may take place at J,K ~ 20.3, 17.8 mag in the Pleiades (absolute values of M_J ~ 14.7 and M_K ~ 12.2 mag). We find evidence of likely low-gravity atmospheres based on the presence of triangular-shape H-band fluxes and the high flux ratio K/H (compatible with red H-K colors) of Calar Pleiades 20, 21, and 22, which is a feature also seen in field low-gravity dwarfs. Weak KI absorption lines at around 1.25 micron are probably seen in two targets. These observations add support to the cluster membership of all seven objects in the Pleiades. The trend delineated by the spectroscopic sequence of Pleiades late-M and L dwarfs resembles that of the field. With masses estimated at 0.012-0.015 Msol (solar metallicity and 120 Myr), Calar Pleiades 20 (L6 +/- 1), 21 (L7 +/- 1), and 22 (L/T) may become the coolest and least massive Pleiades members that are corroborated with photometry, astrometry, and spectroscopy. Calar Pleiades 25 (<0.012 Msol) is a firm free-floating planetary-mass candidate in the Pleiades.

preprint2014arXiv

The QUIJOTE CMB Experiment: status and first results with the multi-frequency instrument

The QUIJOTE (Q-U-I JOint Tenerife) CMB Experiment is designed to observe the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range of 10-40 GHz. The first of the two QUIJOTE telescopes and the multi-frequency (10-20 GHz) instrument have been in operation since November 2012. In 2014 a second telescope and a new instrument at 30GHz will be ready for commissioning, and an additional instrument at 40 GHz is in its final design stages. After three years of effective observations, the data obtained by these telescopes and instruments will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r = 0.05. At the moment, we have completed half of the wide Galactic survey with the multi-frequency instrument covering 18 000 square degrees of the Northern hemisphere. When we finish this survey in early 2014, we shall have reached approximately 14μK per one degree beam at 11, 13, 17 and 19 GHz, in both Q and U.

preprint2014arXiv

Trigonometric parallaxes of young field L dwarfs

(Abridged) We aim to determine the trigonometric parallaxes and proper motions of a sample of ten field L0-L5 dwarfs with spectroscopic evidence for low-gravity atmospheres. We obtained J and Ks imaging data using 2-4-m class telescopes with a typical cadence of one image per month between 2010 January and 2012 December. We also obtained low resolution optical spectra (R~300, 500-1100 nm) using the 10-m GTCs to assess the presence of lithium absorption in four targets and confirm their young age. Trigonometric parallaxes and proper motions were derived to typical accuracies of 1 mas and +/-10 mas/yr. All ten L dwarfs have large motions, and are located at distances between 9 and 47 pc. They lie above and on the sequence of field dwarfs in the absolute J and K_s magnitude versus spectral type and luminosity versus Teff diagrams, implying ages similar to or smaller than those typical of the field. The detection of atomic lithium in the atmosphere of 2MASS J00452143+1634446 is reported for the first time. Three dwarfs have locations in the HR diagram indicative of old ages and high masses consistent with the observed lithium depletion previously published. We did not find evidence for the presence of astrometric companions with minimum detectable masses typically >=25 Mjup and face-on, circular orbits with periods between 60-90 d and 3 yr around eight targets. The astrometric and spectroscopic data indicate that about 60-70% of the field L-type dwarfs in our sample with evidence for low-gravity atmospheres are indeed young-to-intermediate-age brown dwarfs of the solar neighborhood with expected ages and masses in the intervals 10-500 Myr and 11-45 Mjup. The peaked-shape of the H-band spectra of L dwarfs, a signpost of youth, appears to be present up to ages of 120-500 Myr and intermediate-to-high gravities.

preprint2013arXiv

A characterization of the diffuse Galactic emissions at large angular scales using the Tenerife data

The Anomalous microwave emission (AME) has been proved to be an important component of the Galactic diffuse emission in the range from 20 to 60 GHz. To discriminate between different models of AME low frequency microwave data from 10 to 20 GHz are needed. We present here a re-analysis of published and un-published Tenerife data from 10 to 33 GHz at large angular scales (from 5 to 15 degrees). We cross-correlate the Tenerife data to templates of the main Galactic diffuse emissions: synchrotron, free-free and thermal dust. We find evidence of dust correlated emission in the Tenerife data that could be explained as spinning dust grain emission.

preprint2013arXiv

A Community Science Case for E-ELT HIRES

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

preprint2013arXiv

A Frequency Comb calibrated Solar Atlas

The solar spectrum is a primary reference for the study of physical processes in stars and their variation during activity cycles. In Nov 2010 an experiment with a prototype of a Laser Frequency Comb (LFC) calibration system was performed with the HARPS spectrograph of the 3.6m ESO telescope at La Silla during which high signal-to-noise spectra of the Moon were obtained. We exploit those Echelle spectra to study the optical integrated solar spectrum . The DAOSPEC program is used to measure solar line positions through gaussian fitting in an automatic way. We first apply the LFC solar spectrum to characterize the CCDs of the HARPS spectrograph. The comparison of the LFC and Th-Ar calibrated spectra reveals S-type distortions on each order along the whole spectral range with an amplitude of +/-40 m/s. This confirms the pattern found by Wilken et al. (2010) on a single order and extends the detection of the distortions to the whole analyzed region revealing that the precise shape varies with wavelength. A new data reduction is implemented to deal with CCD pixel inequalities to obtain a wavelength corrected solar spectrum. By using this spectrum we provide a new LFC calibrated solar atlas with 400 line positions in the range of 476-530, and 175 lines in the 534-585 nm range. The new LFC atlas improves the accuracy of individual lines by a significant factor reaching a mean value of about 10 m/s. The LFC--based solar line wavelengths are essentially free of major instrumental effects and provide a reference for absolute solar line positions. We suggest that future LFC observations could be used to trace small radial velocity changes of the whole solar photospheric spectrum in connection with the solar cycle and for direct comparison with the predicted line positions of 3D radiative hydrodynamical models of the solar photosphere.

preprint2013arXiv

Fast orbital decays of black hole X-ray binaries: XTE J1118+480 and A0620--00

We present new 10.4m-GTC/OSIRIS spectroscopic observations of the black hole X-ray binary \mbox{XTE J1118+480} that confirm the orbital period decay at $\dot P=-1.90\pm0.57$ ms yr$^{-1}$. This corresponds to a period change of $-0.88\pm0.27$ $μ$s per orbital cycle. We have also collected observations of the black hole X-ray binary A0620--00 to derive an orbital period derivative of $\dot P=-0.60\pm0.08$ ms yr$^{-1}$ ($-0.53\pm0.07$ $μ$s/cycle). Angular momentum losses due to gravitational radiation are unable to explain these large orbital decays in these two short-period black hole binaries. The orbital period decay measured in A0620--00 is very marginally consistent with the predictions of conventional models including magnetic braking, although significant mass loss ($\dot M_{\rm BH} / \dot M_2 \le 20\%$) from the system is required. The fast spiral-in of the star in XTE J1118+480, however, does not fit any standard model and may be driven by magnetic braking under extremely high magnetic fields and/or may require an unknown process or non-standard theories of gravity. This result may suggest an evolutionary sequence in which the orbital period decay begins to speed up as the orbital period decreases. This scenario may have an impact on the evolution and lifetime of black hole X-ray binaries.

preprint2013arXiv

GTC OSIRIS z-band imaging of Y dwarfs

The aim of the project is to contribute to the characterisation of the spectral energy distribution of the coolest brown dwarfs discovered to date, the Y dwarfs. We obtained z-band far-red imaging for six Y dwarfs and a T9+Y0 binary with the OSIRIS (Optical System for Imaging and low Resolution Integrated Spectroscopy) instrument on the 10.4-m Gran Telescopio de Canarias (GTC). We detect five of the seven known Y dwarfs in the $z$-band, infer their optical-to-infrared colours, and measure their proper motions. We find a higher dispersion in the z-J and z-H colours of Y0 dwarfs than in T dwarfs. This dispersion is found to be correlated with H-w2. The high dispersion in the optical-to-infrared colours of Y dwarfs and the possible turn-over towards bluer colours may be a consequence of the presence of sulfide clouds with different thicknesses, the depletion of alcalines, and/or gravity effects.

preprint2013arXiv

Li depletion in solar analogues with exoplanets: Extending the sample

We want to study the effects of the formation of planets and planetary systems on the atmospheric Li abundance of planet host stars. In this work we present new determinations of lithium abundances for 326 Main Sequence stars with and without planets in the T$_\mathrm{eff}$ range 5600-5900 K. 277 stars come from the HARPS sample, the remaining targets have been observed with a variety of high resolution spectrographs. We confirm significant differences in the Li distribution of solar twins (T$_\mathrm{eff}$ = T$_{\odot} \pm$ 80 K, log g = log g$_{\odot}$ $\pm$ 0.2 and [Fe/H] = [Fe/H]$_{\odot} \pm$ 0.2): the full sample of planet host stars (22) shows Li average values lower than "single" stars with no detected planets (60). If we focus in subsamples with narrower ranges in metallicity and age, we observe indications of a similar result though it is not so clear for some of the studied subsamples. Furthermore, we compare the observed spectra of several couples of stars with very similar parameters which show different Li abundances up to 1.6 dex. Therefore we show that neither age, nor mass nor metallicity of a parent star is the only responsible for enhanced Li depletion in solar analogues. We conclude that another variable must account for that difference and suggest that this could be the presence of planets which causes additional rotationally induced mixing in the external layers of planet host stars. Moreover, we find indications that the amount of depletion of Li in planet host solar-type stars is higher when the planets are more massive than Jupiter.

preprint2013arXiv

Planck 2013 results. XII. Component separation

Planck has produced detailed all-sky observations over nine frequency bands between 30 and 857 GHz. These observations allow robust reconstruction of the primordial cosmic microwave background (CMB) temperature fluctuations over nearly the full sky, as well as new constraints on Galactic foregrounds. This paper describes the component separation framework adopted by Planck. We test four foreground-cleaned CMB maps derived using qualitatively different component separation algorithms. The quality of our reconstructions is evaluated through detailed simulations and internal comparisons, and shown through various tests to be internally consistent and robust for CMB power spectrum and cosmological parameter estimation up to l = 2000. The parameter constraints on LambdaCDM cosmologies derived from these maps are consistent with those presented in the cross-spectrum based Planck likelihood analysis. We choose two of the CMB maps for specific scientific goals. We also present maps and frequency spectra of the Galactic low-frequency, CO, and thermal dust emission. The component maps are found to provide a faithful representation of the sky, as evaluated by simulations. For the low-frequency component, the spectral index varies widely over the sky, ranging from about beta = -4 to -2. Considering both morphology and prior knowledge of the low frequency components, the index map allows us to associate a steep spectral index (beta < -3.2) with strong anomalous microwave emission, corresponding to a spinning dust spectrum peaking below 20 GHz, a flat index of beta > -2.3 with strong free-free emission, and intermediate values with synchrotron emission.

preprint2013arXiv

Planck 2013 results. XIII. Galactic CO emission

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

preprint2013arXiv

Planck 2013 results. XIX. The integrated Sachs-Wolfe effect

Based on CMB maps from the 2013 Planck Mission data release, this paper presents the detection of the ISW effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4 sigma, depending on which method is used. We investigate three separate approaches, which cover essentially all previous studies, as well as breaking new ground. (i) Correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection is made using the lensing-induced bispectrum; the correlation between lensing and the ISW effect has a significance close to 2.5 sigma. (ii) Cross-correlation with tracers of LSS, yielding around 3 sigma significance, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) Aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yields a 4 sigma signal when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate results, ranging from negligible to 2.5 sigma at most, but with a more consistent scale and amplitude, the latter being still slightly above what is expected from numerical simulations within LCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe.

preprint2013arXiv

Planck 2013 results. XV. CMB power spectra and likelihood

We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)

preprint2013arXiv

Planck 2013 results. XVIII. Gravitational lensing-infrared background correlation

The multi-frequency capability of the Planck satellite provides information both on the integrated history of star formation (via the cosmic infrared background, or CIB) and on the distribution of dark matter (via the lensing effect on the cosmic microwave background, or CMB). The conjunction of these two unique probes allows us to measure directly the connection between dark and luminous matter in the high redshift (1 < z <3) Universe. We use a three-point statistic optimized to detect the correlation between these two tracers. Following a thorough discussion of possible contaminants and a suite of consistency tests, using lens reconstructions at 100, 143 and 217 GHz and CIB measurements at 100-857 GHz, we report the first detection of the correlation between the CIB and CMB lensing. The well matched redshift distribution of these two signals leads to a detection significance with a peak value of 42 σat 545 GHz and a correlation as high as 80% across these two tracers. Our full set of multi-frequency measurements (both CIB auto- and CIB-lensing cross-spectra) are consistent with a simple halo-based model, with a characteristic mass scale for the halos hosting CIB sources of log_{10}(M/M_sun) = 10.5 \pm 0.6. Leveraging the frequency dependence of our signal, we isolate the high redshift contribution to the CIB, and constrain the star formation rate (SFR) density at z>1. We measure directly the SFR density with around 2 sigma significance for three redshift bins between z=1 and 7, thus opening a new window into the study of the formation of stars at early times.

preprint2013arXiv

Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity

The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result fNL^local= 2.7+/-5.8, fNL^equil= -42+/-75, and fNL^ortho= -25+\-39 (68% CL statistical). NG is detected in the data; using skew-C_l statistics we find a nonzero bispectrum from residual point sources, and the ISW-lensing bispectrum at a level expected in the LambdaCDM scenario. The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are confirmed by skew-C_l, wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, 3-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, c_s \geq 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction r_D \geq 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the 4-point function in the local model tauNL < 2800 (95% CL). These constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure.

preprint2013arXiv

Planck 2013 results. XXV. Searches for cosmic strings and other topological defects

Planck data have been used to provide stringent new constraints on cosmic strings and other defects. We describe forecasts of the CMB power spectrum induced by cosmic strings, calculating these from network models and simulations using line-of-sight Boltzmann solvers. We have studied Nambu-Goto cosmic strings, as well as field theory strings for which radiative effects are important, thus spanning the range of theoretical uncertainty in strings models. We have added the angular power spectrum from strings to that for a simple adiabatic model, with the extra fraction defined as $f_{10}$ at multipole $\ell=10$. This parameter has been added to the standard six parameter fit using COSMOMC with flat priors. For the Nambu-Goto string model, we have obtained a constraint on the string tension of $Gμ/c^2 < 1.5 x 10^{-7}$ and $f_{10} < 0.015$ at 95% confidence that can be improved to $Gμ/c^2 < 1.3 x 10^{-7}$ and $f_{10} < 0.010$ on inclusion of high-$\ell$ CMB data. For the abelian-Higgs field theory model we find, $Gμ_{AH}/c^2 < 3.2 x 10^{-7}$ and $f_{10} < 0.028$. The marginalized likelihoods for $f_{10}$ and in the $f_{10}$--$Ω_b h^2$ plane are also presented. We have also obtained constraints on $f_{10}$ for models with semi-local strings and global textures for which $Gμ/c^2 < 1.1 x 10^{-6}$. We have made complementarity searches for the specific non-Gaussian signatures of cosmic strings, calibrating with all-sky Planck resolution CMB maps generated from networks of post-recombination strings. We have obtained upper limits on the string tension at 95% confidence of $Gμ/c^2 < 8.8 x 10^{-7}$ using modal bispectrum estimation and $Gμ/c^2 < 7.8 x 10^{-7}$ for real space searches with Minkowski functionals. These are conservative upper bounds because only post-recombination string contributions have been included in the non-Gaussian analysis.

preprint2013arXiv

Planck 2013 results. XXVI. Background geometry and topology of the Universe

Planck CMB temperature maps allow detection of large-scale departures from homogeneity and isotropy. We search for topology with a fundamental domain nearly intersecting the last scattering surface (comoving distance $χ_r$). For most topologies studied the likelihood maximized over orientation shows some preference for multi-connected models just larger than $χ_r$. This effect is also present in simulated realizations of isotropic maps and we interpret it as the alignment of mild anisotropic correlations with chance features in a single realization; such a feature can also exist, in milder form, when the likelihood is marginalized over orientations. Thus marginalized, the limits on the radius $R_i$ of the largest sphere inscribed in a topological domain (at log-likelihood-ratio -5) are: in a flat Universe, $R_i>0.9χ_r$ for the cubic torus (cf. $R_i>0.9χ_r$ at 99% CL for a matched-circles search); $R_i>0.7χ_r$ for the chimney; $R_i>0.5χ_r$ for the slab; in a positively curved Universe, $R_i>1.0χ_r$ for the dodecahedron; $R_i>1.0χ_r$ for the truncated cube; $R_i>0.9χ_r$ for the octahedron. Similar limits apply to alternate topologies. We perform a Bayesian search for an anisotropic Bianchi VII$_h$ geometry. In a non-physical setting where the Bianchi parameters are decoupled from cosmology, Planck data favour a Bianchi component with a Bayes factor of at least 1.5 units of log-evidence: a Bianchi pattern is efficient at accounting for some large-scale anomalies in Planck data. However, the cosmological parameters are in strong disagreement with those found from CMB anisotropy data alone. In the physically motivated setting where the Bianchi parameters are fitted simultaneously with standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models: $(ω/H)_0<8\times10^{-10}$ (95% CL). [Abridged]

preprint2013arXiv

Planck 2013 results. XXVIII. The Planck Catalogue of Compact Sources

The Planck Catalogue of Compact Sources (PCCS) is the catalogue of sources detected in the first 15 months of Planck operations, the "nominal" mission. It consists of nine single-frequency catalogues of compact sources, both Galactic and extragalactic, detected over the entire sky. The PCCS covers the frequency range 30--857\,GHz with higher sensitivity (it is 90% complete at 180 mJy in the best channel) and better angular resolution (from ~33' to ~5') than previous all-sky surveys in this frequency band. By construction its reliability is >80% and more than 65% of the sources have been detected at least in two contiguous Planck channels. In this paper we present the construction and validation of the PCCS, its contents and its statistical characterization.

preprint2013arXiv

Planck 2013 results. XXX. Cosmic infrared background measurements and implications for star formation

We present new measurements of CIB anisotropies using Planck. Combining HFI data with IRAS, the angular auto- and cross frequency power spectrum is measured from 143 to 3000 GHz, and the auto-bispectrum from 217 to 545 GHz. The total areas used to compute the CIB power spectrum and bispectrum are about 2240 and 4400 deg^2, respectively. After careful removal of the contaminants, and a complete study of systematics, the CIB power spectrum and bispectrum are measured with unprecedented signal to noise ratio from angular multipoles ell~150 to 2500, and ell~130 to 1100, respectively. Two approaches are developed for modelling CIB power spectrum anisotropies. The first approach takes advantage of the unique measurements by Planck at large angular scales, and models only the linear part of the power spectrum, with a mean bias of dark matter halos hosting dusty galaxies at a given redshift weighted by their contribution to the emissivities. The second approach is based on a model that associates star-forming galaxies with dark matter halos and their subhalos, using a parametrized relation between the dust-processed infrared luminosity and (sub-)halo mass. The two approaches simultaneously fit all auto- and cross- power spectra very well. We find that the star formation history is well constrained up to z~2. However, at higher redshift, the accuracy of the star formation history measurement is strongly degraded by the uncertainty in the spectral energy distribution of CIB galaxies. We also find that CIB galaxies have warmer temperatures as redshift increases. The CIB bispectrum is steeper than that expected from the power spectrum, although well fitted by a power law; this gives some information about the contribution of massive halos to the CIB bispectrum.

preprint2013arXiv

Planck intermediate results. VII. Statistical properties of infrared and radio extragalactic sources from the Planck Early Release Compact Source Catalogue at frequencies between 100 and 857 GHz

(abridged for arXiv) We make use of the Planck all-sky survey to derive number counts and spectral indices of extragalactic sources -- infrared and radio sources -- from the Planck Early Catalogue (ERCSC) at 100 to 857GHz. Our sample contains, after the 80% completeness cut, between 122 and 452 and sources, with flux densities above 0.3 and 1.9Jy at 100 and 857GHz, over about 31 to 40% of the sky. Using Planck HFI, all the sources have been classified as either dust-dominated or synchrotron-dominated on the basis of their spectral energy distributions (SED). Our sample is thus complete, flux-limited and color-selected to differentiate between the two populations. We find an approximately equal number of synchrotron and dusty sources between 217 and 353GHz; at 353GHz or higher (or 217GHz and lower) frequencies, the number is dominated by dusty (synchrotron) sources, as expected. For most of the sources, the spectral indices are also derived. We provide for the first time counts of bright sources from 353 to 857GHz and the contributions from dusty and synchrotron sources at all HFI frequencies in the key spectral range where these spectra are crossing. The observed counts are in the Euclidean regime. The number counts are compared to previously published data (earlier Planck, Herschel, BLAST, SCUBA, LABOCA, SPT, and ACT) and models taking into account both radio or infrared galaxies. We derive the multi-frequency Euclidean level and compare it to WMAP, Spitzer and IRAS results. The submillimetre number counts are not well reproduced by current evolution models of dusty galaxies, whereas the millimetre part appears reasonably well fitted by the most recent model for synchrotron-dominated sources. Finally we provide estimates of the local luminosity density of dusty galaxies, providing the first such measurements at 545 and 857GHz.

preprint2013arXiv

Planck intermediate results. XII: Diffuse Galactic components in the Gould Belt System

We perform an analysis of the diffuse low-frequency Galactic components in the Southern part of the Gould Belt system (130^\circ\leq l\leq 230^\circ and -50^\circ\leq b\leq -10^\circ). Strong ultra-violet (UV) flux coming from the Gould Belt super-association is responsible for bright diffuse foregrounds that we observe from our position inside the system and that can help us improve our knowledge of the Galactic emission. Free-free emission and anomalous microwave emission (AME) are the dominant components at low frequencies (ν< 40 GHz), while synchrotron emission is very smooth and faint. We separate diffuse free-free emission and AME from synchrotron emission and thermal dust emission by using Planck data, complemented by ancillary data, using the "Correlated Component Analysis" (CCA) component separation method and we compare with the results of cross-correlation of foreground templates with the frequency maps. We estimate the electron temperature T_e from H$α$ and free-free emission using two methods (temperature-temperature plot and cross-correlation) and we obtain T_e ranging from 3100 to 5200 K, for an effective fraction of absorbing dust along the line of sight of 30% (f_d=0.3). We estimate the frequency spectrum of the diffuse AME and we recover a peak frequency (in flux density units) of 25.5 \pm 1.5 GHz. We verify the reliability of this result with realistic simulations that include the presence of biases in the spectral model for the AME and in the free-free template. By combining physical models for vibrational and rotational dust emission and adding the constraints from the thermal dust spectrum from Planck and IRAS we are able to get a good description of the frequency spectrum of the AME for plausible values of the local density and radiation field.

preprint2013arXiv

Planck intermediate results. XVI. Profile likelihoods for cosmological parameters

We explore the 2013 Planck likelihood function with a high-precision multi-dimensional minimizer (Minuit). This allows a refinement of the Lambda-cdm best-fit solution with respect to previously-released results, and the construction of frequentist confidence intervals using profile likelihoods. The agreement with the cosmological results from the Bayesian framework is excellent, demonstrating the robustness of the Planck results to the statistical methodology. We investigate the inclusion of neutrino masses, where more significant differences may appear due to the non-Gaussian nature of the posterior mass distribution. By applying the Feldman--Cousins prescription, we again obtain results very similar to those of the Bayesian methodology. However, the profile-likelihood analysis of the CMB combination (Planck+WP+highL) reveals a minimum well within the unphysical negative-mass region. We show that inclusion of the Planck CMB-lensing information regularizes this issue, and provide a robust frequentist upper limit $M_ν< 0.26 eV$ ($95%$ confidence) from the CMB+lensing+BAO data combination.

preprint2012arXiv

A new L-dwarf member of the moderately metal-poor triple system HD 221356

We report on the discovery of a fourth component in the HD 221356 star system, previously known to be formed by an F8V, slightly metal-poor primary ([Fe/H]=-0.26), and a distant M8V+L3V pair. In our ongoing common proper motion search based on VISTA Hemisphere Survey (VHS) and 2MASS catalogues, we have detected a faint (J=13.76+/-0.04 mag) co-moving companion of the F8 star located at angular separation of 12.13+/-0.18 arcsec (position angle of 221.8+/-1.7), corresponding to a projected distance of ~312 AU at 26 pc. Near-infrared spectroscopy of the new companion, covering the 1.5-2.4 micron wavelength range with a resolving power of R~600, indicates an L1+/-1 spectral type. Using evolutionary models the mass of the new companion is estimated at ~0.08 solar masses, which places the object close to the stellar-substellar borderline. This multiple system provides an interesting example of objects with masses slightly above and below the hydrogen burning mass limit. The low mass companions of HD 221356 have slightly bluer colours than field dwarfs with similar spectral type, which is likely a consequence of the sub-solar metallicity of the system.

preprint2012arXiv

First T dwarfs in the VISTA Hemisphere Survey

The aim of the project is to improve our current knowledge of the density of T dwarfs and the shape of the substellar initial mass function by identifying a magnitude-limited sample of T dwarfs in the full southern sky. We present the results of a photometric search aimed at discovering cool brown dwarfs in the Southern sky imaged at infrared wavelengths by the Visible and Infrared Survey Telescope for Astronomy (VISTA) and the Wide Infrared Survey Explorer (WISE) satellite mission. We combined the first data release (DR1) of the VISTA Hemisphere Survey (VHS) and the WISE preliminary data release to extract candidates with red mid-infrared colours and near- to mid-infrared colours characteristics of cool brown dwarfs. The VHS DR1 vs WISE search returned tens of T dwarf candidates, 13 of which are presented here, including two previously published in the literature and five new ones confirmed spectroscopically with spectral types between T4.5 and T8. We estimate that the two T6 dwarfs lie within 16 pc and the T4.5 within 25 pc. The remaining three are 30-50 pc distant. The only T7 dwarf in our sample is the faintest of its spectral class with J=19.28 mag. The other six T dwarf candidates remain without spectroscopic follow-up. We also improve our knowledge on the proper motion accuracy for three bright T dwarfs by combining multi-epoch data from public databases (DENIS, 2MASS, VHS, WISE, Spitzer).

preprint2012arXiv

Holographic Imaging of Crowded Fields: High Angular Resolution Imaging with Excellent Quality at Very Low Cost

We present a method for speckle holography that is optimised for crowded fields. Its two key features are an iterativ improvement of the instantaneous Point Spread Functions (PSFs) extracted from each speckle frame and the (optional) simultaneous use of multiple reference stars. In this way, high signal-to-noise and accuracy can be achieved on the PSF for each short exposure, which results in sensitive, high-Strehl re- constructed images. We have tested our method with different instruments, on a range of targets, and from the N- to the I-band. In terms of PSF cosmetics, stability and Strehl ratio, holographic imaging can be equal, and even superior, to the capabilities of currently available Adaptive Optics (AO) systems, particularly at short near-infrared to optical wavelengths. It outperforms lucky imaging because it makes use of the entire PSF and reduces the need for frame selection, thus leading to higher Strehl and improved sensitivity. Image reconstruction a posteriori, the possibility to use multiple reference stars and the fact that these reference stars can be rather faint means that holographic imaging offers a simple way to image large, dense stellar fields near the diffraction limit of large telescopes, similar to, but much less technologically demanding than, the capabilities of a multi-conjugate adaptive optics system. The method can be used with a large range of already existing imaging instruments and can also be combined with AO imaging when the corrected PSF is unstable.

preprint2012arXiv

Planck Early Results. II. The thermal performance of Planck

The performance of the Planck instruments in space is enabled by their low operating temperatures, 20K for LFI and 0.1K for HFI, achieved through a combination of passive radiative cooling and three active mechanical coolers. The scientific requirement for very broad frequency coverage led to two detector technologies with widely different temperature and cooling needs. Active coolers could satisfy these needs; a helium cryostat, as used by previous cryogenic space missions (IRAS, COBE, ISO, Spitzer, AKARI), could not. Radiative cooling is provided by three V-groove radiators and a large telescope baffle. The active coolers are a hydrogen sorption cooler (<20K), a 4He Joule-Thomson cooler (4.7K), and a 3He-4He dilution cooler (1.4K and 0.1K). The flight system was at ambient temperature at launch and cooled in space to operating conditions. The HFI bolometer plate reached 93mK on 3 July 2009, 50 days after launch. The solar panel always faces the Sun, shadowing the rest of Planck, andoperates at a mean temperature of 384K. At the other end of the spacecraft, the telescope baffle operates at 42.3K and the telescope primary mirror operates at 35.9K. The temperatures of key parts of the instruments are stabilized by both active and passive methods. Temperature fluctuations are driven by changes in the distance from the Sun, sorption cooler cycling and fluctuations in gas-liquid flow, and fluctuations in cosmic ray flux on the dilution and bolometer plates. These fluctuations do not compromise the science data.

preprint2012arXiv

Planck Intermediate Results II: Comparison of Sunyaev-Zeldovich measurements from Planck and from the Arcminute Microkelvin Imager for 11 galaxy clusters

A comparison is presented of Sunyaev-Zeldovich measurements for 11 galaxy clusters as obtained by Planck and by the ground-based interferometer, the Arcminute Microkelvin Imager. Assuming a universal spherically-symmetric Generalised Navarro, Frenk & White (GNFW) model for the cluster gas pressure profile, we jointly constrain the integrated Compton-Y parameter (Y_500) and the scale radius (theta_500) of each cluster. Our resulting constraints in the Y_500-theta_500 2D parameter space derived from the two instruments overlap significantly for eight of the clusters, although, overall, there is a tendency for AMI to find the Sunyaev-Zeldovich signal to be smaller in angular size and fainter than Planck. Significant discrepancies exist for the three remaining clusters in the sample, namely A1413, A1914, and the newly-discovered Planck cluster PLCKESZ G139.59+24.18. The robustness of the analysis of both the Planck and AMI data is demonstrated through the use of detailed simulations, which also discount confusion from residual point (radio) sources and from diffuse astrophysical foregrounds as possible explanations for the discrepancies found. For a subset of our cluster sample, we have investigated the dependence of our results on the assumed pressure profile by repeating the analysis adopting the best-fitting GNFW profile shape which best matches X-ray observations. Adopting the best-fitting profile shape from the X-ray data does not, in general, resolve the discrepancies found in this subset of five clusters. Though based on a small sample, our results suggest that the adopted GNFW model may not be sufficiently flexible to describe clusters universally.

preprint2012arXiv

Planck Intermediate Results. I. Further validation of new Planck clusters with XMM-Newton

We present further results from the ongoing XMM-Newton validation follow-up of Planck cluster candidates, detailing X-ray observations of eleven candidates detected at a signal-to-noise ratio of 4.5<S/N<5.3 in the same 10-month survey maps used in the construction of the Early SZ sample. The sample was selected in order to test internal SZ quality flags, and the pertinence of these flags is discussed in light of the validation results. Ten of the candidates are found to be bona fide clusters lying below the RASS flux limit. Redshift estimates are available for all confirmed systems via X-ray Fe-line spectroscopy. They lie in the redshift range 0.19<z<0.94, demonstrating Planck's capability to detect clusters up to high z. The X-ray properties of the new clusters appear to be similar to previous new detections by Planck at lower z and higher SZ flux: the majority are X-ray underluminous for their mass, estimated using Y_X as mass proxy, and many have a disturbed morphology. We find tentative indication for Malmquist bias in the Y_SZ-Y_X relation, with a turnover at Y_SZ \sim 4 e-4 arcmin^2. We present additional new optical redshift determinations with ENO and ESO telescopes of candidates previously confirmed with XMM-Newton. The X-ray and optical redshifts for a total of 20 clusters are found to be in excellent agreement. We also show that useful lower limits can be put on cluster redshifts using X-ray data only via the use of the Y_X vs. Y_SZ and X-ray flux F_X vs. Y_SZ relations.

preprint2012arXiv

Planck intermediate results. III. The relation between galaxy cluster mass and Sunyaev-Zeldovich signal

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.

preprint2012arXiv

Planck Intermediate Results. IX. Detection of the Galactic haze with Planck

Using precise full-sky observations from Planck, and applying several methods of component separation, we identify and characterize the emission from the Galactic "haze" at microwave wavelengths. The haze is a distinct component of diffuse Galactic emission, roughly centered on the Galactic centre, and extends to |b| ~35 deg in Galactic latitude and |l| ~15 deg in longitude. By combining the Planck data with observations from the WMAP we are able to determine the spectrum of this emission to high accuracy, unhindered by the large systematic biases present in previous analyses. The derived spectrum is consistent with power-law emission with a spectral index of -2.55 +/- 0.05, thus excluding free-free emission as the source and instead favouring hard-spectrum synchrotron radiation from an electron population with a spectrum (number density per energy) dN/dE ~ E^-2.1. At Galactic latitudes |b|<30 deg, the microwave haze morphology is consistent with that of the Fermi gamma-ray "haze" or "bubbles," indicating that we have a multi-wavelength view of a distinct component of our Galaxy. Given both the very hard spectrum and the extended nature of the emission, it is highly unlikely that the haze electrons result from supernova shocks in the Galactic disk. Instead, a new mechanism for cosmic-ray acceleration in the centre of our Galaxy is implied.

preprint2012arXiv

Planck Intermediate Results. V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect

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 3 x R500, 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-Newton 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 slightly flatter than most predictions from numerical simulations. 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 3 x R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

preprint2012arXiv

Planck intermediate results. VI: The dynamical structure of PLCKG214.6+37.0, a Planck discovered triple system of galaxy clusters

The survey of galaxy clusters performed by Planck through the Sunyaev-Zeldovich effect has already discovered many interesting objects, thanks to the whole coverage of the sky. One of the SZ candidates detected in the early months of the mission near to the signal to noise threshold, PLCKG214.6+37.0, was later revealed by XMM-Newton to be a triple system of galaxy clusters. We have further investigated this puzzling system with a multi-wavelength approach and we present here the results from a deep XMM-Newton re-observation. The characterisation of the physical properties of the three components has allowed us to build a template model to extract the total SZ signal of this system with Planck data. We partly reconciled the discrepancy between the expected SZ signal from X-rays and the observed one, which are now consistent at less than 1.2 sigma. We measured the redshift of the three components with the iron lines in the X-ray spectrum, and confirmed that the three clumps are likely part of the same supercluster structure. The analysis of the dynamical state of the three components, as well as the absence of detectable excess X-ray emission, suggest that we are witnessing the formation of a massive cluster at an early phase of interaction.

preprint2012arXiv

Planck intermediate results. VIII. Filaments between interacting clusters

About half of the baryons of the Universe are expected to be in the form of filaments of hot and low density intergalactic medium. Most of these baryons remain undetected even by the most advanced X-ray observatories which are limited in sensitivity to the diffuse low density medium. The Planck satellite has provided hundreds of detections of the hot gas in clusters of galaxies via the thermal Sunyaev-Zel'dovich (tSZ) effect and is an ideal instrument for studying extended low density media through the tSZ effect. In this paper we use the Planck data to search for signatures of a fraction of these missing baryons between pairs of galaxy clusters. Cluster pairs are good candidates for searching for the hotter and denser phase of the intergalactic medium (which is more easily observed through the SZ effect). Using an X-ray catalogue of clusters and the Planck data, we select physical pairs of clusters as candidates. Using the Planck data we construct a local map of the tSZ effect centered on each pair of galaxy clusters. ROSAT data is used to construct X-ray maps of these pairs. After having modelled and subtracted the tSZ effect and X-ray emission for each cluster in the pair we study the residuals on both the SZ and X-ray maps. For the merging cluster pair A399-A401 we observe a significant tSZ effect signal in the intercluster region beyond the virial radii of the clusters. A joint X-ray SZ analysis allows us to constrain the temperature and density of this intercluster medium. We obtain a temperature of kT = 7.1 +- 0.9, keV (consistent with previous estimates) and a baryon density of (3.7 +- 0.2)x10^-4, cm^-3. The Planck satellite mission has provided the first SZ detection of the hot and diffuse intercluster gas.

preprint2012arXiv

Planck intermediate results. X. Physics of the hot gas in the Coma cluster

We present an analysis of Planck satellite data on the Coma Cluster observed via the Sunyaev-Zeldovich effect. Planck is able, for the first time, to detect SZ emission up to r ~ 3 X R_500. We test previously proposed models for the pressure distribution in clusters against the azimuthally averaged data. We find that the Arnaud et al. universal pressure profile does not fit Coma, and that their pressure profile for merging systems provides a good fit of the data only at r<R_500: by r=2XR_500 it underestimates the observed y profile by a factor of ~2. This may indicate that at these larger radii either i) the cluster SZ emission is contaminated by unresolved SZ sources along the line of sight or ii) the pressure profile of Coma is higher at r>R_500 than the mean pressure profile predicted by the simulations. The Planck image shows significant local steepening of the y profile in two regions about half a degree to the west and to the south-east of the cluster centre. These features are consistent with the presence of shock fronts at these radii, and indeed the western feature was previously noticed in the ROSAT PSPC mosaic as well as in the radio. Using Planck y profiles extracted from corresponding sectors we find pressure jumps of 4.5+0.4-0.2 and 5.0+1.3-0.1 in the west and southeast, respectively. Assuming Rankine-Hugoniot pressure jump conditions, we deduce that the shock waves should propagate with Mach number M_w=2.03+0.09-0.04 and M_se=2.05+0.25-0.02 in the West and Southeast, respectively. Finally, we find that the y and radio-synchrotron signals are quasi-linearly correlated on Mpc scales with small intrinsic scatter. This implies either that the energy density of cosmic-ray electrons is relatively constant throughout the cluster, or that the magnetic fields fall off much more slowly with radius than previously thought.

preprint2012arXiv

Planck Intermediate Results. XI: The gas content of dark matter halos: the Sunyaev-Zeldovich-stellar mass relation for locally brightest galaxies

We present the scaling relation between Sunyaev-Zeldovich (SZ) signal and stellar mass for almost 260,000 locally brightest galaxies (LBGs) selected from the Sloan Digital Sky Survey (SDSS). These are predominantly the central galaxies of their dark matter halos. We calibrate the stellar-to-halo mass conversion using realistic mock catalogues based on the Millennium Simulation. Applying a multi-frequency matched filter to the Planck data for each LBG, and averaging the results in bins of stellar mass, we measure the mean SZ signal down to $M_\ast\sim 2\times 10^{11} \Msolar$, with a clear indication of signal at even lower stellar mass. We derive the scaling relation between SZ signal and halo mass by assigning halo properties from our mock catalogues to the real LBGs and simulating the Planck observation process. This relation shows no evidence for deviation from a power law over a halo mass range extending from rich clusters down to $M_{500}\sim 2\times 10^{13} \Msolar$, and there is a clear indication of signal down to $M_{500}\sim 4\times 10^{12} \Msolar$. Planck's SZ detections in such low-mass halos imply that about a quarter of all baryons have now been seen in the form of hot halo gas, and that this gas must be less concentrated than the dark matter in such halos in order to remain consistent with X-ray observations. At the high-mass end, the measured SZ signal is 20% lower than found from observations of X-ray clusters, a difference consistent with Malmquist bias effects in the X-ray sample.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: constraints on the time variation of fundamental constants from the large-scale two-point correlation function

We obtain constraints on the variation of the fundamental constants from the full shape of the redshift-space correlation function of a sample of luminous galaxies drawn from the Data Release 9 of the Baryonic Oscillations Spectroscopic Survey. We combine this information with data from recent CMB, BAO and H_0 measurements. We focus on possible variations of the fine structure constant αand the electron mass m_e in the early universe, and study the degeneracies between these constants and other cosmological parameters, such as the dark energy equation of state parameter w_DE, the massive neutrinos fraction f_ν, the effective number of relativistic species N_eff, and the primordial helium abundance Y_He. When only one of the fundamental constants is varied, our final bounds are α/ α_0 = 0.9957_{-0.0042}^{+0.0041} and m_e /(m_e)_0 = 1.006_{-0.013}^{+0.014}. For their joint variation, our results are α/ α_0 = 0.9901_{-0.0054}^{+0.0055} and m_e /(m_e)_0 = 1.028 +/- 0.019. Although when m_e is allowed to vary our constraints on w_DE are consistent with a cosmological constant, when αis treated as a free parameter we find w_DE = -1.20 +/- 0.13; more than 1 σaway from its standard value. When f_νand αare allowed to vary simultaneously, we find f_ν< 0.043 (95% CL), implying a limit of \sum m_ν< 0.46 eV (95% CL), while for m_e variation, we obtain f_nu < 0.086 (95% CL), which implies \sum m_ν< 1.1 eV (95% CL). When N_eff or Y_He are considered as free parameters, their simultaneous variation with αprovides constraints close to their standard values (when the H_0 prior is not included in the analysis), while when m_e is allowed to vary, their preferred values are significantly higher. In all cases, our results are consistent with no variations of αor m_e at the 1 or 2 σlevel.

preprint2011arXiv

A Planet of an A-Star: HD15082b

Most of the known transiting extrasolar planets orbit slowly rotating F, G or K stars. In here we report on the detection of a transiting planet orbiting the bright, rapidly rotating A5 star HD15082 (WASP-33b, V=8.3, m sin i = 86 km s-1), recently made by SuperWASP. Time resolved spectroscopic observations taken during transit show a hump caused by the planet crossing the line profile. From the analysis of the spectra, we derive the radius of the planet and find that it is orbiting retrograde in respect to the spin of the star. Because of its small distance from an A5 star (the orbital period of only 1.22 days), the equilibrium temperature of the planet is estimated to be 2700 K. The planet thus is one of the hottest planets known, which makes it relatively easy to detect it in the IR. We thus tried to detect it using the TNG but did not succeed.

preprint2011arXiv

Be abundances in cool main-sequence stars with exoplanets

We present new UVES spectra of a sample of 15 cool unevolved stars with and without detected planetary companions. Together with previous determinations, we study Be depletion and possible differences in Be abundances between both groups of stars. We obtain a final sample of 89 and 40 stars with and without planets, respectively, which covers a wide range of effective temperatures, from 4700 K to 6400 K, and includes several cool dwarf stars for the first time. We determine Be abundances for these stars and find that for most of them (the coolest ones) the BeII resonance lines are often undetectable, implying significant Be depletion. While for hot stars Be abundances are aproximately constant, with a slight fall as Teff decreases and the Li-Be gap around 6300 K, we find a steep drop of Be content as Teff decreases for Teff < 5500 K, confirming the results of previous papers. Therefore, for these stars there is an unknown mechanism destroying Be that is not reflected in current models of Be depletion. Moreover, this strong Be depletion in cool objects takes place for all the stars regardless of the presence of planets, thus, the effect of extra Li depletion in solar-type stars with planets when compared with stars without detected planets does not seem to be present for Be, although the number of stars at those temperatures is still small to reach a final conclusion.

preprint2011arXiv

Beryllium abundances in stars with planets:Extending the sample

Context: Chemical abundances of light elements as beryllium in planet-host stars allow us to study the planet formation scenarios and/or investigate possible surface pollution processes. Aims: We present here an extension of previous beryllium abundance studies. The complete sample consists of 70 stars hosting planets and 30 stars without known planetary companions. The aim of this paper is to further assess the trends found in previous studies with less number of objects. This will provide more information on the processes of depletion and mixing of light elements in the interior of late type stars, and will provide possible explanations for the abundance differences between stars that host planets and "single" stars. Methods: Using high resolution UVES spectra, we measure beryllium abundances of 26 stars that host planets and 1 "single" star mainly using the lambda 3131.065 A Be II line, by fitting synthetic spectra to the observational data. We also compile beryllium abundance measurements of 44 stars hosting planets and 29 "single" stars from the literature, resulting in a final sample of 100 objects. Results: We confirm that the beryllium content is roughly the same in stars hosting planets and in "single" stars at temperatures Teff > 5700K. The sample is still small for Teff < 5500K, but it seems that the scatter in Be abundances of dwarf stars is slightly higher at these cooler temperatures. Conclusions: We search for distinctive characteristics of planet hosts through correlations of Be abundance versus Li abundance, age, metallicity and oxygen abundance. These could provide some insight in the formation and evolution of planetary systems, but we did not find any clear correlation.

preprint2011arXiv

Detection of Anomalous Microwave Emission in the Pleiades Reflection Nebula with WMAP and the COSMOSOMAS Experiment

We present evidence for anomalous microwave emission (AME) in the Pleiades reflection nebula, using data from the seven-year release of the Wilkinson Microwave Anisotropy Probe (WMAP) and from the COSMOSOMAS experiment. The flux integrated in a 1-degree radius around R.A.=56.24^{\circ}, Dec.=23.78^{\circ} (J2000) is 2.15 +/- 0.12 Jy at 22.8 GHz, where AME is dominant. COSMOSOMAS data show no significant emission, but allow to set upper limits of 0.94 and 1.58 Jy (99.7% C.L.) respectively at 10.9 and 14.7 GHz, which are crucial to pin down the AME spectrum at these frequencies, and to discard any other emission mechanisms which could have an important contribution to the signal detected at 22.8 GHz. We estimate the expected level of free-free emission from an extinction-corrected H-alpha template, while the thermal dust emission is characterized from infrared DIRBE data and extrapolated to microwave frequencies. When we deduct the contribution from these two components at 22.8 GHz the residual flux, associated with AME, is 2.12 +/- 0.12 Jy (17.7-sigma). The spectral energy distribution from 10 to 60 GHz can be accurately fitted with a model of electric dipole emission from small spinning dust grains distributed in two separated phases of molecular and atomic gas, respectively. The dust emissivity, calculated by correlating the 22.8 GHz data with 100-micron data, is found to be 4.36+/-0.17 muK/MJy/sr, a value that is rather low compared with typical values in dust clouds. The physical properties of the Pleiades nebula indicate that this is indeed a much less opaque object than others were AME has usually been detected. This fact, together with the broad knowledge of the stellar content of this region, provides an excellent testbed for AME characterization in physical conditions different from those generally explored up to now.

preprint2011arXiv

Lucky Imaging Adaptive Optics of the brown dwarf binary GJ569Bab

The potential of combining Adaptive Optics (AO) and Lucky Imaging (LI) to achieve high precision astrometry and differential photometry in the optical is investigated by conducting observations of the close 0\farcs1 brown dwarf binary GJ569Bab. We took 50000 $I$-band images with our LI instrument FastCam attached to NAOMI, the 4.2-m William Herschel Telescope (WHT) AO facility. In order to extract the most of the astrometry and photometry of the GJ569Bab system we have resorted to a PSF fitting technique using the primary star GJ569A as a suitable PSF reference which exhibits an $I$-band magnitude of $7.78\pm0.03$. The AO+LI observations at WHT were able to resolve the binary system GJ569Bab located at $4\farcs 92 \pm 0\farcs05$ from GJ569A. We measure a separation of $98.4 \pm 1.1$ mas and $I$-band magnitudes of $13.86 \pm 0.03$ and $14.48 \pm 0.03$ and $I-J$ colors of 2.72$\pm$0.08 and 2.83$\pm$0.08 for the Ba and Bb components, respectively. Our study rules out the presence of any other companion to GJ569A down to magnitude I$\sim$ 17 at distances larger than 1\arcsec. The $I-J$ colors measured are consistent with M8.5-M9 spectral types for the Ba and Bb components. The available dynamical, photometric and spectroscopic data are consistent with a binary system with Ba being slightly (10-20%) more massive than Bb. We obtain new orbital parameters which are in good agreement with those in the literature.

preprint2011arXiv

Near-infrared Linear Polarization of Ultracool Dwarfs

We report on near-infrared J- and H-band linear polarimetric photometry of eight ultracool dwarfs (two late-M, five L0-L7.5, and one T2.5) with known evidence for photometric variability due to dust clouds, anomalous red infrared colors, or low-gravity atmospheres. The polarimetric data were acquired with the LIRIS instrument on the William Herschel Telescope. We also provide mid-infrared photometry in the interval 3.4-24 micron for some targets obtained with Spitzer and WISE, which has allowed us to confirm the peculiar red colors of five sources in the sample. We can impose modest upper limits of 0.9% and 1.8% on the linear polarization degree for seven targets with a confidence of 99%. Only one source, 2MAS, J02411151-0326587 (L0), appears to be strongly polarized (P ~ 3%) in the J-band with a significance level of P/sigma_P ~ 10. The likely origin of its linearly polarized light and rather red infrared colors may reside in a surrounding disk with an asymmetric distribution of grains. Given its proximity (66 +/- 8 pc), this object becomes an excellent target for the direct detection of the disk.

preprint2011arXiv

Planck Early Results VIII: The all-sky Early Sunyaev-Zeldovich cluster sample

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

preprint2011arXiv

Planck Early Results XI: Calibration of the local galaxy cluster Sunyaev-Zeldovich scaling relations

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters.

preprint2011arXiv

Planck early results XIV: ERCSC validation and extreme radio sources

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

preprint2011arXiv

Planck Early Results XXVI: Detection with Planck and confirmation by XMM-Newton of PLCK G266.6-27.3, an exceptionally X-ray luminous and massive galaxy cluster at z~1

We present first results on PLCK G266.6-27.3, a galaxy cluster candidate detected at a signal-to-noise ratio of 5 in the Planck All Sky survey. An XMM-Newton validation observation has allowed us to confirm that the candidate is a bona fide galaxy cluster. With these X-ray data we measure an accurate redshift, z = 0.94 +/- 0.02, and estimate the cluster mass to be M_500 = (7.8 +/- 0.8)e+14 solar masses. PLCK G266.6-27.3 is an exceptional system: its luminosity of L_X(0.5-2.0 keV)=(1.4 +/- 0.05)e+45 erg/s, equals that of the two most luminous known clusters in the z > 0.5 universe, and it is one of the most massive clusters at z~1. Moreover, unlike the majority of high-redshift clusters, PLCK G266.6-27.3 appears to be highly relaxed. This observation confirms Planck's capability of detecting high-redshift, high-mass clusters, and opens the way to the systematic study of population evolution in the exponential tail of the mass function.

preprint2011arXiv

Planck early results: Cluster Sunyaev-Zeldovich optical scaling relations

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters.

preprint2011arXiv

Planck Early Results: Origin of the submm excess dust emission in the Magellanic Clouds

The integrated Spectral Energy Distributions of the Large and Small Magellanic Cloud appear significantly flatter than expected from dust models based on their FIR and radio emission. The origin of this millimetre excess is still unexplained, and is here investigated using the Planck data. The background CMB contribution is subtracted using an ILC method performed locally around the galaxies. The foreground emission from the Milky Way is subtracted. After subtraction, the emission of both galaxies correlates closely with the gas emission of the LMC and SMC. The millimetre excess in the LMC can be explained by CMB fluctuations, but a significant excess is still present in the SMC SED. The Planck and IRIS data at 100 micron are combined to produce thermal dust temperature and optical depth maps of the two galaxies. The LMC temperature map shows the presence of a warm inner arm already found with the Spitzer data, but also shows the existence of a previously unidentified cold outer arm. Several cold regions are found along this arm, some of which are associated with known molecular clouds. The average emissivity spectral index is found to be consistent with beta=1.5 and beta=1.2 below 500 microns for the LMC and SMC respectively, significantly flatter than the values observed in the Milky Way. Furthermore, there is evidence in the SMC for a further flattening of the SED in the sub-mm. The spatial distribution of the millimetre dust excess in the SMC follows the gas and thermal dust distribution. Different models are explored in order to fit the dust emission in the SMC. It is concluded that the millimetre excess is unlikely to be caused by very cold dust emission and that it could be due to a combination of spinning dust. emission and thermal dust emission by more amorphous dust grains than those present in our Galaxy.

preprint2011arXiv

Planck Early Results: Statistical properties of extragalactic radio sources in the Planck Early Release Compact Source Catalogue

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.

preprint2011arXiv

Planck Early Results: The Galactic Cold Core Population revealed by the first all-sky survey

We present the statistical properties of the first version of the Cold Core Catalogue of Planck Objects (C3PO), in terms of their spatial distribution, temperature, distance, mass, and morphology. We also describe the statistics of the Early Cold Core Catalogue (ECC, delivered with the Early Release Compact Source Catalogue, ERCSC) that is the subset of the 915 most reliable detections of the complete catalogue. We have used the CoCoCoDeT algorithm to extract 10783 cold sources. Temperature and dust emission spectral index β values are derived using the fluxes in the IRAS 100 \mum band and the three highest frequency Planck bands. Temperature spans from 7K to 17K, and peaks around 13K. Data are not consistent with a constant value of β over the all temperature range. β ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including β(T) and the effect of multiple T components folded into the measurements. For one third of the objects the distances are obtained. Most of the detections are within 2 kpc in the Solar neighbourhood, but a few are at distances greater than 4 kpc. The cores are distributed from the deep Galactic plane, despite the confusion, to high latitudes (>30$^{\circle}$). The associated mass estimates range from 1 to $10^5$ solar masses. Using their physical properties these cold sources are shown to be cold clumps, defined as the intermediate cold sub-structures between clouds and cores. These cold clumps are not isolated but mostly organized in filaments associated with molecular clouds. The Cold Core Catalogue of Planck Objects (C3PO) is the first unbiased all-sky catalogue of cold objects. It gives an unprecedented statistical view to the properties of these potential pre-stellar clumps and offers a unique possibility for their classification in terms of their intrinsic properties and environment.

preprint2011arXiv

Planck Early Results: The Planck mission

The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission.

preprint2011arXiv

Planck Early Results: The Planck View of Nearby Galaxies

The all-sky coverage of the Planck Early Release Compact Source Catalogue (ERCSC) provides an unsurpassed survey of galaxies at submillimetre (submm) wavelengths, representing a major improvement in the numbers of galaxies detected, as well as the range of far-IR/submm wavelengths over which they have been observed. We here present the first results on the properties of nearby galaxies using these data. We match the ERCSC catalogue to IRAS-detected galaxies in the Imperial IRAS Faint Source Redshift Catalogue (IIFSCz), so that we can measure the spectral energy distributions (SEDs) of these objects from 60 to 850 microns. This produces a list of 1717 galaxies with reliable associations between Planck and IRAS, from which we select a subset of 468 for SED studies, namely those with strong detections in the three highest frequency Planck bands and no evidence of cirrus contamination. The SEDs are fitted using parametric dust models to determine the range of dust temperatures and emissivities. We find evidence for colder dust than has previously been found in external galaxies, with T<20K. Such cold temperatures are found using both the standard single temperature dust model with variable emissivity beta, or a two dust temperature model with beta fixed at 2. We also compare our results to studies of distant submm galaxies (SMGs) which have been claimed to contain cooler dust than their local counterparts. We find that including our sample of 468 galaxies significantly reduces the distinction between the two populations. Fits to SEDs of selected objects using more sophisticated templates derived from radiative transfer models confirm the presence of the colder dust found through parameteric fitting. We thus conclude that cold (T<20K) dust is a significant and largely unexplored component of many nearby galaxies.

preprint2011arXiv

Planck Early Results: The submillimetre properties of a sample of Galactic cold clumps

(abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the substructures at scales below the Planck beam size.

preprint2011arXiv

Planck Early Results: Thermal dust in Nearby Molecular Clouds

Planck allows unbiased mapping of Galactic sub-millimetre and millimetre emission from the most diffuse regions to the densest parts of molecular clouds. We present an early analysis of the Taurus molecular complex, on line-of-sight-averaged data and without component separation. The emission spectrum measured by Planck and IRAS can be fitted pixel by pixel using a single modified blackbody. Some systematic residuals are detected at 353 GHz and 143 GHz, with amplitudes around -7 % and +13 %, respectively, indicating that the measured spectra are likely more complex than a simple modified blackbody. Significant positive residuals are also detected in the molecular regions and in the 217 GHz and 100 GHz bands, mainly caused by to the contribution of the J=2-1 and J=1-0 12CO and 13CO emission lines. We derive maps of the dust temperature T, the dust spectral emissivity index beta, and the dust optical depth at 250 microns tau. The temperature map illustrates the cooling of the dust particles in thermal equilibrium with the incident radiation field, from 16-17 K in the diffuse regions to 13-14 K in the dense parts. The distribution of spectral indices is centred at 1.78, with a standard deviation of 0.08 and a systematic error of 0.07. We detect a significant T-beta anti-correlation. The dust optical depth map reveals the spatial distribution of the column density of the molecular complex from the densest molecular regions to the faint diffuse regions. We use near-infrared extinction and HI data at 21-cm to perform a quantitative analysis of the spatial variations of the measured dust optical depth at 250 microns per hydrogen atom tau/NH. We report an increase of tau/NH by a factor of about 2 between the atomic phase and the molecular phase, which has a strong impact on the equilibrium temperature of the dust particles.

preprint2011arXiv

Planck early results. III. First assessment of the Low Frequency Instrument in-flight performance

The scientific performance of the Planck Low Frequency Instrument (LFI) after one year of in-orbit operation is presented. We describe the main optical parameters and discuss photometric calibration, white noise sensitivity, and noise properties. A preliminary evaluation of the impact of the main systematic effects is presented. For each of the performance parameters, we outline the methods used to obtain them from the flight data and provide a comparison with pre-launch ground assessments, which are essentially confirmed in flight.

preprint2011arXiv

Planck early results. IX. XMM-Newton follow-up for validation of Planck cluster candidates

We present the XMM-Newton follow-up for confirmation of Planck cluster candidates. Twenty-five candidates have been observed to date using snapshot (~10 ksec) exposures, ten as part of a pilot programme to sample a low range of signal-to-noise ratios (4<S/N<6), and a further 15 in a programme to observe a sample of S/N>5 candidates. The sensitivity and spatial resolution of XMM-Newton allows unambiguous discrimination between clusters and false candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates are confirmed as extended X-ray sources. Seventeen are single clusters, the majority of which are found to have highly irregular and disturbed morphologies (about ~70%). The remaining four sources are multiple systems, including the unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to derive a redshift estimate from the X-ray Fe K line (albeit of variable quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a median redshift of z~0.37. A first determination is made of their X-ray properties including the characteristic size, which is used to improve the estimate of the SZ Compton parameter, Y_SZ. The follow-up validation programme has helped to optimise the Planck candidate selection process. It has also provided a preview of the X-ray properties of these newly-discovered clusters, allowing comparison with their SZ properties, and to the X-ray and SZ properties of known clusters observed in the Planck survey. Our results suggest that Planck may have started to reveal a non-negligible population of massive dynamically perturbed objects that is under-represented in X-ray surveys. However, despite their particular properties, these new clusters appear to follow the Y_SZ-Y_X relation established for X-ray selected objects, where Y_X is the product of the gas mass and temperature.

preprint2011arXiv

Planck Early Results. V. The Low Frequency Instrument data processing

We describe the processing of data from the Low Frequency Instrument (LFI) used in production of the Planck Early Release Compact Source Catalogue (ERCSC). In particular, we discuss the steps involved in reducing the data from telemetry packets to cleaned, calibrated, time-ordered data (TOD) and frequency maps. Data are continuously calibrated using the modulation of the temperature of the cosmic microwave background radiation induced by the motion of the spacecraft. Noise properties are estimated from TOD from which the sky signal has been removed using a generalized least square map-making algorithm. Measured 1/f noise knee-frequencies range from 100mHz at 30GHz to a few tens of mHz at 70GHz. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices required to compute statistical uncertainties on LFI and Planck products are also produced. Main beams are estimated down to the approx -10dB level using Jupiter transits, which are also used for geometrical calibration of the focal plane.

preprint2011arXiv

Planck Early Results. VII. The Early Release Compact Source Catalog

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

preprint2011arXiv

Planck Early Results. X. Statistical analysis of Sunyaev-Zeldovich scaling relations for X-ray galaxy clusters

All-sky data from the Planck survey and the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC) are combined to investigate the relationship between the thermal Sunyaev-Zeldovich (SZ) signal and X-ray luminosity. The sample comprises ~ 1600 X-ray clusters with redshifts up to ~ 1 and spans a wide range in X-ray luminosity. The SZ signal is extracted for each object individually, and the statistical significance of the measurement is maximised by averaging the SZ signal in bins of X-ray luminosity, total mass, or redshift. The SZ signal is detected at very high significance over more than two decades in X-ray luminosity (10^43 erg/s < L_500 E(z)^-7/3 < 2 X 10^45 erg/s). The relation between intrinsic SZ signal and X-ray luminosity is investigated and the measured SZ signal is compared to values predicted from X-ray data. Planck measurements and X-ray based predictions are found to be in excellent agreement over the whole explored luminosity range. No significant deviation from standard evolution of the scaling relations is detected. For the first time the intrinsic scatter in the scaling relation between SZ signal and X-ray luminosity is measured and found to be consistent with the one in the luminosity -- mass relation from X-ray studies. There is no evidence of any deficit in SZ signal strength in Planck data relative to expectations from the X-ray properties of clusters, underlining the robustness and consistency of our overall view of intra-cluster medium properties.

preprint2011arXiv

Planck early results. XV. Spectral energy distributions and radio continuum spectra of northern extragalactic radio sources

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

preprint2011arXiv

Planck early results. XX. New light on anomalous microwave emission from spinning dust grains

Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ~10-60 GHz. Using Planck maps and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and Rho Ophiuchi molecular clouds. The spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the first time. The spectra have a peak in the range 20-40 GHz and are detected at high significances of 17.1 sigma for Perseus and 8.4 sigma for Rho Ophiuchi. In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. In Rho Ophiuchi, the ~30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50-100 GHz, which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and thermal dust. We present spectra for two of the candidates; S140 and S235 are bright HII regions that show evidence for AME, and are well fitted by spinning dust models.

preprint2011arXiv

Planck early results. XXI. Properties of the interstellar medium in the Galactic plane

(abridged) Planck has observed the entire sky from 30 GHz to 857GHz. The observed foreground emission contains contributions from different phases of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular and ionised) in each of a number of Galactocentric rings. Templates are created for various Galactocentric radii using velocity information from atomic (neutral hydrogen) and molecular (12CO) observations. The ionised template is assumed to be traced by free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the above templates, namely "ark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as other ancillary data. The emission per column density of the gas templates allows us to create distinct spectral energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4 GHz to 25 THz (12μm). Apart from the thermal dust and free-free emission, we have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout the Galactic disk. The derived dust propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for $(25\pm5)%$ (statistical) of the total emission at 30 GHz.

preprint2011arXiv

Search and characterization of T-type planetary mass candidates in the sigma Orionis cluster

(Abridged) We aim to: i) confirm the presence of methane absorption in S Ori 73 (a T-type member candidate of the sig Orionis cluster, 3 Myr, 352 pc) through methane imaging; ii) study S Ori 70 and 73 cluster membership via photometric colors and accurate proper motion analysis; iii) perform a new search to identify additional T-type sig Orionis member candidates with likely masses below 7 Mjup. We obtained HAWK-I (VLT) J, H, and CH4off photometry of an area of 119.15 sq. arcmin in sig Orionis down to Jcomp = 21.7 and Hcomp = 21 mag. Near-infrared data were complemented with optical photometry using images acquired with OSIRIS (GTC) and VISTA as part of the VISTA Orion survey. We derived proper motions by comparison of the new HAWK-I and VISTA images with published near-infrared data taken 3.4 - 7.9 yr ago. S Ori 73 has a red H-CH4off color indicating methane absorption in the H-band and a spectral type of T4 +/- 1. S Ori 70 displays a redder methane color than S Ori 73 in agreement with its latter spectral classification. Our proper motion measurements are larger than the motion of sig Orionis, rendering S Ori 70 and 73 cluster membership uncertain. We identified one new photometric candidate with J = 21.69 +/- 0.12 mag and methane color consistent with spectral type greater than T8. S Ori 73 has colors similar to those of T3-T5 field dwarfs, which in addition to its high proper motion suggests that it is probably a field dwarf located at 170-200 pc. The origin of S Ori 70 remains unclear: it can be a field, foreground mid- to late-T free-floating dwarf with peculiar colors, or an orphan planet ejected through strong dynamical interactions from sig Orionis or from a nearby star-forming region in Orion.

preprint2011arXiv

The substellar population of Sigma Orionis: A deep wide survey

We present a deep I,Z photometric survey covering a total area of 1.12 deg^{2} of the Sigma Orionis cluster (Icompl=22 and Zcompl=21.5mag). From I, I-Z color-magnitude diagrams we have selected 153 candidates that fit the previously known sequence of the cluster. Using J-band photometry, we find that 124 of the 151 candidates follow the previously known infrared photometric sequence of the cluster and are probably members. We have studied the spatial distribution of these candidates and found that there are objects located at distances greater than 30 arcmin to the north and west of Sigma Orionis that probably belong to different populations of the Orion's Belt. For the 102 bona fide Sigma Orionis cluster member candidates, we find that the radial surface density can be represented by a decreasing exponential function (sigma = sigma_0 e^{-r/r_0}) with a central density of sigma_0=0.23+/-0.03 object/arcmin^{2} and a characteristic radius of r_0=9.5+/-0.7 arcmin. From a statistical comparison with Monte Carlo simulations, we conclude that the spatial distribution of the cluster member candidates is compatible with a Poissonian distribution and, hence, they are not mainly forming aggregations or sub-clustering. Using near-infrared JHK-band data from 2MASS and UKIDSS and mid-infrared data from IRAC/Spitzer, we find that 5-9 % of the brown dwarf candidates in the Sigma Orionis cluster have K-band excesses and 31+/-7 % of them show mid-infrared excesses at wavelengths longer than 5.8 microns, which are probably related to the presence of disks. We have also calculated the initial mass spectrum (dN/dm) of Sigma Orionis from very low mass stars (0.10 Msol) to the deuterium-burning mass limit (0.012-0.013 Msol). This is a rising function toward lower masses and can be represented by a power-law distribution (dN/dm = m^{-alpha}) with an exponent alpha of 0.7+/-0.3 for an age of 3 Myr.

preprint2010arXiv

A 33 GHz VSA survey of the Galactic plane from 27 to 46 degrees

The Very Small Array (VSA) has been used to survey the l = 27 to 46 deg, |b|<4 deg region of the Galactic plane at a resolution of 13 arcmin. The survey consists of 44 pointings of the VSA, each with a r.m.s. sensitivity of ~90 mJy/beam. These data are combined in a mosaic to produce a map of the area. The majority of the sources within the map are HII regions. We investigated anomalous radio emission from the warm dust in 9 HII regions of the survey by making spectra extending from GHz frequencies to the FIR IRAS frequencies. Acillary radio data at 1.4, 2.7, 4.85, 8.35, 10.55, 14.35 and 94 GHz in addition to the 100, 60, 25 and 12 micron IRAS bands were used to construct the spectra. From each spectrum the free-free, thermal dust and anomalous dust emission were determined for each HII region. The mean ratio of 33 GHz anomalous flux density to FIR 100 micron flux density for the 9 selected HII regions was 1.10 +/-0.21x10^(-4). When combined with 6 HII regions previously observed with the VSA and the CBI, the anomalous emission from warm dust in HII regions is detected with a 33 GHz emissivity of 4.65 +/- 0.4 micro K/ (MJy/sr) at 11.5σ. The anomalous radio emission in HII regions is on average 41+/-10 per cent of the radio continuum at 33 GHz.

preprint2010arXiv

A search for interstellar anthracene toward the Perseus anomalous microwave emission region

We report the discovery of a new broad interstellar (or circumstellar) band at 7088.8 +- 2.0 Åcoincident to within the measurement uncertainties with the strongest band of the anthracene cation (C$_{14}$H$_{10}$$^+$) as measured in gas-phase laboratory spectroscopy at low temperatures (Sukhorukov et al.2004). The band is detected in the line of sight of star Cernis 52, a likely member of the very young star cluster IC 348, and is probably associated with cold absorbing material in a intervening molecular cloud of the Perseus star forming region where various experiments have recently detected anomalous microwave emission. From the measured intensity and available oscillator strength we find a column density of N$_{an^+}$= 1.1(+-0.4) x 10$^{13}$ cm$^{-2}$ implying that ~0.008% of the carbon in the cloud could be in the form of C$_{14}$H$_{10}$$^+$. A similar abundance has been recently claimed for the naphthalene cation (Iglesias-Groth et al. 2008) in this cloud. This is the first location outside the Solar System where specific PAHs are identified. We report observations of interstellar lines of CH and CH$^+$ that support a rather high column density for these species and for molecular hydrogen. The strength ratio of the two prominent diffuse interstellar bands at 5780 and 5797 Åsuggests the presence of a ``zeta'' type cloud in the line of sight (consistent with steep far-UV extinction and high molecular content). The presence of PAH cations and other related hydrogenated carbon molecules which are likely to occur in this type of clouds reinforce the suggestion that electric dipole radiation from fast spinning PAHs is responsible of the anomalous microwave emission detected toward Perseus.

preprint2010arXiv

Constraints on the Polarization of the Anomalous Microwave Emission in the Perseus Molecular Complex from 7-year WMAP data

We have used the seven year Wilkinson Microwave Anisotropy Probe (WMAP) data in order to update the measurements of the intensity signal in the G159.6-18.5 region within the Perseus Molecular Complex, and to set constraints on the polarization level of the anomalous microwave emission in the frequency range where this emission is dominant. At 23, 33 and 41 GHz, we obtain upper limits on the fractional linear polarization of 1.0, 1.8 and 2.7% respectively (with a 95 per cent confidence level). These measurements rule out a significant number of models based on magnetic dipole emission of grains that consist of a simple domain (Draine & Lazarian 1999) as responsible of the anomalous emission. When combining our results with the measurement obtained with the COSMOSOMAS experiment at 11 GHz (Battistelli et al. 2006), we find consistency with the predictions of the electric dipole and resonance relaxation theory (Lazarian & Draine 2000) at this frequency range.

preprint2010arXiv

High contrast optical imaging of companions: the case of the brown dwarf binary HD-130948BC

High contrast imaging at optical wavelengths is limited by the modest correction of conventional near-IR optimized AO systems.We take advantage of new fast and low-readout-noise detectors to explore the potential of fast imaging coupled to post-processing techniques to detect faint companions to stars at small separations. We have focused on I-band direct imaging of the previously detected brown dwarf binary HD130948BC,attempting to spatially resolve the L2+L2 benchmark system. We used the Lucky-Imaging instrument FastCam at the 2.5-m Nordic Telescope to obtain quasi diffraction-limited images of HD130948 with ~0.1" resolution.In order to improve the detectability of the faint binary in the vicinity of a bright (I=5.19 \pm 0.03) solar-type star,we implemented a post-processing technique based on wavelet transform filtering of the image which allows us to strongly enhance the presence of point-like sources in regions where the primary halo dominates. We detect for the first time the BD binary HD130948BC in the optical band I with a SNR~9 at 2.561"\pm 0.007" (46.5 AU) from HD130948A and confirm in two independent dataset that the object is real,as opposed to time-varying residual speckles.We do not resolve the binary, which can be explained by astrometric results posterior to our observations that predict a separation below the NOT resolution.We reach at this distance a contrast of dI = 11.30 \pm 0.11, and estimate a combined magnitude for this binary to I = 16.49 \pm 0.11 and a I-J colour 3.29 \pm 0.13. At 1", we reach a detectability 10.5 mag fainter than the primary after image post-processing. We obtain on-sky validation of a technique based on speckle imaging and wavelet-transform processing,which improves the high contrast capabilities of speckle imaging.The I-J colour measured for the BD companion is slightly bluer, but still consistent with what typically found for L2 dwarfs(~3.4-3.6).

preprint2010arXiv

High spatial resolution and high contrast optical speckle imaging with FASTCAM at the ORM

In this paper, we present an original observational approach, which combines, for the first time, traditional speckle imaging with image post-processing to obtain in the optical domain diffraction-limited images with high contrast (1e-5) within 0.5 to 2 arcseconds around a bright star. The post-processing step is based on wavelet filtering an has analogy with edge enhancement and high-pass filtering. Our I-band on-sky results with the 2.5-m Nordic Telescope (NOT) and the lucky imaging instrument FASTCAM show that we are able to detect L-type brown dwarf companions around a solar-type star with a contrast DI~12 at 2" and with no use of any coronographic capability, which greatly simplifies the instrumental and hardware approach. This object has been detected from the ground in J and H bands so far only with AO-assisted 8-10 m class telescopes (Gemini, Keck), although more recently detected with small-class telescopes in the K band. Discussing the advantage and disadvantage of the optical regime for the detection of faint intrinsic fluxes close to bright stars, we develop some perspectives for other fields, including the study of dense cores in globular clusters. To the best of our knowledge this is the first time that high contrast considerations are included in optical speckle imaging approach.

preprint2010arXiv

Infrared and kinematic properties of the substellar object G 196-3B

We report unusual near- and mid-infrared photometric properties of G 196-3 B, the young substellar companion at 16 arcsec from the active M2.5-type star G 196-3 A, using data taken with the IRAC and MIPS instruments onboard Spitzer. G 196-3 B shows markedly redder colors at all wavelengths from 1.6 up to 24 micron than expected for its spectral type, which is determined at L3 from optical and near-infrared spectra. We discuss various physical scenarios to account for its reddish nature, and conclude that a low-gravity atmosphere with enshrouded upper atmospheric layers and/or a warm dusty disk/envelope provides the most likely explanations, the two of them consistent with an age in the interval 20-300 Myr. We also present new and accurate separate proper motion measurements for G 196-3 A and B confirming that both objects are gravitationally linked and share the same motion within a few mas/yr. After integration of the combined spectrophotometric spectral energy distributions, we obtain that the difference in the bolometric magnitudes of G 196-3 A and B is 6.15 +/- 0.10 mag. Kinematic consideration of the Galactic space motions of the system for distances in the interval 15-30 pc suggests that the pair is a likely member of the Local Association, and that it lay near the past positions of young star clusters like alpha Persei less than 85 Myr ago, where the binary might have originated. At these young ages, the mass of G 196-3 B would be in the range 12-25 Mjup, close to the frontier between planets and brown dwarfs.

preprint2010arXiv

Near-infrared low-resolution spectroscopy of Pleiades L-type brown dwarfs

The fundamental properties of brown dwarfs evolve with age. Models describing the evolution of luminosities and effective temperatures, among other physical parameters, can be empirically constrained using brown dwarfs of various masses in star clusters of well determined age and metallicity. We aim to carry out a spectroscopic and photometric characterization of low-mass brown dwarfs of the ~120 Myr old Pleiades open cluster. We obtained low-resolution near-infrared spectra of the J=17.4-18.8 mag candidate L-type brown dwarfs PLIZ 28 and 35, BRB 17, 21, 23, and 29, which are Pleiades members by photometry and proper motion. We also obtained spectra of the well-known J=15.4-16.1 mag late M-type cluster members PPl 1, Teide 1, and Calar 3. We find that the former six objects have early- to mid-L spectral types and confirm previously reported M-types for the other three objects. The spectra of the L0-type BRB 17 and PLIZ 28 present a triangular H-band continuum shape, indicating that this peculiar spectral feature persists until at least the age of the Pleiades. We add to our sample 36 reported M5-L0-type cluster members, collecting their I_c - and UKIDSS ZYJHK-band photometry. We confirm a possible interleaving of the Pleiades and field L-type sequences in the JHK absolute magnitude versus spectral type diagrams, and quantify marginally redder Pleiades J-K colours, by 0.12+-0.20 mag, possibly related to both reddening and youth. Using field dwarf bolometric correction - and effective temperature - spectral type relations, we obtain the Hertzsprung-Russell diagram of the Pleiades sample. Theoretical models reproduce well the spectral sequence at M5.5-9, but appear to overestimate the luminosity or underestimate the effective temperature at L0-5. For the Pleiades early- to mid L-type brown dwarfs, we estimate theoretical masses to be in the range 0.025-0.035 M_Sol.

preprint2010arXiv

Planck pre-launch status: Design and description of the Low Frequency Instrument

In this paper we present the Low Frequency Instrument (LFI), designed and developed as part of the Planck space mission, the ESA program dedicated to precision imaging of the cosmic microwave background (CMB). Planck-LFI will observe the full sky in intensity and polarisation in three frequency bands centred at 30, 44 and 70 GHz, while higher frequencies (100-850 GHz) will be covered by the HFI instrument. The LFI is an array of microwave radiometers based on state-of-the-art Indium Phosphide cryogenic HEMT amplifiers implemented in a differential system using blackbody loads as reference signals. The front-end is cooled to 20K for optimal sensitivity and the reference loads are cooled to 4K to minimise low frequency noise. We provide an overview of the LFI, discuss the leading scientific requirements and describe the design solutions adopted for the various hardware subsystems. The main drivers of the radiometric, optical and thermal design are discussed, including the stringent requirements on sensitivity, stability, and rejection of systematic effects. Further details on the key instrument units and the results of ground calibration are provided in a set of companion papers.

preprint2010arXiv

Planck pre-launch status: Low Frequency Instrument calibration and expected scientific performance

We give the calibration and scientific performance parameters of the Planck Low Frequency Instrument (LFI) measured during the ground cryogenic test campaign. These parameters characterise the instrument response and constitute our best pre-launch knowledge of the LFI scientific performance. The LFI shows excellent $1/f$ stability and rejection of instrumental systematic effects; measured noise performance shows that LFI is the most sensitive instrument of its kind. The set of measured calibration parameters will be updated during flight operations through the end of the mission.

preprint2010arXiv

Planck pre-launch status: the Planck-LFI programme

This paper provides an overview of the Low Frequency Instrument (LFI) programme within the ESA Planck mission. The LFI instrument has been developed to produce high precision maps of the microwave sky at frequencies in the range 27-77 GHz, below the peak of the cosmic microwave background (CMB) radiation spectrum. The scientific goals are described, ranging from fundamental cosmology to Galactic and extragalactic astrophysics. The instrument design and development are outlined, together with the model philosophy and testing strategy. The instrument is presented in the context of the Planck mission. The LFI approach to ground and inflight calibration is described. We also describe the LFI ground segment. We present the results of a number of tests demonstrating the capability of the LFI data processing centre (DPC) to properly reduce and analyse LFI flight data, from telemetry information to calibrated and cleaned time ordered data, sky maps at each frequency (in temperature and polarization), component emission maps (CMB and diffuse foregrounds), catalogs for various classes of sources (the Early Release Compact Source Catalogue and the Final Compact Source Catalogue). The organization of the LFI consortium is briefly presented as well as the role of the core team in data analysis and scientific exploitation. All tests carried out on the LFI flight model demonstrate the excellent performance of the instrument and its various subunits. The data analysis pipeline has been tested and its main steps verified. In the first three months after launch, the commissioning, calibration, performance, and verification phases will be completed, after which Planck will begin its operational life, in which LFI will have an integral part.

preprint2010arXiv

The Planck/LFI Radiometer Electronics Box Assembly

The Radiometer Electronics Box Assembly (REBA) is the control and data processing on board computer of the Low Frequency Instrument (LFI) of the Planck mission (ESA). The REBA was designed and built incorporating state of the art processors, communication interfaces and real time operating system software in order to meet the scientific performance of the LFI. We present a technical summary of the REBA, including a physical, functional, electrical, mechanical and thermal description. Aspects of the design and development, the assembly, the integration and the verification of the equipment are provided. A brief description of the LFI on board software is given including the Low-Level Software and the main functionalities and architecture of the Application Software. The compressor module, which has been developed as an independent product, later integrated in the application, is also described in this paper. Two identical engineering models EM and AVM, the engineering qualification model EQM, the flight model FM and flight spare have been manufactured and tested. Low-level and Application software have been developed. Verification activities demonstrated that the REBA hardware and software fulfil all the specifications and perform as required for flight operation.

preprint2009arXiv

A census of very-low-mass stars and brown dwarfs in the sigma Orionis cluster

(ABRIDGED) We have analysed the near-infrared photometric data from the Fourth Data Release (DR4) of the UKIRT Infrared Deep Sky Suvey (UKIDSS) Galactic Clusters Survey (GCS) to derive the cluster luminosity and mass functions, evaluate the extent of the cluster, and study the distribution and variability of low-mass stars and brown dwarfs down to the deuterium-burning limit. We have recovered most of the previously published members and found a total of 287 candidate members within the central 30 arcmin in the 0.5-0.009 Msun mass range, including new objects not previously reported in the literature. This new catalogue represents a homogeneous dataset of brown dwarf member candidates over the central 30 arcmin of the cluster. The expected photometric contamination by field objects with similar magnitudes and colours to sigma Orionis members is ~15%. We present evidence of variability at the 99.5% confidence level over ~yearly timescales in 10 member candidates that exhibit signs of youth and the presence of disks. The level of variability is low (<0.3 mag) and does not impact the derivation of the cluster luminosity and mass functions. Furthermore, we find a possible dearth of brown dwarfs within the central five arcmin of the cluster, which is not caused by a lower level of photometric sensitivity around the massive, O-type multiple star sigma Ori in the GCS survey. Using state-of-the-art theoretical models, we derived the luminosity and mass functions within the central 30 arcmin from the cluster centre, with completeness down to J = 19 mag, corresponding to masses ranging from 0.5 Msun down to the deuterium-burning mass boundary (~0.013 Msun). The mass function of sigma Orionis in this mass interval shows a power law index alpha = 0.5+/-0.2.

preprint2009arXiv

Candidate free-floating super-Jupiters in the young sigma Orionis open cluster

Free-floating substellar candidates with estimated theoretical masses of as low as ~5 Jupiter masses have been found in the ~3 Myr old sigma Orionis open cluster. As the overlap with the planetary mass domain increases, the question of how these objects form becomes important. The determination of their number density and whether a mass cut-off limit exists is crucial to understanding their formation. We propose to search for objects of yet lower masses in the cluster and determine the shape of the mass function at low mass. Using new- and (re-analysed) published IZJHKs[3.6]-[8.0]-band data of an area of 840 arcmin2, we performed a search for LT-type cluster member candidates in the magnitude range J=19.5-21.5 mag, based on their expected magnitudes and colours. Besides recovering the T type object S Ori 70 and two other known objects, we find three new cluster member candidates, S Ori 72-74, with J=21 mag and within 12 arcmin of the cluster centre. They have theoretical masses of 4 (-2,+3) M_Jup and are among the least massive free-floating objects detected by direct imaging outside the Solar System. The photometry in archival Spitzer [3.6]-[5.8]-band images infers that S Ori 72 is an L/T transition candidate and S Ori 73 a T-type candidate, following the expected cluster sequence in the mid-infrared. Finally, the L-type candidate S Ori 74 with lower quality photometry is located at 11.8 arcsec (~4250 AU) of a stellar member of sigma Orionis and could be a companion. After contaminant correction in the area complete to J=21.1 mag, we estimate that there remain between zero and two cluster members in the mass interval 6-4 M_Jup. Our result suggests a possible turnover in the substellar mass spectrum below ~6 Jupiter masses, which could be investigated further by wider and deeper photometric surveys.

preprint2009arXiv

On the mass of the neutron star in Cyg X-2

We present new high resolution spectroscopy of the low mass X-ray binary Cyg X-2 which enables us to refine the orbital solution and rotational broadening of the donor star. In contrast with Elebert et al (2009) we find a good agreement with results reported in Casares et al. (1998). We measure $P=9.84450\pm0.00019$ day, $K_2=86.5\pm1.2$ km s$^{-1}$ and $V \sin i=33.7\pm0.9$ km s$^{-1}$. These values imply $q=M_{2}/M_{1}=0.34 \pm 0.02$ and $M_{1}=1.71\pm 0.21$ M$_{\odot}$ (for $i=62.5 \pm 4^{\circ}$). Therefore, the neutron star in Cyg X-2 can be more massive than canonical. We also find no evidence for irradiation effects in our radial velocity curve which could explain the discrepancy between Elebert et al's and our $K_2$ values.

preprint2009arXiv

The Chemical Composition of Cernis 52 (BD+31 640)

We present an abundance analysis of the star Cernis 52 in whose spectrum we recently reported the napthalene cation in absorption at 6707.4 Å. This star is on a line of sight to the Perseus molecular complex. The analysis of high-resolution spectra using a chi^2-minimization procedure and a grid of synthetic spectra provides the stellar parameters and the abundances of O, Mg, Si, S, Ca, and Fe. The stellar parameters of this star are found to be T_{eff} = 8350 +- 200 K, logg= 4.2 +- 0.4 dex. We derived a metallicity of [Fe/H] = -0.01 +- 0.15. These stellar parameters are consistent with a star of $\sim 2$ \Msun in a pre-main-sequence evolutionary stage. The stellar spectrum is significantly veiled in the spectral range 5150-6730 Å up to almost 55 per cent of the total flux at 5150 Å and decreasing towards longer wavelengths. Using Johnson-Cousins and 2MASS photometric data, we determine a distance to Cernis 52 of 231$^{+135}_{-85}$ pc considering the error bars of the stellar parameters. This determination places the star at a similar distance to the young cluster IC 348. This together with its radial velocity, v_r=13.7+-1 km/s, its proper motion and probable young age support Cernis 52 as a likely member of IC 348. We determine a rotational velocity of v\sin i=65 +- 5 km/s for this star. We confirm that the stellar resonance line of \ion{Li}{1} at 6707.8 Å is unable to fit the broad feature at 6707.4 Å. This feature should have a interstellar origin and could possibly form in the dark cloud L1470 surrounding all the cluster IC 348 at about the same distance.

preprint2009arXiv

The spatial distribution of galaxies within the CMB cold spot in the Corona Borealis supercluster

We study the spatial distribution and colours of galaxies within the region covered by the cold spot in the cosmic microwave background (CMB) recently detected by the Very Small Array (VSA; Genova-Santos et al. 2005, 2008) towards the Corona Borealis supercluster (CrB-SC). The spot is in the northern part of a region with a radius ~1 degree (~5 Mpc at the redshift of CrB-SC) enclosing the clusters Abell 2056, 2065, 2059 and 2073, and where the density of galaxies, excluding the contribution from those clusters, is ~2 times higher than the mean value in typical intercluster regions of the CrB-SC. Two of such clusters (Abell 2056 and 2065) are members of the CrB-SC, while the other two are in the background. This high density intercluster region is quite inhomogeneous, being the most remarkable feature a large concentration of galaxies in a narrow filament running from Abell 2065 with a length of ~35 arcmin (~3 Mpc at the redshift of CrB-SC) in the SW-NE direction. This intercluster population of galaxies probably results from the interaction of clusters Abell 2065 and 2056. The area subtended by the VSA cold spot shows an excess of faint (21<r<22) and red (1.1<r-i<1.3) galaxies as compared with typical values within the CrB-SC intercluster regions. This overdensity of galaxies shows a radial dependence and extends out to ~15 arcmin. This could be signature of a previously unnoticed cluster in the background.

preprint2009arXiv

The Sunyaev-Zeldovich effect in superclusters of galaxies using gasdynamical simulations: the case of Corona Borealis

[Abridged] We study the thermal and kinetic Sunyaev-Zel'dovich (SZ) effect associated with superclusters of galaxies using the MareNostrum Universe SPH simulation. We consider superclusters similar to the Corona Borealis Supercluster (CrB-SC). This paper is motivated by the detection at 33GHz of a strong temperature decrement in the CMB towards the core of this supercluster. Multifrequency observations with VSA and MITO suggest the existence of a thermal SZ effect component in the spectrum of this cold spot, which would account for roughly 25% of the total observed decrement. We identify nine regions containing superclusters similar to CrB-SC, obtain the associated SZ maps and calculate the probability of finding such SZ signals arising from hot gas within the supercluster. Our results show that WHIM produces a thermal SZ effect much smaller than the observed value. Neither can summing the contribution of small clusters and galaxy groups in the region explain the amplitude of the SZ signal. When we take into account the actual posterior distribution from the observations, the probability that WHIM can cause a thermal SZ signal like the one observed is <1%, rising up to a 3.2% when the contribution of small clusters and galaxy groups is included. If the simulations provide a suitable description of the gas physics, then we conclude that the thermal SZ component of the CrB spot most probably arises from an unknown galaxy cluster along the line of sight. The simulations also show that the kinetic SZ signal associated with the supercluster cannot provide an explanation for the remaining 75% of the observed cold spot in CrB.

preprint2008arXiv

Chemical abundances of late-type pre-main sequence stars in the $σ$-Orionis cluster

The young $σ$-Orionis cluster is an important location for understanding the formation and evolution of stars, brown dwarfs, and planetary-mass objects. Its metallicity, although being a fundamental parameter, has not been well determined yet. We present the first determination of the metallicity of nine young late-type stars in $σ$-Orionis. Using the optical and near-infrared broadband photometry available in the literature we derive the effective temperatures for these nine cluster stars, which lie in the interval 4300--6500 K (1--3 \Msuno). These parameters are employed to compute a grid of synthetic spectra based on the code MOOG and Kurucz model atmospheres. We employ a $χ^2$-minimization procedure to derive the stellar surface gravity and atmospheric abundances of Al, Ca, Si, Fe, Ni and Li, using multi-object optical spectroscopy taken with WYFFOS+AF2 at at the William Herschel Telescope ($λ/δλ\sim7500$). The average metallicity of the $σ$-Orionis cluster is [Fe/H] $ = -0.02\pm0.09\pm0.13$ (random and systematic errors). The abundances of the other elements, except lithium, seem to be consistent with solar values. Lithium abundances are in agreement with the "cosmic" $^7$Li abundance, except for two stars which show a $\log ε(\mathrm{Li})$ in the range 3.6--3.7 (although almost consistent within the error bars). There are also other two stars with $\log ε(\mathrm{Li})\sim 2.75$. We derived an average radial velocity of the $σ$-Orionis cluster of $28\pm4$km/s. The $σ$-Orionis metallicity is roughly solar.

preprint2008arXiv

The Quijote CMB Experiment

We present the current status of the QUIJOTE (Q-U-I JOint TEnerife) CMB Experiment, a new instrument which will start operations early 2009 at Teide Observatory, with the aim of characterizing the polarization of the CMB and other processes of galactic and extragalactic emission in the frequency range 10-30 GHz and at large angular scales. QUIJOTE will be a valuable complement at low frequencies for the PLANCK mission, and will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

preprint2007arXiv

The CODEX-ESPRESSO experiment: cosmic dynamics, fundamental physics, planets and much more..

CODEX, a high resolution, super-stable spectrograph to be fed by the E-ELT, the most powerful telescope ever conceived, will for the first time provide the possibility of directly measuring the change of the expansion rate of the Universe with time and much more, from the variability of fundamental constants to the search for other earths. A study for the implementation at the VLT of a precursor of CODEX, dubbed ESPRESSO, is presently carried out by a collaboration including ESO, IAC, INAF, IoA Cambridge and Observatoire de Geneve. The present talk is focused on the cosmological aspects of the experiment.

preprint2007arXiv

The substellar mass function in sigma Orionis. II. Optical, near-infrared and IRAC/Spitzer photometry of young cluster brown dwarfs and planetary-mass objects

We investigate the mass function in the substellar domain down to a few Jupiter masses in the young sigma Orionis open cluster (3+/-2 Ma, d = 360^+70_-60 pc). We have performed a deep IJ-band search, covering an area of 790 arcmin^2 close to the cluster centre. This survey was complemented with an infrared follow-up in the HKs- and Spitzer 3.6-8.0 mum-bands. Using colour-magnitude diagrams, we have selected 49 candidate cluster members in the magnitude interval 16.1 mag < I < 23.0 mag. Accounting for flux excesses at 8.0 mum and previously known spectral features of youth, 30 objects are bona fide cluster members. Four are first identified from our optical-near infrared data. Eleven have most probable masses below the deuterium burning limit and are classified as planetary-mass object candidates. The slope of the substellar mass spectrum (Delta N / Delta M = a M^-alpha) in the mass interval 0.11 Msol M < 0.006 Msol is alpha = +0.6+/-0.2. Any opacity mass-limit, if these objects form via fragmentation, may lie below 0.006 Msol. The frequency of sigma Orionis brown dwarfs with circumsubstellar discs is 47+/-15 %. The continuity in the mass function and in the frequency of discs suggests that very low-mass stars and substellar objects, even below the deuterium-burning mass limit, may share the same formation mechanism.

preprint2005arXiv

S Ori J053825.4-024241: A Classical T Tauri-like object at the substellar boundary

We present a spectrophotometric analysis of S Ori J053825.4-024241, a candidate member close to the substellar boundary of the young (1-8 Myr), nearby (~360 pc) sigma Orionis star cluster. Our optical and near-infrared photometry and low-resolution spectroscopy indicate that S Ori J053825.4-024241 is a likely cluster member with a mass estimated from evolutionary models at 0.06+0.07-0.02 Msol, which makes the object a probable brown dwarf. The radial velocity of S Ori J053825.4-024241 is similar to the cluster systemic velocity. This target, which we have classified as an M 6.0+-1.0 low-gravity object, shows excessemission in the near-infrared and anomalously strong photometric variability for its type (from the blue to the J band), suggesting the presence of a surrounding disc. The optical spectroscopic observations show a continuum excess at short wavelengths and a persistent and resolved Halpha emission (pseudo-equivalent width of ~-250 AA) in addition to the presence of other forbidden and permitted emission lines, which we interpret as indicating accretion from the disc and possibly mass loss. We conclude that despite the low mass of S Ori J053825.4-024241, this object exhibits some of the properties typical of active classical T Tauri stars.

preprint2004arXiv

Optical and infrared photometry of new very low-mass stars and brown dwarfs in the Sigma Orionis cluster

We present an RI photometric survey covering an area of 430 arcmin^2 around the multiple star Sigma Orionis. The observations were conducted with the 0.8 m IAC-80 Telescope at the Teide Observatory. The survey limiting R and I magnitudes are 22.5 and 21, and completeness magnitudes 21 and 20, respectively. We have selected 53 candidates from the I vs. R-I colour-magnitude diagram (I=14-20) that follow the previously known photometric sequence of the cluster. Adopting an age of 2-4 Myr for the cluster, we find that these objects span a mass range from 0.35 Msol to 0.015 Msol. We have performed J-band photometry of 52 candidates and Ks photometry for 12 of them, with the result that 50 follow the expected infrared sequence for the cluster, thus confirming with great confidence that the majority of the candidates are bona fide members. JHKs photometry from the Two Micron All Sky Survey (2MASS) is available for 50 of the candidates and are in good agreement with our data. Out of 48 candidates, which have photometric accuracies better than 0.1 mag in all bands, only three appear to show near-infrared excesses.

preprint2001arXiv

Optical spectroscopy of isolated planetary mass objects in the sigma Orionis cluster

We have obtained low resolution optical spectra of 15 isolated planetary mass objects (IPMOs) in the sigma Orionis cluster, and derived spectral types by comparison with nearby M and L dwarfs. The spectral types are in the range late M - mid L, in agreement with our expectations based on colors and magnitudes for bona fide members. Therefore, most of these objects have masses below the deuterium burning limit. About 2/3 show H(alpha) in emission at our spectral resolution. From our spectroscopic and photometric data, we infer that three IPMOs in this sample may be binaries with components of similar masses. These results confirm that the substellar mass function of the sigma Orionis cluster, in the form dN/dM, keeps rising in the planetary domain.

preprint2001arXiv

Oxygen in the Very Early Galaxy

Oxygen abundances in a sample of ultra-metal-poor subdwarfs have been derived from measurements of the oxygen triplet at 7771--5 A and OH lines in the near UV performed in high-resolution and high signal-to-noise spectra obtained with WHT/UES, KeckI/HIRES, and VLT/UVES. Our Fe abundances were derived in LTE and then corrected for NLTE effects following Thevenin and Idiart (1999). The new oxygen abundances confirm previous findings for a progressive linear rise in the oxygen-to-iron ratio with a slope -0.33+-0.02 from solar metallicity to [Fe/H] -3. A slightly higher slope would be obtained if the Fe NLTE corrections were not considered. Below [Fe/H]= -2.5 our stars show [O/Fe] ratios as high as ~ 1.17 (G64-12), which can be interpreted as evidence for oxygen overproduction in the very early epoch of the formation of the halo, possibly associated with supernova events with very massive progenitor stars. We show that the arguments against this linear trend given by Fulbright and Kraft (1999), based on the LTE Fe analysis of two metal-poor stars cannot be sustained when an NLTE analysis is performed. Using 1-D models our analysis of three oxygen indicators available for BD +23 3130 gives consistent abundances within 0.16 dex and average [O/Fe] ratio of 0.91.

preprint2001arXiv

The Substellar Mass Function in sigma Orionis

We combine results from imaging searches for substellar objects in the sigma Orionis cluster and follow-up photometric and spectroscopic observations to derive a census of the brown dwarf population in a region of 847 arcmin^2. We identify 64 very low-mass cluster member candidates in this region. We have available three color (IZJ) photometry for all of them, spectra for 9 objects, and K photometry for 27% of our sample. These data provide a well defined sequence in the I vs I-J, I-K color magnitude diagrams, and indicate that the cluster is affected by little reddening despite its young age (~5 Myr). Using state-of-the-art evolutionary models, we derive a mass function from the low-mass stars (0.2 Msol) across the complete brown dwarf domain (0.075 Msol to 0.013 Msol), and into the realm of free-floating planetary-mass objects (<0.013 Msol). We find that the mass spectrum (dN/dm ~ m^{-alpha}) increases toward lower masses with an exponent alpha = 0.8+/-0.4. Our results suggest that planetary-mass isolated objects could be as common as brown dwarfs; both kinds of objects together would be as numerous as stars in the cluster. If the distribution of stellar and substellar masses in sigma Orionis is representative of the Galactic disk, older and much lower luminosity free-floating planetary-mass objects with masses down to about 0.005 Msol should be abundant in the solar vicinity, with a density similar to M-type stars.

preprint1997arXiv

Studies of CMB structure at Dec=40. II: Analysis and cosmological interpretation

We present a detailed analysis of the cosmic microwave background structure in the Tenerife Dec=+40 degrees data. The effect of local atmospheric contributions on the derived fluctuation amplitude is considered, resulting in an improved separation of the intrinsic CMB signal from noise. Our analysis demonstrates the existence of common structure in independent data scans at 15 and 33 GHz. For the case of fluctuations described by a Gaussian auto-correlation function, a likelihood analysis of our combined results at 15 and 33 GHz implies an intrinsic rms fluctuation level of 48^{+21}_{-15} uK on a coherence scale of 4 degrees; the equivalent analysis for a Harrison-Zel'dovitch model gives a power spectrum normalisation of Q_{rms-ps} = 22^{+10}_{-6} uK. The fluctuation amplitude is seen to be consistent at the 68% confidence level with that reported for the COBE two-year data for primordial fluctuations described by a power law model with a spectral index in the range 1.0 \le n \le 1.6. This limit favours the large scale CMB anisotropy being dominated by scalar fluctuations rather than tensor modes from a gravitational wave background. The large scale Tenerife and COBE results are considered in conjunction with observational results from medium scale experiments in order to place improved limits on the fluctuation spectral index; we find n=1.10 +/- 0.10 assuming standard CDM with H_{0}=50 kms^{-1}Mpc^{-1}.

preprint1996arXiv

Spectroscopy of New Substellar Candidates in the Pleiades: Towards a Spectral Sequence for Young Brown Dwarfs

We present optical and near-IR spectroscopy (600--1000 nm) of 8 faint (I$>$18) very red (R--I$>$2.2) objects discovered in a deep CCD survey of the Pleiades (Zapatero-Osorio et al. 1996). We compare them with reliable cluster members like PPl 15 and Teide 1, and with several field very late-type dwarfs (M4--M9.5), which were observed with similar instrumental configurations. Using pseudocontinuum ratios we classify the new substellar candidates in a spectal sequence defined with reference to field stars of known spectral types. We also reclassify PPl 15 and Teide 1 in a self-consistent way. The likelihood of membership for the new candidates is assesed via the study of their photospheric features, Ha emission, radial velocity, and consistency of their spectral types and I magnitudes with known cluster members. Four of the new substellar candidates are as late or later than PPl~15 (M6.5), but only one, Calar 3 (M8), clearly meets all our membership criteria. It is indeed an object very similar to the brown dwarf Teide 1. Calar 3, together with Teide 1, allows one to compare the spectroscopic characteristics of Pleiades brown dwarfs with those of old very cool dwarfs. The overall spectral properties are similar, but there are slight differences in the NaI doublet (818.3 nm, 819.5 nm), VO molecular band (740 nm), and some spectral ratios, which are probably related to lower surface gravity in the young Pleiades brown dwarfs than in field stars. Finally, we propose a way of improving future CCD-based brown dwarf surveys by using narrow-band near-IR pseudocontinuum filters.

preprint1994arXiv

Constraints to the Masses of Brown Dwarf Candidates from the Lithium Test

We present intermediate dispersion (0.7-2.2 Å~pix$^{-1}$) optical spectroscopic observations aimed at applying the ``Lithium Test'' to a sample of ten brown dwarf candidates located in the general field, two in young open clusters, and two in close binaries. We find evidence for strong Li depletion in all of them, and thus infer lower mass limits of 0.065~M$_\odot$, depending only slightly ($\pm$0.005~M$_\odot$) on the interior models. None of the field brown dwarf candidates in our sample appears to be a very young (age $<$~10$^8$~yr) substellar object. For one of the faintest proper motion Pleiades members known (V=20.7) the Li test implies a mass greater than $\sim$0.08~M$_\odot$, and therefore it is not a brown dwarf. From our spectra we estimate spectral types for some objects and present measurements of Halpha emission strengths and radial velocities. Finally, we compare the positions in the H-R diagram of our sample of brown dwarf candidates with the theoretical region where Li is expected to be preserved (Substellar Lithium Region). We find that certain combinations of temperature calibrations and evolutionary tracks are consistent with the constraints imposed by the observed Li depletion in brown dwarf candidates, while others are not.