Source author record

D. Barrado

D. Barrado 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

86works
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

86 published item(s)

preprint2025arXiv

Improved characterization of the TOI-2141 system: a dense sub-Neptune with non-transiting inner and outer companions

We aim to refine the fundamental parameters of the TOI-2141 planetary system, which includes a transiting sub-Neptune orbiting a Sun-like star in a relatively long orbit of 18.26 days, by combining new photometric and spectroscopic observations. We analyze new space-based photometry from TESS and CHEOPS as well as 61 radial velocity measurements from HARPS-N. We perform individual and joint photometric and RV analyses using several modeling tools within a Bayesian model comparison framework. We refine the radius and mass of the transiting planet TOI-2141 b to 3.15 $\pm$ 0.04 $R_\oplus$ and 20.1 $\pm$ 1.6 $M_\oplus$, respectively, five and two times more precise than the previously reported values. Our radial velocity analysis reveals two additional non-transiting companions with orbital periods of 5.46 and 60.45 days. Despite the innermost planet's high geometric transit probability, we find no evidence for transits in the photometric data. The bulk properties of TOI-2141 b suggest a significant volatile envelope atop an Earth-like core, with modeling indicating a hydrogen-rich atmosphere that may have experienced mild photoevaporation over the system's history. Planets b and c must exhibit a modest mutual inclination of at least 2.4 degrees.

preprint2025arXiv

JWST/MIRI observations of the young TWA 27 system: hydrocarbon disk chemistry, silicate clouds, evidence for a CPD

The Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) enables the characterisation of young self-luminous gas giants at previously inaccessible wavelengths, revealing physical processes in gas, dust, and clouds. We characterise the young planetary system TWA 27 (2M1207) in the mid-infrared (MIR), studying the atmosphere and disk spectra of the M9 brown dwarf TWA 27A and its L6 planetary-mass companion TWA 27b. We obtained data with the MIRI Medium Resolution Spectrometer (MRS) from 4.9 to 20 um, and MIRI Imaging in the F1000W and F1500W filters. We applied high-contrast imaging methods to extract the companion's spectral energy distribution up to 15 um at 0.78 arcsec separation and a contrast of 60. Combining these with published JWST/NIRSpec spectra, we analysed the 1-20 um range using self-consistent atmospheric grids and 0D slab models for molecular disk emission. The atmosphere of TWA 27A is well fitted by a BT-SETTL model with Teff 2780 K, log g 4.3, plus a 740 K blackbody for the inner disk rim. The disk shows at least 11 organic molecules, with no water or silicate dust emission detected. The atmosphere of TWA 27b is matched by a Teff 1400 K low-gravity model with extinction, best fit by the ExoREM grid. MIRI spectra and photometry for TWA 27b reveal a silicate cloud absorption feature between 8-10 um and significant (>5 sigma) infrared excess at 15 um consistent with circumplanetary disk emission. These MIRI observations provide new insights into TWA 27, revealing diverse features to study the formation and evolution of circumplanetary disks and young dusty atmospheres.

preprint2022arXiv

$\tt{KOBEsim}$: a Bayesian observing strategy algorithm for planet detection in radial velocity blind-search surveys

Ground-based observing time is precious in the era of exoplanet follow-up and characterization, especially in high-precision radial velocity instruments. Blind-search radial velocity surveys thus require a dedicated observational strategy in order to optimize the observing time, which is particularly crucial for the detection of small rocky worlds at large orbital periods. We develop an algorithm with the purpose of improving the efficiency of radial velocity observations in the context of exoplanet searches, and we apply it to the K-dwarfs Orbited By habitable Exoplanets (KOBE) experiment. We aim at accelerating exoplanet confirmations or, alternatively, rejecting false signals as early as possible in order to save telescope time and increase the efficiency of both blind-search surveys and follow-up of transiting candidates. Once a minimum initial number of radial velocity datapoints is reached in such a way that a periodicity starts to emerge according to generalized Lomb-Scargle (GLS) periodograms, that period is targeted with the proposed algorithm, named $\texttt{KOBEsim}$. The algorithm selects the next observing date that maximizes the Bayesian evidence for such periodicity in comparison with a model with no Keplerian orbits. By means of simulated data, we prove that the algorithm accelerates the exoplanet detection, needing $29 - 33\,\%$ less observations and $41 - 47\,\%$ less timespan of the full dataset for low-mass planets ($m_{\rm p}\,<\,10\,M_{\oplus}$) in comparison with a conventional monotonic cadence strategy. The enhancement in the number of datapoints for $20\,M_{\oplus}$ planets is also appreciable, $16\,\%$. We also test $\texttt{KOBEsim}$ with real data for a particular KOBE target, and for the confirmed planet $HD~102365\,b$. Both of them demonstrate that the strategy is capable of speeding up the detection up to a factor of $2$.

preprint2022arXiv

55 Cancri e's occultation captured with CHEOPS

Past occultation and phase-curve observations of the ultra-short period super-Earth 55 Cnc e obtained at visible and infrared wavelengths have been challenging to reconcile with a planetary reflection and emission model. In this study, we analyse a set of 41 occultations obtained over a two-year timespan with the CHEOPS satellite. We report the detection of 55 Cnc e's occultation with an average depth of $12\pm3$ ppm. We derive a corresponding 2-$σ$ upper limit on the geometric albedo of $A_g < 0.55$ once decontaminated from the thermal emission measured by Spitzer at 4.5$μ$m. CHEOPS's photometric performance enables, for the first time, the detection of individual occultations of this super-Earth in the visible and identifies short-timescale photometric corrugations likely induced by stellar granulation. We also find a clear 47.3-day sinusoidal pattern in the time-dependent occultation depths that we are unable to relate to stellar noise, nor instrumental systematics, but whose planetary origin could be tested with upcoming JWST occultation observations of this iconic super-Earth.

preprint2022arXiv

A stellar occultation by the transneptunian object (50000) Quaoar observed by CHEOPS

Stellar occultation is a powerful technique that allows the determination of some physical parameters of the occulting object. The result depends on the photometric accuracy, the temporal resolution, and the number of chords obtained. Space telescopes can achieve high photometric accuracy as they are not affected by atmospheric scintillation. Using ESA's CHEOPS space telescope, we observed a stellar occultation by the Transneptunian object (50000) Quaoar. We compare the obtained chord with previous occultations by this object and determine its astrometry with sub-milliarcsecond precision. Also, we determine upper limits to the presence of a global methane atmosphere on the occulting body. We predicted and observed a stellar occultation by Quaoar using the CHEOPS space telescope. We measured the occultation light curve from this data-set and determined the dis- and re-appearance of the star behind the occulting body. Furthermore, a ground-based telescope in Australia was used to constrain Quaoar's limb. Combined with results from previous works, these measurements allow us to obtain a precise position of Quaoar at the occultation time. We present results obtained from the first stellar occultation by a Transneptunian object (TNO) using space telescope orbiting Earth. It was the occultation by Quaoar observed on 2020 June 11. We used the CHEOPS light curve to obtain a surface pressure upper limit of 85 nbar for the detection of a global methane atmosphere. Also, combining this observation with a ground-based observation we fit Quaoar's limb to determine its astrometric position with an uncertainty below 1.0 mas. This observation is a first of its kind, and it shall be considered as a proof of concept of stellar occultation observations of Transneptunian objects with space telescopes orbiting Earth. Moreover, it shows significant prospects for the James Webb Space Telescope.

preprint2022arXiv

Analysis of Early Science observations with the CHaracterising ExOPlanets Satellite (CHEOPS) using pycheops

CHEOPS(CHaracterising ExOPlanet Satellite) is an ESA S-class mission that observes bright stars at high cadence from low-Earth orbit. The main aim of the mission is to characterize exoplanets that transit nearby stars using ultrahigh precision photometry. Here we report the analysis of transits observed by CHEOPS during its Early Science observing programme for four well-known exoplanets: GJ436b, HD106315b, HD97658b and GJ1132b. The analysis is done using pycheops, an open-source software package we have developed to easily and efficiently analyse CHEOPS light curve data using state-of-the-art techniques that are fully described herein. We show that the precision of the transit parameters measured using CHEOPS is comparable to that from larger space telescopes such as Spitzer Space Telescope and Kepler. We use the updated planet parameters from our analysis to derive new constraints on the internal structure of these four exoplanets.

preprint2022arXiv

Infrared Spectroscopy of free-floating planet candidates in Upper Scorpius and Ophiuchus

Context: A rich population of low-mass brown dwarfs and isolated planetary mass objects has been reported recently in the Upper Scorpius and Ophiuchus star forming complex. Aims: We investigate the membership, nature and properties of 17 of these isolated planetary mass candidates using low-resolution near-infrared spectra. Methods: We investigate the membership by looking for evidences of youth using four diagnostics: the slope of the continuum between the J and Ks band, the Hcont and TLI-g gravity sensitive indices, and by comparing the spectra to young and field (old) M and L-dwarf standards. Results: All the targets but one are confirmed as young ultracool objects, with spectral types between L0 and L6 and masses in the range 0.004-0.013 M according to evolutionary models. The status of the last target is unclear at this point. Conclusions: Only one possible contaminant has been identified among the 17 targets, suggesting that the contamination level of the original sample must be lower than 6%

preprint2022arXiv

Lithium depletion boundary, stellar associations, and Gaia

Stellar ages are key to improving our understanding of different astrophysical phenomena. However, many techniques to estimate stellar ages are highly model-dependent. The lithium depletion boundary (LDB), based on the presence or absence of lithium in low-mass stars, can be used to derive ages in stellar associations of between 20 and 500~Ma. The purpose of this work is to revise former LDB ages in stellar associations in a consistent way, taking advantage of the homogeneous \textit{Gaia} parallaxes as well as bolometric luminosity estimations that do not rely on monochromatic bolometric corrections. We studied nine open clusters and three moving groups characterised by a previous determination of the LDB age. We gathered all the available information from our data and the literature: membership, distances, photometric data, reddening, metallicity, and surface gravity. We re-assigned membership and calculated bolometric luminosities and effective temperatures using distances derived from Gaia DR2 and multi-wavelength photometry for individual objects around the former LDB. We located the LDB using a homogeneous method for all the stellar associations. Finally, we estimated the age by comparing it with different evolutionary models. We located the LDB for the twelve stellar associations and derived their ages using several theoretical evolutionary models. We compared the LDB ages among them, along with data obtained with other techniques, such as isochrone fitting, ultimately finding some discrepancies among the various approaches. Finally, we remark that the 32 Ori MG is likely to be composed of at least two populations of different ages.

preprint2022arXiv

Spi-OPS: Spitzer and CHEOPS confirm the near-polar orbit of MASCARA-1 b and reveal a hint of dayside reflection

The light curves of tidally locked hot Jupiters transiting fast-rotating, early-type stars are a rich source of information about both the planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin-orbit angle $Ψ$, a notoriously difficult parameter to measure, from any transit asymmetry resulting from gravity darkening induced by the stellar rotation, the correlations that exist between the transit parameters have led to large disagreements in published values of $Ψ$ for some systems. We aimed to study these phenomena in the light curves of the ultra-hot Jupiter MASCARA-1 b. We obtained optical CHEOPS transit and occultation light curves of MASCARA-1 b, and analysed them jointly with a Spitzer/IRAC 4.5 $μ$m full-phase curve. When fitting the CHEOPS and Spitzer transits together, the degeneracies are greatly diminished and return results consistent with previously published Doppler tomography. Placing priors informed by the tomography achieves even better precision, allowing a determination of $Ψ=72.1^{+2.5}_{-2.4}$ deg. From the occultations and phase variations, we derived dayside and nightside temperatures of $3062^{+66}_{-68}$ K and $1720\pm330$ K, respectively. In addition, we could separately derive geometric albedo $A_g=0.171^{+0.066}_{-0.068}$ and spherical albedo $A_s=0.266^{+0.097}_{-0.100}$ from the CHEOPS data, and Bond albedo $A_B=0.057^{+0.083}_{-0.101}$ from the Spitzer phase curve. Where possible, priors informed by Doppler tomography should be used when fitting transits of fast-rotating stars, though multi-colour photometry may also unlock an accurate measurement of $Ψ$. Our approach to modelling the phase variations at different wavelengths provides a template for how to separate thermal emission from reflected light in spectrally resolved JWST phase curves.

preprint2022arXiv

The young HD 73583 (TOI-560) planetary system: Two 10-M$_\oplus$ mini-Neptunes transiting a 500-Myr-old, bright, and active K dwarf

We present the discovery and characterisation of two transiting planets observed by \textit{TESS} in the light curves of the young and bright (V=9.67) star HD73583 (TOI-560). We perform an intensive spectroscopic and photometric space- and ground-based follow-up in order to confirm and characterise the system. We found that HD73583 is a young ($\sim 500$~Myr) active star with a rotational period of $12.08 \pm 0.11 $\,d, and a mass and radius of $ 0.73 \pm 0.02 M_\odot$ and $0.65 \pm 0.02 R_\odot$, respectively. HD73583 b ($P_b=6.3980420 _{ - 0.0000062 }^{+0.0000067}$ d) has a mass and radius of $10.2 _{-3.1}^{+3.4} M_\oplus$ and$2.79 \pm 0.10 R_\oplus$, respectively, that gives a density of $2.58 _{-0.81}^{ 0.95} {\rm g\,cm^{-3}}$. HD73583 c ($P_c= 18.87974 _{-0.00074 }^{+0.00086}$) has a mass and radius of $9.7_{-1.7} ^ {+1.8} M_\oplus$ and $2.39_{-0.09}^{+0.10} R_\oplus$, respectively, this translates to a density of $3.88 _{-0.80}^{+0.91} {\rm g\,cm^{-3}}$. Both planets are consistent with worlds made of a solid core surrounded by a volatile envelope. Because of their youth and host star brightness, they both are excellent candidates to perform transmission spectroscopy studies. We expect ongoing atmospheric mass-loss for both planets caused by stellar irradiation. We estimate that the detection of evaporating signatures on H and He would be challenging, but doable with present and future instruments.

preprint2021arXiv

A hot mini-Neptune in the radius valley orbiting solar analogue HD 110113

We report the discovery of HD 110113 b (TOI-755.01), a transiting mini-Neptune exoplanet on a 2.5-day orbit around the solar-analogue HD 110113 (Teff = 5730K). Using TESS photometry and HARPS radial velocities gathered by the NCORES program, we find HD 110113 b has a radius of $2.05\pm0.12$ $R_\oplus$ and a mass of $4.55\pm0.62$ $M_\oplus$. The resulting density of $2.90^{+0.75}_{-0.59}$ g cm^{-3} is significantly lower than would be expected from a pure-rock world; therefore, HD 110113 b must be a mini-Neptune with a significant volatile atmosphere. The high incident flux places it within the so-called radius valley; however, HD 110113 b was able to hold onto a substantial (0.1-1\%) H-He atmosphere over its $\sim4$ Gyr lifetime. Through a novel simultaneous gaussian process fit to multiple activity indicators, we were also able to fit for the strong stellar rotation signal with period $20.8\pm1.2$ d from the RVs and confirm an additional non-transiting planet with a mass of $10.5\pm1.2$ $M_\oplus$ and a period of $6.744^{+0.008}_{-0.009}$ d.

preprint2021arXiv

CHEOPS observations of the HD 108236 planetary system: A fifth planet, improved ephemerides, and planetary radii

The detection of a super-Earth and three mini-Neptunes transiting the bright ($V$ = 9.2 mag) star HD 108236 (also known as TOI-1233) was recently reported on the basis of TESS and ground-based light curves. We perform a first characterisation of the HD 108236 planetary system through high-precision CHEOPS photometry and improve the transit ephemerides and system parameters. We characterise the host star through spectroscopic analysis and derive the radius with the infrared flux method. We constrain the stellar mass and age by combining the results obtained from two sets of stellar evolutionary tracks. We analyse the available TESS light curves and one CHEOPS transit light curve for each known planet in the system. We find that HD 108236 is a Sun-like star with $R_{\star}=0.877\pm0.008 R_{\odot}$, $M_{\star}=0.869^{+0.050}_{-0.048} M_{\odot}$, and an age of $6.7_{-5.1}^{+4.0}$ Gyr. We report the serendipitous detection of an additional planet, HD 108236 f, in one of the CHEOPS light curves. For this planet, the combined analysis of the TESS and CHEOPS light curves leads to a tentative orbital period of about 29.5 days. From the light curve analysis, we obtain radii of $1.615\pm0.051$, $2.071\pm0.052$, $2.539_{-0.065}^{+0.062}$, $3.083\pm0.052$, and $2.017_{-0.057}^{+0.052}$ $R_{\oplus}$ for planets HD 108236 b to HD 108236 f, respectively. These values are in agreement with previous TESS-based estimates, but with an improved precision of about a factor of two. We perform a stability analysis of the system, concluding that the planetary orbits most likely have eccentricities smaller than 0.1. We also employ a planetary atmospheric evolution framework to constrain the masses of the five planets, concluding that HD 108236 b and HD 108236 c should have an Earth-like density, while the outer planets should host a low mean molecular weight envelope.

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 HD 137496 system: A dense, hot super-Mercury and a cold Jupiter

Most of the currently known planets are small worlds with radii between that of the Earth and that of Neptune. The characterization of planets in this regime shows a large diversity in compositions and system architectures, with distributions hinting at a multitude of formation and evolution scenarios. Using photometry from the K2 satellite and radial velocities measured with the HARPS and CORALIE spectrographs, we searched for planets around the bright and slightly evolved Sun-like star HD 137496. We precisely estimated the stellar parameters, $M_*$ = 1.035 +/- 0.022 $M_\odot$, $R_*$ = 1.587 +/- 0.028 $R_\odot$, $T_\text{eff}$ = 5799 +/- 61 K, together with the chemical composition of the slightly evolved star. We detect two planets orbiting HD 137496. The inner planet, HD 137496 b, is a super-Mercury (an Earth-sized planet with the density of Mercury) with a mass of $M_b$ = 4.04 +/- 0.55 $M_\oplus$, a radius of $R_b = 1.31_{-0.05}^{+0.06} R_\oplus,$ and a density of $ρ_b = 10.49_{-1.82}^{+2.08}$ $\mathrm{g cm^{-3}}$. With an interior modeling analysis, we find that the planet is composed mainly of iron, with the core representing over 70% of the planet's mass ($M_{core}/M_{total} = 0.73^{+0.11}_{-0.12}$). The outer planet, HD 137496 c, is an eccentric ($e$ = 0.477 +/- 0.004), long period ($P$ = $479.9_{-1.1}^{+1.0}$ days) giant planet ($M_c\sin i_c$ = 7.66 +/- 0.11 $M_{Jup}$) for which we do not detect a transit. HD 137496 b is one of the few super-Mercuries detected to date. The accurate characterization reported here enhances its role as a key target to better understand the formation and evolution of planetary systems. The detection of an eccentric long period giant companion also reinforces the link between the presence of small transiting inner planets and long period gas giants.

preprint2020arXiv

Chamaeleon DANCe. Revisiting the stellar populations of Chamaeleon I and Chamaeleon II with Gaia-DR2 data

Context: Chamaeleon is the southernmost low-mass star-forming complex within 200 pc from the Sun. Its stellar population has been extensively studied in the past, but the current census of the stellar content is not complete yet and deserves further investigation. Aims: We take advantage of the second data release of the \textit{Gaia} space mission to expand the census of stars in Chamaeleon and to revisit the properties of the stellar populations associated to the Chamaeleon I (Cha I) and Chamaeleon II (Cha II) dark clouds. Methods: We perform a membership analysis of the sources in the \textit{Gaia} catalogue over a field of 100 deg$^{2}$ encompassing the Chamaeleon clouds, and use this new census of cluster members to investigate the 6D structure of the complex. Results: We identify 188 and 41 high-probability members of the stellar populations in Cha I and Cha II, respectively, including 19 and 7 new members. Our sample covers the magnitude range from $G=6$ to $G=20$ mag in Cha I, and from $G=12$ to $G=18$ mag in Cha II. We confirm that the northern and southern subgroups of Cha I are located at different distances ($191.4^{+0.8}_{-0.8}$ pc and $186.7^{+1.0}_{-1.0}$ pc), but they exhibit the same space motion within the reported uncertainties. Cha II is located at a distance of $197.5^{+1.0}_{-0.9}$ pc and exhibits a space motion that is consistent with Cha I within the admittedly large uncertainties on the spatial velocities of the stars that come from radial velocity data. The median age of the stars derived from the Hertzsprung-Russell diagram (HRD) and stellar models is about 1-2 Myr, suggesting that they are somewhat younger than previously thought. We do not detect significant age differences between the Chamaeleon subgroups, but we show that Cha II exhibits a higher fraction of disc-bearing stars compared to Cha I.

preprint2020arXiv

Corona-Australis DANCe. I. Revisiting the census of stars with Gaia-DR2 data

Corona-Australis is one of the nearest regions to the Sun with recent and ongoing star formation, but the current picture of its stellar (and substellar) content is not complete yet. We take advantage of the second data release of the Gaia space mission to revisit the stellar census and search for additional members of the young stellar association in Corona-Australis. We applied a probabilistic method to infer membership probabilities based on a multidimensional astrometric and photometric data set over a field of 128 deg$^{2}$ around the dark clouds of the region. We identify 313 high-probability candidate members to the Corona-Australis association, 262 of which had never been reported as members before. Our sample of members covers the magnitude range between $G\gtrsim5$ mag and $G\lesssim20$ mag, and it reveals the existence of two kinematically and spatially distinct subgroups. There is a distributed `off-cloud' population of stars located in the north of the dark clouds that is twice as numerous as the historically known `on-cloud' population that is concentrated around the densest cores. By comparing the location of the stars in the HR-diagram with evolutionary models, we show that these two populations are younger than 10 Myr. Based on their infrared excess emission, we identify 28 Class II and 215 Class III stars among the sources with available infrared photometry, and we conclude that the frequency of Class II stars (i.e. `disc-bearing' stars) in the on-cloud region is twice as large as compared to the off-cloud population. The distance derived for the Corona-Australis region based on this updated census is $d=149.4^{+0.4}_{-0.4}$ pc, which exceeds previous estimates by about 20 pc.In this paper we provide the most complete census of stars in Corona-Australis available to date that can be confirmed with Gaia data.

preprint2020arXiv

Even More Rapidly Rotating Pre-Main Sequence M Dwarfs with Highly Structured Light Curves: An Initial Survey in the Lower Centaurus-Crux and Upper Centaurus-Lupus Associations

Using K2, we recently discovered a new type of periodic photometric variability while analysing the light curves of members of Upper Sco (Stauffer \etal\ 2017). The 23 exemplars of this new variability type are all mid-M dwarfs, with short rotation periods. Their phased light curves have one or more broad flux dips or multiple arcuate structures which are not explicable by photospheric spots or eclipses by solid bodies. Now, using TESS data, we have searched for this type of variability in the other major sections of Sco-Cen, Upper Centaurus-Lupus (UCL) and Lower Centaurus-Crux (LCC). We identify 28 stars with the same light curve morphologies. We find no obvious difference between the Upper Sco and the UCL/LCC representatives of this class in terms of their light curve morphologies, periods or variability amplitudes. The physical mechanism behind this variability is unknown, but as a possible clue we show that the rapidly rotating mid-M dwarfs in UCL/LCC have slightly different colors from the slowly rotating M dwarfs - they either have a blue excess (hot spots?) or a red excess (warm dust?). One of the newly identified stars (TIC242407571) has a very striking light curve morphology. At about every 0.05 in phase are features that resemble icicles, The "icicles" arise because there is a second periodic system whose main feature is a broad flux dip. Using a toy model, we show that the observed light curve morphology results only if the ratio of the two periods and the flux dip width are carefully arranged.

preprint2020arXiv

Masses for the seven planets in K2-32 and K2-233. Four diverse planets in resonant chain and the first young rocky worlds

High-precision planetary densities are key to derive robust atmospheric properties for extrasolar planets. Measuring precise masses is the most challenging part, especially in multi-planetary systems. We measure the masses and densities of a four-planet near resonant chain system (K2-32), and a young ($\sim400$ Myr old) planetary system consisting of three close-in small planets (K2-233). We obtained 199 new HARPS observations for K2-32 and 124 for K2-233 covering a more than three year baseline. We find that K2-32 is a compact scaled-down version of the Solar System's architecture, with a small rocky inner planet (M$_e=2.1^{+1.3}_{-1.1}$~M$_{\oplus}$, P$_e\sim4.35$~days) followed by an inflated Neptune-mass planet (M$_b=15.0^{+1.8}_{-1.7}$~M$_{\oplus}$, P$_b\sim8.99$~days) and two external sub-Neptunes (M$_c=8.1\pm2.4$~M$_{\oplus}$, P$_c\sim20.66$~days; M$_d=6.7\pm2.5$~M$_{\oplus}$, P$_d\sim31.72$~days). K2-32 becomes one of the few multi-planetary systems with four or more planets known with measured masses and radii. Additionally, we constrain the masses of the three planets in K2-233. For the two inner Earth-size planets we constrain their masses to be smaller than M$_b<11.3$ M$_{\oplus}$ (P$_b\sim2.47$~days), M$_c<12.8$ M$_{\oplus}$ (P$_c\sim7.06$~days). The outer planet is a sub-Neptune size planet with an inferred mass of M$_d=8.3^{+5.2}_{-4.7}$ M$_{\oplus}$ (M$_d<21.1$ M$_{\oplus}$, P$_d\sim24.36$~days). Our observations of these two planetary systems confirm for the first time the rocky nature of two planets orbiting a young star, with relatively short orbital periods ($<7$ days). They provide key information for planet formation and evolution models of telluric planets. Additionally, the Neptune-like derived masses of the three planets K2-32 b, c, d puts them in a relatively unexplored regime of incident flux and planet mass, key for transmission spectroscopy studies.

preprint2020arXiv

Precise mass and radius of a transiting super-Earth planet orbiting the M dwarf TOI-1235: a planet in the radius gap?

We report the confirmation of a transiting planet around the bright, inactive M0.5 V star TOI-1235 (TYC 4384-1735-1, V = 11.5 mag), whose transit signal was detected in the photometric time series of Sectors 14, 20, and 21 of the TESS space mission. We confirm the planetary nature of the transit signal, which has a period of 3.44 d, by using precise radial velocity measurements with CARMENES and HARPS-N spectrographs. A comparison of the properties derived for TOI-1235 b's with theoretical models reveals that the planet has a rocky composition, with a bulk density slightly higher than Earth's. In particular, we measure a mass of M_p = 5.9+/-0.6 M_Earth and a radius of R_p = 1.69+/-0.08 R_Earth, which together result in a density of rho_p = 6.7+1.3-1.1 g/cm3. When compared with other well-characterized exoplanetary systems, the particular combination of planetary radius and mass puts our discovery in the radius gap, a transition region between rocky planets and planets with significant atmospheric envelopes, with few known members. While the exact location of the radius gap for M dwarfs is still a matter of debate, our results constrain it to be located at around 1.7 R_Earth or larger at the insolation levels received by TOI-1235 b (~60 S_Earth), which makes it an extremely interesting object for further studies of planet formation and atmospheric evolution.

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. Characterization of the nearby ultra-compact multiplanetary system YZ Ceti

The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters.

preprint2020arXiv

The lithium-rotation connection in the newly discovered young stellar stream Psc-Eri (Meingast 1)

Context. As a fragile element, lithium is a sensitive probe of physical processes occurring in stellar interiors. Aims. We aim at investigating the relationship between lithium abundance and rotation rate in low-mass members of the newly discovered 125~Myr-old Psc-Eri stellar stream. Methods. We obtained high resolution optical spectra and measure the equivalent width of the 607.8 nm LiI line for 40 members of the Psc-Eri stream, whose rotational periods have been derived by arXiv:1905.10588. Results. We show that a tight correlation exists between lithium content and rotation rate among the late-G to early K-type stars of the Psc-Eri stream. Fast rotators are systematically Li-rich, while slow rotators are Li-depleted. This trend mimics the one previously reported for the similar age Pleiades cluster. Conclusions. The lithium-rotation connection thus seems to be universal over a restricted effective temperature range for low-mass stars at or close to the zero-age main sequence, and does not depend on environmental conditions.

preprint2020arXiv

Unveiling the power spectra of $δ$ Scuti stars with TESS. The temperature, gravity, and frequency scaling relation

Thanks to high-precision photometric data legacy from space telescopes like CoRoT and Kepler, the scientific community could detect and characterize the power spectra of hundreds of thousands of stars. Using the scaling relations, it is possible to estimate masses and radii for solar-type pulsators. However, these stars are not the only kind of stellar objects that follow these rules: $δ$ Scuti stars seem to be characterized with seismic indexes such as the large separation ($Δν$). Thanks to long-duration high-cadence TESS light curves, we analysed more than two thousand of this kind of classical pulsators. In that way, we propose the frequency at maximum power ($ν_{\rm max}$) as a proper seismic index since it is directly related with the intrinsic temperature, mass and radius of the star. This parameter seems not to be affected by rotation, inclination, extinction or resonances, with the exception of the evolution of the stellar parameters. Furthermore, we can constrain rotation and inclination using the departure of temperature produced by the gravity-darkening effect. This is especially feasible for fast rotators as most of $δ$ Scuti stars seem to be.

preprint2016arXiv

Close-in planets around giant stars. Lack of hot-Jupiters and prevalence of multi-planetary systems

Extrasolar planets abound in almost any possible configuration. However, until five years ago, there was a lack of planets orbiting closer than 0.5 au to giant or subgiant stars. Since then, recent detections have started to populated this regime by confirming 13 planetary systems. We discuss the properties of these systems in terms of their formation and evolution off the main sequence. Interestingly, we find that $70.0\pm6.6$ % of the planets in this regime are inner components of multiplanetary systems. This value is 4.2$σ$ higher than for main-sequence hosts, which we find to be $42.4\pm0.1$ %. The properties of the known planets seem to indicate that the closest-in planets (a < 0.06 au) to main-sequence stars are massive (i.e., hot Jupiters) and isolated and that they are subsequently engulfed by their host as it evolves to the red giant branch, leaving only the predominant population of multiplanetary systems in orbits 0.06 < a < 0.5 au. We discuss the implications of this emerging observational trend in the context of formation and evolution of hot Jupiters.

preprint2016arXiv

EPIC211089792 b: an aligned and inflated hot jupiter in a young visual binary

In the present paper we report the discovery of a new hot Jupiter, EPIC211089792 b, first detected by the Super-WASP observatory and then by the K2 space mission during its campaign 4. The planet has a period of 3.25d, a mass of 0.73 +/- 0.04 Mjup, and a radius of 1.19 +/- 0.02 Rjup. The host star is a relatively bright (V=12.5) G7 dwarf with a nearby K5V companion. Based on stellar rotation and the abundance of Lithium, we find that the system might be as young as about 450 Myr. The observation of the Rossiter-McLaughlin effect shows the planet is aligned with respect to the stellar spin. Given the deep transit (20mmag), the magnitude of the star and the presence of a nearby stellar companion, the planet is a good target for both space- and ground-based transmission spectroscopy, in particular in the near-infrared where the both stars are relatively bright.

preprint2016arXiv

K2-30b and K2-34b: two inflated hot-Jupiters around Solar-type stars

We report the discovery of the two hot-Jupiters K2-30b and K2-34b. The two planets were detected transiting their main-sequence star with periods $\sim$ 4.099 and $\sim$ 2.996 days, in campaigns 4 and 5 of the extension of the Kepler mission, K2. Subsequent ground-based radial velocity follow-up with SOPHIE, HARPS-N and CAFE, established the planetary nature of the transiting objects. We analyzed the transit signal, radial velocity and spectral energy distributions of the two systems to characterize their properties. Both planets (K2-30b and K2-34b) are bloated hot-Jupiters (1.25 $R_{\rm Jup}$ and 1.33 $R_{\rm Jup}$) around relatively bright (V =13.5 and V=11.5), slow rotating main-sequence (G8 and F9) stars. Thus, these systems are good candidates for detecting the Rossiter-MacLaughlin effect to measure their obliquity and for atmospheric studies.

preprint2016arXiv

Kepler-539: a young extrasolar system with two giant planets on wide orbits and in gravitational interaction

We confirm the planetary nature of Kepler-539b (aka Kepler object of interest K00372.01), a giant transiting exoplanet orbiting a solar-analogue G2 V star. The mass of Kepler-539b was accurately derived thanks to a series of precise radial velocity measurements obtained with the CAFE spectrograph mounted on the CAHA 2.2m telescope. A simultaneous fit of the radial-velocity data and Kepler photometry revealed that Kepler-539b is a dense Jupiter-like planet with a mass of Mp = 0.97 Mjup and a radius of Rp = 0.747 Rjup, making a complete circular revolution around its parent star in 125.6 days. The semi-major axis of the orbit is roughly 0.5 au, implying that the planet is at roughly 0.45 au from the habitable zone. By analysing the mid-transit times of the 12 transit events of Kepler-539b recorded by the Kepler spacecraft, we found a clear modulated transit time variation (TTV), which is attributable to the presence of a planet c in a wider orbit. The few timings available do not allow us to precisely estimate the properties of Kepler-539c and our analysis suggests that it has a mass between 1.2 and 3.6 Mjup, revolving on a very eccentric orbit (0.4<e<0.6) with a period larger than 1000 days. The high eccentricity of planet c is the probable cause of the TTV modulation of planet b. The analysis of the CAFE spectra revealed a relatively high photospheric lithium content, A(Li)=2.48 dex, which, together with both a gyrochronological and isochronal analysis, suggests that the parent star is relatively young.

preprint2016arXiv

Physical parameters of late M-type members of Chamaleon I and TW Hydrae Association: Dust settling, age dispersion and activity

Although mid-to-late type M dwarfs are the most common stars in our stellar neighborhood, our knowledge of these objects is still limited. Open questions include the evolution of their angular momentum, internal structures, dust settling in their atmospheres, age dispersion within populations. In addition, at young ages, late-type Ms have masses below the hydrogen burning limit and therefore are key objects in the debate on the brown dwarf mechanism of formation. In this work we determine and study in detail the physical parameters of two samples of young, late M-type sources belonging to either the Chamaeleon I Dark Cloud or the TW Hydrae Association and compare them with the results obtained in the literature for other young clusters and also for older, field, dwarfs. We used multi-wavelength photometry to construct and analyze SEDs to determine general properties of the photosphere and disk presence. We also used low resolution optical and near-infrared spectroscopy to study activity, accretion, gravity and effective temperature sensitive indicators. We propose a VO-based spectral index that is both temperature and age sensitive. We derived physical parameters using independent techniques confirming the already common feature/problem of the age/luminosity spread. In particular, we highlight two brown dwarfs showing very similar temperatures but clearly different surface gravity (explained invoking extreme early accretion). We also show how, despite large improvement in the dust treatment in theoretical models, there is still room for further progress in the simultaneous reproduction of the optical and near-infrared features of these cold young objects.

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

Search for light curve modulations among Kepler candidates. Three very low-mass transiting companions

Light curve modulations in the sample of Kepler planet candidates allows the disentangling of the nature of the transiting object by photometrically measuring its mass. This is possible by detecting the effects of the gravitational pull of the companion (ellipsoidal modulations) and in some cases, the photometric imprints of the Doppler effect when observing in a broad band (Doppler beaming). We aim to photometrically unveil the nature of some transiting objects showing clear modulations in the phase-folded Kepler light curve. We selected a subsample among the large crop of Kepler objects of interest (KOIs) based on their chances to show detectable light curve modulations, i.e., close ($a<12~R_{\star}$) and large (in terms of radius) candidates. We modeled their phase-folded light curves with consistent equations for the three effects, namely, reflection, ellipsoidal and beaming (known as REB modulations). We provide detailed general equations for the fit of the REB modulations for the case of eccentric orbits. These equations are accurate to the photometric precisions achievable by current and forthcoming instruments and space missions. By using this mathematical apparatus, we find three close-in {very low-mass companions (two of them in the brown dwarf mass domain)} orbiting main-sequence stars (KOI-554, KOI-1074, and KOI-3728), and reject the planetary nature of the transiting objects (thus classifying them as false positives). In contrast, the detection of the REB modulations and transit/eclipse signal allows the measurement of their mass and radius that can provide important constraints for modeling their interiors since just a few cases of low-mass eclipsing binaries are known. Additionally, these new systems can help to constrain the similarities in the formation process of the more massive and close-in planets (hot Jupiters), brown dwarfs, and very low-mass companions.

preprint2016arXiv

The Gaia-ESO Survey: Lithium depletion in the Gamma Velorum cluster and inflated radii in low-mass pre-main-sequence stars

We show that non-magnetic models for the evolution of pre-main-sequence (PMS) stars *cannot* simultaneously describe the colour-magnitude diagram (CMD) and the pattern of lithium depletion seen in the cluster of young, low-mass stars surrounding $γ^2$ Velorum. The age of 7.5+/-1 Myr inferred from the CMD is much younger than that implied by the strong Li depletion seen in the cluster M-dwarfs and the Li depletion occurs at much redder colours than predicted. The epoch at which a star of a given mass depletes its Li and the surface temperature of that star are both dependent on its radius. We demonstrate that if the low-mass stars have radii ~10 per cent larger at a given mass and age, then both the CMD and Li depletion pattern of the Gamma Vel cluster are explained at a common age of 18-21 Myr. This radius inflation could be produced by some combination of magnetic suppression of convection and extensive cool starspots. Models that incorporate radius inflation suggest that PMS stars similar to those in the Gamma Vel cluster, in the range 0.2<M/Msun<0.7, are at least a factor of two older and ~7 per cent cooler than previously thought and that their masses are much larger (by >30 per cent) than inferred from conventional, non-magnetic models in the Hertzsprung-Russell diagram. Systematic changes of this size may be of great importance in understanding the evolution of young stars, disc lifetimes and the formation of planetary systems.

preprint2016arXiv

The seven sisters DANCe.II. Proper motions and the lithium-rotation-activity connection for G and K Pleiads

Stellar clusters are open windows to understand stellar evolution. Specifically, the change with time and the dependence on mass of different stellar properties. As such, they are our laboratories where different theories can be tested. We try to understand the origin of the connection between lithium depletion in F, G and K stars, rotation and activity, in particular in the Pleiades open cluster. We have collected all the relevant data in the literature, including information regarding rotation period, binarity and activity, and cross-matched with proper motions, multi-wavelength photometry and membership probability from the DANCe database. In order to avoid biases, only Pleiades single members with probabilities larger than p=0.75 have been included in the discussion. Results. The analysis confirms that there is a strong link between activity, rotation and the lithium equivalent width excess, specially for the range Lum(bol) = 0.5-0.2 Lsun (about K2-K7 spectral types or 0.75-0.95 Msun). It is not possible to disentangle these effects but we cannot exclude that the observed lithium overabundance is partially an observational effect due to enhanced activity, due to a large coverage by stellar spots induced by high rotation rates. Since a bona fide lithium enhancement is present in young, fast rotators, both activity and rotation should play a role in the lithium problem.

preprint2015arXiv

A submillimeter search for pre- and proto-brown dwarfs in Chamaeleon II

Context. Chamaeleon II molecular cloud is an active star forming region that offers an excellent opportunity for studying the formation of brown dwarfs in the southern hemisphere. Aims. Our aims are to identify a population of pre- and proto- brown dwarfs (5 sigma mass limit threshold of ~0.015 Msun) and provide information on the formation mechanisms of substellar objects. Methods. We performed high sensitivity observations at 870 microns using the LABOCA bolometer at the APEX telescope towards an active star forming region in Chamaeleon II. The data are complemented with an extensive multiwavelength catalogue of sources from the optical to the far-infrared to study the nature of the LABOCA detections. Results. We detect fifteen cores at 870 microns, and eleven of them show masses in the substellar regime. The most intense objects in the surveyed field correspond to the submillimeter counterparts of the well known young stellar objects DK Cha and IRAS 12500-7658. We identify a possible proto-brown dwarf candidate (ChaII-APEX-L) with IRAC emission at 3.6 and 4.5 microns. Conclusions. Our analysis indicates that most of the spatially resolved cores are transient, and that the point-like starless cores in the sub-stellar regime (with masses between 0.016 Msun and 0.066 Msun) could be pre-brown dwarfs cores gravitationally unstable if they have radii smaller than 220 AU to 907 AU (1.2" to 5" at 178 pc) respectively for different masses. ALMA observations will be the key to reveal the energetic state of these pre-brown dwarfs candidates.

preprint2015arXiv

Eclipsing binaries and fast rotators in the Kepler sample. Characterization via radial velocity analysis from Calar Alto

The Kepler mission has provided high-accurate photometric data in a long time span for more than two hundred thousands stars, looking for planetary transits. Among the detected candidates, the planetary nature of around 15% has been established or validated by different techniques. But additional data is needed to characterize the rest of the candidates and reject other possible configurations. We started a follow-up program to validate, confirm, and characterize some of the planet candidates. In this paper we present the radial velocity analysis (RV) of those presenting large variations, compatible with being eclipsing binaries. We also study those showing large rotational velocities, which prevents us from obtaining the necessary precision to detect planetary-like objects. We present new RV results for 13 Kepler objects of interest (KOIs) obtained with the CAFE spectrograph at the Calar Alto Observatory, and analyze their high-spatial resolution images and the Kepler light curves of some interesting cases. We have found five spectroscopic and eclipsing binaries. Among them, the case of KOI-3853 is of particular interest. This system is a new example of the so-called heartbeat stars, showing dynamic tidal distortions in the Kepler light curve. We have also detected duration and depth variations of the eclipse. We suggest possible scenarios to explain such effect, including the presence of a third substellar body possibly detected in our RV analysis. We also provide upper mass limits to the transiting companions of other six KOIs with large rotational velocities. This property prevents the RV method to obtain the necessary precision to detect planetary-like masses. Finally, we analyze the large RV variations of other two KOIs, incompatible with the presence of planetary-mass objects. These objects are likely to be stellar binaries but a longer timespan is still needed.

preprint2015arXiv

First detection of thermal radio jets in a sample of proto-brown dwarf candidates

We observed with the JVLA at 3.6 and 1.3 cm a sample of 11 proto-brown dwarf candidates in Taurus in a search for thermal radio jets driven by the most embedded brown dwarfs. We detected for the first time four thermal radio jets in proto-brown dwarf candidates. We compiled data from UKIDSS, 2MASS, Spitzer, WISE and Herschel to build the Spectral Energy Distribution (SED) of the objects in our sample, which are similar to typical Class~I SEDs of Young Stellar Objects (YSOs). The four proto-brown dwarf candidates driving thermal radio jets also roughly follow the well-known trend of centimeter luminosity against bolometric luminosity determined for YSOs, assuming they belong to Taurus, although they present some excess of radio emission compared to the known relation for YSOs. Nonetheless, we are able to reproduce the flux densities of the radio jets modeling the centimeter emission of the thermal radio jets using the same type of models applied to YSOs, but with corresponding smaller stellar wind velocities and mass-loss rates, and exploring different possible geometries of the wind or outflow from the star. Moreover, we also find that the modeled mass outflow rates for the bolometric luminosities of our objects agree reasonably well with the trends found between the mass outflow rates and bolometric luminosities of YSOs, which indicates that, despite the "excess" centimeter emission, the intrinsic properties of proto-brown dwarfs are consistent with a continuation of those of very low mass stars to a lower mass range. Overall, our study favors the formation of brown dwarfs as a scaled-down version of low-mass stars.

preprint2015arXiv

Kepler-447b: a hot-Jupiter with an extremely grazing transit

We present the radial velocity confirmation of the extrasolar planet Kepler-447b, initially detected as a candidate by the Kepler mission. In this work, we analyze its transit signal and the radial velocity data obtained with the Calar Alto Fiber-fed Echelle spectrograph (CAFE). By simultaneously modeling both datasets, we obtain the orbital and physical properties of the system. According to our results, Kepler-447b is a Jupiter-mass planet ($M_p=1.37^{+0.48}_{-0.46}~M_{\rm Jup}$), with an estimated radius of $R_p=1.65^{+0.59}_{-0.56}~R_{\rm Jup}$ (uncertainties provided in this work are $3σ$ unless specified). This translates into a sub-Jupiter density. The planet revolves every $\sim7.8$ days in a slightly eccentric orbit ($e=0.123^{+0.037}_{-0.036}$) around a G8V star with detected activity in the Kepler light curve. Kepler-447b transits its host with a large impact parameter ($b=1.076^{+0.112}_{-0.086}$), being one of the few planetary grazing transits confirmed so far and the first in the Kepler large crop of exoplanets. We estimate that only around 20% of the projected planet disk occults the stellar disk. The relatively large uncertainties in the planet radius are due to the large impact parameter and short duration of the transit. Planets with such an extremely large impact parameter can be used to detect and analyze interesting configurations such as additional perturbing bodies, stellar pulsations, rotation of a non-spherical planet, or polar spot-crossing events. All these scenarios would periodically modify the transit properties (depth, duration, and time of mid-transit), what could be detectable with sufficient accurate photometry. Short-cadence photometric data (at the 1 minute level) would help in the search for these exotic configurations in grazing planetary transits like that of Kepler-447b.

preprint2015arXiv

Messier 35 (NGC2168) DANCe I. Membership, proper motions and multi-wavelength photometry

We aim at identifying the cluster's members by deriving membership probabilities for the sources within 1 degree of the cluster's center, going further away than equivalent previous studies. We measure accurate proper motions and multi-wavelength (optical and near-infrared) photometry using ground based archival images of the cluster. We use these measurements to compute membership probabilities. The list of candidate members from Barrado+2001 is used as training set to identify the cluster's locus in a multi-dimensional space made of proper motions, luminosities and colors. The final catalog includes 338892 sources with multi-wavelength photometry. Approximately half (194452) were detected at more than two epochs and we measured their proper motion and used it to derive membership probability. A total of 4349 candidate members with membership probabilities greater than 50% are found in this sample in the luminosity range between 10 and 22mag. The slow proper motion of the cluster and the overlap of its sequence with the field and background sequences in almost all color-magnitude and color-color diagrams complicate the analysis and the contamination level is expected to be significant. Our study nevertheless provides a coherent and quantitative membership analysis of Messier 35 based on a large fraction of the best ground-based data sets obtained over the past 18 years. As such, it represents a valuable input for follow-up studies using in particular the Kepler K2 photometric time series.

preprint2015arXiv

Photo-dynamical mass determination of the multi-planetary system K2-19

K2-19 is the second multi-planetary system discovered with K2 observations. The system is composed of two Neptune size planets close to the 3:2 mean-motion resonance. To better characterise the system we obtained two additional transit observations of K2-19b and five additional radial velocity observations. These were combined with K2 data and fitted simultaneously with the system dynamics (photo-dynamical model) which increases the precision of the transit time measurements. The higher transit time precision allows us to detect the chopping signal of the dynamic interaction of the planets that in turn permits to uniquely characterise the system. Although the reflex motion of the star was not detected, dynamic modelling of the system allowed us to derive planetary masses of $M_b= 44 \pm 12\, M_{\oplus}$ and $M_c = 15.9 \pm 7.0\, M_{\oplus}$ for the inner and the outer planets respectively, leading to densities close to Uranus. We also show that our method allows the derivation of mass ratios using only the 80 days of observations during the first campaign of K2.

preprint2015arXiv

Search for associations containing young stars (SACY). VI. Is multiplicity universal? Stellar multiplicity in the range 3-1000 au from adaptive-optics observations

Context. Young loose nearby associations are unique samples of close (<150 pc), young (approx 5-100 Myr) pre-main sequence (PMS) stars. A significant number of members of these associations have been identified in the SACY collaboration. We can use the proximity and youth of these members to investigate key ingredients in star formation processes, such as multiplicity. Aims. We present the statistics of identified multiple systems from 113 confirmed SACY members. We derive multiplicity frequencies, mass-ratio, and physical separation distributions in a consistent parameter space, and compare our results to other PMS populations and the field. Methods. We have obtained adaptive-optics assisted near-infrared observations with NACO (ESO/VLT) and IRCAL (Lick Observatory) for at least one epoch of all 113 SACY members. We have identified multiple systems using co-moving proper-motion analysis and using contamination estimates. We have explored ranges in projected separation and mass-ratio of a [3-1000 au], and q [0.1-1], respectively. Results. We have identified 31 multiple systems (28 binaries and 3 triples). We derive a multiplicity frequency (MF) of MF_(3-1000au)=28.4 +4.7, -3.9% and a triple frequency (TF) of TF_(3-1000au)=2.8 +2.5, -0.8% in the separation range of 3-1000 au. We do not find any evidence for an increase in the MF with primary mass. The estimated mass-ratio of our statistical sample (with power-law index gamma=-0.04 +/- 0.14) is consistent with a flat distribution (gamma = 0). Conclusions. We show further similarities (but also hints of discrepancies) between SACY and the Taurus region: flat mass-ratio distributions and statistically similar MF and TF values. We also compared the SACY sample to the field (in the separation range of 19-100 au), finding that the two distributions are indistinguishable, suggesting a similar formation mechanism.

preprint2015arXiv

Spectro-astrometry of LkCa 15 with X-Shooter: Searching for emission from LkCa 15b

Planet formation is one explanation for the partial clearing of dust observed in the disks of some T Tauri stars. Indeed studies using state-of-the-art high angular resolution techniques have very recently begun to observe planetary companions in these so-called transitional disks. The goal of this work is to use spectra of the transitional disk object LkCa 15 obtained with X-Shooter on the Very Large Telescope to investigate the possibility of using spectro-astrometry to detect planetary companions to T Tauri stars. It is argued that an accreting planet should contribute to the total emission of accretion tracers such as H$α$ and therefore planetary companions could be detected with spectro-astrometry in the same way as it has been used to detect stellar companions to young stars. A probable planetary-mass companion was recently detected in the disk of LkCa 15. Therefore, it is an ideal target for this pilot study. We studied several key accretion lines in the wavelength range 300 nm to 2.2 $μ$m with spectro-astrometry. While no spectro-astrometric signal is measured for any emission lines the accuracy achieved in the technique is used to place an upper limit on the contribution of the planet to the flux of the H$α$, Pa$γ$, and Pa$β$ lines. The derived upper limits on the flux allows an upper limit of the mass accretion rate, log($\dot{M}_{acc}$) = -8.9 to -9.3 for the mass of the companion between 6 M$_{Jup}$ and 15 M$_{Jup}$, respectively, to be estimated (with some assumptions).

preprint2015arXiv

The EChO science case

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

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 Seven Sisters DANCe. I. Empirical isochrones, Luminosity and Mass Functions of the Pleiades cluster

The DANCe survey provides photometric and astrometric (position and proper motion) measurements for approximately 2 millions unique sources in a region encompassing $\approx$80deg$^{2}$ centered around the Pleiades cluster. We aim at deriving a complete census of the Pleiades, and measure the mass and luminosity function of the cluster. Using the probabilistic selection method described in Sarro+2014, we identify high probability members in the DANCe ($i\ge$14mag) and Tycho-2 ($V\lesssim$12mag) catalogues, and study the properties of the cluster over the corresponding luminosity range. We find a total of 2109 high probability members, of which 812 are new, making it the most extensive and complete census of the cluster to date. The luminosity and mass functions of the cluster are computed from the most massive members down to $\approx$0.025M$_{\odot}$. The size, sensitivity and quality of the sample result in the most precise luminosity and mass functions observed to date for a cluster. Our census supersedes previous studies of the Pleiades cluster populations, both in terms of sensitivity and accuracy.

preprint2014arXiv

Cluster membership probabilities from proper motions and multiwavelength photometric catalogues: I. Method and application to the Pleiades cluster

We present a new technique designed to take full advantage of the high dimensionality (photometric, astrometric, temporal) of the DANCe survey to derive self-consistent and robust membership probabilities of the Pleiades cluster. We aim at developing a methodology to infer membership probabilities to the Pleiades cluster from the DANCe multidimensional astro-photometric data set in a consistent way throughout the entire derivation. The determination of the membership probabilities has to be applicable to censored data and must incorporate the measurement uncertainties into the inference procedure. We use Bayes' theorem and a curvilinear forward model for the likelihood of the measurements of cluster members in the colour-magnitude space, to infer posterior membership probabilities. The distribution of the cluster members proper motions and the distribution of contaminants in the full multidimensional astro-photometric space is modelled with a mixture-of-Gaussians likelihood. We analyse several representation spaces composed of the proper motions plus a subset of the available magnitudes and colour indices. We select two prominent representation spaces composed of variables selected using feature relevance determination techniques based in Random Forests, and analyse the resulting samples of high probability candidates. We consistently find lists of high probability (p > 0.9975) candidates with $\approx$ 1000 sources, 4 to 5 times more than obtained in the most recent astro-photometric studies of the cluster. The methodology presented here is ready for application in data sets that include more dimensions, such as radial and/or rotational velocities, spectral indices and variability.

preprint2014arXiv

Gas and dust in the Beta Pictoris Moving Group as seen by the Herschel Space Observatory

Context. Debris discs are thought to be formed through the collisional grinding of planetesimals, and can be considered as the outcome of planet formation. Understanding the properties of gas and dust in debris discs can help us to comprehend the architecture of extrasolar planetary systems. Herschel Space Observatory far-infrared (IR) photometry and spectroscopy have provided a valuable dataset for the study of debris discs gas and dust composition. This paper is part of a series of papers devoted to the study of Herschel PACS observations of young stellar associations. Aims. This work aims at studying the properties of discs in the Beta Pictoris Moving Group (BPMG) through far-IR PACS observations of dust and gas. Methods. We obtained Herschel-PACS far-IR photometric observations at 70, 100 and 160 microns of 19 BPMG members, together with spectroscopic observations of four of them. Spectroscopic observations were centred at 63.18 microns and 157 microns, aiming to detect [OI] and [CII] emission. We incorporated the new far-IR observations in the SED of BPMG members and fitted modified blackbody models to better characterise the dust content. Results. We have detected far-IR excess emission toward nine BPMG members, including the first detection of an IR excess toward HD 29391.The star HD 172555, shows [OI] emission, while HD 181296, shows [CII] emission, expanding the short list of debris discs with a gas detection. No debris disc in BPMG is detected in both [OI] and [CII]. The discs show dust temperatures in the range 55 to 264 K, with low dust masses (6.6*10^{-5} MEarth to 0.2 MEarth) and radii from blackbody models in the range 3 to 82 AU. All the objects with a gas detection are early spectral type stars with a hot dust component.

preprint2014arXiv

Herschel -PACS observations of [OI] and H2O in Cha II

Gas plays a major role in the dynamical evolution of protoplanetary discs. Its coupling with the dust is the key to our understanding planetary formation. Studying the gas content is therefore a crucial step towards understanding protoplanetary discs evolution. Such a study can be made through spectroscopic observations of emission lines in the far-infrared, where some of the most important gas coolants emit, such as the [OI] 3P1-3 P2 transition at 63.18 microns. We aim at characterising the gas content of protoplanetary discs in the intermediate-aged Chamaeleon II (Cha II) star forming region. We also aim at characterising the gaseous detection fractions within this age range, which is an essential step tracing gas evolution with age in different star forming regions. We obtained Herschel-PACS line scan spectroscopic observations at 63 microns of 19 Cha II Class I and II stars. The observations were used to trace [OI] and o-H2O at 63 microns. The analysis of the spatial distribution of [OI], when extended, can be used to understand the origin of the emission. We have detected [OI] emission toward seven out of the nineteen systems observed, and o-H2O emission at 63.32 microns in just one of them, Sz 61. Cha II members show a correlation between [OI] line fluxes and the continuum at 70 microns, similar to what is observed in Taurus. We analyse the extended [OI] emission towards the star DK Cha and study its dynamical footprints in the PACS Integral Field Unit (IFU). We conclude that there is a high velocity component from a jet combined with a low velocity component with an origin that may be a combination of disc, envelope and wind emission. The stacking of spectra of objects not detected individually in [OI] leads to a marginal 2.6sigma detection that may indicate the presence of gas just below our detection limits for some, if not all, of them.

preprint2014arXiv

High-resolution imaging of $Kepler$ planet host candidates. A comprehensive comparison of different techniques

The Kepler mission has discovered thousands of planet candidates. Currently, some of them have already been discarded; more than 200 have been confirmed by follow-up observations, and several hundreds have been validated. However, most of them are still awaiting for confirmation. Thus, priorities (in terms of the probability of the candidate being a real planet) must be established for subsequent observations. The motivation of this work is to provide a set of isolated (good) host candidates to be further tested by other techniques. We identify close companions of the candidates that could have contaminated the light curve of the planet host. We used the AstraLux North instrument located at the 2.2 m telescope in the Calar Alto Observatory to obtain diffraction-limited images of 174 Kepler objects of interest. The lucky-imaging technique used in this work is compared to other AO and speckle imaging observations of Kepler planet host candidates. We define a new parameter, the blended source confidence level (BSC), to assess the probability of an object to have blended non-detected eclipsing binaries capable of producing the detected transit. We find that 67.2% of the observed Kepler hosts are isolated within our detectability limits, and 32.8% have at least one visual companion at angular separations below 6 arcsec. We find close companions (below 3 arcsec) for the 17.2% of the sample. The planet properties of this sample of non-isolated hosts are revised. We report one possible S-type binary (KOI-3158). We also report three possible false positives (KOIs 1230.01, 3649.01, and 3886.01) due to the presence of close companions. The BSC parameter is calculated for all the isolated targets and compared to both the value prior to any high-resolution image and, when possible, to observations from previous high-spatial resolution surveys in the Kepler sample.

preprint2014arXiv

Orion revisited. II. The foreground population to Orion A

Following the recent discovery of a large population of young stars in front of the Orion Nebula, we carried out an observational campaign with the DECam wide-field camera covering ~10~deg^2 centered on NGC 1980 to confirm, probe the extent of, and characterize this foreground population of pre-main-sequence stars. We confirm the presence of a large foreground population towards the Orion A cloud. This population contains several distinct subgroups, including NGC1980 and NGC1981, and stretches across several degrees in front of the Orion A cloud. By comparing the location of their sequence in various color-magnitude diagrams with other clusters, we found a distance and an age of 380pc and 5~10Myr, in good agreement with previous estimates. Our final sample includes 2123 candidate members and is complete from below the hydrogen-burning limit to about 0.3Msun, where the data start to be limited by saturation. Extrapolating the mass function to the high masses, we estimate a total number of ~2600 members in the surveyed region. We confirm the presence of a rich, contiguous, and essentially coeval population of about 2600 foreground stars in front of the Orion A cloud, loosely clustered around NGC1980, NGC1981, and a new group in the foreground of the OMC-2/3. For the area of the cloud surveyed, this result implies that there are more young stars in the foreground population than young stars inside the cloud. Assuming a normal initial mass function, we estimate that between one to a few supernovae must have exploded in the foreground population in the past few million years, close to the surface of Orion A, which might be responsible, together with stellar winds, for the structure and star formation activity in these clouds. This long-overlooked foreground stellar population is of great significance, calling for a revision of the star formation history in this region of the Galaxy.

preprint2014arXiv

Physical parameter determinations of young Ms. Taking advantage of the Virtual Observatory to compare methodologies

One of the very first steps astronomers working in stellar physics perform to advance in their studies, is to determine the most common/relevant physical parameters of the objects of study (effective temperature, bolometric luminosity, surface gravity, etc.). Different methodologies exist depending on the nature of the data, intrinsic properties of the objects, etc. One common approach is to compare the observational data with theoretical models passed through some simulator that will leave in the synthetic data the same imprint than the observational data carries, and see what set of parameters reproduce the observations best. Even in this case, depending on the kind of data the astronomer has, the methodology changes slightly. After parameters are published, the community tend to quote, praise and criticize them, sometimes paying little attention on whether the possible discrepancies come from the theoretical models, the data themselves or just the methodology used in the analysis. In this work we perform the simple, yet interesting, exercise of comparing the effective temperatures obtained via SED and more detailed spectral fittings (to the same grid of models), of a sample of well known and characterized young M-type objects members to different star forming regions and show how differences in temperature of up to 350 K can be expected just from the difference in methodology/data used. On the other hand we show how these differences are smaller for colder objects even when the complexity of the fit increases like for example introducing differential extinction. To perform this exercise we benefit greatly from the framework offered by the Virtual Observatory.

preprint2014arXiv

Radial velocity confirmation of Kepler-91 b. Additional evidence of its planetary nature using the Calar Alto/CAFE instrument

The object transiting the star Kepler-91 was recently assessed as being of planetary nature. The confirmation was achieved by analysing the light-curve modulations observed in the Kepler data. However, quasi-simultaneous studies claimed a self-luminous nature for this object, thus rejecting it as a planet. In this work, we apply an {independent} approach to confirm the planetary mass of Kepler-91b by using multi-epoch high-resolution spectroscopy obtained with the Calar Alto Fiber-fed Echelle spectrograph (CAFE). We obtain the physical and orbital parameters with the radial velocity technique. In particular, we derive a value of $1.09 \pm 0.20\,M_{\mathrm{Jup}}$ for the mass of Kepler-91b, in excellent agreement with our previous estimate that was based on the orbital brightness modulation.

preprint2014arXiv

The first planet detected in the WTS: an inflated hot-Jupiter in a 3.35 day orbit around a late F-star [ERRATUM]

We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope (UKIRT) in August 2007. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J ~ 16 were constructed for ~60000 stars and searched for periodic transit signals. For one of the most promising transiting candidates, high-resolution spectra taken at the Hobby-Eberly Telescope (HET) allowed us to estimate the spectroscopic parameters of the host star, a late-F main sequence dwarf (V=16.13) with possibly slightly subsolar metallicity, and to measure its radial velocity variations. The combined analysis of the light curves and spectroscopic data resulted in an orbital period of the substellar companion of 3.35 days, a planetary mass of 4.01 +- 0.35 Mj and a planetary radius of 1.49+0.16-0.18 Rj. WTS-1b has one of the largest radius anomalies among the known hot Jupiters in the mass range 3-5 Mj. The high irradiation from the host star ranks the planet in the pM class.

preprint2014arXiv

Towards an unbiased stellar census in open clusters using multi-wavelength photometry

We look for very low-mass members in open clusters with different ages placed at different distances. The main goal is to produce a reliable census of low-mass stars for each of the clusters and derive the Initial Mass Function. To achieve this, we combine deep optical and infrared photometry from our own observing runs and from different public databases. We also characterise the individual stellar parameters of our targets.

preprint2014arXiv

WTS-2 b: a hot Jupiter orbiting near its tidal destruction radius around a K-dwarf

We report the discovery of WTS-2 b, an unusually close-in 1.02-day hot Jupiter (Mp=1.12MJ, Rp=1.363RJ) orbiting a K2V star, which has a possible gravitationally-bound M-dwarf companion at 0.6 arcsec separation contributing ~20 percent of the total flux in the observed J-band light curve. The planet is only 1.5 times the separation from its host star at which it would be destroyed by Roche lobe overflow, and has a predicted remaining lifetime of just ~40 Myr, assuming a tidal dissipation quality factor of Q'*=10^6. Q'* is a key factor in determining how frictional processes within a host star affect the orbital evolution of its companion giant planets, but it is currently poorly constrained by observations. We calculate that the orbital decay of WTS-2 b would correspond to a shift in its transit arrival time of T_shift~17 seconds after 15 years assuming Q'*=10^6. A shift less than this would place a direct observational constraint on the lower limit of Q'* in this system. We also report a correction to the previously published expected T_shift for WASP-18 b, finding that T_shift=356 seconds after 10 years for Q'*=10^6, which is much larger than the estimated 28 seconds quoted in WASP-18 b discovery paper. We attempted to constrain Q'* via a study of the entire population of known transiting hot Jupiters, but our results were inconclusive, requiring a more detailed treatment of transit survey sensitivities at long periods. We conclude that the most informative and straight-forward constraints on Q'* will be obtained by direct observational measurements of the shift in transit arrival times in individual hot Jupiter systems. We show that this is achievable across the mass spectrum of exoplanet host stars within a decade, and will directly probe the effects of stellar interior structure on tidal dissipation.

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 sensitivity analysis of the WFCAM Transit Survey for short-period giant planets around M dwarfs

The WFCAM Transit Survey (WTS) is a near-infrared transit survey running on the United Kingdom Infrared Telescope (UKIRT), designed to discover planets around M dwarfs. The WTS acts as a poor-seeing backup programme for the telescope, and represents the first dedicated wide-field near-infrared transit survey. In this paper we describe the observing strategy of the WTS and the processing of the data to generate lightcurves. We describe the basic properties of our photometric data, and measure our sensitivity based on 950 observations. We show that the photometry reaches a precision of ~4mmag for the brightest unsaturated stars in lightcurves spanning almost 3 years. Optical (SDSS griz) and near-infrared (UKIRT ZYJHK) photometry is used to classify the target sample of 4600 M dwarfs with J magnitudes in the range 11-17. Most have spectral-types in the range M0-M2. We conduct Monte Carlo transit injection and detection simulations for short period (<10 day) Jupiter- and Neptune-sized planets to characterize the sensitivity of the survey. We investigate the recovery rate as a function of period and magnitude for 4 hypothetical star-planet cases: M0-2+Jupiter, M2-4+Jupiter, M0-2+Neptune, M2-4+Neptune. We find that the WTS lightcurves are very sensitive to the presence of Jupiter-sized short-period transiting planets around M dwarfs. Hot Neptunes produce a much weaker signal and suffer a correspondingly smaller recovery fraction. Neptunes can only be reliably recovered with the correct period around the rather small sample (~100) of the latest M dwarfs (M4-M9) in the WTS. The non-detection of a hot-Jupiter around an M dwarf by the WFCAM Transit Survey allows us to place an upper limit of 1.7-2.0 per cent (at 95 per cent confidence) on the planet occurrence rate.

preprint2013arXiv

CAFE: Calar Alto Fiber-fed Echelle spectrograph

We present here CAFE, the Calar Alto Fiber-fed Echelle spectrograph, a new instrument built at the Centro Astronomico Hispano Alemán (CAHA). CAFE is a single fiber, high-resolution ($R\sim$70000) spectrograph, covering the wavelength range between 3650-9800Å. It was built on the basis of the common design for Echelle spectrographs. Its main aim is to measure radial velocities of stellar objects up to $V\sim$13-14 mag with a precision as good as a few tens of $m s^{-1}$. To achieve this goal the design was simplified at maximum, removing all possible movable components, the central wavelength is fixed, so the wavelentgth coverage; no filter wheel, one slit and so on, with a particular care taken in the thermal and mechanical stability. The instrument is fully operational and publically accessible at the 2.2m telescope of the Calar Alto Observatory. In this article we describe (i) the design, summarizing its manufacturing phase; (ii) characterize the main properties of the instrument; (iii) describe the reduction pipeline; and (iv) show the results from the first light and commissioning runs. The preliminar results indicate that the instrument fulfill the specifications and it can achieve the foreseen goals. In particular, they show that the instrument is more efficient than anticipated, reaching a $S/N\sim$20 for a stellar object as faint as $V\sim$14.5 mag in $\sim$2700s integration time. The instrument is a wonderful machine for exoplanetary research (by studying large samples of possible systems cotaining massive planets), galactic dynamics (high precise radial velocities in moving groups or stellar associations) or astrochemistry.

preprint2013arXiv

CALIFA, the Calar Alto Legacy Integral Field Area survey. II. First public data release

We present the first public data release of the CALIFA survey. It consists of science-grade optical datacubes for the first 100 of eventually 600 nearby (0.005<z<0.03) galaxies, obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto observatory. The galaxies in DR1 already cover a wide range of properties in color-magnitude space, morphological type, stellar mass, and gas ionization conditions. This offers the potential to tackle a variety of open questions in galaxy evolution using spatially resolved spectroscopy. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the nominal wavelength range 3745-7500A with a spectral resolution of 6.0A (FWHM), and (ii) a medium-resolution V1200 setup covering the nominal wavelength range 3650-4840A with a spectral resolution of 2.3A (FWHM). We present the characteristics and data structure of the CALIFA datasets that should be taken into account for scientific exploitation of the data, in particular the effects of vignetting, bad pixels and spatially correlated noise. The data quality test for all 100 galaxies showed that we reach a median limiting continuum sensitivity of 1.0x10^-18erg/s/cm^2/A/arcsec^2 at 5635A and 2.2x10^-18erg/s/cm^2/A/arcsec^2 at 4500A for the V500 and V1200 setup respectively, which corresponds to limiting r and g band surface brightnesses of 23.6mag/arcsec^2 and 23.4mag/arcsec^2, or an unresolved emission-line flux detection limit of roughly 1x10^-17erg/s/cm^2/arcsec^2 and 0.6x10^-17erg/s/cm^2/arcsec^2, respectively. The median spatial resolution is 3.7", and the absolute spectrophotometric calibration is better than 15% (1sigma). We also describe the available interfaces and tools that allow easy access to this first public CALIFA data.

preprint2013arXiv

DUst Around NEarby Stars. The survey observational results

Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system's counterparts are the asteroid and Edgeworth-Kuiper belts. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. We used {\it Herschel}/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 $μ$m were obtained, complemented in some cases with observations at 70 $μ$m, and at 250, 350 and 500 $μ$m using SPIRE. The observing strategy was to integrate as deep as possible at 100 $μ$m to detect the stellar photosphere. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of $\sim$ 12.1% $\pm$ 5% before \emph{Herschel} to $\sim$ 20.2% $\pm$ 2%. A significant fraction ($\sim$ 52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70-160$μ$m range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

preprint2013arXiv

Dynamical Analysis of Nearby ClustErs. Automated astrometry from the ground: precision proper motions over wide field

The kinematic properties of the different classes of objects in a given association hold important clues about its member's history, and offer a unique opportunity to test the predictions of the various models of stellar formation and evolution. DANCe (standing for Dynamical Analysis of Nearby ClustErs) is a survey program aimed at deriving a comprehensive and homogeneous census of the stellar and substellar content of a number of nearby (<1kpc) young (<500Myr) associations. Whenever possible, members will be identified based on their kinematics properties, ensuring little contamination from background and foreground sources. Otherwise, the dynamics of previously confirmed members will be studied using the proper motion measurements. We present here the method used to derive precise proper motion measurements, using the Pleiades cluster as a test bench. Combining deep wide field multi-epoch panchromatic images obtained at various obervatories over up to 14 years, we derive accurate proper motions for the sources present in the field of the survey. The datasets cover ~80 square degrees, centered around the Seven Sisters. Using new tools, we have computed a catalog of 6116907 unique sources, including proper motion measurements for 3577478 of them. The catalogue covers the magnitude range between i=12~24mag, achieving a proper motion accuracy <1mas/yr for sources as faint as i=22.5mag. We estimate that our final accuracy reaches 0.3mas/yr in the best cases, depending on magnitude, observing history, and the presence of reference extragalactic sources for the anchoring onto the ICRS.

preprint2013arXiv

Gas and dust in the TW Hydrae Association as seen by the Herschel Space Observatory

Context. Circumstellar discs are the places where planets form, therefore knowledge of their evolution is crucial for our understanding of planet formation. The Herschel Space Observatory is providing valuable data for studying disc systems, thanks to its sensitivity and wavelength coverage. This paper is one of several devoted to analysing and modelling Herschel-PACS observations of various young stellar associations from the GASPS Open Time Key Programme. Aims. The aim of this paper is to elucidate the gas and dust properties of circumstellar discs in the 10 Myr TW Hya Association (TWA) using new far-infrared (IR) imaging and spectroscopy from Herschel-PACS. Methods. We obtained far-IR photometric data at 70, 100, and 160 microns of 14 TWA members; spectroscopic observations centred on the [OI] line at 63.18 microns were also obtained for 9 of the 14. The new photometry for each star was incorporated into its full spectral energy distribution (SED). Results. We detected excess IR emission that is characteristic of circumstellar discs from five TWA members, and computed upper limits for another nine. Two TWA members (TWA 01 and TWA 04B) also show [OI] emission at 63.18 microns. Discs in the TWA association display a variety of properties, with a wide range of dust masses and inner radii, based on modified blackbody modelling. Both transitional and debris discs are found in the sample. Models for sources with a detected IR excess give dust masses in the range from 0.15 Msun to 63 Msun.

preprint2013arXiv

Kepler-91b: a planet at the end of its life. Planet and giant host star properties via light-curve variations

The evolution of planetary systems is intimately linked to the evolution of their host star. Our understanding of the whole planetary evolution process is based on the large planet diversity observed so far. To date, only few tens of planets have been discovered orbiting stars ascending the Red Giant Branch. Although several theories have been proposed, the question of how planets die remains open due to the small number statistics. In this work we study the giant star Kepler-91 (KOI-2133) in order to determine the nature of a transiting companion. This system was detected by the Kepler Space Telescope. However, its planetary confirmation is needed. We confirm the planetary nature of the object transiting the star Kepler-91 by deriving a mass of $ M_p=0.88^{+0.17}_{-0.33} ~M_{\rm Jup}$ and a planetary radius of $R_p=1.384^{+0.011}_{-0.054} ~R_{\rm Jup}$. Asteroseismic analysis produces a stellar radius of $R_{\star}=6.30\pm 0.16 ~R_{\odot}$ and a mass of $M_{\star}=1.31\pm 0.10 ~ M_{\odot} $. We find that its eccentric orbit ($e=0.066^{+0.013}_{-0.017}$) is just $1.32^{+0.07}_{-0.22} ~ R_{\star}$ away from the stellar atmosphere at the pericenter. Kepler-91b could be the previous stage of the planet engulfment, recently detected for BD+48 740. Our estimations show that Kepler-91b will be swallowed by its host star in less than 55 Myr. Among the confirmed planets around giant stars, this is the planetary-mass body closest to its host star. At pericenter passage, the star subtends an angle of $48^{\circ}$, covering around 10% of the sky as seen from the planet. The planetary atmosphere seems to be inflated probably due to the high stellar irradiation.

preprint2012arXiv

A search for pre-substellar cores and proto-brown dwarf candidates in Taurus: multiwavelength analysis in the B213-L1495 clouds

In an attempt to study whether the formation of brown dwarfs (BDs) takes place as a scaled-down version of low-mass stars, we conducted IRAM30m/MAMBO-II observations at 1.2 mm in a sample of 12 proto-BD candidates selected from Spitzer/IRAC data in the B213-L1495 clouds in Taurus. Subsequent observations with the CSO at 350 micron, VLA at 3.6 and 6 cm, and IRAM30m/EMIR in the 12CO(1-0), 13CO(1-0), and N2H+(1-0) transitions were carried out toward the two most promising Spitzer/IRAC source(s), J042118 and J041757. J042118 is associated with a compact (<10 arcsec or <1400 AU) and faint source at 350 micron, while J041757 is associated with a partially resolved (~16 arcsec or ~2000 AU) and stronger source emitting at centimetre wavelengths with a flat spectral index. The corresponding masses of the dust condensations are ~1 and ~5 Mjup for J042118 and J041757, respectively. In addition, about 40 arcsec to the northeast of J041757 we detect a strong and extended submillimetre source, J041757-NE, which is not associated with NIR/FIR emission down to our detection limits, but is clearly detected in 13CO and N2H+ at ~7 km/s, and for which we estimated a total mass of ~100 Mjup, close to the mass required to be gravitationally bound. In summary, our observational strategy has allowed us to find in B213-L1495 two proto-BD candidates and one pre-substellar core candidate, whose properties seem to be consistent with a scaled-down version of low-mass stars.

preprint2012arXiv

A tidal disruption-like X-ray flare from the quiescent galaxy SDSS J120136.02+300305.5

SDSS J120136.02+300305.5 was detected in an XMM-Newton slew from June 2010 with a flux 56 times higher than an upper limit from ROSAT, corresponding to Lx~3x10^44 ergs/s. It has the optical spectrum of a quiescent galaxy (z=0.146). Overall the X-ray flux has evolved consistently with the canonical t^-5/3 model, expected for returning stellar debris from a tidal disruption event, fading by a factor ~300 over 300 days. In detail the source is very variable and became invisible to Swift between 27 and 48 days after discovery, perhaps due to self-absorption. The X-ray spectrum is soft but is not the expected tail of optically thick thermal emission. It may be fit with a Bremsstrahlung or double-power-law model and is seen to soften with time and declining flux. Optical spectra taken 12 days and 11 months after discovery indicate a deficit of material in the broad line and coronal line regions of this galaxy, while a deep radio non-detection implies that a jet was not launched during this event.

preprint2012arXiv

Discovery and characterisation of detached M-dwarf eclipsing binaries in the WFCAM Transit Survey

We report the discovery of 16 detached M-dwarf eclipsing binaries with J<16 mag and provide a detailed characterisation of three of them, using high-precision infrared light curves from the WFCAM Transit Survey (WTS). Such systems provide the most accurate and model-independent method for measuring the fundamental parameters of these poorly understood yet numerous stars, which currently lack sufficient observations to precisely calibrate stellar evolution models. We fully solve for the masses and radii of three of the systems, finding orbital periods in the range 1.5<P<4.9 days, with masses spanning 0.35-0.50 Msun and radii between 0.38-0.50 Rsun, with uncertainties of ~3.5-6.4% in mass and ~2.7-5.5% in radius. Close-companions in short-period binaries are expected to be tidally-locked into fast rotational velocities, resulting in high levels of magnetic activity. This is predicted to inflate their radii by inhibiting convective flow and increasing star spot coverage. The radii of the WTS systems are inflated above model predictions by ~3-12%, in agreement with the observed trend, despite an expected lower systematic contribution from star spots signals at infrared wavelengths. We searched for correlation between the orbital period and radius inflation by combining our results with all existing M-dwarf radius measurements of comparable precision, but we found no statistically significant evidence for a decrease in radius inflation for longer period, less active systems. Radius inflation continues to exists in non-synchronised systems indicating that the problem remains even for very low activity M-dwarfs. Resolving this issue is vital not only for understanding the most populous stars in the Universe, but also for characterising their planetary companions, which hold the best prospects for finding Earth-like planets in the traditional habitable zone.

preprint2012arXiv

Four ultra-short period eclipsing M-dwarf binaries in the WFCAM Transit Survey

We report on the discovery of four ultra-short period (P<0.18 days) eclipsing M-dwarf binaries in the WFCAM Transit Survey. Their orbital periods are significantly shorter than of any other known main-sequence binary system, and are all significantly below the sharp period cut-off at P~0.22 days as seen in binaries of earlier type stars. The shortest-period binary consists of two M4 type stars in a P=0.112 day orbit. The binaries are discovered as part of an extensive search for short-period eclipsing systems in over 260,000 stellar lightcurves, including over 10,000 M-dwarfs down to J=18 mag, yielding 25 binaries with P<0.23 days. In a popular paradigm, the evolution of short period binaries of cool main-sequence stars is driven by loss of angular momentum through magnetised winds. In this scheme, the observed P~0.22 day period cut-off is explained as being due to timescales that are too long for lower-mass binaries to decay into tighter orbits. Our discovery of low-mass binaries with significantly shorter orbits implies that either these timescales have been overestimated for M-dwarfs, e.g. due to a higher effective magnetic activity, or that the mechanism for forming these tight M-dwarf binaries is different from that of earlier type main-sequence stars.

preprint2012arXiv

HD 172555: detection of 63 microns [OI] emission in a debris disc

Context. HD 172555 is a young A7 star belonging to the Beta Pictoris Moving Group that harbours a debris disc. The Spitzer IRS spectrum of the source showed mid-IR features such as silicates and glassy silica species, indicating the presence of a warm dust component with small grains, which places HD 172555 among the small group of debris discs with such properties. The IRS spectrum also shows a possible emission of SiO gas. Aims. We aim to study the dust distribution in the circumstellar disc of HD 172555 and to asses the presence of gas in the debris disc. Methods. As part of the GASPS Open Time Key Programme, we obtained Herschel-PACS photometric and spectroscopic observations of the source. We analysed PACS observations of HD 172555 and modelled the Spectral Energy Distribution (SED) with a modified blackbody and the gas emission with a two-level population model with no collisional de-excitation. Results. We report for the first time the detection of OI atomic gas emission at 63.18 microns in the HD 172555 circumstellar disc. We detect excesses due to circumstellar dust toward HD 172555 in the three photometric bands of PACS (70, 100, and 160 microns). We derive a large dust particle mass of 4.8e-4 Earth masses and an atomic oxygen mass of 2.5e-2*R^2 Earth masses, where R in AU is the separation between the star and the inner disc. Thus, most of the detected mass of the disc is in the gaseous phase.

preprint2012arXiv

Multiplicity in transiting planet host stars. A Lucky Imaging study of Kepler candidates

In the exoplanetary era, the Kepler spacecraft is causing a revolution by discovering thousands of new planet candidates. However, a follow up program is needed in order to reject false candidates and to fully characterize the bona-fide exoplanets. Our main aims are: 1./ Detect and analyze close companions inside the typical Kepler PSF to study if they are the responsible of the dim in the Kepler light curves, 2./ Study the change in the stellar and planetary parameters due to the presence of an unresolved object, 3./ Help to validate those Kepler Objects of Interest that do not present any object inside the Kepler PSF and 4./ Study the multiplicity rate in planet host candidates. Such a large sample of observed planet host candidates allows us to do statistics about the presence of close (visual or bounded) companions to the harboring star. We present here Lucky Imaging observations for a total amount of 98 Kepler Objects of Interest. This technique is based on the acquisition of thousands of very short exposure time images. Then, a selection and combination of a small amount of the best quality frames provides a high resolution image with objects having a 0.1 arcsec PSF. We applied this technique to carry out observations in the Sloan i and Sloan z filters of our Kepler candidates. We find blended objects inside the Kepler PSF for a significant percentage of KOIs. On one hand, only 58.2% of the hosts do not present any object within 6 arcsec. On the other hand, we have found 19 companions closer than 3 arcsec in 17 KOIs. According to their magnitudes and i-z color, 8 of them could be physically bounded to the host star.

preprint2012arXiv

The first planet detected in the WTS: an inflated hot-Jupiter in a 3.35 day orbit around a late F-star

We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J=16 were constructed for 60000 stars and searched for periodic transit signals. For one of the most promising transiting candidates, high-resolution spectra taken at the Hobby-Eberly Telescope allowed us to estimate the spectroscopic parameters of the host star, a late-F main sequence dwarf (V=16.13) with possibly slightly subsolar metallicity, and to measure its radial velocity variations. The combined analysis of the light curves and spectroscopic data resulted in an orbital period of the substellar companion of 3.35 days, a planetary mass of 4.01+-0.35 Mj and a planetary radius of 1.49+-0.17 Rj. WTS-1b has one of the largest radius anomalies among the known hot Jupiters in the mass range 3-5 Mj.

preprint2011arXiv

An icy Kuiper-Belt around the young solar-type star HD 181327

HD 181327 is a young F5/F6V star belonging to the Beta Pictoris moving group (12 Myr). It harbors an optically thin belt of circumstellar material at 90 AU. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. We obtained far-IR observations with the Herschel/PACS instrument, and 3.2 mm observations with the ATCA array. The geometry of the belt is constrained with newly reduced HST/NICMOS images that break the degeneracy between the disk geometry and the dust properties. We use the radiative transfer code GRaTer to compute a large grid of models, and we identify the grain models that best reproduce the Spectral Energy Distribution through a Bayesian analysis. We attempt to detect the [OI] and [CII] lines with PACS spectroscopy, providing observables to our photochemical code ProDiMo. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected in the far-IR with PACS and the disk is resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the line fluxes. We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a total dust mass of 0.05 M\oplus (up to 1 mm). We discuss evidences that the grains are fluffy aggregates. The upper limits on the gas lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than 17 M\oplus. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

preprint2011arXiv

Detection of warm water vapour in Taurus protoplanetary discs by Herschel

Line spectra of 68 Taurus T Tauri stars were obtained with the Herschel-PACS (Photodetector Array Camera & Spectrometer) instrument as part of the GASPS (Gas Evolution in Protoplanetary Systems) survey of protoplanetary discs. A careful examination of the line scans centred on the [OI] 63.18 microns fine-structure line unveiled a line at 63.32 micron in some of these spectra. We identify this line with a transition of ortho-water. It is detected confidently (i.e., >3 sigma) in eight sources, i.e., 24% of the sub-sample with gas-rich discs. Several statistical tests were used to search for correlations with other disc and stellar parameters such as line fluxes of [OI] 6300 Armstrong and 63.18 microns; X-ray luminosity and continuum levels at 63 microns and 850 microns. Correlations are found between the water line fluxes and the [OI] 63.18 microns line luminosity, the dust continuum, and possibly with the stellar X-ray luminosity. This is the first time that this line of warm water vapour has been detected in protoplanetary discs. We discuss its origins, in particular whether it comes from the inner disc and/or disc surface or from shocks in outflows and jets. Our analysis favours a disc origin, with the observed water vapour line produced within 2-3AU from the central stars, where the gas temperature is of the order of 500-600 K.

preprint2011arXiv

Spectroscopy of Very Low Mass Stars and Brown Dwarfs in the Lambda Orionis Star Forming Region. I. Enlarging the census down to the planetary mass domain in Collinder 69

Whilst there is a generally accepted evolutionary scheme for the formation of low-mass stars, the analogous processes when moving down in mass to the brown dwarfs regime are not yet well understood. In this first paper we try to build the most complete and unbiased spectroscopically confirmed census of the population of Collinder 69, the central cluster of the Lambda Orionis star forming region, as a first step in addressing the question of how brown dwarfs and planetary mass objects form. We have studied age dependent features in optical and near-infrared spectra of candidate members to the cluster (such as alkali lines and accretion associated indicators). In addition, we have complemented that study with the analysis of other youth indicators like X-ray emission or mid-infrared excess. We have confirmed the membership to Collinder 69 of \sim90 photometric candidate members. As a byproduct we have determined a temperature scale for young M, very low-mass stars and brown dwarfs. We have assembled one of the most complete Initial Mass Functions from 0.016 to 20 M\odot. And, finally, we have studied the implications of the spatial distribution of the confirmed members on the proposed mechanisms of brown dwarfs formation.

preprint2011arXiv

The substellar population of Sigma Orionis: A deep wide survey

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

preprint2011arXiv

YSOVAR: the first sensitive, wide-area, mid-IR photometric monitoring of the ONC

We present initial results from time series imaging at infrared wavelengths of 0.9 sq. degrees in the Orion Nebula Cluster (ONC). During Fall 2009 we obtained 81 epochs of Spitzer 3.6 and 4.5 micron data over 40 consecutive days. We extracted light curves with ~3% photometric accuracy for ~2000 ONC members ranging from several solar masses down to well below the hydrogen burning mass limit. For many of the stars, we also have time-series photometry obtained at optical (Ic) and/or near-infrared (JKs) wavelengths. Our data set can be mined to determine stellar rotation periods, identify new pre-main-sequence (PMS) eclipsing binaries, search for new substellar Orion members, and help better determine the frequency of circumstellar disks as a function of stellar mass in the ONC. Our primary focus is the unique ability of 3.6 & 4.5 micron variability information to improve our understanding of inner disk processes and structure in the Class I and II young stellar objects (YSOs). In this paper, we provide a brief overview of the YSOVAR Orion data obtained in Fall 2009, and we highlight our light curves for AA-Tau analogs - YSOs with narrow dips in flux, most probably due to disk density structures passing through our line of sight. Detailed follow-up observations are needed in order to better quantify the nature of the obscuring bodies and what this implies for the structure of the inner disks of YSOs.

preprint2010arXiv

Age determination of the HR8799 planetary system using asteroseismology

Discovery of the first planetary system by direct imaging around HR8799 has made the age determination of the host star a very important task. This determination is the key to derive accurate masses of the planets and to study the dynamical stability of the system. The age of this star has been estimated using different procedures. In this work we show that some of these procedures have problems and large uncertainties, and the real age of this star is still unknown, needing more observational constraints. Therefore, we have developed a comprehensive modeling of HR8799, and taking advantage of its gamma Doradus-type pulsations, we have estimated the age of the star using asteroseismology. The accuracy in the age determination depends on the rotation velocity of the star, and therefore an accurate value of the inclination angle is required to solve the problem. Nevertheless, we find that the age estimate for this star previously published in the literature ([30,160] Myr) is unlikely, and a more accurate value might be closer to the Gyr. This determination has deep implications on the value of the mass of the objects orbiting HR8799. An age around $\approx$ 1 Gyr implies that these objects are brown dwarfs.

preprint2010arXiv

Asteroseismic determination of the physical characteristics of the planetary system host HR\,8799 ($λ$ Bootis nature and age)

HR\,8799 is a $λ$ Bootis, $γ$ Doradus star hosting a planetary system and a debris disk with two rings. This makes this system a very interesting target for asteroseismic studies. In particular, this work is devoted to the determination of the internal metallicity of this star, linked with its $λ$ Bootis nature, and its age, taking the advantage of its $γ$ Doradus-type pulsations. To do so we have used the equilibrium code CESAM and the non-adiabatic pulsational code GraCo. We have applied the Frequency Ratio Method and the Time Dependent Convection theory to estimate the mode identification, the Brunt-Vaïsälä frequency integral and the mode instability, making a selection of the possible models fulfilling all observational constraints. Using the position of the star in the HR diagram, the solar metallicity models is discarded. This result contradicts one of the main assumptions of the most accepted hypothesis explaining the $λ$ Bootis nature, the accretion/diffusion of gas from a star with solar metallicity. Therefore, in sight of these new results, a revision of this hypothesis is needed. The inclusion of accurate internal chemical mixing is necessary. The use of the asteroseismological constraints provides a very accurate determination of the physical characteristics of HR\,8799: an age in the ranges [1123, 1625] and [26, 430] Myr, and a mass in the ranges [1.32, 1.33] and [1.44, 1.45] $M_{\odot}$, respectively, depending on the visual angle $i$. The determination of this angle and more accurate multicolor photometric observations can definitively fix the mass, metallicity and age of this star. In particular, an accurate age estimation is needed for a correct understanding of the planetary system. In our study we have found that the age widely used for modelling the system is unlike.

preprint2010arXiv

Cold dust around nearby stars (DUNES). First results: A resolved exo-Kuiper belt around the solar-like star zeta^2 Ret

We present the first far-IR observations of the solar-type stars delta Pav, HR 8501, 51 Peg and zeta^2 Ret, taken within the context of the DUNES Herschel Open Time Key Programme (OTKP). This project uses the PACS and SPIRE instruments with the objective of studying infrared excesses due to exo-Kuiper belts around nearby solar-type stars. The observed 100 um fluxes from delta Pav, HR 8501, and 51 Peg agree with the predicted photospheric fluxes, excluding debris disks brighter than Ldust/Lstar ~ 5 x 10^-7 (1 sigma level) around those stars. A flattened, disk-like structure with a semi-major axis of ~ 100 AU in size is detected around zeta^2 Ret. The resolved structure suggests the presence of an eccentric dust ring, which we interpret as an exo-Kuiper belt with Ldust/Lstar ~ 10^-5.

preprint2010arXiv

Gas in Protoplanetary Discs (GASPS) 1. First results

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

preprint2010arXiv

Gas in the protoplanetary disc of HD 169142: Herschel's view

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

preprint2010arXiv

Herschel-PACS observation of the 10 Myr old T Tauri disk TW Hya: Constraining the disk gas mass

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

preprint2010arXiv

Multi-wavelength study of the disk around the very low-mass star Par-Lup3-4

Par-Lup3-4 is a very low-mass star (spectral type M5) in the Lupus III star-forming region. The object is underluminous by ~4 mag when compared to objects of similar mass in the same association. To better understand the origin of its underluminosity, we have analyzed high angular resolution near-IR imaging data and mid-IR spectroscopy. We have also compared the SED of the target (from the optical to the sub-millimeter regime) to a grid of radiative transfer models of circumstellar disks. The diffraction-limited infrared observations do not show obvious extended emission, allowing us to put an upper limit to the disk outer radius of ~20AU. The lack of extended emission, together with the non detection of a strong 9.8 microns silicate in absorption indicates that Par-Lup3-4 is probably in a Class II (rather than Class I) evolutionary stage. The SED of Par-Lup3-4 resembles that of objects with edge-on disks seen in scattered light, that is, a double peaked-SED and a dip at ~10 microns. We can fit the whole SED with a single disk model with an inclination of 81+/-6 degrees which provides a natural explanation for the under-luminosity of the target. Our analysis allows to put constraints on the disk inner radius, Rin < 0.05 AU, which is very close to the dust sublimation radius, and the maximum size of the dust grains, a_max > 10 microns, which indicates that dust processing has already taken place in Par-Lup3-4.

preprint2010arXiv

Resolving the cold debris disc around a planet-hosting star: PACS photometric imaging observations of q1 Eri (HD10647, HR506)

About two dozen exo-solar debris systems have been spatially resolved. These debris discs commonly display a variety of structural features such as clumps, rings, belts, eccentric distributions and spiral patterns. In most cases, these features are believed to be formed, shaped and maintained by the dynamical influence of planets orbiting the host stars. In very few cases has the presence of the dynamically important planet(s) been inferred from direct observation. The solar-type star q1 Eri is known to be surrounded by debris, extended on scales of < 30''. The star is known to host at least one planet, albeit on an orbit far too small to make it responsible for structures at distances of tens to hundreds of AU. The aim of the present investigation is twofold: to determine the optical and material properties of the debris and to infer the spatial distribution of the dust, which may hint at the presence of additional planets. The photodetector array camera and spectrometer (PACS) aboard the Herschel Space Observatory allows imaging observations in the far infrared at unprecedented resolution, i.e. at better than 6'' to 12'' over the wavelength range of 60 μm to 210 μm. Together with the results from ground-based observations, these spatially resolved data can be modelled to determine the nature of the debris and its evolution more reliably than would be possible from unresolved data alone. For the first time has the q1 Eri disc been resolved at far infrared wavelengths. The PACS observations at 70, 100 and 160 μm reveal an oval image showing a disc-like structure in all bands, the size of which increases with wavelength. Assuming a circular shape yields the inclination of its equatorial plane with respect to that of the sky, i > 53deg. The results of image de-convolution indicate that i likely is larger than 63deg, where 90deg corresponds to an edge-on disc. {abridged}

preprint2010arXiv

The Calar Alto Observatory: current status and future instrumentation

The Calar Alto Observatory, located at 2168m height above the sea level in continental Europe, holds a significant number of astronomical telescopes and experiments, covering a large range of the electromagnetic domain, from gamma-ray to near-infrared. It is a very well characterized site, with excellent logistics. Its main telescopes includes a large suite of instruments. At the present time, new instruments, namely CAFE, PANIC and Carmenes, are under development. We are also planning a new operational scheme in order to optimize the observatory resources.

preprint2010arXiv

The Herschel view of Gas in Protoplanetary Systems (GASPS). First comparisons with a large grid of models

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

preprint2010arXiv

The planetary system host HR\,8799: On its $λ$ Bootis nature

HR\,8799 is a $λ$ Bootis, $γ$ Doradus star hosting a planetary system and a debris disk with two rings. This makes this system a very interesting target for asteroseismic studies. This work is devoted to the determination of the internal metallicity of this star, linked with its $λ$ Bootis nature (i.e., solar surface abundances of light elements, and subsolar surface abundances of heavy elements), taking advantage of its $γ$ Doradus pulsations. This is the most accurate way to obtain this information, and this is the first time such a study is performed for a planetary-system-host star. We have used the equilibrium code CESAM and the non-adiabatic pulsational code GraCo. We have applied the Frequency Ratio Method (FRM) and the Time Dependent Convection theory (TDC) to estimate the mode identification, the Brunt-Vaïsälä frequency integral and the mode instability, making the selection of the possible models. When the non-seismological constraints (i.e its position in the HR diagram) are used, the solar abundance models are discarded. This result contradicts one of the main hypothesis for explaining the $λ$ Bootis nature, namely the accretion/diffusion of gas by a star with solar abundance. Therefore, according to these results, a revision of this hypothesis is needed. The inclusion of accurate internal chemical mixing processes seems to be necessary to explain the peculiar abundances observed in the surface of stars with internal subsolar metallicities. The use of the asteroseismological constraints, like those provided by the FRM or the instability analysis, provides a very accurate determination of the physical characteristics of HR 8799. However, a dependence of the results on the inclination angle $i$ still remains. The determination of this angle, more accurate multicolour photometric observations, and high resolution spectroscopy can definitively fix the mass and metallicity of this star.

preprint2010arXiv

XMM-Newton investigations of the Lambda Orionis star-forming region (XILO). I. The young cluster Collinder 69

This is the first paper of a series devoted to the Lambda Orionis star-forming region, from the X-ray perspective, which will provide a comprehensive view of this complex region. In this paper we focus in uncovering the population of the central, young cluster Collinder 69 (C69), and in particular those diskless members not identified by previous near- and mid-infrared surveys, and to establish the X-ray luminosity function for the association. We have combined two exposures taken with the XMM-Newton satellite with an exhaustive data set of optical, near- and mid-infrared photometry to assess the membership of the X-ray sources based on color-color and color-magnitude diagrams, as well as other properties, such as effective temperatures, masses and bolometric luminosities. We detected a total of 164 X-ray sources, of which 66 are probable and possible cluster members. A total of 16 are newly identified probable members. The two XMM-Newton pointings east and west of the cluster center have allowed us to verify the heterogeneous spatial distribution of young stars, probably related to the large scale structure of the region. The disk fraction of the X-ray detected cluster sample is very low, close to 10%, in remarkable contrast to the low-mass stellar and substellar population (mostly undetected in X-rays) where the disk fraction reaches about 50%. The X-ray luminosity function of C69 provides support for an age of several Myr when compared with other well known young associations. With our improved cluster census we confirm previous reports on the untypically low disk fraction compared to other clusters of several Myr age. The different disk fractions of X-ray detected (essentially solar-like) and undetected (mostly low-mass stars and brown dwarfs) members can be understood as a consequence of a mass-dependence of the time-scale for disk evolution.