• ### The seven sisters DANCe IV. Bayesian hierarchical model(1804.08360)

April 23, 2018 astro-ph.SR
Aims. We develop, test and characterise of a new statistical tool (intelligent system) for the sifting and analysis of nearby young open cluster (NYOC) populations. Methods. Using a Bayesian formalism, this statistical tool is able to obtain the posterior distributions of parameters governing the cluster model. It also uses hierarchical bayesian models to establish weakly informative priors, and incorporates the treatment of missing values and non-homogeneous (heteroscedastic) observational uncertainties. Results. From simulations, we estimate that this statistical tool renders kinematic (proper motion) and photometric (luminosity) distributions of the cluster population with a contamination rate of $5.8 \pm 0.2$ %. The luminosity distributions and present day mass function agree with the ones found by Bouy et al. (2015b) on the completeness interval of the survey. At the probability threshold of maximum accuracy, the classifier recovers $\sim$ 90% of Bouy et al. (2015b) candidate members and finds 10% of new ones. Conclusions. A new statistical tool for the analysis of NYOC is introduced, tested and characterised. Its comprehensive modelling of the data properties allows it to get rid of the biases present in previous works. In particular, those resulting from the use of only completely observed (non-missing) data and the assumption of homoskedastic uncertainties. Also, its Bayesian framework allows it to properly propagate observational uncertainties into membership probabilities and cluster velocity and luminosity distributions. Our results are in a general agreement with those from the literature, although we provide the most up-to-date and extended list of candidate members of the Pleiades cluster.
• ### The Seven Sisters DANCe III: Projected spatial distribution(1711.11037)

April 23, 2018 astro-ph.GA, astro-ph.SR
Methods. We compute Bayesian evidences and Bayes Factors for a set of variations of the classical radial models by King (1962), Elson et al. (1987) and Lauer et al. (1995). The variations incorporate different degrees of model freedom and complexity, amongst which we include biaxial (elliptical) symmetry, and luminosity segregation. As a by-product of the model comparison, we obtain posterior distributions and maximum a posteriori estimates for each set of model parameters. Results. We find that the model comparison results depend on the spatial extent of the region used for the analysis. For a circle of 11.5 parsecs around the cluster centre (the most homogeneous and complete region), we find no compelling reason to abandon Kings model, although the Generalised King model, introduced in this work, has slightly better fitting properties. Furthermore, we find strong evidence against radially symmetric models when compared to the elliptic extensions. Finally, we find that including mass segregation in the form of luminosity segregation in the J band, is strongly supported in all our models. Conclusions. We have put the question of the projected spatial distribution of the Pleiades cluster on a solid probabilistic framework, and inferred its properties using the most exhaustive and least contaminated list of Pleiades candidate members available to date. Our results suggest however that this sample may still lack about 20% of the expected number of cluster members. Therefore, this study should be revised when the completeness and homogeneity of the data can be extended beyond the 11.5 parsecs limit. Such study will allow a more precise determination of the Pleiades spatial distribution, its tidal radius, ellipticity, number of objects and total mass.
• ### Searching for H$_{\alpha}$ emitting sources around MWC758: SPHERE/ZIMPOL high-contrast imaging(1803.09264)

March 25, 2018 astro-ph.SR, astro-ph.EP
MWC758 is a young star surrounded by a transitional disk. Recently, a protoplanet candidate has been detected around MWC758 through high-resolution $L'$-band observations. The candidate is located inside the disk cavity at a separation of $\sim$111 mas from the central star, and at an average position angle of $\sim$165.5 degrees. We have performed simultaneous adaptive optics observations of MWC758 in the H$_{\alpha}$ line and the adjacent continuum using SPHERE/ZIMPOL at the Very Large Telescope (VLT). We aim at detecting accreting protoplanet candidates through spectral angular differential imaging observations. The data analysis does not reveal any H$_{\alpha}$ signal around the target. The derived contrast curve in the B_Ha filter allows us to derive a 5$\sigma$ upper limit of $\sim$7.6 mag at 111 mas, the separation of the previously detected planet candidate. This contrast translates into a H$_{\alpha}$ line luminosity of $L_{\rm H_{\alpha}}\lesssim$ 5$\times$10$^{-5}$ $L_{\odot}$ at 111 mas, and an accretion luminosity of $L_{acc} <$3.7$\times$10$^{-4}\,L_{\odot}$. For the predicted mass range of MWC758b, 0.5-5 $M_{\rm Jup}$, this implies accretion rates of $\dot M \lesssim$ 3.4$\times$(10$^{-8}$-10$^{-9})\,M_{\odot}/yr$, for an average planet radius of 1.1 $R_{\rm Jup}$. Therefore, our estimates are consistent with the predictions of accreting circumplanetary accretion models for $R_{\rm in} = 1 R_{\rm Jup}$. In any case, the non-detection of any H$_{\alpha}$ emitting source in the ZIMPOL images does not allow us to unveil the true nature of the $L'$ detected source.
• ### Early phases in the stellar and substellar formation and evolution: Infrared and submillimeter data in the Barnard 30 dark cloud(1712.06399)

Dec. 18, 2017 astro-ph.GA, astro-ph.SR
The early evolutionary stage of brown dwarfs are not very well characterized, specially during the embedded phase. To gain insight into the dominant formation mechanism of very low-mass objects and brown dwarfs, we conducted deep observations at 870$\mu$m with the LABOCA bolometer at the APEX telescope. Our goal was to identify young sub-mm sources in the Barnard 30 dark cloud. We complemented these data with multi-wavelength observations from the optical to the far-IR and. As a result, we have identified 34 submm sources and a substantial number of possible and probable Barnard 30 members within each individual APEX/LABOCA beam. They can be classified in three distinct groups. First, 15 out of these 34 have a clear optical or IR counterpart to the submm peak and nine of them are potential proto-BDs candidates. Moreover, a substantial number of them could be multiple systems. A second group of 13 sources comprises candidate members with significant infrared excesses located away from the central submm emission. All of them include brown dwarf candidates, some displaying IR excess, but their association with submm emission is unclear. In addition, we have found six starless cores and, based on the total dust mass estimate, three might be pre-substellar (or pre-BDs) cores. Finally, the complete characterization of our APEX/LABOCA sources, focusing on those detected at 24 and/or 70 $\mu$m, indicates that in our sample of 34 submm sources there are, at least: two WTTs, four CTTs, five young stellar objects (YSOs), eight proto-BD candidates (with another three dubious cases), and one Very Low Luminosity object (VeLLO).
• ### The Gaia-ESO Survey: open clusters in Gaia-DR1 - a way forward to stellar age calibration(1711.07699)

Nov. 21, 2017 astro-ph.GA, astro-ph.SR
We describe the methodologies that, taking advantage of Gaia-DR1 and the Gaia-ESO Survey data, enable the comparison of observed open star cluster sequences with stellar evolutionary models. The final, long-term goal is the exploitation of open clusters as age calibrators. We perform a homogeneous analysis of eight open clusters using the Gaia-DR1 TGAS catalogue for bright members, and information from the Gaia-ESO Survey for fainter stars. Cluster membership probabilities for the Gaia-ESO Survey targets are derived based on several spectroscopic tracers. The Gaia-ESO Survey also provides the cluster chemical composition. We obtain cluster parallaxes using two methods. The first one relies on the astrometric selection of a sample of bona fide members, while the other one fits the parallax distribution of a larger sample of TGAS sources. Ages and reddening values are recovered through a Bayesian analysis using the 2MASS magnitudes and three sets of standard models. Lithium depletion boundary (LDB) ages are also determined using literature observations and the same models employed for the Bayesian analysis. For all but one cluster, parallaxes derived by us agree with those presented in Gaia Collaboration et al. (2017), while a discrepancy is found for NGC 2516; we provide evidence supporting our own determination. Inferred cluster ages are robust against models and are generally consistent with literature values. The systematic parallax errors inherent in the Gaia DR1 data presently limit the precision of our results. Nevertheless, we have been able to place these eight clusters onto the same age scale for the first time, with good agreement between isochronal and LDB ages where there is overlap. Our approach appears promising and demonstrates the potential of combining Gaia and ground-based spectroscopic datasets.
• ### Mining the Kilo-Degree Survey for solar system objects(1711.02780)

Nov. 8, 2017 astro-ph.EP, astro-ph.IM
The search for minor bodies in the solar system promises insights into its formation history. Wide imaging surveys offer the opportunity to serendipitously discover and identify these traces of planetary formation and evolution. We aim to present a method to acquire position, photometry, and proper motion measurements of solar system objects in surveys using dithered image sequences. The application of this method on the Kilo-Degree Survey is demonstrated. Optical images of 346 square degree fields of the sky are searched in up to four filters using the AstrOmatic software suite to reduce the pixel to catalog data. The solar system objects within the acquired sources are selected based on a set of criteria depending on their number of observation, motion, and size. The Virtual Observatory SkyBoT tool is used to identify known objects. We observed 20,221 SSO candidates, with an estimated false-positive content of less than 0.05%. Of these SSO candidates, 53.4% are identified by SkyBoT. KiDS can detect previously unknown SSOs because of its depth and coverage at high ecliptic latitude, including parts of the Southern Hemisphere. Thus we expect the large fraction of the 46.6% of unidentified objects to be truly new SSOs. Our method is applicable to a variety of dithered surveys such as DES, LSST, and Euclid. It offers a quick and easy-to-implement search for solar system objects. SkyBoT can then be used to estimate the completeness of the recovered sample.
• ### 280 one-opposition near-Earth asteroids recovered by the EURONEAR with the Isaac Newton Telescope(1711.00709)

Nov. 3, 2017 astro-ph.EP, astro-ph.IM
One-opposition near-Earth asteroids (NEAs) are growing in number, and they must be recovered to prevent loss and mismatch risk, and to improve their orbits, as they are likely to be too faint for detection in shallow surveys at future apparitions. We aimed to recover more than half of the one-opposition NEAs recommended for observations by the Minor Planet Center (MPC) using the Isaac Newton Telescope (INT) in soft-override mode and some fractions of available D-nights. During about 130 hours in total between 2013 and 2016, we targeted 368 NEAs, among which 56 potentially hazardous asteroids (PHAs), observing 437 INT Wide Field Camera (WFC) fields and recovering 280 NEAs (76% of all targets). Engaging a core team of about ten students and amateurs, we used the THELI, Astrometrica, and the Find_Orb software to identify all moving objects using the blink and track-and-stack method for the faintest targets and plotting the positional uncertainty ellipse from NEODyS. Most targets and recovered objects had apparent magnitudes centered around V~22.8 mag, with some becoming as faint as V~24 mag. One hundred and three objects (representing 28% of all targets) were recovered by EURONEAR alone by Aug 2017. Orbital arcs were prolonged typically from a few weeks to a few years; our oldest recoveries reach 16 years. The O-C residuals for our 1,854 NEA astrometric positions show that most measurements cluster closely around the origin. In addition to the recovered NEAs, 22,000 positions of about 3,500 known minor planets and another 10,000 observations of about 1,500 unknown objects (mostly main-belt objects) were promptly reported to the MPC by our team. Four new NEAs were discovered serendipitously in the analyzed fields, increasing the counting to nine NEAs discovered by the EURONEAR in 2014 and 2015.
• Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.
• Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.
• ### The Gaia-ESO Survey: Structural and dynamical properties of the young cluster Chamaeleon I(1701.03741)

Jan. 13, 2017 astro-ph.GA, astro-ph.SR
The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $\gamma$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density \Sigma and the level of mass segregation $\Lambda_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured.
• ### A revised moving cluster distance to the Pleiades open cluster(1610.05641)

Oct. 18, 2016 astro-ph.SR
Aims: We provide a new distance estimate for the Pleiades based on the moving cluster method, which will be useful to further discuss the so-called Pleiades distance controversy and compare it with the very precise parallaxes from the Gaia space mission. Methods: We apply a refurbished implementation of the convergent point search method to an updated census of Pleiades stars to calculate the convergent point position of the cluster from stellar proper motions. Then, we derive individual parallaxes for 64 cluster members using radial velocities compiled from the literature, and approximate parallaxes for another 1146 stars based on the spatial velocity of the cluster. This represents the largest sample of Pleiades stars with individual distances to date. Results: The parallaxes derived in this work are in good agreement with previous results obtained in different studies (excluding Hipparcos) for individual stars in the cluster. We report a mean parallax of $7.44\pm 0.08$~mas and distance of $134.4^{+2.9}_{-2.8}$pc that is consistent with the weighted mean of $135.0\pm 0.6$pc obtained from the non-Hipparcos results in the literature. Conclusions: Our result for the distance to the Pleiades open cluster is not consistent with the Hipparcos catalog, but favors the recent and more precise distance determination of $136.2\pm 1.2$pc obtained from Very Long Baseline Interferometry observations. It is also in good agreement with the mean distance of $133\pm 5$pc obtained from the first trigonometric parallaxes delivered by the Gaia satellite for the brightest cluster members in common with our sample.
• ### Orion revisited III. The Orion Belt population(1609.04948)

Sept. 16, 2016 astro-ph.GA, astro-ph.SR
This paper continues our study of the foreground population to the Orion molecular clouds. The goal is to characterize the foreground population north of NGC 1981 and to investigate the star formation history in the large Orion star-forming region. We focus on a region covering about 25 square degrees, centered on the $\epsilon$ Orionis supergiant (HD 37128, B0\,Ia) and covering the Orion Belt asterism. We used a combination of optical (SDSS) and near-infrared (2MASS) data, informed by X-ray (\textit{XMM-Newton}) and mid-infrared (WISE) data, to construct a suite of color-color and color-magnitude diagrams for all available sources. We then applied a new statistical multiband technique to isolate a previously unknown stellar population in this region. We identify a rich and well-defined stellar population in the surveyed region that has about 2\,000 objects that are mostly M stars. We infer the age for this new population to be at least 5\, Myr and likely $\sim10$\,Myr and estimate a total of about 2\,500 members, assuming a normal IMF. This new population, which we call the Orion Belt population, is essentially extinction-free, disk-free, and its spatial distribution is roughly centered near $\epsilon$ Ori, although substructure is clearly present. The Orion Belt population is likely the low-mass counterpart to the Ori OB Ib subgroup. Although our results do not rule out Blaauw's sequential star formation scenario for Orion, we argue that the recently proposed blue streams scenario provides a better framework on which one can explain the Orion star formation region as a whole. We speculate that the Orion Belt population could represent the evolved counterpart of an Orion nebula-like cluster.
• ### Proper motion survey and kinematic analysis of the Rho Ophiuchi embedded cluster(1609.04963)

Sept. 16, 2016 astro-ph.GA, astro-ph.SR
We aim at performing a kinematic census of young stellar objects (YSOs) in the Rho Ophiuchi F core and partially in the E core of the L1688 dark cloud. We run a proper motion program at the ESO New Technology Telescope (NTT) with the Son of ISAAC (SOFI) instrument over nine years in the near-infrared. We complemented these observations with various public image databases to enlarge the time base of observations and the field of investigation to 0.5 deg X 0.5 deg. We derived positions and proper motions for 2213 objects. From these, 607 proper motions were derived from SOFI observations with a ~1.8 mas/yr accuracy while the remaining objects were measured only from auxiliary data with a mean precision of about ~3 mas/yr. We performed a kinematic analysis of the most accurate proper motions derived in this work, which allowed us to separate cluster members from field stars and to derive the mean properties of the cluster. From the kinematic analysis we derived a list of 68 members and 14 candidate members, comprising 26 new objects with a high membership probability. These new members are generally fainter than the known ones. We measured a mean proper motion of (mu_RA*, mu_DEC)=(-8.2, -24.3)+/-0.8 mas/yr for the L1688 dark cloud. A supervised classification was applied to photometric data of members to allocate a spectral energy distribution (SED) classification to the unclassified members.} We kinematically confirmed that the 56 members that were known from previous studies of the Rho Ophiuchi F cluster and that were also part of our survey are members of the cluster, and we added 26 new members. We defined the evolutionary status of the unclassified members of the cluster. We showed that a large part (23) of these new members are probably brown dwarfs, which multiplies the number of known substellar objects in the cluster by a factor of 3.3.
• ### The seven sisters DANCe.II. Proper motions and the lithium-rotation-activity connection for G and K Pleiads(1608.02751)

Aug. 9, 2016 astro-ph.SR
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.
• ### Search for light curve modulations among Kepler candidates. Three very low-mass transiting companions(1606.02398)

June 8, 2016 astro-ph.SR, astro-ph.EP
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.
• ### Rotation in the Pleiades With K2: III. Speculations on Origins and Evolution(1606.00057)

June 2, 2016 astro-ph.SR
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.
• ### Rotation in the Pleiades with K2: I. Data and First Results(1606.00052)

May 31, 2016 astro-ph.SR
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.
• ### Rotation in the Pleiades with K2: II. Multi-Period Stars(1606.00055)

May 31, 2016 astro-ph.SR
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.
• ### Constraining the properties of transitional disks in Chamaeleon I with Herschel(1602.03183)

Transitional disks are protoplanetary disks with opacity gaps/cavities in their dust distribution, a feature that may be linked to planet formation. We perform Bayesian modeling of the three transitional disks SZ Cha, CS Cha and T25 including photometry from the Herschel Space Observatory to quantify the improvements added by these new data. We find disk dust masses between 2x10^-5 and 4x10^-4 Msun, and gap radii in the range of 7-18 AU, with uncertainties of ~ one order of magnitude and ~ 4 AU, respectively. Our results show that adding Herschel data can significantly improve these estimates with respect to mid-infrared data alone, which have roughly twice as large uncertainties on both disk mass and gap radius. We also find weak evidence for different density profiles with respect to full disks. These results open exciting new possibilities to study the distribution of disk masses for large samples of disks.
• ### A submillimeter search for pre- and proto-brown dwarfs in Chamaeleon II(1512.00418)

Dec. 1, 2015 astro-ph.SR
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.
• ### Herschel-Planck dust optical depth and column density maps - II. Perseus(1511.08503)

Nov. 26, 2015 astro-ph.SR
We present optical depth and temperature maps of the Perseus molecular cloud, obtained combining dust emission data from the Herschel and Planck satellites and 2MASS/NIR dust extinction maps. The maps have a resolution of 36 arcsec in the Herschel regions, and of 5 arcmin elsewhere. The dynamic range of the optical depth map ranges from $1\times10^{-2}\, \mathrm{mag}$ up to $20 \,\mathrm{mag}$ in the equivalent K band extinction. We also evaluate the ratio between the $2.2 \,\mathrm{\mu m}$ extinction coefficient and the $850 \,\mathrm{\mu m}$ opacity. The value we obtain is close to the one found in the Orion B molecular cloud. We show that the cumulative and the differential area function of the data (which is proportional to the probability distribution function of the cloud column density) follow power laws with index respectively $\simeq -2$, and $\simeq -3$. We use WISE data to improve current YSO catalogues based mostly on \emph{Spitzer} data and we build an up-to-date selection of Class~I/0 objects. Using this selection, we evaluate the local Schmidt law, $\Sigma_{\mathrm{YSO}} \propto \Sigma_{\mathrm{gas}}^{\beta}$, showing that $\beta=2.4 \pm 0.6$. Finally, we show that the area-extinction relation is important for determining the star formation rate in the cloud, which is in agreement with other recent works.
• ### Parallax of the L4.5 dwarf 2M1821$+$14 from high-precision astrometry with OSIRIS at GTC(1509.08718)

Sept. 29, 2015 astro-ph.SR, astro-ph.EP
We used the OSIRIS camera at the 10.4 m Gran Telescopio Canarias (GTC) to monitor the astrometric motion of the L4.5 dwarf 2M1821$+$14 over 17 months. The astrometric residuals of eleven epochs have a r.m.s. dispersion of 0.4 mas, which is larger than the average precision of 0.23 mas per epoch and hints towards an additional signal or excess noise. Comparison of the point-spread-functions in OSIRIS and FORS2/VLT images reveals no differences critical for high-precision astrometry, despite the GTC's segmented primary mirror. We attribute the excess noise to an unknown effect that may be uncovered with additional data. For 2M1821$+$14, we measured a relative parallax of $106.15 \pm 0.18$ mas and determined a correction of $0.50\pm0.05$ mas to absolute parallax, leading to a distance of $9.38 \pm0.03$ pc. We excluded at 3-$\sigma$ confidence the presence of a companion to 2M1821$+$14 down to a mass ratio of 0.1 ($\approx 5\, M_\mathrm{Jupiter}$) with a period of 50--1000 days and a separation of 0.1--0.7 au. The accurate parallax allowed us to estimate the age and mass of 2M1821$+$14 of 120--700 Myr and 0.049$^{+0.014}_{-0.024}$ M$_\odot$, thus confirming its intermediate age and substellar mass. We complement our study with a parallax and proper motion catalogue of 587 stars ($i'\simeq15.5-22$) close to 2M1821$+$14, used as astrometric references. This study demonstrates sub-mas astrometry with the GTC, a capability applicable for a variety of science cases including the search for extrasolar planets and relevant for future astrometric observations with E-ELT and TMT.
• ### YSOVAR: Mid-Infrared Variability in NGC 1333(1508.04705)

Aug. 19, 2015 astro-ph.SR
As part of the Young Stellar Object VARiability (YSOVAR) program, we monitored NGC 1333 for ~35 days at 3.6 and 4.5 um using the Spitzer Space Telescope. We report here on the mid-infrared variability of the point sources in the ~10x~20arcmin area centered on 03:29:06, +31:19:30 (J2000). Out of 701 light curves in either channel, we find 78 variables over the YSOVAR campaign. About half of the members are variable. The variable fraction for the most embedded SEDs (Class I, flat) is higher than that for less embedded SEDs (Class II), which is in turn higher than the star-like SEDs (Class III). A few objects have amplitudes (10-90th percentile brightness) in [3.6] or [4.5]>0.2 mag; a more typical amplitude is 0.1-0.15 mag. The largest color change is >0.2 mag. There are 24 periodic objects, with 40% of them being flat SED class. This may mean that the periodic signal is primarily from the disk, not the photosphere, in those cases. We find 9 variables likely to be 'dippers', where texture in the disk occults the central star, and 11 likely to be 'bursters', where accretion instabilities create brightness bursts. There are 39 objects that have significant trends in [3.6]-[4.5] color over the campaign, about evenly divided between redder-when-fainter (consistent with extinction variations) and bluer-when-fainter. About a third of the 17 Class 0 and/or jet-driving sources from the literature are variable over the YSOVAR campaign, and a larger fraction (~half) are variable between the YSOVAR campaign and the cryogenic-era Spitzer observations (6-7 years), perhaps because it takes time for the envelope to respond to changes in the central source. The NGC 1333 brown dwarfs do not stand out from the stellar light curves in any way except there is a much larger fraction of periodic objects (~60% of variable brown dwarfs are periodic, compared to ~30% of the variables overall).
• ### First detection of thermal radio jets in a sample of proto-brown dwarf candidates(1505.07586)

May 28, 2015 astro-ph.GA
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.
• ### Search for associations containing young stars (SACY). VI. Is multiplicity universal? Stellar multiplicity in the range 3-1000 au from adaptive-optics observations(1505.07837)

May 28, 2015 astro-ph.SR
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.