In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 \mu$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

In this paper, we address two issues related to primordial disk evolution in three clusters (NGC 1333, IC 348, and Orion A) observed by the INfrared Spectra of Young Nebulous Clusters (IN-SYNC) project. First, in each cluster, averaged over the spread of age, we investigate how disk lifetime is dependent on stellar mass. The general relation in IC 348 and Orion A is that primordial disks around intermediate mass stars (2--5$M_{\odot}$) evolve faster than those around loss mass stars (0.1--1$M_{\odot}$), which is consistent with previous results. However, considering only low mass stars, we do not find a significant dependence of disk frequency on stellar mass. These results can help to better constrain theories on gas giant planet formation timescales. Secondly, in the Orion A molecular cloud, in the mass range of 0.35--0.7$M_{\odot}$, we provide the most robust evidence to date for disk evolution within a single cluster exhibiting modest age spread. By using surface gravity as an age indicator and employing 4.5 $\mu m$ excess as a primordial disk diagnostic, we observe a trend of decreasing disk frequency for older stars. The detection of intra-cluster disk evolution in NGC 1333 and IC 348 is tentative, since the slight decrease of disk frequency for older stars is a less than 1-$\sigma$ effect.

We present an analysis of spectrosopic and astrometric data from APOGEE-2 and Gaia DR2 to identify structures towards the Orion Complex. By applying a hierarchical clustering algorithm to the 6-dimensional stellar data, we identify spatially and/or kinematically distinct groups of young stellar objects with ages ranging from 1 to 12 Myr. We also investigate the star forming history within the Orion Complex, and identify peculiar sub-clusters. With this method we reconstruct the older populations in the regions that are presently largely devoid of molecular gas, such as Orion C (which includes the $\sigma$ Ori cluster), and Orion D (the population that traces Ori OB1a, OB1b, and Orion X). We report on the distances, kinematics, and ages of the groups within the Complex. The Orion D groups is in the process of expanding. On the other hand, Orion B is still in the process of contraction. In $\lambda$ Ori the proper motions are consistent with a radial expansion due to an explosion from a supernova; the traceback age from the expansion exceeds the age of the youngest stars formed near the outer edges of the region, and their formation would have been triggered when they were half-way from the cluster center to their current positions. We also present a comparison between the parallax and proper motion solutions obtained by Gaia DR2, and those obtained towards star-forming regions by Very Long Baseline Array.

The Orion Star Forming Complex (OSFC) is a central target for the APOGEE-2 Young Cluster Survey. Existing membership catalogs span limited portions of the OSFC, reflecting the difficulty of selecting targets homogeneously across this extended, highly structured region. We have used data from wide field photometric surveys to produce a less biased parent sample of young stellar objects (YSOs) with infrared (IR) excesses indicative of warm circumstellar material or photometric variability at optical wavelengths across the full 420 square degrees extent of the OSFC. When restricted to YSO candidates with H < 12.4, to ensure S/N ~100 for a six visit source, this uniformly selected sample includes 1307 IR excess sources selected using criteria vetted by Koenig & Liesawitz and 990 optical variables identified in the Pan-STARRS1 3$\pi$ survey: 319 sources exhibit both optical variability and evidence of circumstellar disks through IR excess. Objects from this uniformly selected sample received the highest priority for targeting, but required fewer than half of the fibers on each APOGEE-2 plate. We fill the remaining fibers with previously confirmed and new color-magnitude selected candidate OSFC members. Radial velocity measurements from APOGEE-1 and new APOGEE-2 observations taken in the survey's first year indicate that ~90% of the uniformly selected targets have radial velocities consistent with Orion membership.The APOGEE-2 Orion survey will include >1100 bona fide YSOs whose uniform selection function will provide a robust sample for comparative analyses of the stellar populations and properties across all sub-regions of Orion.

We study the core mass function (CMF) of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array via 1.3~mm continuum emission at a resolution of 1.0\arcsec\ (2500~au). We have mapped a field of 5.3\arcmin$\times$5.3\arcmin\ centered on the protocluster clump. We measure the CMF in the central region, exploring various core detection algorithms, which give source numbers ranging from 60 to 125, depending on parameter selection. We estimate completeness corrections due to imperfect flux recovery and core identification via artificial core insertion experiments. For masses $M\gtrsim1\:M_\odot$, the fiducial dendrogram-identified CMF can be fit with a power law of the form ${\rm{d}}N/{\rm{d}}{\rm{log}}M\propto{M}^{-\alpha}$ with $\alpha \simeq1.24\pm0.17$, slightly shallower than, but still consistent with, the index of the Salpeter stellar initial mass function of 1.35. Clumpfind-identified CMFs are significantly shallower with $\alpha\simeq0.64\pm0.13$. While raw CMFs show a peak near $1\:M_\odot$, completeness-corrected CMFs are consistent with a single power law extending down to $\sim 0.5\:M_\odot$, with only a tentative indication of a shallowing of the slope around $\sim1\:M_\odot$. We discuss the implications of these results for star and star cluster formation theories.

We study the occurrence of spectroscopic binaries in young star-forming regions using the INfrared Spectroscopy of Young Nebulous Clusters(IN-SYNC) survey, carried out in SDSS-III with the APOGEE spectrograph. Multi-epoch observations of thousands of low-mass stars in Orion A, NGC 2264, NGC 1333, IC 348, and the Pleiades have been carried out, yielding H-band spectra with a nominal resolution of R=22,500 for sources with H $\le$ 12 mag. Radial velocity precisions of $\sim$0.3 $km\:s^{-1}$ were achieved, which we use to identify radial velocity variations indicative of undetected companions. We use Monte Carlo simulations to assess the types of spectroscopic binaries to which we are sensitive, finding sensitivity to binaries with orbital periods $< 10^{4}$ d, for stars with $2500 {\rm K} \le T_\mathrm{eff} \le 6000 {\rm K}$ and $\it{v} \sin \it{i}$ $\le$ 100 $km\:s^{-1}$. Using Bayesian inference, we find evidence for a decline in the spectroscopic binary fraction, by a factor of 3-4 from the age of our pre-main-sequence sample to the Pleiades age . The significance of this decline is weakened if spot-induced radial-velocity jitter is strong in the sample, and is only marginally significant when comparing any one of the pre-main-sequence clusters against the Pleiades. However, the same decline in both sense and magnitude is found for each of the five pre-main-sequence clusters, and the decline reaches statistical significance of greater than 95% confidence when considering the pre-main-sequence clusters jointly. Our results suggest that dynamical processes disrupt the widest spectroscopic binaries ($P_{\rm orb} \approx 10^3 - 10^4$ d) as clusters age, indicating that this occurs early in the stars' evolution, while they still reside within their nascent clusters.

The kinematics and dynamics of young stellar populations enable us to test theories of star formation. With this aim, we continue our analysis of the SDSS-III/APOGEE IN-SYNC survey, a high resolution near infrared spectroscopic survey of young clusters. We focus on the Orion A star-forming region, for which IN-SYNC obtained spectra of $\sim2700$ stars. In Paper IV we used these data to study the young stellar population. Here we study the kinematic properties through radial velocities ($v_r$). The young stellar population remains kinematically associated with the molecular gas, following a $\sim10\:{\rm{km\:s}}^{-1}$ gradient along filament. However, near the center of the region, the $v_r$ distribution is slightly blueshifted and asymmetric; we suggest that this population, which is older, is slightly in foreground. We find evidence for kinematic subclustering, detecting statistically significant groupings of co-located stars with coherent motions. These are mostly in the lower-density regions of the cloud, while the ONC radial velocities are smoothly distributed, consistent with it being an older, more dynamically evolved cluster. The velocity dispersion $\sigma_v$ varies along the filament. The ONC appears virialized, or just slightly supervirial, consistent with an old dynamical age. Here there is also some evidence for on-going expansion, from a $v_r$--extinction correlation. In the southern filament, $\sigma_v$ is $\sim2$--$3$ times larger than virial in the L1641N region, where we infer a superposition along the line of sight of stellar sub-populations, detached from the gas. On the contrary, $\sigma_v$ decreases towards L1641S, where the population is again in agreement with a virial state.

We present radial velocity measurements for 70 high confidence, and 34 potential binary systems in fields containing the Perseus Molecular Cloud, Pleiades, NGC 2264, and the Orion A star forming region. 18 of these systems have been previously identified as binaries in the literature. Candidate double-lined spectroscopic binaries (SB2s) are identified by analyzing the cross-correlation functions (CCFs) computed during the reduction of each APOGEE spectrum. We identify sources whose CCFs are well fit as the sum of two Lorentzians as likely binaries, and provide an initial characterization of the system based on the radial velocities indicated by that dual fit. For systems observed over several epochs, we present mass ratios and systemic velocities; for two systems with observations on eight or more epochs, and which meet our criteria for robust orbital coverage, we derive initial orbital parameters. The distribution of mass ratios for multi-epoch sources in our sample peaks at q=1, but with a significant tail toward lower q values. Tables reporting radial velocities, systemic velocities, and mass ratios are provided online. We discuss future improvements to the radial velocity extraction method we employ, as well as limitations imposed by the number of epochs currently available in the APOGEE database. The Appendix contains brief notes from the literature on each system in the sample, and more extensive notes for select sources of interest.

We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accurate stellar parameters (T_eff, log g, v sin i) and extinctions, and placed the sources in the Hertzsprung-Russel Diagram (HRD). We have also extracted radial velocities for the kinematic characterization of the population. We compare our measurements with literature results for a sub-sample of targets in order to assess the performances and accuracy of the survey. Source extinction shows evidence for dust grains that are larger than those in the diffuse interstellar medium (ISM): we estimate an average R_V=5.5 in the region. Importantly, we find a clear correlation between HRD inferred ages and spectroscopic surface-gravity inferred ages. This clearly indicates a real spread of stellar radii at fixed temperature, and together with additional correlations with extinction and with disk presence, strongly suggests a real spread of ages large than a few Myr. Focussing on the young population around NGC1980 iota Ori, which has previously been suggested to be a separate, foreground, older cluster, we confirm its older (5Myr) age and low A_V, but considering that its radial velocity distribution is indistinguishable from the Orion A's population, we suggest that NGC1980 is part of Orion A's star formation activity. Based on their stellar parameters and kinematic properties, we identify 383 new candidate members of Orion A, most of which are diskless sources in areas of the region poorly studied by previous works.

Age and mass determinations for isolated stellar objects remain model-dependent. While stellar interior and atmospheric theoretical models are rapidly evolving, we need a powerful tool to test them. Open clusters are good candidates for this role. We complement previous studies on the Hyades multiplicity by Lucky Imaging observations with the AstraLux Norte camera. This allows us to exclude possible binary and multiple systems with companions outside 2--7 AU separation and to create a "single-star sequence" for the Hyades. The sequence encompasses 250 main-sequence stars ranging from A5V to M6V. Using the "Tool for Astrophysical Data Analysis" (TA-DA), we create various theoretical isochrones applying different combinations of interior and atmospheric models. We compare the isochrones with the observed Hyades single-star sequence on J vs. J - K_s, J vs. J - H and K_s vs. H - K_s color-magnitude diagrams. As a reference we also compute absolute fluxes and magnitudes for all stars from X-ray to mid-infrared based on photometric measurements available in the literature(ROSAT X-ray, GALEX UV, APASS gri, 2MASS JHK_s, and WISE W1 to W).We find that combinations of both PISA and DARTMOUTH stellar interior models with BT-Settl 2010 atmospheric models describe the observed sequence well. The full sequence covers the mass range 0.13 to 2.3 Msun, and effective temperatures between 3060 K and 8200 K. Within the measurement uncertainties, the current generation of models agree well with the single-star sequence. The primary limitations are the uncertainties in the measurement of the distance to individual Hyades members, and uncertainties in the photometry. Additionally, a small (~0.05 mag) systematic offset can be noted on J vs. J - K and K vs. H - K diagrams - the observed sequence is shifted to redder colors from the theoretical predictions.

Most field stars will have encountered the highest stellar density and hence the largest number of interactions in their birth environment. Yet the stellar dynamics during this crucial phase are poorly understood. Here we analyze the radial velocities measured for 152 out of 380 observed stars in the 2-6 Myr old star cluster IC 348 as part of the SDSS-III APOGEE. The radial velocity distribution of these stars is fitted with one or two Gaussians, convolved with the measurement uncertainties including binary orbital motions. Including a second Gaussian improves the fit; the high-velocity outliers that are best fit by this second component may either (1) be contaminants from the nearby Perseus OB2 association, (2) be a halo of ejected or dispersing stars from IC 348, or (3) reflect that IC 348 has not relaxed to a Gaussian velocity distribution. We measure a velocity dispersion for IC 348 of $0.72 \pm 0.07$ km s$^{-1}$ (or $0.64 \pm 0.08$ km s$^{-1}$ if two Gaussians are fitted), which implies a supervirial state, unless the gas contributes more to the gravitational potential than expected. No evidence is found for a dependence of this velocity dispersion on distance from the cluster center or stellar mass. We also find that stars with lower extinction (in the front of the cloud) tend to be redshifted compared with stars with somewhat higher extinction (towards the back of the cloud). This data suggests that the stars in IC 348 are converging along the line of sight. We show that this correlation between radial velocity and extinction is unlikely to be spuriously caused by the small cluster rotation of $0.024 \pm 0.013$ km s$^{-1}$ arcmin$^{-1}$ or by correlations between the radial velocities of neighboring stars. This signature, if confirmed, will be the first detection of line-of-sight convergence in a star cluster(...)

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial velocity of the region and the diffuse gas velocity dispersion, but significantly larger than the velocity dispersion of the dense, star-forming cores, which have a sub-virial velocity dispersion of 0.5 km/s. Since the NGC 1333 cluster is dynamically young and deeply embedded, this measurement provides a strong constraint on the initial velocity dispersion of newly-formed stars. We propose that the difference in velocity dispersion between stars and dense cores may be due to the influence of a 70 micro-Gauss magnetic field acting on the dense cores, or be the signature of a cluster with initial sub-structure undergoing global collapse.

We study the degree of angular substructure in the stellar position distribution of young members of Galactic star-forming regions, looking for correlations with distance from cluster center, surface number density of stars, and local dynamical age. To this end we adopt the catalog of members in 18 young ($\sim$1-3 Myr) clusters from the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Survey and the statistical analysis of the Angular Dispersion Parameter, $\delta_{\rm ADP}$. We find statistically significant correlation between $\delta_{\rm ADP}$ and physical projected distance from the center of the clusters, with the centers appearing smoother than the outskirts, consistent with more rapid dynamical processing on local dynamical, free-fall or orbital timescales. Similarly, smoother distributions are seen in regions of higher surface density, or older dynamical ages. These results indicate that dynamical processing that erases substructure is already well-advanced in young, sometimes still-forming, clusters. Such observations of the dissipation of substructure have the potential to constrain theoretical models of the dynamical evolution of young and forming clusters.

The spatial morphology and dynamical status of a young, still-forming stellar cluster provide valuable clues on the conditions during the star formation event and the processes that regulated it. We analyze the Orion Nebula Cluster (ONC), utilizing the latest censuses of its stellar content and membership estimates over a large wavelength range. We determine the center of mass of the ONC, and study the radial dependence of angular substructure. The core appears rounder and smoother than the outskirts, consistent with a higher degree of dynamical processing. At larger distances the departure from circular symmetry is mostly driven by the elongation of the system, with very little additional substructure, indicating a somewhat evolved spatial morphology or an expanding halo. We determine the mass density profile of the cluster, which is well fitted by a power law that is slightly steeper than a singular isothermal sphere. Together with the ISM density, estimated from average stellar extinction, the mass content of the ONC is insufficient by a factor $\sim 1.8$ to reproduce the observed velocity dispersion from virialized motions, in agreement with previous assessments that the ONC is moderately supervirial. This may indicate recent gas dispersal. Based on the latest estimates for the age spread in the system and our density profiles, we find that, at the half-mass radius, 90% of the stellar population formed within $\sim 5$-$8$ free-fall times ($t_{\rm ff}$). This implies a star formation efficiency per $t_{\rm ff}$ of $\epsilon_{\rm ff}\sim 0.04$-$0.07$, i.e., relatively slow and inefficient star formation rates during star cluster formation.

Over two years 8,859 high-resolution H-band spectra of 3,493 young (1 - 10 Myr) stars were gathered by the multi-object spectrograph of the APOGEE project as part of the IN-SYNC ancillary program of that SDSS-III survey. Here we present the forward modeling approach used to derive effective temperatures, surface gravities, radial velocities, rotational velocities, and H-band veiling from these near-infrared spectra. We discuss in detail the statistical and systematic uncertainties in these stellar parameters. In addition we present accurate extinctions by measuring the E(J-H) of these young stars with respect to the single-star photometric locus in the Pleiades. Finally we identify an intrinsic stellar radius spread of about 25% for late-type stars in IC 348 using three (nearly) independent measures of stellar radius, namely the extinction-corrected J-band magnitude, the surface gravity and the $R \sin i$ from the rotational velocities and literature rotation periods. We exclude that this spread is caused by uncertainties in the stellar parameters by showing that the three estimators of stellar radius are correlated, so that brighter stars tend to have lower surface gravities and larger $R \sin i$ than fainter stars at the same effective temperature. Tables providing the spectral and photometric parameters for the Pleiades and IC 348 have been provided online.

Calibration and benchmarking of evolutionary and atmospheric models is essential for the study of low-mass stars and brown dwarfs, as for isolated objects these models are the only way to determine basic parameters like mass and age. The Hyades star cluster with an age of around 625 Myr and located at the distance of ~ 45 pc, is the most accessible cluster in the solar neighborhood. R\"oser et al. (2011) establish a list of 724 likely Hyades members. Using available literature data (HST, HIPPARCOS, WDS, Patience et al. 1998, Mermilliod et al. 2009, Morzinski 2011) and our own lucky imaging observations with AstraLux at the 2.2m telescope in Calar Alto, we establish a single-star sequence containing 255 Hyades members spanning the mass range ~ 0.2--1.5 solar masses . This sequence is used for testing and calibration of various existing stellar and atmospheric models (PADOVA, DARTMOUTH, BCAH, BT-Settl)

The temporal decay of mass accretion in young stars is a fundamental tracer of the early evolution of circumstellar disks. Through population syntheses, we study how correlated uncertainties between the estimated parameters of young stars (luminosity, temperature, mass, age) and mass accretion rates Mdot, as well as observational selection effects, can bias the temporal decay of mass accretion rates (Mdot t^-eta) inferred from a comparison of measured Mdot with isochronal ages in young stellar clusters. We find that the presence of realistic uncertainties reduces the measured value of eta by up to a factor of 3, leading to the inference of shallower decays than the true value. This suggests a much faster temporal decay of Mdot than generally assumed. When considering the minimum uncertainties in ages affecting the Orion Nebula Cluster, the observed value eta~1.4, typical of Galactic star forming regions, can only be reproduced if the real decay exponent is eta>4. This effect becomes more severe if one assumes that observational uncertainties are larger, as required by some fast star formation scenarios. Our analysis shows that while selection effects due to sample incompleteness do bias eta, they can not alter this main result and strengthen it in many cases. A remaining uncertainty in our work is that it applies to the most commonly used and simple relationship between Mdot, the accretion luminosity and the stellar parameters. We briefly explore how a more complex interplay between these quantities might change the results.

Current planet formation models are largely based on the observational constraint that protoplanetary disks have lifetime 3Myr. Recent studies, however, report the existence of PMS stars with signatures of accretion (strictly connected with the presence of circumstellar disks)and photometrically determined ages of 30 Myr, or more. Here we present a spectroscopic study of two major age outliers in the ONC. We use broad band, intermediate resolution VLT/X-Shooter spectra combined with an accurate method to determine the stellar parameters and the related age of the targets to confirm their peculiar age estimates and the presence of ongoing accretion.The analysis is based on a multi-component fitting technique, which derives simultaneously SpT, extinction, and accretion properties of the objects. With this method we confirm and quantify the ongoing accretion. From the photospheric parameters of the stars we derive their position on the HRD, and the age given by evolutionary models. Together with other age indicators like the lithium equivalent width we estimate with high accuracy the age of the objects. Our study shows that the two objects analyzed are not older than the typical population of the ONC. Our results show that, while photometric determination of the photospheric parameters are an accurate method to estimate the parameters of the bulk of young stellar populations, those of individual objects with high accretion rates and extinction may be affected by large uncertainties. Broad band spectroscopic determinations should thus be used to confirm the nature of individual objects. Our analysis shows that this method allows us to obtain an accurate determination of the photospheric parameters of accreting YSOs in any nearby star-forming region. We suggest that our detailed, broad- band spectroscopy method should be used to derive accurate properties of candidate old and accreting YSOs.

The Hubble Space Telescope (HST) Treasury Program on the Orion Nebula Cluster has used 104 orbits of HST time to image the Great Orion Nebula region with the Advanced Camera for Surveys (ACS), the Wide-Field/Planetary Camera 2 (WFPC2) and the Near Infrared Camera and Multi Object Spectrograph (NICMOS) instruments in 11 filters ranging from the U-band to the H-band equivalent of HST. The program has been intended to perform the definitive study of the stellar component of the ONC at visible wavelengths, addressing key questions like the cluster IMF, age spread, mass accretion, binarity and cirumstellar disk evolution. The scanning pattern allowed to cover a contiguous field of approximately 600 square arcminutes with both ACS and WFPC2, with a typical exposure time of approximately 11 minutes per ACS filter, corresponding to a point source depth AB(F435W) = 25.8 and AB(F775W)=25.2 with 0.2 magnitudes of photometric error. We describe the observations, data reduction and data products, including images, source catalogs and tools for quick look preview. In particular, we provide ACS photometry for 3399 stars, most of them detected at multiple epochs, WFPC2 photometry for 1643 stars, 1021 of them detected in the U-band, and NICMOS JH photometry for 2116 stars. We summarize the early science results that have been presented in a number of papers. The final set of images and the photometric catalogs are publicly available through the archive as High Level Science Products at the STScI Multimission Archive hosted by the Space Telescope Science Institute.

We present TA-DA, a new software aimed at greatly simplify and improve the analysis of stellar photometric data in comparison with theoretical models, and allow the derivation of stellar parameters from multi-band photometry. Its flexibility allows one to address a number of such problems: from the interpolation of stellar models, or sets of stellar physical parameters in general, to the computation of synthetic photometry in arbitrary filters or units; from the analysis of observed color-magnitude diagrams, to a Bayesian derivation of stellar parameters (and extinction) based on multi-band data. TA-DA is available as a pre-compiled IDL widget-based application; its graphical user interface makes it considerably user-friendly. In this paper we describe the software and its functionalities.

Atmospheric and evolutionary models for low-mass stars rely on approximate assumptions on the physics of the stellar structure and the atmospheric radiative transfer. This leads to biased theoretical predictions on the photospheric Spectral Energy Distributions of Pre-Main Sequence (PMS) stars, and affects the derivation of stellar parameters from photometric data. Our goal is to correct the biases present in the theoretical predictions for the near-IR photometry of low-mass PMS stars. Using empirical intrinsic IR colors, we assess the accuracy of current synthetic spectral libraries and evolutionary models. We consider a sample of ~300 PMS stars in the Orion Nebula Cluster (age 1 Myr) with measured luminosities, temperatures and photospheric JHKs photometry. By analyzing the photospheric colors of our sample of stars, we find that the synthetic JHKs photometry provided by theoretical spectral templates for late spectral types (>K6) are accurate at the level of ~0.2 mag, while colors are accurate at ~0.1 mag. We tabulate the intrinsic photospheric colors, appropriate for the Orion Nebula Cluster, in the range K6-M8.5. They can be conveniently used as templates for the intrinsic colors of other young (age<5 Myr) stellar clusters. An empirical correction of the atmospheric templates can fix the discrepancies between expected and observed colors. Still, other biases in the evolutionary models prevent a more robust comparison between observations and theoretical absolute magnitudes. In particular, PMS evolutionary models seem to consistently underestimate the intrinsic near-infrared flux at the very late spectral types, and this may introduce spurious features in the low-mass end of the photometrically-determined Initial Mass Function of young clusters.

The present observational understanding of the evolution of the mass accretion rates (Macc) in pre-main sequence stars is limited by the lack of accurate measurements of Macc over homogeneous and large statistical samples of young stars. Such observational effort is needed to properly constrain the theory of star formation and disk evolution. Based on HST/WFPC2 observations, we present a study of Macc for a sample of \sim 700 sources in the Orion Nebula Cluster, ranging from the Hydrogen-burning limit to M\ast \sim 2M\odot. We derive Macc from both the U-band excess and the H{\alpha} luminosity (LH{\alpha}), after determining empirically both the shape of the typical accretion spectrum across the Balmer jump and the relation between the accretion luminosity (Lacc) and LH{\alpha}, that is Lacc/L\odot = (1.31\pm0.03)\cdotLH{\alpha}/L\odot + (2.63\pm 0.13). Given our large statistical sample, we are able to accurately investigate relations between Macc and the parameters of the central star such as mass and age. We clearly find Macc to increase with stellar mass, and decrease over evolutionary time, but we also find strong evidence that the decay of Macc with stellar age occurs over longer timescales for more massive PMS stars. Our best fit relation between these parameters is given by: log(Macc/M\odot\cdotyr)=(-5.12 \pm 0.86) -(0.46 \pm 0.13) \cdot log(t/yr) -(5.75 \pm 1.47)\cdot log(M\ast/M\odot) + (1.17 \pm 0.23)\cdot log(t/yr) \cdot log(M\ast/M\odot). These results also suggest that the similarity solution model could be revised for sources with M\ast > 0.5M\odot. Finally, we do not find a clear trend indicating environmental effects on the accretion properties of the sources.

We present new HST/WFPC2 observations of 20 fields centered around T Tauri stars in the Chamaeleon I star forming region. Images have been obtained in the F631N ([OI]6300A), F656N (Ha) and F673N ([SII]6716A+6731A) narrow-band filters, plus the Johnson V-band equivalent F547M filter. We detect 31 T Tauri stars falling within our fields. We discuss the optical morphology of 10 sources showing evidence of either binarity, circumstellar material, or mass loss. We supplement our photometry with a compilation of optical, infrared and sub-millimeter data from the literature, together with new sub-mm data for three objects, to build the Spectral Energy Distributions (SED) of 19 single sources. Using an SED model fitting tool, we self-consistently estimate a number of stellar and disk parameters, while mass accretion rates are directly derived from our Ha photometry. We find that bolometric luminosities derived from dereddened optical data tend to be underestimated in systems with high alpha(2-24} IR spectral index, suggesting that disks seen nearly edge-on may occasionally be interpreted as low luminosity (and therefore more evolved) sources. On the other hand, the same alpha(2-24) spectral index, a tracer of the amount of dust in the warmer layers of the circumstellar disks, and the mass accretion rate appear to decay with the isocronal stellar age, suggesting that the observed age spread (~0.5-5 Myr) within the cluster is real. Our sample contains a few outliers that may have dissipated their circumstellar disks on shorter time-scale.

We present a new method to determine the age spread of resolved stellar populations in a starburst cluster. The method relies on a two-step process. In the first step, kinematic members of the cluster are identified based on multi-epoch astrometric monitoring. In the second step, a Bayesian analysis is carried out, comparing the observed photometric sequence of cluster members with sets of theoretical isochrones. When applying this methodology to optical and near-infrared high angular resolution Hubble Space Telescope (HST) and adaptive optics observations of the ~5 Myr old starburst cluster Westerlund 1 and ~2 Myr old starburst cluster NGC 3603 YC, we derive upper limits for the age spreads of 0.4 and 0.1 Myr, respectively. The results strongly suggest that star formation in these starburst clusters happened almost instantaneously.

We present a spectroscopic study of the most massive stars in the young (4 Myr old) stellar cluster LH 95 in the Large Magellanic Cloud. This analysis allows us to complete the census of the stellar population of the system, previously investigated by us down to 0.4 solar masses with deep HST Advanced Camera for Surveys photometry. We perform spectral classification of the five stars in our sample, based on high resolution optical spectroscopy obtained with 2.2m MPG/ESO FEROS. We use complementary ground-based photometry, previously performed by us, to place these stars in the Hertzsprung-Russel diagram. We derive their masses and ages by interpolation from evolutionary models. The average ages and age spread of the most massive stars are found to be in general comparable with those previously derived for the cluster from its low mass PMS stars. We use the masses of the 5 sample stars to extend to the high-mass end the stellar initial mass function of LH 95 previously established by us. We find that the initial mass function follows a Salpeter relation down to the intermediate-mass regime at 2 Msun. The second most massive star in LH 95 shows broad Balmer line emission and infrared excess, which are compatible with a classical Be star. The existence of such a star in the system adds a constrain to the age of the cluster, which is well covered by our age and age spread determinations. The most massive star, a 60-70 Msun O2 giant is found to be younger (<1 Myr) than the rest of the population. Its mass in relation to the total mass of the system does not follow the empirical relation of the maximum stellar mass versus the hosting cluster mass, making LH 95 an exception to the average trend.

We present a new census of the Orion Nebula Cluster (ONC) over a large field of view (>30'x30'), significantly increasing the known population of stellar and substellar cluster members with precisely determined properties. We develop and exploit a technique to determine stellar effective temperatures from optical colors, nearly doubling the previously available number of objects with effective temperature determinations in this benchmark cluster. Our technique utilizes colors from deep photometry in the I-band and in two medium-band filters at lambda~753 and 770nm, which accurately measure the depth of a molecular feature present in the spectra of cool stars. From these colors we can derive effective temperatures with a precision corresponding to better than one-half spectral subtype, and importantly this precision is independent of the extinction to the individual stars. Also, because this technique utilizes only photometry redward of 750nm, the results are only mildly sensitive to optical veiling produced by accretion. Completing our census with previously available data, we place some 1750 sources in the Hertzsprung-Russel diagram and assign masses and ages down to 0.02 solar masses. At faint luminosities, we detect a large population of background sources which is easily separated in our photometry from the bona fide cluster members. The resulting initial mass function of the cluster has good completeness well into the substellar mass range, and we find that it declines steeply with decreasing mass. This suggests a deficiency of newly formed brown dwarfs in the cluster compared to the Galactic disk population.

Aims. We present a study of the distribution of stellar ages in the Orion Nebula Cluster (ONC) based on accurate HST photometry taken from the HST Treasury Program observations of the ONC utilizing the most recent estimate of the cluster's distance (Menten et al. 2007). We investigate the presence of an intrinsic age spread in the region and a possible trend of age with the spatial distribution. Methods. We estimate the extinction and accretion luminosity towards each source by performing synthetic photometry on an empirical calibration of atmospheric models (Da Rio et al. 2010) using the package Chorizos (Maiz-Apellaniz 2004). The position of the sources in the HR-diagram is compared with different theoretical isochrones to estimate the mean cluster age and age dispersion. Through Monte Carlo simulations we quantify the amount of intrinsic age spread in the region, taking into account uncertainties on the distance, spectral type, extinction, unresolved binaries, accretion and photometric variability. Results. According to Siess et al. (2000) evolutionary models the mean age of the Cluster is 2.2 Myr with a scatter of few Myrs. With Monte Carlo simulations we find that the observed age spread is inconsistent with a coeval stellar population, but is in agreement with a star formation activity between 1.5 and 3.5 Myrs. We also observe light evidence for a trend of ages with spatial distribution.

MAD is the first multi-conjugated adaptive optics system at the VLT. We present H and Ks observations of the young massive cluster Trumpler 14 revealing the power of MCAO systems by providing a homogeneous Strehl ratio over a large field of view. Mean Strehl ratios of 6.0 and 5.9 per cent with maximum Strehl ratios of 9.8 and 12.6 per cent in H and Ks, respectively, show significant improvement of the spatial PSF stability compared to single-conjugated adaptive optics systems. Photometry of our observations cover a dynamic range of ~10 mag including 2-3 times more sources than comparable seeing-limited observations. The colour-magnitude diagram reveals that the massive cluster originated in a recent starburst-like event 1+/-0.5 Myr ago. We tentatively detect hints for an older population of 3 Myr suggesting that low intensity star formation might have been going on in the HII region for a few Myr. We derive the luminosity function and mass function between 0.1 M_sun and 3.2 M_sun and identify a change of the power law slope of the mass function at m_c~0.53(+0.12/-0.10) M_sun. The MF appears shallow with power law slopes of Gamma1=-0.50+/-0.11 above m_c and Gamma2=0.63+/-0.32 below m_c.

We present our investigation of pre--main-sequence (PMS) stellar populations in the Large Magellanic Cloud (LMC) from imaging with Hubble Space Telescope WFPC2 camera. Our targets of interest are four star-forming regions located at the periphery of the super-giant shell LMC 4 (Shapley Constellation III). The PMS stellar content of the regions is revealed through the differential Hess diagrams and the observed color-magnitude diagrams (CMDs). Further statistical analysis of stellar distributions along cross-sections of the faint part of the CMDs allowed the quantitative assessment of the PMS stars census, and the isolation of faint PMS stars as the true low-mass stellar members of the regions. These distributions are found to be well represented by a double Gaussian function, the first component of which represents the main-sequence field stars and the second the native PMS stars of each region. Based on this result, a cluster membership probability was assigned to each PMS star according to its CMD position. The higher extinction in the region LH 88 did not allow the unambiguous identification of its native stellar population. The CMD distributions of the PMS stars with the highest membership probability in the regions LH 60, LH 63 and LH 72 exhibit an extraordinary similarity among the regions, suggesting that these stars share common characteristics, as well as common recent star formation history. Considering that the regions are located at different areas of the edge of LMC 4, this finding suggests that star formation along the super-giant shell may have occurred almost simultaneously.

We present a new method for the evaluation of the age and age-spread among pre-main-sequence (PMS) stars in star-forming regions in the Magellanic Clouds, accounting simultaneously for photometric errors, unresolved binarity, differential extinction, stellar variability, accretion and crowding. The application of the method is performed with the statistical construction of synthetic color-magnitude diagrams using PMS evolutionary models. We convert each isochrone into 2D probability distributions of artificial PMS stars in the CMD by applying the aforementioned biases that dislocate these stars from their original CMD positions. A maximum-likelihood technique is then applied to derive the probability for each observed star to have a certain age, as well as the best age for the entire cluster. We apply our method to the photometric catalog of ~2000 PMS stars in the young association LH 95 in the LMC, based on the deepest HST/ACS imaging ever performed toward this galaxy, with a detection limit of V~28, corresponding to M~0.2 Msun. Our treatment shows that the age determination is very sensitive to the considered grid of evolutionary models and the assumed binary fraction. The age of LH 95 is found to vary from 2.8 Myr to 4.4 Myr, depending on these factors. Our analysis allows us to disentangle a real age-spread from the apparent CMD-broadening caused by the physical and observational biases. We find that LH 95 hosts an age-spread well represented by a gaussian distribution with a FWHM of the order of 2.8 Myr to 4.2 Myr depending on the model and binary fraction. We detect a dependence of the average age of the system with stellar mass. This dependence does not appear to have any physical meaning, being rather due to imperfections of the PMS evolutionary models, which tend to predict lower ages for the intermediate masses, and higher ages for low-mass stars.

We present a new analysis of the stellar population of the Orion Nebula Cluster (ONC) based on multi-band optical photometry and spectroscopy. We study the color-color diagrams in BVI, plus a narrow-band filter centered at 6200A, finding evidences that intrinsic color scales valid for main-sequence dwarfs are incompatible with the ONC, while a better agreement is found employing synthetic intrinsic colors obtained constraining the typical lower surface gravity of young stars. We refine these model colors even further, empirically, by comparison with a sample of ONC stars with no accretion and no extinction. We consider the stars with known spectral types from the literature, and add 65 newly classified stars from slit spectroscopy and 182 M-type from narrow-band photometry; in this way we isolate a sample of about 1000 stars with known spectral type. We introduce a new method to self-consistently derive reddening and accretion excess from the location of each star in the BVI color-color diagram. This enables us to accurately determine the extinction of the ONC members. We adopt a lower distance for the ONC than previously assumed, based on recent parallax measurements. With a careful choice also of the spectral type-temperature transformation, we produce the new H-R diagram of the ONC population, more populated than previous works. With respect to previous works, we find higher luminosity for late-type stars and a lower luminosity for early types. We determine the age distribution of the population, peaking at 2-3 Myr, a higher age than previously estimated. We study the distribution of the members in the mass-age plane, and find that taking into account selection effects due to incompleteness removes an apparent correlation between mass and age. We derive the IMF for low- and intermediate-mass members of the ONC, which turns out to be model-dependent, and shows a turn-over at ~<0.2Msun.

We have analyzed two epochs of HST/WFPC2 observations of the young Galactic starburst cluster in NGC 3603 with the aim to study its internal dynamics and stellar population. Relative proper motions measured over 10.15 yrs of more than 800 stars enable us to distinguish cluster members from field stars. The best-fitting isochrone yields Av=4.6-4.7 mag, a distance of 6.6-6.9 kpc, and an age of 1 Myr for NGC 3603 Young Cluster (NYC). We identify pre-main-sequence/main-sequence transition stars located in the short-lived radiative-convective gap, which in the NYC occurs in the mass range 3.5-3.8 Msun. We also identify a sparse population of stars with an age of 4 Myr, which appear to be the lower mass counterparts to previously discovered blue supergiants located in the giant HII region NGC 3603. For the first time, we are able to measure the internal velocity dispersion of a starburst cluster from 234 stars with I < 18.5 mag to {\sigma}_pm1D=141+/-27 {\mu}as/yr (4.5+/-0.8 km/s at a distance of 6.75 kpc). As stars with masses between 1.7 and 9 Msun all exhibit the same velocity dispersion, the cluster stars have not yet reached equipartition of kinetic energy (i.e., the cluster is not in virial equilibrium). The results highlight the power of combining high-precision astrometry and photometry, and emphasize the role of NYC as a benchmark object for testing stellar evolution models and dynamical models for young clusters and as a template for extragalactic starburst clusters.

We present J, H and K photometry of the Orion Nebula Cluster obtained at the CTIO/Blanco 4 m telescope in Cerro Tololo with the ISPI imager. From the observations we have assembled a catalog of about 7800 sources distributed over an area of approximately 30'x40', the largest of any survey deeper than 2MASS in this region. The catalog provides absolute coordinates accurate to about 0.15 arcseconds and 3sigma photometry in the 2MASS system down to J 19.5mag, H 18.0mag, K 18.5mag, enough to detect planetary size objects 1 Myr old under Av 10mag of extinction at the distance of the Orion Nebula. We present a preliminary analysis of the catalog, done comparing the (J-H, H-K) color-color diagram, the (H, J-H) and (K, H-K) color-magnitude diagrams and the JHK luminosity functions of three regions at increasing projected distance from the Trapezium. Sources in the inner region typically show IR colors compatible with reddened T Tauri stars, whereas the outer fields are dominated by field stars seen through an amount of extinction which decreases with the distance from the center. The color-magnitude diagrams make it possible to clearly distinguish between the main ONC population, spread across the full field, and background sources. The luminosity functions of the inner region, corrected for completeness, remain relatively flat in the sub-stellar regime regardless of the strategy adopted to remove background contamination.

We present U, B, V, I broad-band, 6200A TiO medium-band and Halpha photometry of the Orion Nebula Cluster obtained with the WFI imager at the ESO/MPI 2.2 telescope. The nearly-simultaneous observations cover the entire ONC in a field of about 34x34 arcmin. They enable us to determine stellar colors avoiding the additional scatter in the photometry induced by stellar variability typical of pre-main sequence stars. We identify 2,612 point-like sources in I band, 58%, 43% and 17% of them detected also in V, B and U, respectively. 1040 sources are identified in Halpha band. In this paper we present the observations, the calibration techniques, and the resulting catalog. We show the derived CMD of the population and discuss the completeness of our photometry. We define a spectro-photometric TiO index from the fluxes in V, I, and TiO-band. We find a correlation between the index and the spectral type valid for M-type stars, that is accurate to better than 1 spectral sub-class for M3-M6 types and better than 2 spectral subclasses for M0-M2 types. This allows us to newly classify 217 stars. We subtract from our Halpha photometry the photospheric continuum at its wavelength, deriving calibrated line excess for the full sample. This represents the largest Halpha star catalog obtained to date on the ONC. This data set enables a full re-analysis of the properties of the Pre-Main Sequence population in the Orion Nebula Cluster to be presented, in an accompanying paper.

In this photometric study of the stellar association LH 95 in the Large Magellanic Cloud (LMC) we focus on the pre-main Sequence (PMS) population in order to construct, for the first time, the sub-solar initial mass function (IMF) in the LMC. We use the deepest photometry ever performed in the LMC with the Advanced Camera for Surveys (ACS) on-board the Hubble Space Telescope(HST). We carry out a Monte Carlo technique to subtract the contribution of the general field of LMC and we isolate the central region in the observed area of the association. We study the mass function of its field-subtracted population. For this purpose, we introduce a new set of evolutionary models, derived from the calculations on the evolution of PMS stars by Siess et al. We use these models with our observations of LH 95 to derive the IMF of the system, which is reliably constructed down to 0.43 M_sun, the lowest mass ever observed within reasonable completeness in the Magellanic Clouds. Consequently, its construction offers an outstanding improvement in our understanding of the low-mass star formation in the LMC. The system IMF of LH 95 shows a definite change in its slope at 1 M_sun, where it becomes more shallow. In general, the shape of this IMF agrees very well with the average Galactic IMF, down to the sub-solar regime. As far as the slope of this system IMF is concerned, it is found to be somewhat more shallow than the corresponding classical Galactic IMF in the sub-solar regime, probably due to unresolved binarity, while for stars with M>1M_sun it becomes slightly steeper. We do not find significant differences in the shape of the overall IMF of LH 95 from that of each of the three individual sub-clusters, suggesting that the IMF of LH 95 is not subject to local variability.

We have found a photoevaporated disk in the Orion Nebula that includes a wide binary. HST/ACS observations of the proplyd 124-132 show two point-like sources separated by 0".15, or about 60 AU at the distance of Orion. The two sources have nearly identical I and z magnitudes. We analyze the brightest component, Source N, comparing the observed magnitudes with those predicted using a 1 Myr Baraffe/NEXTGEN isochrone with different accretion luminosities and extinctions. We find that a low mass (\simeq 0.04 M_\odot) brown dwarf ~1 Myr old with mass accretion rate \log\dot{M}\simeq -10.3, typical for objects of this mass, and about 2 magnitudes of visual extinction provides the best fit to the data. This is the first observation of a circumbinary disk undergoing photoevaporation and, if confirmed by spectroscopic observations, the first direct detection of a wide substellar pair still accreting and enshrouded in its circumbinary disk.