Source author record

G. Micela

G. Micela 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

121works
6topics
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

121 published item(s)

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

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

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

Investigating the architecture and internal structure of the TOI-561 system planets with CHEOPS, HARPS-N and TESS

We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of $62$ HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b ($P = 0.45$ d, $R = 1.42$ R$_\oplus$, $M = 2.0$ M$_\oplus$), c ($P = 10.78$ d, $R = 2.91$ R$_\oplus$, $M = 5.4$ M$_\oplus$), d ($P = 25.7$ d, $R = 2.82$ R$_\oplus$, $M = 13.2$ M$_\oplus$) and e ($P = 77$ d, $R = 2.55$ R$_\oplus$, $M = 12.6$ M$_\oplus$). Moreover, we identify an additional, long-period signal ($>450$ d) in the RVs, which could be due to either an external planetary companion or to stellar magnetic activity. The precise masses and radii obtained for the four planets allowed us to conduct interior structure and atmospheric escape modelling. TOI-561 b is confirmed to be the lowest density ($ρ_{\rm b} = 3.8 \pm 0.5$ g cm$^{-3}$) ultra-short period (USP) planet known to date, and the low metallicity of the host star makes it consistent with the general bulk density-stellar metallicity trend. According to our interior structure modelling, planet b has basically no gas envelope, and it could host a certain amount of water. In contrast, TOI-561 c, d, and e likely retained an H/He envelope, in addition to a possibly large water layer. The inferred planetary compositions suggest different atmospheric evolutionary paths, with planets b and c having experienced significant gas loss, and planets d and e showing an atmospheric content consistent with the original one. The uniqueness of the USP planet, the presence of the long-period planet TOI-561 e, and the complex architecture make this system an appealing target for follow-up studies.

preprint2022arXiv

Low mass young stars in the Milky Way unveiled by DBSCAN and Gaia EDR3. Mapping the star forming regions within 1.5 Kpc

With an unprecedented astrometric and photometric data precision, Gaia EDR3 gives us, for the first time, the opportunity to systematically detect and map in the optical bands, the low mass populations of the star forming regions (SFRs) in the Milky Way. We provide a catalogue of the Gaia EDR3 data (photometry, proper motions and parallaxes) of the young stellar objects (YSOs) identified in the Galactic Plane (|b|<30 deg) within about 1.5 kpc. The catalogue of the SFRs to which they belong is also provided to study the properties of the very young clusters and put them in the context of the Galaxy structure. We applied the machine learning unsupervised clustering algorithm DBSCAN on a sample of Gaia EDR3 data photometrically selected on the region where very young stars (t<10 Myr) are expected to be found, with the aim to identify co-moving and spatially consistent stellar clusters. A subsample of 52 clusters, selected among the 7323 found with DBSCAN, has been used as template data set, to identify very young clusters from the pattern of the observed color-absolute magnitude diagrams through a pattern match process. We find 124440 candidate YSOs clustered in 354 SFRs and stellar clusters younger than 10 Myr and within about 1.5 Kpc. In addition, 65863 low mass members of 322 stellar clusters located within about 500 pc and with ages 10 Myr<t<100 Myr were also found. The selected YSOs are spatially correlated with the well known SFRs. Most of them are associated with well concentrated regions or complex structures of the Galaxy and a substantial number of them have been recognized for the first time. The massive SFRs, such as, for example, Orion, Sco-Cen and Vela, located within 600-700 pc trace a very complex three-dimensional pattern, while the farthest ones seem to follow a more regular pattern along the Galactic Plane.

preprint2022arXiv

The GAPS Programme at TNG XXXIII. HARPS-N detects multiple atomic species in emission from the dayside of KELT-20b

The detection of lines in emission in planetary atmospheres provides direct evidence of temperature inversion. We confirm the trend of ultra-hot Jupiters orbiting A-type stars showing temperature inversions on their daysides, by detecting metals emission lines in the dayside of KELT-20b. We first detect the planetary emission by using the G2 stellar mask of the HARPS-N pipeline, which is mainly composed of neutral iron lines, as a template. Using neutral iron templates, we perform a retrieval of the atmospheric temperature-pressure profile of the planet, confirming a thermal inversion. Then we create models of planetary emission of different species using the retrieved inverted temperature-pressure profile. By using the cross-correlation technique, we detect FeI, FeII and CrI at signal-to-noise ratio levels of 7.1, 3.9 and 3.6, respectively. The latter is detected for the first time in emission in the atmosphere of an exoplanet. Contrary to FeI, FeII and CrI are detected only after the occultation and not before, hinting for different atmospheric properties in view on the pre- and post- occultation orbital phases. A further retrieval of the temperature-pressure profile performed independently on the pre- and post- occultation phases, while not highly significant, points to a steeper thermal inversion in the post-occultation.

preprint2022arXiv

The GAPS programme at TNG XXXIV. Activity-rotation, flux-flux relationships, and active region evolution through stellar age

Active region evolution plays an important role in the generation and variability of magnetic fields on the surface of lower main-sequence stars. However, determining the lifetime of active region growth and decay as well as their evolution is a complex task. We aim to test whether the lifetime for active region evolution shows any dependency on the stellar parameters. We identify a sample of stars with well-defined ages via their kinematics. We made use of high-resolution spectra to compute rotational velocities, activity levels, and emission excesses. We use these data to revisit the activity-rotation-age relationship. The time-series of the main optical activity indicators were analysed together with the available photometry by using Gaussian processes to model the stellar activity of these stars. Autocorrelation functions of the available photometry were also analysed. We use the derived lifetimes for active region evolution to search for correlations with the stellar age, the spectral type, and the level of activity. We also use the pooled variance technique to characterise the activity behaviour of our targets. Our analysis confirms the decline of activity and rotation as the star ages. We also confirm that the rotation rate decays with age more slowly for cooler stars and that, for a given age, cooler stars show higher levels of activity. We show that F- and G-type young stars also depart from the inactive stars in the flux-flux relationship. The gaussian process analysis of the different activity indicators does not seem to provide any useful information on active region's lifetime and evolution. On the other hand, active region's lifetimes derived from the light-curve analysis might correlate with the stellar age and temperature. Although we caution the small number statistics, our results suggest that active regions seem to live longer on younger, cooler, and more active stars.

preprint2022arXiv

The GAPS Programme at TNG: XXXVI. Measurement of the Rossiter-McLaughlin effect and revising the physical and orbital parameters of the HAT-P-15, HAT-P-17, HAT-P-21, HAT-P-26, HAT-P-29 eccentric planetary systems

Aim: We aim to refine the orbital and physical parameters and determine the sky-projected planet orbital obliquity of five eccentric transiting planetary systems: HAT-P-15, HAT-P-17, HAT-P-21, HAT-P-26, and HAT-P-29. Each of the systems hosts a hot Jupiter, except for HAT-P-26 which hosts a Neptune-mass planet. Methods: We observed transit events of these planets with the HARPS-N spectrograph, obtaining high-precision radial velocity measurements that allow us to measure the Rossiter-McLaughlin effect for each of the target systems. We used these new HARPS-N spectra and archival data, including those from Gaia, to better characterise the stellar atmospheric parameters. The photometric parameters for four of the hot Jupiters were recalculated using 17 new transit light curves, obtained with an array of medium-class telescopes, and data from the TESS space telescope. HATNet time-series photometric data were checked for the signatures of rotation periods of the target stars and their spin axis inclination. Results: From the analysis of the Rossiter-McLaughlin effect, we derived a sky-projected obliquity of 13, -26.3, -0.7, -26 degree for HAT-P-15b, HAT-P-17b, HAT-P-21b and HAT-P-29b, respectively. Due to the quality of the data, we were not able to well constrain the sky-projected obliquity for HAT-P-26b, although a prograde orbit is favoured. The stellar activity of HAT-P-21 indicates a rotation period of 15.88 days, which allowed us to determine its true misalignment angle (25 degree). Our new analysis of the physical parameters of the five exoplanetary systems returned values compatible with those existing in the literature. Using TESS and the available transit light curves, we reviewed the orbital ephemeris for the five systems and confirmed that the HAT-P-26 system shows transit timing variations, which may tentatively be attributed to the presence of a third body.

preprint2022arXiv

The GAPS Programme with HARPS-N at TNG. XXXV. Fundamental properties of transiting exoplanet host stars

Exoplanetary properties depend on stellar properties: to know the planet with accuracy and precision it is necessary to know the star as accurately and precisely as possible. Our immediate aim is to characterize in a homogeneous and accurate way a sample of 27 transiting planet-hosting stars observed within the GAPS program. We determined stellar parameters (effective temperature, surface gravity, rotational velocity) and abundances of 26 elements (Li,C,N,O,Na,Mg,Al,Si,S,Ca,Sc,Ti,V,Cr,Fe,Mn,Co,Ni,Cu,Zn,Y,Zr,Ba,La,Nd,Eu). Our study is based on high-resolution HARPS-N@TNG and FEROS@ESO spectra and uniform techniques. We derived kinematic properties from Gaia data and estimated for the first time in exoplanet host stars ages using elemental ratios as chemical clocks. Teff of our stars is of 4400-6700 K, while [Fe/H] is within -0.3 and 0.4 dex. Lithium is present in 7 stars. [X/H] and [X/Fe] abundances vs [Fe/H] are consistent with the Galactic Chemical Evolution. The dependence of [X/Fe] with the condensation temperature is critically analyzed with respect to stellar and kinematic properties. All targets with measured C and O abundances show C/O<0.8, compatible with Si present in rock-forming minerals. Most of targets show 1.0<Mg/Si<1.5, compatible with Mg distributed between olivine and pyroxene. HAT-P-26, the target hosting the lowest-mass planet, shows the highest Mg/Si ratio. From our chemo-dinamical analysis we find agreement between ages and position within the Galactic disk. We note a tendency for higher density planets to be around metal-rich stars and hints of higher stellar abundances of some volatiles for lower mass planets. We cannot exclude that part of our results could be also related to the location of the stars within the Galactic disk. We trace the planetary migration scenario from the composition of the planets related to the chemical composition of the hosting stars

preprint2022arXiv

The GAPS Programme with HARPS-N at TNG. XXXVII. A precise density measurement of the young ultra-short period planet TOI-1807 b

Great strides have been made in recent years in the understanding of the mechanisms involved in the formation and evolution of planetary systems; despite this, many observational facts still do not have an explanation. A great contribution to the study of planetary formation processes comes from the study of young, low-mass planets, with short orbital periods. In the last years, the TESS satellite has identified many planets of this kind, and their characterization is mandatory to understand how they formed and evolved. Within the framework of the GAPS project, we performed the validation and characterization of the ultra-short period planet (USPP) TOI-1807b, orbiting its young host star BD+39 2643 (~300 Myr) in only 13 hours. This is the youngest USPP discovered so far. Thanks to a joint modeling of the stellar activity and planetary signals in the TESS light curve and in HARPS-N radial-velocity measurements, combined with accurate estimation of stellar parameters, we validated the planetary nature of TOI-1807b and measured its orbital and physical parameters. By using astrometric, photometric, and spectroscopic observations we found that BD+39 2643 is a young, active K dwarf star, member of a 300+/-80 Myr old moving group and that it rotates in Prot=8.8+/-0.1 days. This star hosts an USPP with an orbital period of only P_b=0.54937+/-0.00001 d. Thanks to the exquisite photometric and spectroscopic series, and the accurate information on the stellar activity, we measured both the radius and the mass of TOI-1807b with high precision, obtaining R_b=1.37+/-0.09 R_Earth and M_b=2.57+/-0.50 M_Earth. These planet parameters correspond to a rocky planet with an Earth-like density and no extended H/He envelope. From the analysis of the age-R_P distribution for planets with well measured ages, we inferred that TOI-1807b may have already lost a large part of its atmosphere during its 300 Myr life.

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.

preprint2022arXiv

X-ray variability of HD 189733 across eight years of XMM-Newton observations

The characterization of exoplanets, their formation, evolution, and chemical changes is tightly linked to our knowledge of their host stars. In particular, stellar X-rays and UV emission have a strong impact on the dynamical and chemical evolution of planetary atmospheres. We analyzed 25 XMM-Newton observations encompassing about eight years and totaling about 958 ks in order to study the X-ray emission of HD 189733 A. We find that the corona of HD 189733 A has an average temperature of 0.4 keV and it is only during flares that the mean temperature increases to 0.9 keV. Apart from the flares, there is no significant change in the flux and hardness of the coronal emission on a timescale of several months to years. Thus, we conclude that there is no detectable activity cycle on such timescales. We identified the flares and built their energy distribution. The number of flares observed around the phases of the planetary eclipses is not statistically different from the number of flares during transit phases. However, we do find a hint of a difference in the flare-energy distributions, as the flares observed around the planetary eclipses tend to be more energetic than the flares observed around the primary transits of the planet. We modeled the distribution of the number of flares per day with a power law, showing that it is steeper than the one observed in the Sun and in other Main Sequence stars. The steepness hints at a significant fraction of undetected micro-flares. Altogether, the plasma temperatures below 1 keV observed during the flares, along with the slightly larger fraction of energetic flares seen at the secondary transits highlight the peculiarity of the corona of HD 189733 A and points to star-planet interaction as the plausible origin of part of its X-ray emission. However, more observational and modeling efforts are required to confirm or disprove this scenario.

preprint2021arXiv

New constraints on the future evaporation of the young exoplanets in the V1298 Tau system

Transiting planets at young ages are key targets for improving our understanding of the evolution of exo-atmospheres. We present results of a new X-ray observation of V1298 Tau with XMM-Newton, aimed to determine more accurately the high-energy irradiation of the four planets orbiting this pre-main-sequence star, and the possible variability due to magnetic activity on short and long time scales. Following the first measurements of planetary masses in the V1298 Tau system, we revise early guesses of the current escape rates from the planetary atmospheres, employing our updated atmospheric evaporation models to predict the future evolution of the system. Contrary to previous expectations, we find that the two outer Jupiter-size planets will not be affected by any evaporation on Gyr time scales, and the same occurs for the two smaller inner planets, unless their true masses are lower than $\sim 40$ M$_\oplus$. These results confirm that relatively massive planets can reach their final position in the mass-radius diagram very early in their evolutionary history.

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

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.

preprint2021arXiv

The GAPS Programme at TNG. XXIX. No detection of reflected light from 51 Peg b using optical high-resolution spectroscopy

The analysis of exoplanetary atmospheres by means of high-resolution spectroscopy is an expanding research field which provides information on chemical composition, thermal structure, atmospheric dynamics and orbital velocity of exoplanets. In this work, we aim at the detection of the light reflected by the exoplanet 51~Peg~b employing optical high-resolution spectroscopy. To detect the light reflected by the planetary dayside we use optical HARPS and HARPS-N spectra taken near the superior conjunction of the planet, when the flux contrast between the planet and the star is maximum. To search for the weak planetary signal, we cross-correlate the observed spectra with a high S/N stellar spectrum. We homogeneously analyze the available datasets and derive a $10^{-5}$ upper limit on the planet-to-star flux contrast in the optical. The upper limit on the planet-to-star flux contrast of $10^{-5}$ translates into a low albedo of the planetary atmosphere ($\rm A_g\lesssim0.05-0.15$ for an assumed planetary radius in the range $\rm 1.5-0.9~R_{Jup}$, as estimated from the planet's mass).

preprint2021arXiv

The GAPS Programme at TNG. XXXII. The revealing non-detection of metastable HeI in the atmosphere of the hot Jupiter WASP-80b

The hot Jupiter WASP-80b has been identified as a possible excellent target for detecting and measuring HeI absorption in the upper atmosphere. We observed 4 primary transits of WASP-80b in the optical and near-IR using the HARPS-N and GIANO-B high-resolution spectrographs, focusing on the HeI triplet. We further employed a three-dimensional hydrodynamic aeronomy model to understand the observational results. We did not find any signature of planetary absorption at the position of the HeI triplet with an upper limit of 0.7% (i.e. 1.11 planetary radii; 95% confidence level). We re-estimated the stellar high-energy emission that we combined with a stellar photospheric model to generate the input for the hydrodynamic modelling. We obtained that, assuming a solar He to H abundance ratio, HeI absorption should have been detected. Considering a stellar wind 25 times weaker than solar, we could reproduce the non-detection only assuming a He to H abundance ratio about 16 times smaller than solar. Instead, considering a stellar wind 10 times stronger than solar, we could reproduce the non-detection only with a He to H abundance ratio about 10 times smaller than solar. We attempted to understand this result by collecting all past HeI measurements looking for correlations with stellar high-energy emission and planetary gravity, but without finding any. WASP-80b is not the only planet with a sub-solar estimated He to H abundance ratio, suggesting the presence of efficient physical mechanisms (e.g. phase separation, magnetic fields) capable of significantly modifying the He to H content in the upper atmosphere of hot Jupiters. The planetary macroscopic properties and the shape of the stellar spectral energy distribution are not sufficient for predicting the presence or absence of detectable metastable He in a planetary atmosphere, as also the He abundance appears to play a major role.

preprint2020arXiv

A systematic study of \ce{CO2} planetary atmospheres and their link to the stellar environment

The Milky Way Galaxy is literally teeming with exoplanets; thousands of planets have been discovered, with thousands more planet candidates identified. Terrestrial-like planets are quite common around other stars, and are expected to be detected in large numbers in the future. Such planets are the primary targets in the search for potentially habitable conditions outside the solar system. Determining the atmospheric composition of exoplanets is mandatory to understand their origin and evolution, as atmospheric processes play crucial roles in many aspects of planetary architecture. In this work we construct and exploit a 1D radiative transfer model based on the discrete-ordinates method in plane-parallel geometry. Radiative results are linked to a convective flux that redistributes energy at any altitude producing atmospheric profiles in radiative-convective equilibrium. The model has been applied to a large number (6250) of closely dry synthetic \ce{CO2} atmospheres, and the resulting pressure and thermal profiles have been interpreted in terms of parameter variability. Although less accurate than 3D general circulation models, not properly accounting for e.g., clouds and atmospheric and ocean dynamics, 1D descriptions are computationally inexpensive and retain significant value by allowing multidimensional parameter sweeps with relative ease.

preprint2020arXiv

An ultra-short period rocky super-Earth orbiting the G2-star HD 80653

Ultra-short period (USP) planets are a class of exoplanets with periods shorter than one day. The origin of this sub-population of planets is still unclear, with different formation scenarios highly dependent on the composition of the USP planets. A better understanding of this class of exoplanets will, therefore, require an increase in the sample of such planets that have accurate and precise masses and radii, which also includes estimates of the level of irradiation and information about possible companions. Here we report a detailed characterization of a USP planet around the solar-type star HD 80653 $\equiv$ EP 251279430 using the K2 light curve and 108 precise radial velocities obtained with the HARPS-N spectrograph, installed on the Telescopio Nazionale Galileo. From the K2 C16 data, we found one super-Earth planet ($R_{b}=1.613\pm0.071 R_{\oplus}$) transiting the star on a short-period orbit ($P_{\rm b}=0.719573\pm0.000021$ d). From our radial velocity measurements, we constrained the mass of HD 80653 b to $M_{b}=5.60\pm0.43 M_{\oplus}$. We also detected a clear long-term trend in the radial velocity data. We derived the fundamental stellar parameters and determined a radius of $R_{\star}=1.22\pm0.01 R_{\odot}$ and mass of $M_{\star}=1.18\pm0.04 M_{\odot}$, suggesting that HD 80653, has an age of $2.7\pm1.2$ Gyr. The bulk density ($ρ_{b} = 7.4 \pm 1.1$ g cm$^{-3}$) of the planet is consistent with an Earth-like composition of rock and iron with no thick atmosphere. Our analysis of the K2 photometry also suggests hints of a shallow secondary eclipse with a depth of 8.1$\pm$3.7 ppm. Flux variations along the orbital phase are consistent with zero. The most important contribution might come from the day-side thermal emission from the surface of the planet at $T\sim3480$ K.

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

GJ 357 b: A Super-Earth Orbiting an Extremely Inactive Host Star

In this paper we present a deep X-ray observation of the nearby M dwarf GJ 357 and use it to put constraints on the atmospheric evolution of its planet, GJ 357 b. We also analyse the systematic errors in the stellar parameters of GJ 357 in order to see how they affect the perceived planetary properties. We estimate the age of GJ 357 b by comparing the observed X-ray luminosity of its host star, derived from a recent {\em XMM-Newton} observation {($\log{L_{\rm x}}\,{\rm [erg/s]} = 25.73$), with $L_{\rm x} -$ age relations for M dwarfs. We find that GJ 357 presents one of the lowest X-ray activity levels ever measured for an M dwarf, and we put a lower limit on its age of $5$\,Gyr.} Using this age limit, we perform a backwards reconstruction of the original primordial atmospheric reservoir. Furthermore, by considering the systematic errors in the stellar parameters, we find a range of possible planetary masses, radii, and densities. From the backwards reconstruction of GJ 357 b's irradiation history we find that the upper limit of its initial primordial atmospheric mass is $\sim \rm 38M_{\oplus}$. An initial atmospheric reservoir significantly larger than this may have survived through the X-ray and ultraviolet irradiation history, hence being inconsistent with current observations that suggest a telluric composition. In spite of the unlikelihood of a currently existing primordial envelope, volcanism and outgassing may have contributed to a secondary atmosphere. Under this assumption, we present three different synthetic infrared spectra for GJ 357 b that one might expect, consisting of $100\%~\rm CO_{2}$, $100\%~\rm SO_{2}$, and $75\%~ \rm N_{2}$, $24\%~\rm CO_{2}$ and $1\%~\rm H_{2}O$.

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

Neutral Iron Emission Lines From The Day-side Of KELT-9b -- The GAPS Programme With HARPS-N At TNG XX

We present the first detection of atomic emission lines from the atmosphere of an exoplanet. We detect neutral iron lines from the day-side of KELT-9b (Teq $\sim$ 4, 000 K). We combined thousands of spectrally resolved lines observed during one night with the HARPS-N spectrograph (R $\sim$ 115, 000), mounted at the Telescopio Nazionale Galileo. We introduce a novel statistical approach to extract the planetary parameters from the binary mask cross-correlation analysis. We also adapt the concept of contribution function to the context of high spectral resolution observations, to identify the location in the planetary atmosphere where the detected emission originates. The average planetary line profile intersected by a stellar G2 binary mask was found in emission with a contrast of 84 $\pm$ 14 ppm relative to the planetary plus stellar continuum (40 $\pm$ 5$\%$ relative to the planetary continuum only). This result unambiguously indicates the presence of an atmospheric thermal inversion. Finally, assuming a modelled temperature profile previously published (Lothringer et al. 2018), we show that an iron abundance consistent with a few times the stellar value explains the data well. In this scenario, the iron emission originates at the $10^{-3}$-$10^{-5}$ bar level.

preprint2020arXiv

Principal Component Analysis to correct data systematics. Case study: K2 light curves

Instrumental data are affected by systematic effects that dominate the errors and can be relevant when searching for small signals. This is the case of the K2 mission, a follow up of the Kepler mission, that, after a failure on two reaction wheels, has lost its stability properties rising strongly the systematics in the light curves and reducing its photometric precision. In this work, we have developed a general method to remove time related systematics from a set of light curves, that has been applied to K2 data. The method uses the Principal Component Analysis to retrieve the correlation between the light curves due to the systematics and to remove its effect without knowing any information other than the data itself. We have applied the method to all the K2 campaigns available at the Mikulski Archive for Space Telescopes, and we have tested the effectiveness of the procedure and its capability in preserving the astrophysical signal on a few transits and on eclipsing binaries. One product of this work is the identification of stable sources along the ecliptic plane that can be used as photometric calibrators for the upcoming Atmospheric Remote-sensing Exoplanet Large-survey mission.

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

The GAPS Programme at TNG -- XXIII. HD 164922 d: a close-in super-Earth discovered with HARPS-N in a system with a long-period Saturn mass companion

In the framework of the Global Architecture of Planetary Systems (GAPS) project we collected more than 300 spectra with HARPS-N at the TNG for the bright G9V star HD164922. This target is known to host one gas giant planet in a wide orbit (Pb~1200 days, semi-major axis ~2 au) and a Neptune-mass planet with a period Pc ~76 days. We searched for additional low-mass companions in the inner region of the system. We compared the radial velocities (RV) and the activity indices derived from the HARPS-N time series to measure the rotation period of the star and used a Gaussian process regression to describe the behaviour of the stellar activity. We exploited this information in a combined model of planetary and stellar activity signals in an RV time-series composed of almost 700 high-precision RVs, both from HARPS-N and literature data. We performed a dynamical analysis to evaluate the stability of the system and the allowed regions for additional potential companions. Thanks to the high sensitivity of the HARPS-N dataset, we detect an additional inner super-Earth with an RV semi-amplitude of 1.3+/-0.2 m/s, a minimum mass of ~4+/-1 M_E and a period of 12.458+/-0.003 days. We disentangle the planetary signal from activity and measure a stellar rotation period of ~42 days. The dynamical analysis shows the long term stability of the orbits of the three-planet system and allows us to identify the permitted regions for additional planets in the semi-major axis ranges 0.18-0.21 au and 0.6-1.4 au. The latter partially includes the habitable zone of the system. We did not detect any planet in these regions, down to minimum detectable masses of 5 and 18 M_E, respectively. A larger region of allowed planets is expected beyond the orbit of planet b, where our sampling rules-out bodies with minimum mass > 50 M_E.

preprint2020arXiv

The GAPS Programme at TNG -- XXIV. An eccentric Neptune-mass planet near the inner edge of the BD-11 4672 habitable zone

With the growth of comparative exoplanetology, it is increasingly clear that investigating the relationships between inner and outer planets plays a key role in discriminating between competing formation and evolution models. To do so, it is important to probe the inner region of systems hosting long-period giants in search for undetected lower-mass planetary companions. In this work we present our results on the K-dwarf star BD-11~4672, already known to host a long-period giant planet, as the first output of a subsample of the GAPS programme specifically aimed at assessing the impact of inefficient migration of planets formed beyond the snowline by searching for Neptune-mass and super-Earths planetary companions of known longer-period giants. The high-precision HARPS-N observations of BD-11~4672 are used in conjunction with literature time series in order to search for additional inner planetary signals to be fitted using differential evolution Markov chain Monte Carlo. The long-term stability of the new orbital solutions is tested by N-body dynamical simulations. We report the detection of BD-11~4672~c, a new Neptune-mass planet with an orbital period of $74.20_{-0.08}^{+0.06}\,\rm{d}$, eccentricity $0.40_{-0.15}^{+0.13}$, semimajor axis $0.30\pm0.01\,\rm{au}$ and minimum mass $15.37_{-2.81}^{+2.97}\,M_\oplus$ orbiting slightly outside the inner edge of the optimistic circumstellar habitable zone. In order to assess its impact on the dynamical stability of the habitable zone we compute the angular momentum deficit of the system, showing that planet c has a severe negative impact on the stability of possible additional lower-mass temperate planets. The BD-11~4672 system is notable for its architecture, hosting both a long-period giant planet and an inner lower-mass planet, the latter being also among the most eccentric Neptune-mass planets known at similar periods.

preprint2020arXiv

The GAPS Programme at TNG -- XXV. Stellar atmospheric parameters and chemical composition through GIARPS optical and near-infrared spectra

The detailed chemical composition of stars is important in many astrophysical fields, among which the characterisation of exoplanetary systems. Previous studies seem to indicate an anomalous chemical pattern of the youngest stellar population in the solar vicinity with a sub-solar metal content. This can influence various observational relations linking the properties of exoplanets to the characteristics of the host stars, for example the giant planet-metallicity relation. In this framework, we aim to expand our knowledge of the chemical composition of intermediate-age stars and understand whether these peculiarities are real or related to spectroscopic analysis techniques. We analysed high-resolution optical and near-infrared GIARPS spectra of intermediate-age stars (< 700Myr). To overcome issues related to the young ages of the stars, we applied a new spectroscopic method that uses titanium lines to derive the atmospheric parameters, in particular surface gravities and microturbulence velocity parameter. We also derived abundances of 14 different atomic species. The lack of systematic trends between elemental abundances and effective temperatures validates our method. However, we observed that the coolest (<5400 K) stars in the sample, display higher abundances for the Cr II, and for high-excitation potential C I lines. We found a positive correlation between the higher abundances measured of C I and Cr II and the activity index logR$_{HK}$. Instead, we found no correlations between the C abundances obtained from CH molecular band at 4300Å, and both effective temperatures and activity. Thus, we suggest that these are better estimates for C abundances in young and cool stars. Finally, we found an indication of an increasing abundance ratio [X/H] with the condensation temperature for HD167389, indicating possible episodes of planet engulfment.

preprint2020arXiv

The GAPS Programme at TNG XXI -- A GIARPS case-study of known young planetary candidates: confirmation of HD 285507 b and refutation of AD Leo b

The existence of hot Jupiters is still not well understood. Two main channels are thought to be responsible for their current location: a smooth planet migration through the proto-planetary disk or the circularization of an initial high eccentric orbit by tidal dissipation leading to a strong decrease of the semimajor axis. Different formation scenarios result in different observable effects, such as orbital parameters (obliquity/eccentricity), or frequency of planets at different stellar ages. In the context of the GAPS Young-Objects project, we are carrying out a radial velocity survey with the aim to search and characterize young hot-Jupiter planets. Our purpose is to put constraints on evolutionary models and establish statistical properties, such as the frequency of these planets from a homogeneous sample. Since young stars are in general magnetically very active, we performed multi-band (visible and near-infrared) spectroscopy with simultaneous GIANO-B + HARPS-N (GIARPS) observing mode at TNG. This helps to deal with stellar activity and distinguish the nature of radial velocity variations: stellar activity will introduce a wavelength-dependent radial velocity amplitude, whereas a Keplerian signal is achromatic. As a pilot study, we present here the cases of two already claimed hot Jupiters orbiting young stars: HD285507 b and AD Leo b. Our analysis of simultaneous high-precision GIARPS spectroscopic data confirms the Keplerian nature of HD285507's radial velocities variation and refines the orbital parameters of the hot Jupiter, obtaining an eccentricity consistent with a circular orbit. On the other hand, our analysis does not confirm the signal previously attributed to a planet orbiting AD Leo. This demonstrates the power of the multi-band spectroscopic technique when observing active stars.

preprint2020arXiv

The GAPS programme at TNG XXII. The GIARPS view of the extended helium atmosphere of HD189733 b accounting for stellar activity

Exoplanets orbiting very close to their host star are strongly irradiated. This can lead the upper atmospheric layers to expand and evaporate into space. The metastable helium (HeI) triplet at 1083.3nm has recently been shown to be a powerful diagnostic to probe extended and escaping exoplanetary atmosphere. We perform high-resolution transmission spectroscopy of the transiting hot Jupiter HD189733b with the GIARPS (GIANO-B + HARPS-N) observing mode of the Telescopio Nazionale Galileo, taking advantage of the simultaneous optical+near infrared spectral coverage to detect HeI in the planet's extended atmosphere and to gauge the impact of stellar magnetic activity on the planetary absorption signal. Observations were performed during five transit events of HD189733b. By comparison of the in- and out-of-transit GIANO-B observations we compute high-resolution transmission spectra, on which we perform equivalent width measurements and light-curves analyses to gauge the excess in-transit absorption in the HeI triplet. We detect an absorption signal during all five transits. The mean in-transit absorption depth amounts to 0.75+/-0.03%. We detect night-to-night variations in the HeI absorption signal likely due to the transit events occurring in presence of stellar surface inhomogeneities. We evaluate the impact of stellar-activity pseudo-signals on the true planetary absorption using a comparative analysis of the HeI and the H$α$ lines. We interpret the time-series of the HeI absorption lines in the three nights not affected by stellar contamination -exhibiting a mean in-transit absorption depth of 0.77+/-0.04%- using a 3-d atmospheric code. Our simulations suggest that the helium layers only fill part of the Roche lobe. Observations can be explained with a thermosphere heated to $\sim$12000 K, expanding up to $\sim$1.2 planetary radii, and losing $\sim$1 g/s of metastable helium.

preprint2020arXiv

The spectral impact of magnetic activity on disk-integrated HARPS-N solar observations: exploring new activity indicators

Stellar activity is the major roadblock on the path to finding true Earth-analogue planets with the Doppler technique. Thus, identifying new indicators that better trace magnetic activity (i.e. faculae and spots) is crucial to aid in disentangling these signals from that of a planet's Doppler wobble. In this work, we investigate activity related features as seen in disk-integrated spectra from the HARPS-N solar telescope. We divide high-activity spectral echelle orders by low-activity master templates (as defined using both log R'HK and images from the Solar Dynamics Observatory, SDO), creating "relative spectra". With resolved images of the surface of the Sun (via SDO), the faculae and spot filling factors can be calculated, giving a measure of activity independent of, and in addition to, log R'HK. We find pseudo-emission (and pseudo-absorption) features in the relative spectra that are similar to those reported in our previous work on alpha Cen B. In alpha Cen B, the features are shown to correlate better to changes in faculae filling factor than spot filling factor. In this work we more confidently identify changes in faculae coverage of the visible hemisphere of the Sun as the source of features produced in the relative spectra. Finally, we produce trailed spectra to observe the RV component of the features, which show that the features move in a redward direction as one would expect when tracking active regions rotating on the surface of a star.

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.

preprint2016arXiv

A Chandra X-ray study of the young star cluster NGC 6231: low-mass population and initial mass function

NGC6231 is a massive young star cluster, near the center of the Sco OB1 association. While its OB members are well studied, its low-mass population has received little attention. We present high-spatial resolution Chandra ACIS-I X-ray data, where we detect 1613 point X-ray sources. Our main aim is to clarify global properties of NGC6231 down to low masses through a detailed membership assessment, and to study the cluster stars' spatial distribution, the origin of their X-ray emission, the cluster age and formation history, and initial mass function. We use X-ray data, complemented by optical/IR data, to establish cluster membership. The spatial distribution of different stellar subgroups also provides highly significant constraints on cluster membership, as does the distribution of X-ray hardness. We perform spectral modeling of group-stacked X-ray source spectra. We find a large cluster population down to ~0.3 Msun (complete to ~1 Msun), with minimal non-member contamination, with a definite age spread (1-8 Myrs) for the low-mass PMS stars. We argue that low-mass cluster stars also constitute the majority of the few hundreds unidentified X-ray sources. We find mass segregation for the most massive stars. The fraction of circumstellar-disk bearing members is found to be ~5%. Photoevaporation of disks under the action of massive stars is suggested by the spatial distribution of the IR-excess stars. We also find strong Halpha emission in 9% of cluster PMS stars. The dependence of X-ray properties on mass, stellar structure, and age agrees with extrapolations based on other young clusters. The cluster initial mass function, computed over ~2 dex in mass, has a slope Gamma~-1.14. The total mass of cluster members above 1 Msun is 2280 Msun, and the inferred total mass is 4380 Msun. We also study the peculiar, hard X-ray spectrum of the Wolf-Rayet star WR79.

preprint2016arXiv

Gaia-ESO Survey: Gas dynamics in the Carina Nebula through optical emission lines

We present observations from the Gaia-ESO Survey in the lines of H$α$, [N II], [S II] and He I of nebular emission in the central part of the Carina Nebula. We investigate the properties of the two already known kinematic components (approaching and receding, respectively), which account for the bulk of emission. Moreover, we investigate the features of the much less known low-intensity high-velocity (absolute RV $>$50 km/s) gas emission. We show that gas giving rise to H$α$ and He I emission is dynamically well correlated, but not identical, to gas seen through forbidden-line emission. Gas temperatures are derived from line-width ratios, and densities from [S II] doublet ratios. The spatial variation of N ionization is also studied, and found to differ between the approaching and receding components. The main result is that the bulk of the emission lines in the central part of Carina arises from several distinct shell-like expanding regions, the most evident found around $η$ Car, the Trumpler 14 core, and the star WR25. Such "shells" are non-spherical, and show distortions probably caused by collisions with other shells or colder, higher-density gas. Part of them is also obscured by foreground dust lanes, while only very little dust is found in their interior. Preferential directions, parallel to the dark dust lanes, are found in the shell geometries and physical properties, probably related to strong density gradients in the studied region. We also find evidence that the ionizing flux emerging from $η$ Car and the surrounding Homunculus nebula varies with polar angle. The high-velocity components in the wings of H$α$ are found to arise from expanding dust reflecting the $η$ Car spectrum.

preprint2016arXiv

Photoevaporation and close encounters: how the environment around Cygnus OB2 affects the evolution of protoplanetary disks

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy.

preprint2016arXiv

Rotation in the Pleiades with K2: I. Data and First Results

Young (125 Myr), populous ($>$1000 members), and relatively nearby, the Pleiades has provided an anchor for stellar angular momentum models for both younger and older stars. We used K2 to explore the distribution of rotation periods in the Pleiades. With more than 500 new periods for Pleiades members, we are vastly expanding the number of Pleiads with periods, particularly at the low mass end. About 92\% of the members in our sample have at least one measured spot-modulated rotation period. For the $\sim$8\% of the members without periods, non-astrophysical effects often dominate (saturation, etc.), such that periodic signals might have been detectable, all other things being equal. We now have an unusually complete view of the rotation distribution in the Pleiades. The relationship between $P$ and $(V-K_{\rm s})_0$ follows the overall trends found in other Pleiades studies. There is a slowly rotating sequence for $1.1\lesssim(V-K_{\rm s})_0\lesssim 3.7$, and a primarily rapidly rotating population for $(V-K_{\rm s})_0\gtrsim 5.0$. There is a region in which there seems to be a disorganized relationship between $P$ and $(V-K_{\rm s})_0$ for $3.7 \lesssim(V-K_{\rm s})_0\lesssim 5.0$. Paper II continues the discussion, focusing on multi-period structures, and Paper III speculates about the origin and evolution of the period distribution in the Pleiades.

preprint2016arXiv

Rotation in the Pleiades with K2: II. Multi-Period Stars

We use K2 to continue the exploration of the distribution of rotation periods in Pleiades that we began in Paper I. We have discovered complicated multi-period behavior in Pleiades stars using these K2 data, and we have grouped them into categories, which are the focal part of this paper. About 24% of the sample has multiple, real frequencies in the periodogram, sometimes manifesting as obvious beating in the light curves. Those having complex and/or structured periodogram peaks, unresolved multiple periods, and resolved close multiple periods are likely due to spot/spot group evolution and/or latitudinal differential rotation; these largely compose the slowly rotating sequence in $P$ vs.~$(V-K_{\rm s})_0$ identified in Paper I. The fast sequence in $P$ vs.~$(V-K_{\rm s})_0$ is dominated by single-period stars; these are likely to be rotating as solid bodies. Paper III continues the discussion, speculating about the origin and evolution of the period distribution in the Pleiades.

preprint2016arXiv

Rotation in the Pleiades With K2: III. Speculations on Origins and Evolution

We use high quality K2 light curves for hundreds of stars in the Pleiades to understand better the angular momentum evolution and magnetic dynamos of young, low mass stars. The K2 light curves provide not only rotational periods but also detailed information from the shape of the phased light curve not available in previous studies. A slowly rotating sequence begins at $(V-K_{\rm s})_0\sim$1.1 (spectral type F5) and ends at $(V-K_{\rm s})_0\sim$ 3.7 (spectral type K8), with periods rising from $\sim$2 to $\sim$11 days in that interval. Fifty-two percent of the Pleiades members in that color interval have periods within 30\% of a curve defining the slow sequence; the slowly rotating fraction decreases significantly redward of $(V-K_{\rm s})_0$=2.6. Nearly all of the slow-sequence stars show light curves that evolve significantly on timescales less than the K2 campaign duration. The majority of the FGK Pleiades members identified as photometric binaries are relatively rapidly rotating, perhaps because binarity inhibits star-disk angular momentum loss mechanisms during pre-main sequence evolution. The fully convective, late M dwarf Pleiades members (5.0 $<(V-K_{\rm s})_0<$ 6.0) nearly always show stable light curves, with little spot evolution or evidence of differential rotation. During PMS evolution from $\sim$3 Myr (NGC2264 age) to $\sim$125 Myr (Pleiades age), stars of 0.3 $M_{\odot}$ shed about half their angular momentum, with the fractional change in period between 3 and 125 Myr being nearly independent of mass for fully convective stars. Our data also suggest that very low mass binaries form with rotation periods more similar to each other and faster than would be true if drawn at random from the parent population of single stars.

preprint2016arXiv

Spot modelling of periodic weak-line T Tauri stars observed by CoRoT in NGC 2264

The space telescope CoRoT has provided light curves of T Tauri stars belonging to the star-forming region of NGC 2264 with unprecedented continuity and precision in the framework of a coordinated multi-wavelength observational project. We perform spot modelling of the optical light curves of five weak-line T Tauri stars whose variability is dominated by starspots. We apply two-spot and evolving single spot models in the framework of a Bayesian Monte Carlo Markov Chain approach to derive the a posteriori distribution of the starspot parameters and the inclination of the star rotation axis. We focus on the rotation periods of the spots that can provide evidence for differential rotation in those stars. We find meaningful results in the case of three stars with an inclination higher than 50 deg and a slow variation of the light curve amplitude. The relative difference of the spot rotation periods ranges from 0.02 to 0.05 that is 3-10 times larger than the amplitude of the differential rotation found in similar stars with Doppler imaging techniques. We conclude that the intrinsic starspot evolution, although very slow, has a significant impact on the determination of the differential rotation by means of our spot modelling approach. We estimate typical timescales for the evolution of the starspot pattern between about 20 and 50 rotation periods in our stars.

preprint2016arXiv

Star-disk interaction in classical T Tauri stars revealed using wavelet analysis

The extension of the corona of classical T Tauri stars (CTTS) is under discussion. The standard model of magnetic configuration of CTTS predicts that coronal magnetic flux tubes connect the stellar atmosphere to the inner region of the disk. However, differential rotation may disrupt these long loops. The results from Hydrodynamic modeling of X-ray flares observed in CTTS confirming the star-disk connection hypothesis are still controversial. Some authors suggest the presence of the accretion disk prevent the stellar corona to extent beyond the co-rotation radius, while others simply are not confident with the methods used to derive loop lengths. We use independent procedures to determine the length of flaring loops in stars of the Orion Nebula Cluster previously analyzed using Hydrodynamic models. Our aim is to disentangle between the two scenarios proposed. We present a different approach to determine the length of flaring loops based on the oscillatory nature of the loops after strong flares. We use wavelet tools to reveal oscillations during several flares. The subsequent analysis of such oscillations is settle on the Physics of coronal seismology. Our results likely confirm the large extension of the corona of CTTS and the hypothesis of star-disk magnetic interaction in at least three CTTS of the Orion Nebula Cluster. Analyzing oscillations in flaring events is a powerful tool to determine the physical characteristics of magnetic loops in coronae in stars other than the Sun. The results presented in this work confirm the star-disk magnetic connection in CTTS.

preprint2016arXiv

Supersaturation and activity-rotation relation in PMS stars: the young Cluster h Per

The magnetic activity of late-type MS stars is characterized by different regimes, and their activity levels are well described by Ro, the ratio between P_rot and the convective turnover time. Very young PMS stars show, similarly to MS stars, intense magnetic activity. However they do not show clear activity-rotation trends, and it still debated which stellar parameters determine their magnetic activity levels. To bridge the gap between MS and PMS stars, we studied the activity-rotation relation in the young cluster h Per, a ~13 Myr old cluster, that contains both fast and slow rotators, whose members have ended their accretion phase and have already developed a radiative core. It offers us the opportunity to study the activity level of intermediate-age PMS stars with different rotational velocities, excluding any interactions with the circumstellar environment. We constrained the magnetic activity levels of h Per members measuring their X-ray emission from a Chandra observation, while P_rot were obtained by Moraux et al. (2013). We collected a final catalog of 414 h Per members with known P_rot, T_eff, M_star, with 169 of them having also detected X-ray emission. We found that h Per members, with 1.0 M_sun < M_star < 1.4 M_sun, display different activity regimes: fast rotators show supersaturation, while slower rotators have activity levels compatible to the non-saturated regime. At 13 Myr h Per is therefore the youngest cluster showing activity-rotation regimes analogous to that of MS stars, indicating that, at this age, magnetic field production is likely regulated by the alpha-Omega type dynamo. Moreover we observed that supersaturation is better described by P_rot than Ro, and that the observed patterns are compatible with the hypothesis of centrifugal stripping. In this scenario we inferred that coronae can produce structures as large as ~2 R_star above the stellar surface.

preprint2016arXiv

The Gaia-ESO Survey: A lithium-rotation connection at 5 Myr?

We investigate here the lithium content of young low-mass stars in the 5~Myr-old star forming region NGC~2264 and its relationship with rotation. We combine lithium equivalent width measurements, EW(Li), from the Gaia-ESO Survey with the determination of rotational periods from the CSI~2264 survey. We consider only bona fide non accreting cluster members in order to minimize uncertainties on EW(Li). We report the existence of a relationship between lithium content and rotation in NGC~2264 at an age of 5~Myr. The Li-rotation connection is seen over a restricted temperature range, Teff=3800--4400K, where fast rotators are Li-rich compared to slow ones. This correlation is similar to, albeit of lower amplitude than, the Li-rotation connection previously reported for K dwarfs in the 125 Myr-old Pleiades cluster. We investigate whether the non-standard pre-main sequence models developed so far to explain the Pleiades results, which are based on episodic accretion, pre-main sequence core-envelope decoupling, and/or radius inflation due to enhanced magnetic activity, can account for an early development of the Li-rotation connection. While radius inflation appears to be the most promising possibility, each of these models has issues. We therefore also discuss external causes that might operate during the first few Myr of pre-main sequence evolution, such as planet engulfment and/or steady disk accretion, as possible candidates for the common origin for Li-excess and fast rotation in young low-mass pre-main sequence stars. The emergence of a connection between lithium content and rotation rate at such an early age as 5~Myr suggests a complex link between accretion processes, early angular momentum evolution, and possibly planet formation, which likely impacts early stellar evolution and still is to be fully deciphered.

preprint2016arXiv

The Gaia-ESO Survey: membership and Initial Mass Function of the Gamma Velorum cluster

Understanding the properties of young open clusters, such as the Initial Mass Function (IMF), star formation history and dynamic evolution, is crucial to obtain reliable theoretical predictions of the mechanisms involved in the star formation process. We want to obtain a list, as complete as possible, of confirmed members of the young open cluster Gamma Velorum, with the aim of deriving general cluster properties such as the IMF. We used all available spectroscopic membership indicators within the Gaia-ESO public archive together with literature photometry and X-ray data and, for each method, we derived the most complete list of candidate cluster members. Then, we considered photometry, gravity and radial velocities as necessary conditions to select a subsample of candidates whose membership was confirmed by using the lithium and H$α$ lines and X-rays as youth indicators. We found 242 confirmed and 4 possible cluster members for which we derived masses using very recent stellar evolutionary models. The cluster IMF in the mass range investigated in this study shows a slope of $α=2.6\pm0.5$ for $0.5<M/M_\odot <1.3$ and $α=1.1\pm0.4$ for $0.16<M/M_\odot <0.5$ and is consistent with a standard IMF. The similarity of the IMF of the young population around $γ^2 $Vel to that in other star forming regions and the field suggests it may have formed through very similar processes.

preprint2016arXiv

The GAPS programme with HARPS-N at TNG XI. Pr~0211 in M~44: the first multi-planet system in an open cluster

Open cluster (OC) stars share the same age and metallicity, and, in general, their age and mass can be estimated with higher precision than for field stars. For this reason, OCs are considered an important laboratory to study the relation between the physical properties of the planets and those of their host stars, and the evolution of planetary systems. We started an observational campaign within the GAPS collaboration to search for and characterize planets in OCs We monitored the Praesepe member Pr0211 to improve the eccentricity of the Hot-Jupiter (HJ) already known to orbit this star and search for additional planets. An eccentric orbit for the HJ would support a planet-planet scattering process after its formation. From 2012 to 2015, we collected 70 radial velocity (RV) measurements with HARPS-N and 36 with TRES of Pr0211. Simultaneous photometric observations were carried out with the robotic STELLA telescope in order to characterize the stellar activity. We discovered a long-term trend in the RV residuals that we show to be due to the presence of a second, massive, outer planet. Orbital parameters for the two planets are derived by simultaneously fitting RVs and photometric light curves, with the activity signal modelled as a series of sinusoids at the rotational period of the star and its harmonics. We confirm that Pr0211b has a nearly circular orbit ($e = 0.02 \pm 0.01$), with an improvement of a factor two with respect to the previous determination of its eccentricity, and estimate that Pr0211c has a mass $M_p\sin i = 7.9 \pm 0.2 M_J$, a period $P>$3500 days and a very eccentric orbit ($e>$0.60). Such peculiar systems may be typical of open clusters if the planet-planet scattering phase leading to the formation of HJs is caused by stellar encounters rather than unstable primordial orbits. Pr0211 is the first multi-planet system discovered around an OC star. (abridged)

preprint2015arXiv

Coordinated X-ray and Optical observations of Star-Planet Interaction in HD 17156

The large number of close-in Jupiter-size exoplanets prompts the question whether star-planet interaction (SPI) effects can be detected. We focused our attention on the system HD 17156, having a Jupiter-mass planet in a very eccentric orbit. Here we present results of the XMM-Newton observations and of a five months coordinated optical campaign with the HARPS-N spectrograph. We observed HD 17156 with XMM-Newton when the planet was approaching the apoastron and then at the following periastron passage, quasi simultaneously with HARPS-N. We obtained a clear ($\approx 5.5σ$) X-ray detection only at the periastron visit, accompanied by a significant increase of the $R'_{\rm HK}$ chromospheric index. We discuss two possible scenarios for the activity enhancement: magnetic reconnection and flaring or accretion onto the star of material tidally stripped from the planet. In any case, this is possibly the first evidence of a magnetic SPI effect caught in action.

preprint2015arXiv

FUV variability of HD 189733. Is the star accreting material from its hot Jupiter?

Hot Jupiters are subject to strong irradiation from the host stars and, as a consequence, they do evaporate. They can also interact with the parent stars by means of tides and magnetic fields. Both phenomena have strong implications for the evolution of these systems. Here we present time resolved spectroscopy of HD~189733 observed with the Cosmic Origin Spectrograph (COS) on board to HST. The star has been observed during five consecutive HST orbits, starting at a secondary transit of the planet ($ϕ$ ~0.50-0.63). Two main episodes of variability of ion lines of Si, C, N and O are detected, with an increase of line fluxes. Si IV lines show the highest degree of variability. The FUV variability is a signature of enhanced activity in phase with the planet motion, occurring after the planet egress, as already observed three times in X-rays. With the support of MHD simulations, we propose the following interpretation: a stream of gas evaporating from the planet is actively and almost steadily accreting onto the stellar surface, impacting at $70-90°$ ahead of the sub-planetary point.

preprint2015arXiv

Gaia-ESO Survey: The analysis of pre-main sequence stellar spectra

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

preprint2015arXiv

Large Binocular Telescope view of the atmosphere of GJ1214b

The atmospheric composition and vertical structure of the super-Earth GJ1214b has been a subject of debate since its discovery in 2009. Recent studies have indicated that high-altitude clouds might mask the lower layers. However, some data points that were gathered at different times and facilities do not fit this picture, probably because of a combination of stellar activity and systematic errors. We observed two transits of GJ1214b with the Large Binocular Camera, the dual-channel camera at the Large Binocular Telescope. For the first time, we simultaneously measured the relative planetary radius $k=R_\mathrm{p}/R_\star$ at blue and red optical wavelengths ($B+R$), thus constraining the Rayleigh scattering on GJ1214b after correcting for stellar activity effects. To the same purpose, a long-term photometric follow-up of the host star was carried out with WiFSIP at STELLA, revealing a rotational period that is significantly longer than previously reported. Our new unbiased estimates of $k$ yield a flat transmission spectrum extending to shorter wavelengths, thus confirming the cloudy atmosphere scenario for GJ1214b.

preprint2015arXiv

Revisiting Spitzer transit observations with Independent Component Analysis: new results for the GJ436 system

We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 μm of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to extract conclusively an atmospheric signal; those analyses also reported discrepant results, hence the necessity of this reanalysis. The method we used has been proposed in Morello et al. (2014) to analyze 3.6 μm transit light-curves of the hot Jupiter HD189733b; it performes an Independent Component Analysis (ICA) on a set of pixel-light-curves, i.e. time series read by individual pixels, from the same photometric observation. Our method only assumes the independence of instrumental and astrophysical signals, and therefore guarantees a higher degree of objectivity compared to parametric detrending techniques published in the literature. The datasets we analyzed in this paper represent a more challenging test compared to the previous ones. Contrary to previous results reported in the literature, our results (1) do not support any detectable inter-epoch variations of orbital and stellar parameters, (2) are photometrically stable at the level 10e-4 in the IR, and (3) the transit depth measurements at the two wavelengths are consistent within 1σ. We also (4) detect a possible transit duration variation (TDV) of 80 s (2 σ significance level), that has not been pointed out in the literature, and (5) confirm no transit timing variations (TTVs) >30 s.

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 Gaia-ESO Survey: Discovery of a spatially extended low-mass population in the Vela OB2 association

The nearby (distance~350-400 pc), rich Vela OB2 association, includes $γ^2$ Velorum, one of the most massive binaries in the solar neighbourhood and an excellent laboratory for investigating the formation and early evolution of young clusters. Recent Gaia-ESO survey observations have led to the discovery of two kinematically distinct populations in the young (10-15 Myr) cluster immediately surrounding $γ^2$ Velorum. Here we analyse the results of Gaia-ESO survey observations of NGC 2547, a 35 Myr cluster located two degrees south of $γ^2$ Velorum. The radial velocity distribution of lithium-rich pre-main sequence stars shows a secondary population that is kinematically distinct from and younger than NGC 2547. The radial velocities, lithium absorption lines, and the positions in a colour-magnitude diagram of this secondary population are consistent with those of one of the components discovered around $γ^2$ Velorum. This result shows that there is a young, low-mass stellar population spread over at least several square degrees in the Vela OB2 association. This population could have originally been part of a cluster around $γ^2$ Velorum that expanded after gas expulsion or formed in a less dense environment that is spread over the whole Vela OB2 region.

preprint2015arXiv

The Gaia-ESO Survey: Empirical determination of the precision of stellar radial velocities and projected rotation velocities

The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (v sin i) for representative samples of Galactic stars, that will complement information obtained by the Gaia astrometry satellite. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and v sin i using spectra from repeated exposures of the same stars. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and v sin i, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student's t-distributions than by normal distributions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the v sin i precision for stars in young clusters, as a function of S/N, v sin i and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22-0.26 km/s, dependent on instrumental configuration.

preprint2015arXiv

The Gaia-ESO Survey: N-body modelling of the Gamma Velorum cluster

The Gaia-ESO Survey has recently unveiled the complex kinematic signature of the Gamma Velorum cluster: this cluster is composed of two kinematically distinct populations (hereafter, population A and B), showing two different velocity dispersions and a relative ~2 km s^-1 radial velocity (RV) shift. In this paper, we propose that the two populations of the Gamma Velorum cluster originate from two different sub-clusters, born from the same parent molecular cloud. We investigate this possibility by means of direct-summation N-body simulations. Our scenario is able to reproduce not only the RV shift and the different velocity dispersions, but also the different centroid (~0.5 pc), the different spatial concentration and the different line-of-sight distance (~5 pc) of the two populations. The observed 1-2 Myr age difference between the two populations is also naturally explained by our scenario, in which the two sub-clusters formed in two slightly different star formation episodes. Our simulations suggest that population B is strongly supervirial, while population A is close to virial equilibrium. We discuss the implications of our models for the formation of young star clusters and OB associations in the Milky Way.

preprint2015arXiv

The GAPS Programme with HARPS-N at TNG VIII: Observations of the Rossiter-McLaughlin effect and characterisation of the transiting planetary systems HAT-P-36 and WASP-11/HAT-P-10

We determine the true and the projected obliquity of HAT-P-36 and WASP-11/HAT-P-10 systems, respectively, which are both composed of a relatively cool star and a hot-Jupiter planet. Thanks to the high-resolution spectrograph HARPS-N, we observed the Rossiter-McLaughlin effect for both the systems by acquiring precise radial-velocity measurements during planetary transit events. We also present photometric observations comprising six light curves covering five transit events, obtained using three medium-class telescopes and the telescope-defocussing technique. One transit of WASP-11/HAT-P-10 was followed simultaneously from two observatories. The three transit light curves of HAT-P-36b show anomalies that are attributable to starspot complexes on the surface of the parent star, in agreement with the analysis of its spectra that indicate a moderate activity. By analysing the complete HATNet data set of HAT-P-36, we estimated the stellar rotation period by detecting a periodic photometric modulation in the light curve caused by star spots, obtaining Prot=15.3 days, which implies that the inclination of the stellar rotational axis with respect to the line of sight is 65 degree. We used the new spectroscopic and photometric data to revise the main physical parameters and measure the sky-projected misalignment angle of the two systems. We found λ=-14 degree for HAT-P-36 and λ=7 degree for WASP-11/HAT-P-10, indicating in both cases a good spin-orbit alignment. In the case of HAT-P-36, we also measured its real obliquity, which turned out to be 25 degrees.

preprint2015arXiv

The GAPS Programme with HARPS-N at TNG. VII. Putting exoplanets in the stellar context: magnetic activity and asteroseismology of $τ$ Bootis A

Aims. We observed the $τ$ Boo system with the HARPS-N spectrograph to test a new observational strategy aimed at jointly studying asteroseismology, the planetary orbit, and star-planet magnetic interaction. Methods. We collected high-cadence observations on 11 nearly consecutive nights and for each night averaged the raw FITS files using a dedicated software. In this way we obtained spectra with a high signal-to-noise ratio, used to study the variation of the CaII H&K lines and to have radial velocity values free from stellar oscillations, without losing the oscillations information. We developed a dedicated software to build a new custom mask that we used to refine the radial velocity determination with the HARPS-N pipeline and perform the spectroscopic analysis. Results. We updated the planetary ephemeris and showed the acceleration caused by the stellar binary companion. Our results on the stellar activity variation suggest the presence of a high-latitude plage during the time span of our observations. The correlation between the chromospheric activity and the planetary orbital phase remains unclear. Solar-like oscillations are detected in the radial velocity time series: we estimated asteroseismic quantities and found that they agree well with theoretical predictions. Our stellar model yields an age of $0.9\pm0.5$ Gyr for $τ$ Boo and further constrains the value of the stellar mass to $1.38\pm0.05$ M$_\odot$.

preprint2015arXiv

The GAPS Programme with HARPS-N at TNG. X. Differential abundances in the XO-2 planet hosting binary

Binary stars hosting exoplanets are a unique laboratory where chemical tagging can be performed to measure with high accuracy the elemental abundances of both stellar components, with the aim to investigate the formation of planets and their subsequent evolution. Here, we present a high-precision differential abundance analysis of the XO-2 wide stellar binary based on high resolution HARPS-N@TNG spectra. Both components are very similar K-dwarfs and host planets. Since they formed presumably within the same molecular cloud, we expect they should possess the same initial elemental abundances. We investigate if the presence of planets can cause some chemical imprints in the stellar atmospheric abundances. We measure abundances of 25 elements for both stars with a range of condensation temperature $T_{\rm C}=40-1741$ K, achieving typical precisions of $\sim 0.07$ dex. The North component shows abundances in all elements higher by $+0.067 \pm 0.032$ dex on average, with a mean difference of +0.078 dex for elements with $T_{\rm C} > 800$ K. The significance of the XO-2N abundance difference relative to XO-2S is at the $2σ$ level for almost all elements. We discuss the possibility that this result could be interpreted as the signature of the ingestion of material by XO-2N or depletion in XO-2S due to locking of heavy elements by the planetary companions. We estimate a mass of several tens of $M_{\oplus}$ in heavy elements. The difference in abundances between XO-2N and XO-2S shows a positive correlation with the condensation temperatures of the elements, with a slope of $(4.7 \pm 0.9) \times 10^{-5}$ dex K$^{-1}$, which could mean that both components have not formed terrestrial planets, but that first experienced the accretion of rocky core interior to the subsequent giant planets.

preprint2015arXiv

The GAPS Programme with HARPS-N@TNG IX. The multi-planet system KELT-6: detection of the planet KELT-6 c and measurement of the Rossiter-McLaughlin effect for KELT-6 b

Aims. For more than 1.5 years we monitored spectroscopically the star KELT-6 (BD+312447), known to host the transiting hot Saturn KELT-6b, because a previously observed long-term trend in radial velocity time series suggested the existence of an outer companion. Methods. We collected a total of 93 new spectra with the HARPS-N and TRES spectrographs. A spectroscopic transit of KELT-6b was observed with HARPS-N, and simultaneous photometry was obtained with the IAC-80 telescope. Results. We proved the existence of an outer planet with a mininum mass M$_{\rm p}$sini=3.71$\pm$0.21 M$_{\rm Jup}$ and a moderately eccentric orbit ($e=0.21_{-0.036}^{+0.039}$) of period P$\sim$3.5 years. We improved the orbital solution of KELT-6b and obtained the first measurement of the Rossiter-McLaughlin effect, showing that the planet has a likely circular, prograde, and slightly misaligned orbit, with a projected spin-orbit angle $λ$=$-$36$\pm$11 degrees. We improved the KELT-6b transit ephemeris from photometry, and we provided new measurements of the stellar parameters. KELT-6 appears as an interesting case to study the formation and evolution of multi-planet systems.

preprint2015arXiv

The GAPS Programme with HARPS-N@TNG V. A comprehensive analysis of the XO-2 stellar and planetary systems

XO-2 is the first confirmed wide stellar binary system where the almost twin components XO-2N and XO-2S have planets. This stimulated a detailed characterization study of the stellar and planetary components based on new observations. We collected high-resolution spectra with the HARPS-N spectrograph and multi-band light curves. Spectral analysis led to an accurate determination of the stellar atmospheric parameters and characterization of the stellar activity. We collected 14 transit light curves of XO-2Nb used to improve the transit parameters. Photometry provided accurate magnitude differences between the stars and a measure of their rotation periods. The iron abundance of XO-2N was found to be +0.054 dex greater, within more than 3-sigma, than that of XO-2S. We confirm a long-term variation in the radial velocities of XO-2N, and we detected a turn-over with respect to previous measurements. We suggest the presence of a second massive companion in an outer orbit or the stellar activity cycle as possible causes of the observed acceleration. The latter explanation seems more plausible with the present dataset. We obtained an accurate value of the projected spin-orbit angle for the XO-2N system (lambda=7+/-11 degrees), and estimated the real 3-D spin-orbit angle (psi=27 +12/-27 degrees). We measured the XO-2 rotation periods, and found a value of P=41.6 days in the case of XO-2N, in excellent agreement with the predictions. The period of XO-2S appears shorter, with an ambiguity between 26 and 34.5 days that we cannot solve with the present dataset alone. XO-2N appears to be more active than the companion, and this could be due to the fact that we sampled different phases of their activity cycle, or to an interaction between XO-2N and its hot Jupiter that we could not confirm.

preprint2015arXiv

The GAPS Programme with HARPS-N@TNG VI: The Curious Case of TrES-4b

We revisit the TrES-4 system parameters based on high-precision HARPS-N radial-velocity measurements and new photometric light curves. A combined spectroscopic and photometric analysis allows us to determine a spectroscopic orbit with an amplitude $K=51\pm3$ m s$^{-1}$. The derived mass of TrES-4b is found to be $M_{\rm p} = 0.49\pm0.04 \rm M_{Jup}$, significantly lower than previously reported. Combined with the large radius ($R_{\rm p} = 1.84_{-0.09}^{+0.08} \rm R_{Jup}$) inferred from our analysis, TrES-4b becomes the second-lowest density transiting hot Jupiter known. We discuss several scenarios to explain the puzzling discrepancy in the mass of TrES-4b in the context of the exotic class of highly inflated transiting giant planets.

preprint2015arXiv

The HARPS-N Rocky Planet Search I. HD219134b: A transiting rocky planet in a multi-planet system at 6.5 pc from the Sun

We present here the detection of a system of four low-mass planets around the bright (V=5.5) and close-by (6.5 pc) star HD219134. This is the first result of the Rocky Planet Search program with HARPS-N on the TNG in La Palma. The inner planet orbits the star in 3.0937 +/-0.0004 days, on a quasi-circular orbit with a semi-major axis of 0.0382 +/- 0.0003 AU. Spitzer observations allowed us to detect the transit of the planet in front of the star making HD219134b the nearest known transiting planet to date. From the amplitude of the radial-velocity variation (2.33 +/- 0.24 m/s) and observed depth of the transit (359 +/- 38 ppm), the planet mass and radius are estimated to be 4.46 +/- 0.47 M_{\oplus} and 1.606 +/- 0.086 R_{\oplus} leading to a mean density of 5.89 +/- 1.17 g/cc, suggesting a rocky composition. One additional planet with minimum mass of 2.67 +/- 0.59 M_{\oplus} moves on a close-in, quasi-circular orbit with a period of 6.765 +/- 0.005 days. The third planet in the system has a period of 46.78 +/- 0.16 days and a minimum mass of 8.7 +/- 1.1 M{\oplus}, at 0.234 +/- 0.002 AU from the star. Its eccentricity is 0.32 +/- 0.14. The period of this planet is close to the rotational period of the star estimated from variations of activity indicators (42.3 +/- 0.1 days). The planetary origin of the signal is, however, the preferred solution as no indication of variation at the corresponding frequency is observed for activity-sensitive parameters. Finally, a fourth additional longer-period planet of mass of 62 +/- 6 M_{\oplus} orbits the star in 1190 days, on an eccentric orbit (e=0.27 +/- 0.11) at a distance of 2.14 +/- 0.27 AU.

preprint2014arXiv

Characterization of the Kepler-101 planetary system with HARPS-N. A hot super-Neptune with an Earth-sized low-mass companion

We report on the characterization of the Kepler-101 planetary system, thanks to a combined DE-MCMC analysis of Kepler data and forty radial velocities obtained with the HARPS-N spectrograph. This system was previously validated by Rowe et al. (2014) and is composed of a hot super-Neptune, Kepler-101b, and an Earth-sized planet, Kepler-101c. These two planets orbit the slightly evolved and metal-rich G-type star in 3.49 and 6.03 days, respectively. With mass $M_{\rm p}=51.1_{-4.7}^{+5.1}~M_{\oplus}$, radius $R_{\rm p}=5.77_{-0.79}^{+0.85}~R_{\oplus}$, and density $ρ_{\rm p}=1.45_{-0.48}^{+0.83} \rm g\;cm^{-3}$, Kepler-101b is the first fully-characterized super-Neptune, and its density suggests that heavy elements make up a significant fraction of its interior; more than $60\%$ of its total mass. Kepler-101c has a radius of $1.25_{-0.17}^{+0.19}~R_{\oplus}$, which implies the absence of any H/He envelope, but its mass could not be determined due to the relative faintness of the parent star for highly precise radial-velocity measurements ($K_{\rm p}=13.8$) and the limited number of radial velocities. The $1~σ$ upper limit, $M_{\rm p} < 3.8~M_{\oplus}$, excludes a pure iron composition with a $68.3\%$ probability. The architecture of the Kepler-101 planetary system - containing a close-in giant planet and an outer Earth-sized planet with a period ratio slightly larger than the 3:2 resonance - is certainly of interest for planet formation and evolution scenarios. This system does not follow the trend, seen by Ciardi et al. (2013), that in the majority of Kepler systems of planet pairs with at least one Neptune-size or larger planet, the larger planet has the longer period.

preprint2014arXiv

EChO Payload electronics architecture and SW design

EChO is a three-modules (VNIR, SWIR, MWIR), highly integrated spectrometer, covering the wavelength range from 0.55 $μ$m, to 11.0 $μ$m. The baseline design includes the goal wavelength extension to 0.4 $μ$m while an optional LWIR module extends the range to the goal wavelength of 16.0 $μ$m. An Instrument Control Unit (ICU) is foreseen as the main electronic subsystem interfacing the spacecraft and collecting data from all the payload spectrometers modules. ICU is in charge of two main tasks: the overall payload control (Instrument Control Function) and the housekeepings and scientific data digital processing (Data Processing Function), including the lossless compression prior to store the science data to the Solid State Mass Memory of the Spacecraft. These two main tasks are accomplished thanks to the Payload On Board Software (P-OBSW) running on the ICU CPUs.

preprint2014arXiv

EChO spectra and stellar activity II. The case of dM stars

EChO is a dedicated mission to investigate exoplanetary atmospheres. When extracting the planetary signal, one has to take care of the variability of the hosting star, which introduces spectral distortion that can be mistaken as planetary signal. Magnetic variability is a major deal in particular for M stars. To this purpose, assuming a one spot dominant model for the stellar photosphere, we develop a mixed observational-theoretical tool to extract the spot's parameters from the observed optical spectrum. This method relies on a robust library of spectral M templates, which we derive using the observed spectra of quiet M dwarfs in the SDSS database. Our procedure allows to correct the observed spectra for photospheric activity in most of the analyzed cases, reducing the spectral distortion down to the noise levels. Ongoing refinements of the template library and the algorithm will improve the efficiency of our algorithm.

preprint2014arXiv

Effects of XUV radiation on circumbinary planets

Several circumbinary planets have recently been discovered. The orbit of a planet around a binary stellar system poses several dynamic constraints. The effects that radiation from the host stars may have on the planet atmospheres must be considered. Because of the configuration of a close binary system, these stars have a high rotation rate, which causes a permanent state of high stellar activity and copious XUV radiation. The accumulated effects are stronger than for exoplanets around single stars, and cause a faster evaporation of their atmospheres. We evaluate the effects that stellar radiation has on the evaporation of exoplanets around binary systems and on the survival of these planets. We considered the XUV spectral range to account for the photons that are easily absorbed by a planet atmosphere that is mainly composed of hydrogen. A more complex atmospheric composition is expected to absorb this radiation more efficiently. We used direct X-ray observations to evaluate the energy in the X-rays range and coronal models to calculate the (nondetectable) EUV part of the spectrum. The simulations show that exoplanets in a close orbit will suffer strong photoevaporation that may cause a total loss of atmosphere in a short time. A binary system of two solar-like stars will be highly efficient in evaporating the atmosphere of the planet. These systems will be difficult to find, even if they are dynamically stable. Still, planets may orbit around binary systems of low mass stars for wider orbits. Currently known circumbinary planets are not substantially affected by thermal photoevaporation processes, unless Kepler-47 b has an inflated atmosphere. The distribution of the orbital periods of circumbinary planets is shifted to much longer periods than the average of Kepler planets, which supports a scenario of strong photoevaporation in close-in circumbinary planets.

preprint2014arXiv

Empirical classification of VLT/Giraffe stellar spectra in the wavelength range 6440-6810 A in the gamma Vel cluster, and calibration of spectral indices

We study spectral diagnostics available from optical spectra with R=17000 obtained with the VLT/Giraffe HR15n setup, using observations from the Gaia-ESO Survey, on the gamma Vel young cluster, in order to determine the fundamental parameters of these stars. We define a set of spectroscopic indices, sampling TiO bands, H-alpha core and wings, and many temperature- and gravity-sensitive lines. Combined indices tau (gamma) are also defined as Teff (log g) indicators over a wide spectral-type range. H-alpha emission-line indices are also chromospheric activity or accretion indicators. A metallicity-sensitive index is also defined. These indices enable us to find a clear difference between gravities of main-sequence and pre-main-sequence stars (as well as giant stars): the (gamma,tau) diagram is thus argued to be a promising distance-independent age measurement tool for young clusters. Our indices were quantitatively calibrated by means of photometry and literature reference spectra (from UVES-POP and ELODIE 3.1 Libraries, and other datasets), over a wide range of stellar parameters. The ability of our indices to select peculiar, or other rare classes of stars is also established. Finally, our gravity determinations support the existence of an older pre-main-sequence population in the gamma Vel sky region, in agreement with the lithium depletion pattern of those stars (abridged).

preprint2014arXiv

Gaia-ESO Survey: The analysis of high-resolution UVES spectra of FGK-type stars

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.

preprint2014arXiv

The Gaia-ESO Survey: Chromospheric Emission, Accretion Properties, and Rotation in $γ$ Velorum and Chamaeleon I

We use the fundamental parameters delivered by the GES consortium in the first internal data release to select the members of $γ$ Vel and Cha I among the UVES and GIRAFFE spectroscopic observations. A total of 140 $γ$ Vel members and 74 Cha I members were studied. We calculated stellar luminosities through spectral energy distributions, while stellar masses were derived by comparison with evolutionary tracks. The spectral subtraction of low-activity and slowly rotating templates, which are rotationally broadened to match the $v\sin i$ of the targets, enabled us to measure the equivalent widths (EWs) and the fluxes in the H$α$ and H$β$ lines. The H$α$ line was also used for identifying accreting objects and for evaluating the mass accretion rate ($\dot M_{\rm acc}$). The distribution of $v\sin i$ for the members of $γ$ Vel displays a peak at about 10 km s$^{-1}$ with a tail toward faster rotators. There is also some indication of a different $v\sin i$ distribution for the members of its two kinematical populations. Only a handful of stars in $γ$ Vel display signatures of accretion, while many more accretors were detected in the younger Cha~I. Accreting and active stars occupy two different regions in a $T_{\rm eff}$-flux diagram and we propose a criterion for distinguishing them. We derive $\dot M_{\rm acc}$ in the ranges $10^{-11}$-$10^{-9} M_\odot$yr$^{-1}$ and $10^{-10}$-$10^{-7} M_\odot$yr$^{-1}$ for $γ$ Vel and Cha I accretors, respectively. We find less scatter in the $\dot M_{\rm acc}-M_\star$ relation derived through the H$α$ EWs, when compared to the H$α$ $10\%W$ diagnostics, in agreement with other authors.

preprint2014arXiv

The Gaia-ESO Survey: Kinematic structure in the Gamma Velorum Cluster

Context: A key science goal of the Gaia-ESO survey (GES) is to use the kinematics of low-mass stars in young clusters to probe their dynamical histories and how they populate the field as they become unbound. The clustering of low-mass stars around the massive W-R binary gamma^2 Velorum was one of the first GES targets. Aims: To empirically determine the radial velocity (RV) precision of GES data, construct a kinematically unbiased sample of cluster members and characterise their dynamical state. Methods: Targets were selected from colour-magnitude diagrams and intermediate resolution spectroscopy used to derive RVs and assess membership from the strength of the Li6708A line. The RV distribution was analysed using a maximum likelihood technique that accounts for unresolved binaries. Results: The GES RV precision is about 0.25km/s and sufficient to resolve velocity structure in the low-mass population around gamma^2 Vel. The structure is well fitted by two kinematic components with roughly equal numbers of stars; the first has an intrinsic dispersion of 0.34+/-0.16km/s, consistent with virial equilibrium. The second has a broader dispersion of 1.60+/-0.37km/s and is offset from the first by ~2km/s. The first population is older by 1-2Myr based on a greater level of Li depletion seen among its M-stars and is probably more centrally concentrated around gamma^2 Vel. Conclusions: We consider several formation scenarios, concluding that the two kinematic components are a bound remnant of the original, denser cluster that formed gamma^2 Vel, and a dispersed population from the wider Vela OB2 association, of which gamma^2 Vel is the most massive member. The apparent youth of gamma^2 Vel compared to the older (>=10Myr) low-mass population surrounding it suggests a scenario where the massive binary formed in a clustered environment after the formation of the bulk of the low-mass stars.[abridged]

preprint2014arXiv

The Gaia-ESO Survey: processing of the FLAMES-UVES spectra

The Gaia-ESO Survey is a large public spectroscopic survey that aims to derive radial velocities and fundamental parameters of about 10^5 Milky Way stars in the field and in clusters. Observations are carried out with the multi-object optical spectrograph FLAMES, using simultaneously the medium resolution (R~20,000) GIRAFFE spectrograph and the high resolution (R~47,000) UVES spectrograph. In this paper, we describe the methods and the software used for the data reduction, the derivation of the radial velocities, and the quality control of the FLAMES-UVES spectra. Data reduction has been performed using a workflow specifically developed for this project. This workflow runs the ESO public pipeline optimizing the data reduction for the Gaia-ESO Survey, performs automatically sky subtraction, barycentric correction and normalisation, and calculates radial velocities and a first guess of the rotational velocities. The quality control is performed using the output parameters from the ESO pipeline, by a visual inspection of the spectra and by the analysis of the signal-to-noise ratio of the spectra. Using the observations of the first 18 months, specifically targets observed multiple times at different epochs, stars observed with both GIRAFFE and UVES, and observations of radial velocity standards, we estimated the precision and the accuracy of the radial velocities. The statistical error on the radial velocities is sigma~0.4 km s^-1 and is mainly due to uncertainties in the zero point of the wavelength calibration. However, we found a systematic bias with respect to the GIRAFFE spectra (~0.9 km s^-1) and to the radial velocities of the standard stars (~0.5 kms^-1) retrieved from the literature. This bias will be corrected in the future data releases, when a common zero point for all the setups and instruments used for the survey will be established.

preprint2014arXiv

The Gaia-ESO Survey: radial metallicity gradients and age-metallicity relation of stars in the Milky Way disk

We study the relationship between age, metallicity, and alpha-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpc to 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages > 9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more alpha-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

preprint2014arXiv

The Gaia-ESO Survey: the chemical structure of the Galactic discs from the first internal data release

Most high-resolution spectroscopic studies of the Galactic discs were mostly confined to objects in the solar vicinity. Here we aim at enlarging the volume in which individual chemical abundances are used to characterise both discs, using the first internal data release of the Gaia-ESO survey. We derive and discuss the abundances of eight elements (Mg, Al, Si, Ca, Ti, Fe, Cr, Ni, and Y). The trends of these elemental abundances with iron are very similar to those in the solar neighbourhood. We find a natural division between alpha-rich and alpha-poor stars, best seen in the bimodality of the [Mg/M] distributions in bins of metallicity, which we attribute to thick- and thin-disc sequences, respectively. With the possible exception of Al, the observed dispersion around the trends is well described by the expected errors, leaving little room for astrophysical dispersion. Using previously derived distances from Recio-Blanco et al. (2014b), we further find that the thick-disc is more extended vertically and is more centrally concentrated towards the inner Galaxy than the thin-disc, which indicates a shorter scale-length. We derive the radial and vertical gradients in metallicity, iron, four alpha-element abundances, and Al for the two populations, taking into account the identified correlation between R_GC and |Z|. Radial metallicity gradient is found in the thin disc. The positive radial individual [alpha/M] gradients found are at variance from the gradients observed in the RAVE survey. The thin disc also hosts a negative vertical metallicity gradient, accompanied by positive individual [alpha/M] and [Al/M] gradients. The thick-disc, presents no radial metallicity gradient, a shallower vertical metallicity gradient than the thin-disc, an alpha-elements-to-iron radial gradient in the opposite sense than that of the thin disc, and positive vertical individual [alpha/M] and [Al/M] gradients.

preprint2014arXiv

The Gaia-ESO Survey: the first abundance determination of the pre-main-sequence cluster Gamma Velorum

Knowledge of the abundance distribution of star forming regions and young clusters is critical to investigate a variety of issues, from triggered star formation and chemical enrichment by nearby supernova explosions to the ability to form planetary systems.In spite of this, detailed abundance studies are currently available for relatively few regions. In this context, we present the analysis of the metallicity of the Gamma Velorum cluster, based on the products distributed in the first internal release of the Gaia-ESO Survey. The Gamma Velorum candidate members have been observed with FLAMES, using both UVES and Giraffe, depending on the target brightness and spectral type. In order to derive a solid metallicity determination for the cluster, membership of the observed stars must be first assessed. To this aim, we use several membership criteria including radial velocities, surface gravity estimates, and the detection of the photospheric lithium line. Out of the 80 targets observed with UVES, we identify 14 high-probability members. We find that the metallicity of the cluster is slightly subsolar, with a mean [Fe/H]=-0.057+/-0.018 dex. Although J08095427-4721419 is one of the high-probability members, its metallicity is significantly larger than the cluster average. We speculate about its origin as the result of recent accretion episodes of rocky bodies of ~60 M_Sun hydrogen-depleted material from the circumstellar disk.

preprint2014arXiv

The Gaia-ESO Survey: the Galactic Thick to Thin Disc transition

(Abridged) We have used the atmospheric parameters, [alpha/Fe] abundances and radial velocities, determined from the Gaia-ESO Survey GIRAFFE spectra of FGK-type stars (iDR1), to provide a chemo-kinematical characterisation of the disc stellar populations. We focuss on a subsample of 1016 stars with high quality parameters, covering the volume |Z|<4.5kpc and R in the range 2-13kpc. We have identified a thin to thick disc separation in the [alpha/Fe] vs [M/H] plane, thanks to the presence of a low-density region in the number density distribution. The thick disc stars seem to lie in progressively thinner layers above the Galactic plane, as metallicity increases and [alpha/Fe] decreases. The thin disc population presents a constant value of the mean distance to the plane at all metallicities. Our data confirm the already known correlations between V_phi and [M/H] for the two discs. For the thick disc sequence, a study of the possible contamination by thin disc stars suggests a gradient up to 64km/s/dex. The distributions of V_phi, V_Z, and orbital parameters are analysed for the chemically separated samples. Concerning the gradients with galactocentric radius, we find for the thin disc a flat behaviour of V_phi, a [M/H] gradient of -0.058dex/kpc and a small positive [alpha/Fe] gradient. For the thick disc, flat gradients in [M/H] and [alpha/Fe] are derived. Our chemo-kinematical analysis suggests a picture in which the thick disc seems to have experienced a settling process, during which its rotation increased progressively, and, possibly, the V_phi dispersion decreased. At [M/H]-0.25dex and [alpha/Fe]0.1dex, the mean characteristics of the thick disc in distance to the Galactic plane, V_phi, V_phi dispersion and eccentricity agree with those of the thin disc stars, suggesting a possible connection between these populations at a certain epoch of the disc evolution.

preprint2014arXiv

The GAPS Programme with HARPS-N at TNG. III: The retrograde orbit of HAT-P-18b

The measurement of the Rossiter-McLaughlin effect for transiting exoplanets places constraints on the orientation of the orbital axis with respect to the stellar spin axis, which can shed light on the mechanisms shaping the orbital configuration of planetary systems. Here we present the interesting case of the Saturn-mass planet HAT-P-18b, which orbits one of the coolest stars for which the Rossiter-McLaughlin effect has been measured so far. We acquired a spectroscopic time-series, spanning a full transit, with the HARPS-N spectrograph mounted at the TNG telescope. The very precise radial velocity measurements delivered by the HARPS-N pipeline were used to measure the Rossiter-McLaughlin effect. Complementary new photometric observations of another full transit were also analysed to obtain an independent determination of the star and planet parameters. We find that HAT-P-18b lies on a counter-rotating orbit, the sky-projected angle between the stellar spin axis and the planet orbital axis being lambda=132 +/- 15 deg. By joint modelling of the radial velocity and photometric data we obtain new determinations of the star (M_star = 0.770 +/- 0.027 M_Sun; R_star= 0.717 +/- 0.026 R_Sun; Vsin(I_star) = 1.58 +/- 0.18 km/s) and planet (M_pl = 0.196 +/- 0.008 M_J; R_pl = 0.947 +/- 0.044 R_J) parameters. Our spectra provide for the host star an effective temperature T_eff = 4870 +/- 50 K, a surface gravity of log(g_star) = 4.57 +/- 0.07 cm/s, and an iron abundance of [Fe/H] = 0.10 +/- 0.06. HAT-P-18b is one of the few planets known to transit a star with T_eff < 6250 K on a retrograde orbit. Objects such as HAT-P-18b (low planet mass and/or relatively long orbital period) most likely have a weak tidal coupling with their parent stars, therefore their orbits preserve any original misalignment. As such, they are ideal targets to study the causes of orbital evolution in cool main-sequence stars.

preprint2014arXiv

The GAPS programme with HARPS-N@TNG IV: A planetary system around XO-2S

We performed an intensive radial velocity monitoring of XO-2S, the wide companion of the transiting planet-host XO-2N, using HARPS-N at TNG in the framework of the GAPS programme. The radial velocity measurements indicate the presence of a new planetary system formed by a planet that is slightly more massive than Jupiter at 0.48 au and a Saturn-mass planet at 0.13 au. Both planetary orbits are moderately eccentric and were found to be dynamically stable. There are also indications of a long-term trend in the radial velocities. This is the first confirmed case of a wide binary whose components both host planets, one of which is transiting, which makes the XO-2 system a unique laboratory for understanding the diversity of planetary systems.

preprint2014arXiv

The Visible and Near Infrared module of EChO

The Visible and Near Infrared (VNIR) is one of the modules of EChO, the Exoplanets Characterization Observatory proposed to ESA for an M-class mission. EChO is aimed to observe planets while transiting by their suns. Then the instrument had to be designed to assure a high efficiency over the whole spectral range. In fact, it has to be able to observe stars with an apparent magnitude Mv= 9-12 and to see contrasts of the order of 10-4 - 10-5 necessary to reveal the characteristics of the atmospheres of the exoplanets under investigation. VNIR is a spectrometer in a cross-dispersed configuration, covering the 0.4-2.5 micron spectral range with a resolving power of about 330 and a field of view of 2 arcsec. It is functionally split into two channels respectively working in the 0.4-1 and 1.0-2.5 micron spectral ranges. Such a solution is imposed by the fact the light at short wavelengths has to be shared with the EChO Fine Guiding System (FGS) devoted to the pointing of the stars under observation. The spectrometer makes use of a HgCdTe detector of 512 by 512 pixels, 18 micron pitch and working at a temperature of 45K as the entire VNIR optical bench. The instrument has been interfaced to the telescope optics by two optical fibers, one per channel, to assure an easier coupling and an easier colocation of the instrument inside the EChO optical bench.

preprint2014arXiv

Young Stellar Object Variability (YSOVAR): Long Timescale Variations in the Mid-Infrared

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 & 4.5 um) time-series photometry of the Orion Nebula Cluster plus smaller footprints in eleven other star-forming cores (AFGL490, NGC1333, MonR2, GGD 12-15, NGC2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs about twice per day for ~40d. We also define a "standard sample of members", consisting of all the IR-selected members and X-ray selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data--the Stetson index, a chi^2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of ~6 years, by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data; out of members and field objects combined, at most 0.02% may have transient IR excesses.

preprint2013arXiv

A multiwavelength view of star-disk interaction in NGC 2264

Variability is a signature property of cool young stars, particularly for those surrounded by disks. Traditional single-band time series display complex features associated with accretion, disk structure, and accompanying stellar activity, but these processes are challenging to model. To make progress in connecting observed time domain properties with the underlying physics of young stars and their disks, we have embarked on an unprecedented multiwavelength monitoring campaign: the Coordinated Synoptic Investigation of NGC 2264 ("CSI 2264"). Beginning in December 2011, CSI 2264 has acquired 30 continuous days of mid-infrared time series from Spitzer, simultaneous optical monitoring from CoRoT and MOST, X-ray observations with Chandra, as well as complementary data from a number of ground-based telescopes. The extraordinary photometric precision, cadence, and time baseline of these observations enable detailed correlation of variability properties at different wavelengths, corresponding to locations from the stellar surface to the inner AU of the disk. We present the early results of the program, and discuss the need for further modeling efforts into young stars and their disks.

preprint2013arXiv

An X-rays Survey of the Young Stellar Population of the Lynds 1641 and Iota Orionis Regions

We present an XMM-Newton survey of the part of Orion A cloud south of the Orion Nebula. This survey includes the Lynds 1641 (L1641) dark cloud, a region of the Orion A cloud with very few massive stars and hence a relatively low ambient UV flux, and the region around the O9 III star Iota Orionis. In addition to proprietary data, we used archival XMM data of the Orion Nebula Cluster (ONC) to extend our analysis to a major fraction of the Orion A cloud. We have detected 1060 X-ray sources in L1641 and Iota Ori region. About 94% of the sources have 2MASS & Spitzer counterparts, 204 and 23 being Class II and Class I or protostars objects, respectively. In addition, we have identified 489 X-ray sources as counterparts to Class III candidates, given they are bright in X-rays and appear as normal photospheres at mid-IR wavelengths. The remaining 205 X-ray sources are likely distant AGNs or other galactic sources not related to Orion A. We find that Class III candidates appear more concentrated in two main clusters in L1641. The first cluster of Class III stars is found toward the northern part of L1641, concentrated around Iota Ori. The stars in this cluster are more evolved than those in the Orion Nebula. We estimate a distance of 300-320 pc for this cluster and thus it is closer than the Orion A cloud. Another cluster rich in Class III stars is located in L1641 South and appears to be a slightly older cluster embedded in the Orion A cloud. Furthermore, other evolved Class III stars are found north of the ONC toward NGC 1977.

preprint2013arXiv

Astrometric Detection of Giant Planets Around Nearby M Dwarfs: The Gaia Potential

[abridged] We carry out numerical simulations to gauge the Gaia potential for precision astrometry of exoplanets orbiting a sample of known dM stars within 30 pc from the Sun. (1) It will be possible to accurately determine orbits and masses for Jupiter-mass planets with orbital periods in the range 0.2<P<6.0 yr and with an astrometric signal-to-noise ratio > 10. Given present-day estimates of the planet fraction f_p around M dwarfs, 100 giant planets could be found by Gaia around the sample. Comprehensive screening by Gaia of the reservoir of 4x10^5 M dwarfs within 100 pc could result in 2600 detections and as many as 500 accurate orbit determinations. The value of f_p could then be determined with an accuracy of 2%, an improvement by over an order of magnitude with respect to the most precise values available to-date; (2) in the same period range, inclination angles corresponding to quasi-edge-on configurations will be determined with enough precision (a few percent) so that it will be possible to identify intermediate-separation planets which are potentially transiting within the errors. Gaia could alert us of the existence of 10 such systems. More than 250 candidates could be identified assuming solutions compatible with transit configurations within 10% accuracy, although a large fraction of these (85%) could be false positives; (3) for well-sampled orbits, the uncertainties on planetary ephemerides, separation and position angle, will degrade at typical rates of < 1 mas/yr and < 2 deg/yr, respectively; (4) planetary phases will be measured with typical uncertainties of several degrees, resulting (under the assumption of purely scattering atmospheres) in phase-averaged errors on the phase function of 0.05, and expected uncertainties in the determination of the emergent flux of intermediate-separation (0.3<a<2.0 AU) giant planets of 20%. [abridged]

preprint2013arXiv

Monitoring young associations and open clusters with Kepler in two-wheel mode

We outline a proposal to use the Kepler spacecraft in two-wheel mode to monitor a handful of young associations and open clusters, for a few weeks each. Judging from the experience of similar projects using ground-based telescopes and the CoRoT spacecraft, this program would transform our understanding of early stellar evolution through the study of pulsations, rotation, activity, the detection and characterisation of eclipsing binaries, and the possible detection of transiting exoplanets. Importantly, Kepler's wide field-of-view would enable key spatially extended, nearby regions to be monitored in their entirety for the first time, and the proposed observations would exploit unique synergies with the GAIA ESO spectroscopic survey and, in the longer term, the GAIA mission itself. We also outline possible strategies for optimising the photometric performance of Kepler in two-wheel mode by modelling pixel sensitivity variations and other systematics.

preprint2013arXiv

Spectroscopic observations of blue stars with infrared excesses in NGC 6611

The young open cluster NGC 6611 includes among its candidate members a class of peculiar objects with interesting properties: blue stars with infrared IR excesses. These stars show excesses in IR bands, signature of the presence of a circumstellar disk, but optical colors typical of older field stars. In order to confirm their membership to the cluster, it is therefore important to use new spectroscopic observations, together with previous photometric data. We aim at confirming the membership of these objects and at investigating their physical properties to verify whether the observed colors are intrinsic or altered by the presence of the disk or by the accretion processes. We analyze the intermediate resolution spectroscopic data obtained for a subsample of blue stars in NGC 6611 with FLAMES. In particular, we focus on the study of: 1) the profile of the Ha emission line, to select stars with accretion and outflow activity; 2) the Li absorption line, used as a youth indicator; 3) the radial velocity. Using the spectroscopic analysis, it has been possible to investigate the presence of the Li absorption line as well as to discriminate between stars with inert or active disk. In particular, from the analysis of the Ha emission line we were able to infer the activity due to the accretion and outflow processes and the variability of the emission. We also investigated the binarity of the blue stars and their membership to NGC 6611.

preprint2013arXiv

The GAPS programme with HARPS-N at TNG. I: Observations of the Rossiter-McLaughlin effect and characterisation of the transiting system Qatar-1

A long-term multi-purpose observational programme has started with HARPS-N@TNG aimed to characterise the global architectural properties of exoplanetary systems. In this first paper we fully characterise the transiting system Qatar-1. We exploit HARPS-N high-precision radial velocity measurements obtained during a transit to measure the Rossiter-McLaughlin effect in the Qatar-1 system, and out-of-transit measurements to redetermine the spectroscopic orbit. New photometric transit light-curves are analysed and a spectroscopic characterisation of the host star atmospheric parameters is performed based on various methods (line equivalent width ratios, spectral synthesis, spectral energy distribution). We achieved a significant improvement in the accuracy of the orbital parameters and derived the spin-orbit alignment of the system; this information, combined with the spectroscopic determination of the host star properties, allows us to derive the fundamental physical parameters for star and planet (masses and radii). The orbital solution for the Qatar-1 system is consistent with a circular orbit and the system presents a sky-projected obliquity of lambda = -8.4+-7.1 deg. The planet, with a mass of 1.33+-0.05 M_J, is found to be significantly more massive than previously reported. The host star is confirmed to be metal-rich ([Fe/H]= 0.20+-0.10) and slowly rotating (vsinI = 1.7+-0.3 km/s), though moderately active, as indicated by strong chromospheric emission in the Ca II H&K line cores (logR'_HK about -4.60). The system is well aligned and fits well within the general lambda vs Teff trend. We definitely rule out any significant orbital eccentricity. The evolutionary status of the system is inferred based on gyrochronology, and the present orbital configuration and timescale for orbital decay are discussed in terms of star-planet tidal interactions.

preprint2013arXiv

The GAPS Programme with HARPS-N@TNG II: No giant planets around the metal-poor star HIP 11952

In the context of the program Global Architecture of Planetary Systems (GAPS), we have performed radial velocity monitoring of the metal-poor star HIP 11952 on 35 nights over about 150 days using the newly installed high resolution spectrograph HARPS-N at the TNG and HARPS at ESO 3.6m telescope. The radial velocities show a scatter of 7 m/s, compatible with the measurement errors for such a moderately warm metal-poor star (Teff = 6040+-120K; [Fe/H] =-1.9+-0.1). We then exclude the presence of the two giant planets with periods of 6.95+-0.01 d and 290.0+-16.2 d and radial velocity semi-amplitudes of 100.3+-19.4 m/s and 105.2+-14.7 m/s, respectively, which had recently been announced. This result is important considering that HIP 11952 was thought to be the most metal-poor star hosting a planetary system with giant planets, thus challenging some models of planet formation.

preprint2013arXiv

The Hot and Energetic Universe: Solar system and exoplanets

The high resolution non-dispersive spectroscopy and unprecedented sensitivity of Athena+ will revolutionize solar system observing: the origin of the ions producing Jupiter's X-ray aurorae via charge exchange will be conclusively established, as well as their dynamics, giving clues to their acceleration mechanisms. X-ray aurorae on Saturn will be searched for to a depth unattainable by current Earth-bound observatories. The X-ray Integral Field Unit of Athena+ will map Mars' expanding exosphere, which has a line-rich solar wind charge exchange spectrum, under differing solar wind conditions and through the seasons; relating Mars' X-ray emission to its atmospheric loss will have significant impact also on the study of exoplanet atmospheres. Spectral mapping of cometary comae, which are spectacular X-ray sources with extremely line-rich spectra, will probe solar wind composition and speed at varying distances from the Sun. Athena+ will provide unique contributions also to exoplanetary astrophysics. Athena+ will pioneer the study of ingress/eclipse/egress effects during planetary orbits of hot-Jupiters, and will confirm/improve the evidence of Star-Planet Interactions (SPI) in a wider sample of planetary systems. Finally Athena+ will drastically improve the knowledge of the X-ray incident radiation on exoplanets, a key element for understanding the effects of atmospheric mass loss and of the chemical and physical evolution of planet atmospheres, particularly relevant in the case of young systems.

preprint2013arXiv

The UV and X-ray activity of the M dwarfs within 10pc of the Sun

We systematically study the X-ray and ultraviolet emission of a subsample of M dwarfs from a recent proper-motion survey, selecting all M dwarfs within 10pc to obtain a nearly volume-limited sample (~90% completeness). Archival ROSAT, XMM-Newton and GALEX data are combined with published spectroscopic studies of Halpha emission and rotation to obtain a broad picture of stellar activity on M dwarfs. We make use of synthetic model spectra to determine the relative contributions of photospheric and chromospheric emission to the ultraviolet flux. We also analyse the same diagnostics for a comparison sample of young M dwarfs in the TWHya association (~10Myrs). We find that generally the emission in the GALEX bands is dominated by the chromosphere but the photospheric component is not negligible in early-M field dwarfs. The surface fluxes for the Halpha, near-ultraviolet, far-ultraviolet and X-ray emission are connected via a power law dependence. We present here for the first time such flux-flux relations involving broad-band ultraviolet emission for M dwarfs. For given spectral type the activity indices, defined as flux ratio between the activity diagnostic and the bolometric flux of the star, display a spread of 2-3 dex which is largest for M4 stars. The mean activity index for fast rotators, likely representing the saturation level, decreases from X-rays over the FUV to the NUV band and Halpha, i.e. the fractional radiation output increases with atmospheric height. The comparison to the ultraviolet and X-ray properties of TWHya members shows a drop of nearly three orders of magnitude for the luminosity in these bands between ~10Myr and few Gyrs age. A few young field dwarfs (< 1Gyr) in the 10pc sample bridge the gap indicating that the drop in magnetic activity with age is a continuous process. The slope of the age decay is steeper for the X-ray than for the UV luminosity.

preprint2012arXiv

Chandra/ACIS-I study of the X-ray properties of the NGC 6611 and M16 stellar population

Mechanisms regulating the origin of X-rays in YSOs and the correlation with their evolutionary stage are under debate. Studies of the X-ray properties in young clusters allow to understand these mechanisms. One ideal target for this analysis is the Eagle Nebula (M16), with its central cluster NGC6611. At 1750 pc from the Sun, it harbors 93 OB stars, together with a population of low-mass stars from embedded protostars to disk-less Class III objects, with age <= 3Myrs. We study an archival 78 ksec Chandra/ACIS-I observation of NGC6611, and two new 80ksec observations of the outer region of M16, one centered on the Column V, and one on a region of the molecular cloud with ongoing star-formation. We detect 1755 point sources, with 1183 candidate cluster members (219 disk-bearing and 964 disk-less). We study the global X-ray properties of M16 and compare them with those of the Orion Nebula Cluster. We also compare the level of X-ray emission of Class II and Class III stars, and analyze the X-ray spectral properties of OB stars. Our study supports the lower level of X-ray activity for the disk-bearing stars with respect to the disk-less members. The X-ray Luminosity Function (XLF) of M16 is similar to that of Orion, supporting the universality of the XLF in young clusters. 85% of the O stars of NGC6611 have been detected in X-rays. With only one possible exception, they show soft spectra with no hard component, indicating that mechanisms for the production of hard X-ray emission in O stars are not operating in NGC 6611.

preprint2012arXiv

Multiwavelength flux variations induced by stellar magnetic activity: effects on planetary transits

Stellar magnetic activity is a source of noise in the study of the transits of extrasolar planets. It induces flux variations which affect significantly the transit depth determination and the derivations of planetary and stellar parameters. Furthermore, the colour dependence of stellar activity may significantly influence the characterization of planetary atmospheres. Here we present a systematic approach to quantify the corresponding stellar flux variations as a function of wavelength bands. We consider a star with spots covering a given fraction of its disc and model the variability in the UBVRIJHK photometric system and in the Spitzer/IRAC wavebands for dwarf stars from G to M spectral types. We compare activity-induced flux variations in different passbands with planetary transits and quantify how they affect the determination of the planetary radius and the analysis of the transmission spectroscopy in the study of planetary atmospheres. We suggest that the monitoring of the systems by using broad band photometry, from visible to infrared, helps to constraining activity effects. The ratio of the relative variations of the stellar fluxes in short wavelength optical bands (e.g., U or B) to near infrared ones (e.g., J or K) can be used to distinguish starspot brightness dips from planetary transits in a stellar light curve. In addition to the perturbations in the measurement of the planetary radius, we find that starspots can affect the determination of the relative semimajor axis and the inclination of the planetary orbit which have a significant impact on the derivation of the stellar density from the transit light curves.

preprint2012arXiv

Spectral classification and HR diagram of pre-main sequence stars in NGC6530

Mechanisms involved in the star formation process and in particular the duration of the different phases of the cloud contraction are not yet fully understood. Photometric data alone suggest that objects coexist in the young cluster NGC6530 with ages from ~1 Myr up to 10 Myrs. We want to derive accurate stellar parameters and, in particular, stellar ages to be able to constrain a possible age spread in the star-forming region NGC6530. We used low-resolution spectra taken with VIMOS@VLT and literature spectra of standard stars to derive spectral types of a subsample of 94 candidate members of this cluster. We assign spectral types to 86 of the 88 confirmed cluster members and derive individual reddenings. Our data are better fitted by the anomalous reddening law with R$_{\rm V}$=5. We confirm the presence of strong differential reddening in this region. We derive fundamental stellar parameters, such as effective temperatures, photospheric colors, luminosities, masses, and ages for 78 members, while for the remaining 8 YSOs we cannot determine the interstellar absorption, since they are likely accretors, and their V-I colors are bluer than their intrinsic colors. The cluster members studied in this work have masses between 0.4 and 4 M$_\odot$ and ages between 1-2 Myrs and 6-7 Myrs. We find that the SE region is the most recent site of star formation, while the older YSOs are loosely clustered in the N and W regions. The presence of two distint generations of YSOs with different spatial distribution allows us to conclude that in this region there is an age spread of ~6-7 Myrs. This is consistent with the scenario of sequential star formation suggested in literature.

preprint2012arXiv

X-ray variability of PMS stars - Toward an explanation of the different X-ray properties of CTTS and WTTS

The intense X-ray emission from coronae and accretion shocks in young PMS stars is likely to play an important role in the evolution and dispersal of circumstellar disks. Several aspects of the physics of this X-ray emission remain mysterious, e.g., whether and how much accretion affects coronal emission. We studied the X-ray variability of ~1 Myr old low-mass PMS stars as a function of timescale, stellar rotation, and stellar characteristics, in order to gain insights on the working of PMS coronae, their X-ray emission, and the circumstellar environment in which they are immersed. We have exploited the ~850 ksec long Chandra observation of the Orion Nebula Cluster obtained by the COUP collaboration in Jan. 2003, and statistically analyzed the X-ray lightcurves of low-mass stars in several subsamples. In particular, we characterized the different X-ray behavior of stars with and without circumstellar accretion disks. Accreting stars (Classical T Tauri Stars, CTTSs) are found to be more variable than non accreting ones (Weak-lined T Tauri Stars, WTTSs) at all timescales and in all the X-ray energy bands considered. Variability is seen to increase with time-scale up to ~10 days, i.e. the longest probed. Signatures of rotational modulation are observed for both CTTSs and WTTSs, and most clearly for CTTSs in the soft X-ray band. Lower mass stars are more variable than higher mass ones. We propose that the difference in variability between CTTSs and WTTSs may be explained assuming that the X-ray emission of CTTS is affected by time-variable absorption due circumstellar structures, such as warps in the inner disk and/or accretion streams. This suggestion is appealing because, in the hypothesis that the coronae of CTTSs and WTTSs are similar, it may also explain why CTTSs have lower and more scattered X-ray emission levels with respect to WTTSs.

preprint2011arXiv

Characterising the Atmospheres of Transiting Planets with a Dedicated Space Telescope

Exoplanetary science is among the fastest evolving fields of today's astronomical research. Ground-based planet-hunting surveys alongside dedicated space missions (Kepler, CoRoT) are delivering an ever-increasing number of exoplanets, now numbering at ~690, with ESA's GAIA mission planned to bring this number into the thousands. The next logical step is the characterisation of these worlds: what is their nature? Why are they as they are? The use of the HST and Spitzer Space Telescope to probe the atmospheres of transiting hot, gaseous exoplanets has demonstrated that it is possible with current technology to address this ambitious goal. The measurements have also shown the difficulty of understanding the physics and chemistry of these environments when having to rely on a limited number of observations performed on a handful of objects. To progress substantially in this field, a dedicated facility for exoplanet characterization with an optimised instrument design (detector performance, photometric stability, etc.), able to observe through time and over a broad spectral range a statistically significant number of planets, will be essential. We analyse the performances of a 1.2/1.4m space telescope for exoplanet transit spectroscopy from the visible to the mid IR, and present the SNR ratio as function of integration time and stellar magnitude/spectral type for the acquisition of spectra of planetary atmospheres in a variety of scenarios: hot, warm, and temperate planets, orbiting stars ranging in spectral type from hot F to cool M dwarfs. We include key examples of known planets (e.g. HD 189733b, Cancri 55 e) and simulations of plausible terrestrial and gaseous planets, with a variety of thermodynamical conditions. We conclude that even most challenging targets, such as super-Earths in the habitable-zone of late-type stars, are within reach of a M-class, space-based spectroscopy mission.

preprint2011arXiv

EChO - Exoplanet Characterisation Observatory

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

preprint2011arXiv

Estimation of the XUV radiation onto close planets and their evaporation

Context: The current distribution of planet mass vs. incident stellar X-ray flux supports the idea that photoevaporation of the atmosphere may take place in close-in planets. Integrated effects have to be accounted for. A proper calculation of the mass loss rate due to photoevaporation requires to estimate the total irradiation from the whole XUV range. Aims: The purpose of this paper is to extend the analysis of the photoevaporation in planetary atmospheres from the accessible X-rays to the mostly unobserved EUV range by using the coronal models of stars to calculate the EUV contribution to the stellar spectra. The mass evolution of planets can be traced assuming that thermal losses dominate the mass loss of their atmospheres. Methods: We determine coronal models for 82 stars with exoplanets that have X-ray observations available. Then a synthetic spectrum is produced for the whole XUV range (~1-912 Å). The determination of the EUV stellar flux, calibrated with real EUV data, allows us to calculate the accumulated effects of the XUV irradiation on the planet atmosphere with time, as well as the mass evolution for planets with known density. Results: We calibrate for the first time a relation of the EUV luminosity with stellar age valid for late-type stars. In a sample of 109 exoplanets, few planets with masses larger than ~1.5 Mj receive high XUV flux, suggesting that intense photoevaporation takes place in a short period of time, as previously found in X-rays. The scenario is also consistent with the observed distribution of planet masses with density. The accumulated effects of photoevaporation over time indicate that HD 209458b may have lost 0.2 Mj since an age of 20 Myr. Conclusions: Coronal radiation produces rapid photoevaporation of the atmospheres of planets close to young late-type stars. More complex models are needed to explain fully the observations.

preprint2011arXiv

High resolution spectroscopy of Ne II emission from young stellar objects

Constraining the spatial and thermal structure of the gaseous component of circumstellar disks is crucial to understand star and planet formation. Models predict that the [Ne II] line at 12.81 μm detected in young stellar objects with Spitzer traces disk gas and its response to high energy radiation, but such [Ne II] emission may also originate in shocks within powerful outflows. To distinguish between these potential origins for mid-infrared [Ne II] emission and to constrain disk models, we observed 32 young stellar objects using the high resolution (R~30000) mid-infrared spectrograph VISIR at the VLT. We detected the 12.81 μm [Ne II] line in 12 objects, tripling the number of detections of this line in young stellar objects with high spatial and spectral resolution spectrographs. We obtain the following main results: a) In Class I objects the [Ne II] emission observed from Spitzer is mainly due to gas at a distance of more than 20-40 AU from the star, where neon is, most likely, ionized by shocks due to protostellar outflows. b) In transition and pre-transition disks, most of the emission is confined to the inner disk, within 20-40 AU from the central star. c) Detailed analysis of line profiles indicates that, in transition and pre-transition disks, the line is slightly blue-shifted (2-12 km s{^-1}) with respect to the stellar velocity, and the line width is directly correlated with the disk inclination, as expected if the emission is due to a disk wind. d) Models of EUV/X-ray irradiated disks reproduce well the observed relation between the line width and the disk inclination, but underestimate the blue-shift of the line.

preprint2011arXiv

Planets and X-rays: a radiation diet

According to theory, high energy emission from the coronae of cool stars can severely erode the atmosphere of orbiting planets. To test the long term effects of the erosion we study a large sample of planet-hosting stars observed in X-rays. The results reveal that massive planets (Mp sin i > 1.5 Mj) may survive only if exposed to low accumulated coronal radiation. The planet HD 209458 b might have lost more than 1 Mj already, and other cases, like tau Boo b, could be losing mass at a rate of 3.4 Earth masses per Gyr. The strongest erosive effects would take place during the first stages of the stellar life, when the faster rotation generates more energetic coronal radiation. The planets with higher density seem to resist better the radiation effects, as foreseen by models. Current models need to be improved to explain the observed distribution of planetary masses with the coronal radiation received.

preprint2011arXiv

Star Formation in the Outer Galaxy: Coronal Properties of NGC 1893

We investigate the X-ray properties of NGC 1893, a young cluster (~ 1-2 Myr) in the outer part of the Galaxy (galactic radius \geq 11 kpc) where we expect differences in the disk evolution and in the mass distribution of the stars, to explore the X-ray emission of its members and compare it with that of young stars in star forming regions near to the Sun. We analyze 5 deep Chandra ACIS-I observations with a total exposure time of 450 ks. Source events of the 1021 X-ray sources have been extracted with the IDL-based routine ACIS-Extract. Using spectral fitting and quantile analysis of X-ray spectra, we derive X-ray luminosities and compare the respective properties of Class II and Class III members. We also evaluate the variability of sources using the Kolmogorov-Smirnov test and we identify flares in the lightcurves. The X-ray luminosity of NGC 1893 X-ray members is in the range 10^29.5 - 10^31.5 erg/s. Diskless stars are brighter in X-rays than disk-bearing stars, given the same bolometric luminosity. We found that 34% of the 1021 lightcurves appear variable and that they show 0.16 flare per source, on average. Comparing our results with those relative to the Orion Nebula Cluster, we find that, accounting for observational biases, the X-ray properties of NGC 1893 and the Orion ones are very similar. The X-ray properties in NGC 1893 are not affected by the environment and the stellar population in the outer Galaxy may have the same coronal properties of nearby star forming regions. The X-ray luminosity properties and the X-ray luminosity function appear to be universal and can therefore be used for distance estimations and for determining stellar properties as already suggested by Feigelson and collaborators.

preprint2011arXiv

The radial distribution of dust species in young brown dwarf disks

We present a study of the radial distribution of dust species in young brown dwarf disks. Our work is based on a compositional analysis of the 10 and 20 micron silicate emission features for brown dwarfs in the Taurus-Auriga star-forming region. A fundamental finding of our work is that brown dwarfs exhibit stronger signs of dust processing in the cold component of the disk, compared to the higher mass T Tauri stars in Taurus. For nearly all of our targets, we find a flat disk structure, which is consistent with the stronger signs of dust processing observed in these disks. For the case of one brown dwarf, 2M04230607, we find the forsterite mass fraction to be a factor of ~3 higher in the outer disk compared to the inner disk region. Simple large-scale radial mixing cannot account for this gradient in the dust chemical composition, and some local crystalline formation mechanism may be effective in this disk. The relatively high abundance of crystalline silicates in the outer cold regions of brown dwarf disks provides an interesting analogy to comets. In this context, we have discussed the applicability of the various mechanisms that have been proposed for comets on the formation and the outward transport of high-temperature material. We also present Chandra X-ray observations for two Taurus brown dwarfs, 2M04414825 and CFHT-BD-Tau 9. We find 2M04414825, which has a ~12% crystalline mass fraction, to be more than an order of magnitude brighter in X-ray than CFHT-BD-Tau 9, which has a ~35% crystalline mass fraction. Combining with previous X-ray data, we find the inner disk crystalline mass fractions to be anti-correlated with the X-ray strength.

preprint2011arXiv

X-ray emission from protostellar jet HH 154: the first evidence of a diamond shock?

X-ray emission from about ten protostellar jets has been discovered and it appears as a feature common to the most energetic jets. Although X-ray emission seems to originate from shocks internal to jets, the mechanism forming these shocks remains controversial. One of the best studied X-ray jet is HH 154 that has been observed by Chandra over a time base of about 10 years. We analyze the Chandra observations of HH 154 by investigating the evolution of its X-ray source. We show that the X-ray emission consists of a bright stationary component and a faint elongated component. We interpret the observations by developing a hydrodynamic model describing a protostellar jet originating from a nozzle and compare the X-ray emission synthesized from the model with the X-ray observations. The model takes into account the thermal conduction and radiative losses and shows that the jet/nozzle leads to the formation of a diamond shock at the nozzle exit. The shock is stationary over the period covered by our simulations and generates an X-ray source with luminosity and spectral characteristics in excellent agreement with the observations. We conclude that the X-ray emission from HH 154 is consistent with a diamond shock originating from a nozzle through which the jet is launched into the ambient medium. We suggest that the physical origin of the nozzle could be related to the dense gas in which the HH 154 driving source is embedded and/or to the magnetic field at the jet launching/collimation region.

preprint2011arXiv

X-ray view of IC348 in the light of an updated cluster census

We study the properties of the coronae of the low-mass stars in the young (~2-3Myr), nearby (~310pc) open cluster IC348 combining X-ray and optical/infrared data. The four existing Chandra observations of IC348 are merged, thus providing a deeper and spatially more complete X-ray view than previous X-ray studies of the cluster. We have compiled a comprehensive catalog of IC348 members taking into account recent updates to the cluster census. Our data collection comprises fundamental stellar parameters, infrared excess indicating the presence of disks, Halpha emission as a tracer of chromospheric emission or accretion and mass accretion rates. We have detected 290 X-ray sources in four merged Chandra exposures, of which 187 are associated with known cluster members. Only four of the X-ray sources are brown dwarfs (spectral type M6 and later). The detection rate is highest for diskless Class III stars and increases with stellar mass. This may be explained with higher X-ray luminosities for higher mass and later evolutionary stage that is evident in the X-ray luminosity functions. In particular, we find that for the lowest examined masses (0.1-0.25 Msun) there is a difference between the X-ray luminosity functions of accreting and non-accreting stars (classified on the basis of their Halpha emission strength) as well as those of disk-bearing and diskless stars (classified on the basis of the slope of the spectral energy distribution). These differences disappear for higher masses. This is related to our finding that the L_x/L_bol ratio is non-constant across the mass/luminosity sequence of IC348 with a decrease towards lower luminosity stars. Our analysis of an analogous stellar sample in the Orion Nebula Cluster suggests that the decline of L_x/L_ bol for young stars at the low-mass end of the stellar sequence is likely universal.

preprint2010arXiv

A detailed study of the rise phase of a long duration X-ray flare in the young star TWA 11B

We analyzed a long duration flare observed in a serendipitous XMM-Newton detection of the M star CD-39 7717B (TWA 11B), member of the young stellar association TW Hya (~ 8 Myr). Only the rise phase (with a duration of ~ 35 ks) and possibly the flare peak were observed. We took advantage of the high count-rate of the X-ray source to carry out a detailed analysis of its spectrum during the whole exposure. After a careful analysis, we interpreted the rise phase as resulting from the ignition of a first group of loops (event A) which triggered a subsequent two-ribbon flare (event B). Event A was analyzed using a single-loop model, while a two-ribbon model was applied for event B. Loop semi-lengths of ~ 4 R* were obtained. Such large structures had been previously observed in very young stellar objects (~ 1 - 4 Myr). This is the first time that they have been inferred in a slightly more evolved star. The fluorescent iron emission line at 6.4 keV was detected during event B. Since TWA 11B seems to have no disk, the most plausible explanation found for its presence in the X-ray spectrum of this star is collisional - or photo- ionization. As far as we are concerned, this is only the third clear detection of Fe photospheric fluorescence in stars other than the Sun.

preprint2010arXiv

A scenario of planet erosion by coronal radiation

Context: According to theory, high-energy emission from the coronae of cool stars can severely erode the atmospheres of orbiting planets. No observational tests of the long term effects of erosion have yet been made. Aims: To analyze the current distribution of planetary mass with X-ray irradiation of the atmospheres in order to make an observational assessment of the effects of erosion by coronal radiation. Methods: We study a large sample of planet-hosting stars with XMM-Newton, Chandra and ROSAT; make a careful identification of X-ray counterparts; and fit their spectra to make accurately measurements of the stellar X-ray flux. Results: The distribution of the planetary masses with X-ray flux suggests that erosion has taken place: most surviving massive planets, (M_p sin i >1.5 M_J), have been exposed to lower accumulated irradiation. Heavy erosion during the initial stages of stellar evolution is followed by a phase of much weaker erosion. A line dividing these two phases could be present, showing a strong dependence on planet mass. Although a larger sample will be required to establish a well-defined erosion line, the distribution found is very suggestive. Conclusions: The distribution of planetary mass with X-ray flux is consistent with a scenario in which planet atmospheres have suffered the effects of erosion by coronal X-ray and EUV emission. The erosion line is an observational constraint to models of atmospheric erosion.

preprint2010arXiv

Chronology of star formation and disk evolution in the Eagle Nebula

Massive SFR are characterized by intense ionizing fluxes, strong stellar winds and supernovae explosions, all of which have important effects on the surrounding media, on the star-formation (SF) process and on the evolution of YSOs and their disks. We present a multiband study of the massive young cluster NGC6611 and M16, to study how OB stars affect the early stellar evolution and the SF. We search for evidence of triggered SF by OB stars in NGC6611 on a large spatial scale (~10 pc) and how the efficiency of disks photoevaporation depends on the central stars mass. We assemble a multiband catalog with photometric data, from B band to 8.0micron, and X-ray data obtained with 2 new and 1 archival ACIS-I observation. We select the stars with disks from IR photometry and disk-less from X-ray emission, both in NGC6611 and the outer region of M16. We study induced photoevaporation searching for the spatial variation of disk frequency for distinct stellar mass ranges. The triggering of SF by OB stars has been investigated by deriving the history of SF across the nebula. We find evidence of sequential SF in the Eagle Nebula going from the SE (2.6 Myrs) to the NW (0.3 Myrs), with the median age of ~1 Myear. We observe a drop of the disk frequency close to OB stars (up to an average distance of 1 pc), without effects at larger distances. Furthermore, disks are more frequent around low-mass stars (<1 M(solar)) than in high-mass stars, regardless of the distance from OB stars. The SF chronology in M16 does not support the hypothesis of a large-scale SF triggered by OB stars in NGC6611. Instead, we speculate that it was triggered by the encounter (~3 Myrs ago) with a giant molecular shell created ~6 Myrs ago.

preprint2010arXiv

Correlated optical and X-ray variability in CTTS - Indications of absorption-modulated emission

Optical and X-ray emission from classical T Tauri stars (CTTSs) has long been known to be highly variable. Our long, uninterrupted optical observation of the NGC2264 region with CoRoT allows the optical variability in CTTS to be studied with unprecedented accuracy and time coverage. Two short Chandra observations obtained during the CoRoT pointing with a separation of 16 days allow us to study whether there is a correlation between optical and X-ray variability on this timescale, thus probing the physical mechanisms driving the variability in both bands. We have computed the optical and X-ray fractional variability between the two 30 ks duration windows covered by both the Chandra and CoRoT observations, for a sample of classical and weak line T Tauri stars (WTTSs) in NGC2264. A scatter plot clearly shows that the variability of CTTSs in the optical and soft X-ray (0.5-1.5 keV) bands is correlated, while no correlation is apparent in the hard (1.5-8.0 keV) band. Also, no correlation in either band is present for WTTSs. We show that the correlation between soft X-ray and optical variability of CTTSs can be naturally explained in terms of time-variable shading (absorption) from circumstellar material orbiting the star, in a scenario rather similar to the one invoked to explain the observed phenomenology in the CTTS AA Tau. The slope of the observed correlation implies (in the hypothesis of homogeneous shading) a significant dust depletion in the circumstellar material (with a gas-to-dust ratio approximately 5 times lower than the standard value for interstellar material).

preprint2010arXiv

Debris Disks of Members of the Blanco 1 Open Cluster

We have used the Spitzer Space Telescope to obtain Multiband Imaging Photometer for Spitzer (MIPS) 24 um photometry for 37 members of the ~100 Myr old open cluster Blanco 1. For the brightest 25 of these stars (where we have 3sigma uncertainties less than 15%), we find significant mid-IR excesses for eight stars, corresponding to a debris disk detection frequency of about 32%. The stars with excesses include two A stars, four F dwarfs and two G dwarfs. The most significant linkage between 24 um excess and any other stellar property for our Blanco 1 sample of stars is with binarity. Blanco 1 members that are photometric binaries show few or no detected 24 um excesses whereas a quarter of the apparently single Blanco 1 members do have excesses. We have examined the MIPS data for two other clusters of similar age to Blanco 1 -- NGC 2547 and the Pleiades. The AFGK photometric binary star members of both of these clusters also show a much lower frequency of 24 um excesses compared to stars that lie near the single-star main sequence. We provide a new determination of the relation between V-Ks color and Ks-[24] color for main sequence photospheres based on Hyades members observed with MIPS. As a result of our analysis of the Hyades data, we identify three low mass Hyades members as candidates for having debris disks near the MIPS detection limit.

preprint2010arXiv

Pre-main sequence stars with disks in the Eagle Nebula observed in scattered light

NGC6611 and its parental cloud, the Eagle Nebula (M16), are well-studied star-forming regions, thanks to their large content of both OB stars and stars with disks and the observed ongoing star formation. We identified 834 disk-bearing stars associated with the cloud, after detecting their excesses in NIR bands from J band to 8.0 micron. In this paper, we study in detail the nature of a subsample of disk-bearing stars that show peculiar characteristics. They appear older than the other members in the V vs. V-I diagram, and/or they have one or more IRAC colors at pure photospheric values, despite showing NIR excesses, when optical and infrared colors are compared. We confirm the membership of these stars to M16 by a spectroscopic analysis. The physical properties of these stars with disks are studied by comparing their spectral energy distributions (SEDs) with the SEDs predicted by models of T-Tauri stars with disks and envelopes. We show that the age of these stars estimated from the V vs. V-I diagram is unreliable since their V-I colors are altered by the light scattered by the disk into the line of sight. Only in a few cases their SEDs are compatible with models with excesses in V band caused by optical veiling. Candidate members with disks and photospheric IRAC colors are selected by the used NIR disk diagnostic, which is sensitive to moderate excesses, such as those produced by disks with low masses. In 1/3 of these cases, scattering of stellar flux by the disks can also be invoked. The photospheric light scattered by the disk grains into the line of sight can affect the derivation of physical parameters of ClassII stars from photometric optical and NIR data. Besides, the disks diagnostic we defined are useful for selecting stars with disks, even those with moderate excesses or whose optical colors are altered by veiling or photospheric scattered light.

preprint2010arXiv

Results from DROXO. III. Observation, source list and X-ray properties of sources detected in the "Deep Rho Ophiuchi XMM-Newton Observation"

X-rays from very young stars are powerful probes to investigate the mechanisms at work in the very first stages of the star formation and the origin of X-ray emission in very young stars. We present results from a 500 ks long observation of the Rho Ophiuchi cloud with a XMM-Newton large program named DROXO, aiming at studying the X-ray emission of deeply embedded Young Stellar Objects (YSOs). The data acquired during the DROXO program were reduced with SAS software, and filtered in time and energy to improve the signal to noise of detected sources; light curves and spectra were obtained. We detected 111 sources, 61 of them associated with rho Ophiuchi YSOs as identified from infrared observations with ISOCAM. Specifically, we detected 9 out of 11 Class I, 31 out of 48 Class II and 15 out 16 Class III objects. Six objects out of 21 classified Class III candidates are also detected. At the same time we suggest that 15 Class III candidates that remain undetected at log Lx < 28.3 are not related to the cloud population. The global detection rate is ~64%. We have achieved a flux sensitivity of ~5 x 10^{-15} erg s^{-1} cm^{-2}. The Lx to L_bol ratio shows saturation at a value of ~ -3.5 for stars with T_eff <= 5000 K or 0.7 M_sun as observed in the Orion Nebula. The plasma temperatures and the spectrum absorption show a decline with YSO class, with Class I YSOs being hotter and more absorbed than Class II and III YSOs. In one star (GY 266) with infrared counterpart in 2MASS and Spitzer catalogs we have detected a soft excess in the X-ray spectrum which is best fitted by a cold thermal component less absorbed than the main thermal component of the plasma. Such a soft component hints to the presence of plasma heated by shocks due to jets outside the dense circumstellar material.

preprint2010arXiv

Star formation in the outer Galaxy: membership and fundamental parameters of the young open cluster NGC 1893

Different environmental conditions can play a crucial role in determining final products of the star formation process and in this context, less favorable activities of star formation are expected in the external regions of our Galaxy. We studied the properties of the young open cluster NGC 1893 located about 12 Kpc from the galactic center, to investigate how different physical conditions can affect the process of star formation. By adopting a multiwavelength approach, we compiled a catalog extending from X-rays to NIR data to derive the cluster membership. In addition, optical and NIR photometric properties are used to evaluate the cluster parameters. We find 415 diskless candidate members plus 1061 young stellar objects with a circumstellar disk or class II candidate members, 125 of which are also Halpha emitters. Considering the diskless candidate members, we find that the cluster distance is 3.6$\pm$0.2 kpc and the mean interstellar reddening is E(B-V)=0.6$\pm$0.1 with evidence of differential reddening in the whole surveyed region. NGC 1893 contains a conspicuous population of pre-main sequence stars together with the well studied main sequence cluster population; we found a disk fraction of about 70% similar to that found in clusters of similar age in the solar neighbour and then, despite expected unfavorable conditions for star formation, we conclude that very rich young clusters can form also in the outer regions of our Galaxy.

preprint2010arXiv

The enigmatic young brown dwarf binary FU Tau: accretion and activity

FU Tau belongs to a rare class of young, wide brown dwarf binaries. We have resolved the system in a Chandra X-ray observation and detected only the primary, FU Tau A. Hard X-ray emission, presumably from a corona, is present but, unexpectedly, we detect also a strong and unusually soft component from FU Tau A. Its X-ray properties, so far unique among brown dwarfs, are very similar to those of the T Tauri star TW Hya. The analogy with TW Hya suggests that the dominating soft X-ray component can be explained by emission from accretion shocks. However, the typical free-fall velocities of a brown dwarf are too low for an interpretation of the observed X-ray temperature as post-shock region. On the other hand, velocities in excess of the free-fall speed are derived from archival optical spectroscopy, and independent pieces of evidence for strong accretion in FU Tau A are found in optical photometry. The high X-ray luminosity of FU Tau A coincides with a high bolometric luminosity confirming an unexplained trend among young brown dwarfs. In fact, FU Tau A is overluminous with respect to evolutionary models while FU Tau B is on the 1 Myr isochrone suggesting non-contemporaneous formation of the two components in the binary. The extreme youth of FU Tau A could be responsible for its peculiar X-ray properties, in terms of atypical magnetic activity or accretion. Alternatively, rotation and magnetic field effects may reduce the efficiency of convection which in turn affects the effective temperature and radius of FU Tau A shifting its position in the HR diagram. Although there is no direct prove of this latter scenario so far we present arguments for its plausibility.

preprint2010arXiv

The New Hard X-ray Mission

The Italian New Hard X-ray Mission (NHXM) has been designed to provide a real breakthrough on a number of hot astrophysical issues that includes: black holes census, the physics of accretion, the particle acceleration mechanisms, the effects of radiative transfer in highly magnetized plasmas and strong gravitational fields. NHXM is an evolution of the HEXIT-Sat concept and it combines fine imaging capability up to 80 keV, today available only at E<10 keV, with sensitive photoelectric imaging polarimetry. It consists of four identical mirrors, with a 10 m focal length, achieved after launch by means of a deployable structure. Three of the four telescopes will have at their focus identical spectral-imaging cameras, while X-ray imaging polarimetric cameras will be placed at the focus of the fourth. In order to ensure a low and stable background, NHXM will be placed in a low Earth equatorial orbit. In this paper we provide an overall description of this mission that is currently in phase B.

preprint2009arXiv

Hard X-ray flux from low-mass stars in the Cygnus OB2 Association

Context. The Cygnus OB2 association, the central engine of the Cygnus X star-forming region, is the subject of an extensive INTEGRAL Key Project that will accumulate 6Ms of observations. Analysis of 2Ms of observations by De Becker and co-workers provides the most sensitive limit yet obtained on hard X-ray emission from the cluster. Aims. We investigate the X-ray emission in the 20-40 keV band expected from the flaring low-mass stellar population in Cygnus OB2. We discuss whether such emission needs to be considered in the interpretation of existing and future X-ray observations of the region, and whether such observations might provide insight into the high-energy processes on low-mass pre-main sequence stars. Methods. The total hard X-ray flux from low-mass stars is estimated by assuming the observed soft X-ray emission stems from a superposition of flares. We further assume the ratio of hard X-ray to soft X-ray emission is described by a scaling found for solar flares by Isola and co-workers. Results. We estimate the low-mass stellar hard X-ray flux in the 20-40 keV band to lie in the range ~2x10^31-6x10^32 erg/s and discuss some potential biases that might affect this result. Conclusions. Hard X-ray emission could lie at a level not much below the current observed flux upper limits for Cygnus OB2. If this emission could be detected, it would provide insight into the hard X-ray production of large flares on pre-main sequence stars. We highlight the penetrating power of hard X-rays from low-mass stellar populations as a possible pointer to our Galaxy's hidden star-forming clusters and super-clusters using more sensitive observations from future missions.

preprint2009arXiv

Quantifying the contamination by old main-sequence stars in young moving groups: the case of the Local Association

The associations and moving groups of young stars are excellent laboratories for investigating stellar formation in the solar neighborhood. Previous results have confirmed that a non-negligible fraction of old main-sequence stars is present in the lists of possible members of young stellar kinematic groups. A detailed study of the properties of these samples is needed to separate the young stars from old main-sequence stars with similar space motion, and identify the origin of these structures. We used stars possible members of the young (~ 10 - 650 Myr) moving groups from the literature. To determine the age of the stars, we used several suitable age indicators for young main sequence stars, i.e., X-ray fluxes and other photometric data. We also used spectroscopic data, in particular the equivalent width of the lithium line Li I and Halpha, to constrain the range of ages of the stars. By combining photometric and spectroscopic data, we were able to separate the young stars (10 - 650 Myr) from the old (> 1 Gyr) field ones. We found, in particular, that the Local Association is contaminated by old field stars at the level of ~30%. This value must be considered as the contamination for our particular sample, and not of the entire Local Association. For other young moving groups, it is more difficult to estimate the fraction of old stars among possible members. However, the level of X-ray emission can, at least, help to separate two age populations: stars with <200 Myr and stars older than this. Our results are consistent with a scenario in which the moving groups contain both groups of young stars formed in a recent star-formation episode and old field stars with similar space motion. Only by combining X-ray and optical spectroscopic data is it possible to distinguish between these two age populations.

preprint2009arXiv

Results from DROXO. II. - [Ne II] and X-ray emission from Rho Ophiuchi young stellar objects

The infrared [Ne II] and [Ne III] fine structure lines at 12.81um and 15.55um are predicted to trace the circumstellar disk gas subject to X-ray heating and ionization. We investigate the origin of these lines by comparing observations with models of X-ray irradiated disks and by searching for empirical correlations between the line luminosities and stellar and circumstellar parameters. We measure neon line fluxes and X-ray luminosities for 28 young stellar objects in the Rho Ophiuchi star formation region for which good quality infrared spectra and X-ray data have been obtained, the former with the Spitzer IRS and the latter with the Deep Rho Ophiuchi XMM-Newton Observation. We detect the [Ne II] and the [Ne III] lines in 10 and 1 cases, respectively. Line luminosities show no correlation with X-ray emission. The luminosity of the [Ne II] line for one star, and that of both the [Ne II] and [Ne III] lines for a second star, match the predictions of published models of X-ray irradiated disks; for the remaining 8 objects the [Ne II] emission is 1-3 dex higher than predicted on the basis of their L_X. Class I objects show significantly stronger [Ne II] lines than Class II and Class III ones. A correlation is moreover found between the [Ne II] line emission and the disk mass accretion rates. This might point toward a role of accretion-generated UV emission in the generation of the line or to other mechanisms related to mass inflows from circumstellar disks and envelopes and/or to the associated mass outflows (winds and jets). We conclude that the X-ray luminosity is not the only parameter that determines the [Ne II] emission. Explaining the strong [Ne II] emission of Class I objects likely requires the inclusion in the models of additional physical components such as the envelope, inflows, and outflows.

preprint2009arXiv

The X-ray emission from Z CMa during an FUor-like outburst and the detection of its X-ray jet

Accretion shocks have been recognized as important X-ray emission mechanism for pre-main sequence stars. Yet the X-ray properties of FUor outbursts, events that are caused by violent accretion, have been given little attention. We have observed the FUor object Z CMa during optical outburst and quiescence with Chandra. No significant changes in X-ray brightness and spectral shape are found, suggesting that the X-ray emission is of coronal nature. Due to the binary nature of Z CMa the origin of the X-ray source is ambiguous. However, the moderate hydrogen column density derived from our data makes it unlikely that the embedded primary star is the X-ray source. The secondary star, which is the FUor object, is thus responsible for both the X-ray emission and the presently ongoing accretion outburst, which seem however to be unrelated phenomena. The secondary is also known to drive a large outflow and jet, that we detect here for the first time in X-rays. The distance of the X-ray emitting outflow source to the central star is higher than in jets of low-mass stars.

preprint2008arXiv

The X-ray cycle in the solar-type star HD 81809

(abridged) Our long-term XMM-Newton program of long-term monitoring of a solar-like star with a well-studied chromospheric cycle, HD 81809 aims to study whether an X-ray cycle is present, along with studying its characteristics and its relation to the chromospheric cycle. Regular observations of HD 81809 were performed with XMM-Newton, spaced by 6 months from 2001 to 2007. We studied the variations in the resulting coronal luminosity and temperature, and compared them with the chromospheric CaII variations. We also modeled the observations in terms of a mixture of active regions, using a methodology originally developed to study the solar corona. Our observations show a well-defined cycle with an amplitude exceeding 1 dex and an average luminosity approximately one order of magnitude higher than in the Sun. The behavior of the corona of HD 81809 can be modeled well in terms of varying coverage of solar-like active regions, with a larger coverage than for the Sun, showing it to be compatible with a simple extension of the solar case.

preprint2007arXiv

Coronal Abundances in Orion Nebula Cluster Stars

Following the Chandra Orion Ultradeep Project (COUP) observation, we have studied the chemical composition of the hot plasma in a sample of 146 X-ray bright pre-main sequence stars in the Orion Nebula Cluster. We report measurements of individual element abundances for a subsample of 86 slightly-absorbed and bright X-ray sources, using low resolution X-ray spectra obtained from the Chandra ACIS instrument. The X-ray emission originates from a plasma with temperatures and elemental abundances very similar to those of active coronae in older stars. A clear pattern of abundances vs. First Ionization Potential (FIP) is evident if solar photospheric abundances are assumed as reference. The results are validated by extensive simulations. The observed abundance distributions are compatible with a single pattern of abundances for all stars, although a weak dependence on flare loop size may be present. The abundance of calcium is the only one which appears to vary substantially between stars, but this quantity is affected by relatively large uncertainties. The ensemble properties of the X-ray bright COUP sources confirm that the iron in the emitting plasma is underabundant with respect to both the solar composition and to the average stellar photospheric values. Comparison of the present plasma abundances with those of the stellar photospheres and those of the gaseous component of the nebula, indicates a good agreement for all the other elements with available measurements, and in particular for the high-FIP elements (Ne, Ar, O, and S) and for the low-FIP element Si. We conclude that there is evidence of a significant chemical fractionation effect only for iron, which appears to be depleted by a factor 1.5--3 with respect to the stellar composition.

preprint2004arXiv

Chandra Orion Ultradeep Project: Observations and Source Lists

We present a description of the data reduction methods and the derived catalog of more than 1600 X-ray point sources from the exceptionally deep January 2003 Chandra X-ray Observatory observation of the Orion Nebula Cluster and embedded populations around OMC-1. The observation was obtained with Chandra's Advanced CCD Imaging Spectrometer (ACIS) and has been nicknamed the Chandra Orion Ultradeep Project (COUP). With an 838 ks exposure made over a continuous period of 13.2 days, the COUP observation provides the most uniform and comprehensive dataset on the X-ray emission of normal stars ever obtained in the history of X-ray astronomy.

preprint2004arXiv

X-ray spectroscopy of the unsteady quiescent corona of AD Leo with Chandra

We present the analysis and interpretation of an observation of the flare star AD Leo (dM3e) with Chandra/LETG. The high resolution X-ray spectrum allowed us to infer the plasma emission measure distribution (EMD) vs. temperature, and the abundances of individual elements in the corona of this magnetically-active star, during a typical state characterized by significant variability but no large flaring event. We have also measured plasma densities at various temperatures using spectroscopic diagnostics provided by He-like triplets and Fe XXI lines. We show that the present EMD is similar to those previously obtained from EUVE spectra during quiescent and flaring states, confirming the long-term stability of the corona of AD Leo. At variance with the case of other active stars, the EMD of AD Leo is characterized by a significantly shallower slope, compatible with that predicted by static models of isobaric loops with constant cross-section and uniform heating. We discuss such coronal modeling to infer the average properties of the corona in terms of loop populations, including estimates of the surface filling factor. We also show that the EMD is compatible with the model of a corona continuously heated by flares. The coronal composition is such that the element abundances, relative to solar values, tend to increase with the First Ionization Potential, with few exceptions. The line to continuum ratios suggest a nearly solar coronal metallicity.