The large astrometric and photometric survey performed by the Gaia mission allows for a panoptic view of the Galactic disc and in its stellar cluster population. Hundreds of clusters were only discovered after the latest G data release (DR2) and have yet to be characterised. Here we make use of the deep and homogeneous Gaia photometry down to G=18 to estimate the distance, age, and interstellar reddening for about 2000 clusters identified with Gaia~DR2 astrometry. We use these objects to study the structure and evolution of the Galactic disc. We rely on a set of objects with well-determined parameters in the literature to train an artificial neural network to estimate parameters from the Gaia photometry of cluster members and their mean parallax. We obtain reliable parameters for 1867 clusters. Our new homogeneous catalogue confirms the relative lack of old clusters in the inner disc (with a few notable exceptions). We also quantify and discuss the variation of scale height with cluster age, and detect the Galactic warp in the distribution of old clusters. This work results in a large and homogenous cluster catalogue. However, the present sample is still unable to trace the Outer spiral arm of the Milky Way, which indicates that the outer disc cluster census might still be incomplete.

FU Orionis objects (FUors) are rapidly-accreting, pre-main sequence objects that are known to exhibit large outbursts at optical and near-infrared wavelengths, with post-eruption, small-scale photometric variability superimposed on longer-term trends. In contrast, little is known about the variability of FUors at longer wavelengths. To explore this further, we observed six FUor objects using the NOrthern Extended Millimeter Array (NOEMA) and for a subset of three objects we obtained coordinated observations with NOEMA and the Lowell Discovery Telescope (LDT). In combination with previously published NOEMA observations from 2014, our 2017 observations of V1735 Cyg provide the first detection of variability in an FUor object at 2.7 mm. In the absence of significant optical variability, we discount the possibility that the mm flux density changed as a result of irradiation from the central disk. In addition, a change in the dust mass due to infall is highly unlikely. A plausible explanation for the change in 2.7 mm flux density is variability in free-free emission due to changes in the object's jet/wind. Thus, it may be that free-free emission in some FUor objects is significant at $\sim$3 mm and must be considered when deriving disk masses in order to help constrain the mechanism responsible for triggering FUor outbursts.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We present the results of the hierarchical clustering analysis of the Gaia DR2 data to search for clusters, co-moving groups, and other stellar structures. The current paper builds on the sample from the previous work, extending it in distance from 1 kpc to 3 kpc, increasing the number of identified structures up to 8292. To aid in the analysis of the population properties, we developed a neural network called Auriga to robustly estimate the age, extinction, and distance of a stellar group based on the input photometry and parallaxes of the individual members. We apply Auriga to derive the properties of not only the structures found in this paper, but also previously identified open clusters. Through this work, we examine the temporal structure of the spiral arms. Specifically, we find that the Sagittarius arm has moved by >500 pc in the last 100 Myr, and the Perseus arm has been experiencing a relative lull in star formation activity over the last 25 Myr. We confirm the findings from the previous paper on the transient nature of the spiral arms, with the timescale of transition of a few 100 Myr. Finally, we find a peculiar ~1 Gyr old stream of stars that appears to be heliocentric. It is unclear what is the origin of it.

Machine learning allows efficient extraction of physical properties from stellar spectra that have been obtained by large surveys. The viability of ML approaches has been demonstrated for spectra covering a variety of wavelengths and spectral resolutions, but most often for main sequence or evolved stars, where reliable synthetic spectra provide labels and data for training. Spectral models of young stellar objects (YSOs) and low mass main sequence (MS) stars are less well-matched to their empirical counterparts, however, posing barriers to previous approaches to classify spectra of such stars. In this work we generate labels for YSOs and low mass MS stars through their photometry. We then use these labels to train a deep convolutional neural network to predict log g, Teff, and Fe/H for stars with APOGEE spectra in the DR14 dataset. This "APOGEE Net" has produced reliable predictions of log g for YSOs, with uncertainties of within 0.1 dex and a good agreement with the structure indicated by pre-main sequence evolutionary tracks, and correlate well with independently derived stellar radii. These values will be useful for studying pre-main sequence stellar populations to accurately diagnose membership and ages.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

Aims:We take advantage of the second data release of the Gaia space mission and the state-of-the-art astrometry delivered from very long baseline interferometry observations to revisit the structure and kinematics of the nearby Taurus star-forming region. Methods: We apply a hierarchical clustering algorithm for partitioning the stars in our sample into groups (i.e., clusters) that are associated with the various molecular clouds of the complex, and derive the distance and spatial velocity of individual stars and their corresponding molecular clouds. Results: We show that the molecular clouds are located at different distances and confirm the existence of important depth effects in this region reported in previous studies. For example, we find that the L 1495 molecular cloud is located at $d=129.9^{+0.4}_{-0.3}$ pc, while the filamentary structure connected to it (in the plane of the sky) is at $d=160.0^{+1.2}_{-1.2}$ pc. We report B 215 and L 1558 as the closest ($d=128.5^{+1.6}_{-1.6}$ pc) and most remote ($d=198.1^{+2.5}_{-2.5}$ pc) substructures of the complex, respectively. The median inter-cloud distance is 25 pc and the relative motion of the subgroups is on the order of a few km/s. We find no clear evidence for expansion (or contraction) of the Taurus complex, but signs of the potential effects of a global rotation. Finally, we compare the radial velocity of the stars with the velocity of the underlying $^{13}$CO molecular gas and report a mean difference of $0.04\pm0.12$ km/s (with r.m.s. of 0.63 km/s) confirming that the stars and the gas are tightly coupled.

The star forming region of the Orion Nebula (ONC) is ideal to study the stellar dynamics of young stars in a clustered environment. Using Gaia DR2 we search for the pre-main sequence stars with unusually high proper motions that may be representative of a dynamical ejection from unstable young triple systems or other close three-body encounters. We identify twenty-six candidate stars that are likely to have had such an encounter in the last 1 Myr. Nine of these stars could be traced back to the densest central-most region of the ONC, the Trapezium, while five others have likely interactions with other OB-type stars in the cluster. Seven stars originate from other nearby populations within the Orion Complex that coincidentally scattered towards the ONC. A definitive point of origin cannot be identified for the remaining sources. These observations shed light on the frequency of the ejection events in young clusters.

Gaia DR2 provides unprecedented precision in measurements of the distance and kinematics of stars in the solar neighborhood. Through applying unsupervised machine learning on DR2's 5-dimensional dataset (3d position + 2d velocity), we identify a number of clusters, associations, and co-moving groups within 1 kpc and $|b|<30^\circ$ (many of which have not been previously known). We estimate their ages with the precision of $\sim$0.15 dex. Many of these groups appear to be filamentary or string-like, oriented in parallel to the Galactic plane, and some span hundreds of pc in length. Most of these string lack a central cluster, indicating that their filamentary structure is primordial, rather than the result of tidal stripping or dynamical processing. The youngest strings ($<$100 Myr) are orthogonal to the Local Arm. The older ones appear to be remnants of several other arm-like structures that cannot be presently traced by dust and gas. The velocity dispersion measured from the ensemble of groups and strings increase with age, suggesting a timescale for dynamical heating of $\sim$300 Myr. This timescale is also consistent with the age at which the population of strings begins to decline, while the population in more compact groups continues to increase, suggesting that dynamical processes are disrupting the weakly bound string populations, leaving only individual clusters to be identified at the oldest ages. These data shed a new light on the local galactic structure and a large scale cloud collapse.

Recent analyses of Gaia data have provided direct evidence that most young stellar clusters are in a state of expansion, with velocities of the order of ~0.5 km/s. Traditionally, expanding young clusters have been pictured as entities that became unbound due to the lack of gravitational binding once the gas from the parental cloud that formed the cluster has been expelled by the stellar radiation of the massive stars in the cluster. In the present contribution, we used radiation-magnetohydrodynamic numerical simulations of molecular cloud formation and evolution to understand how stellar clusters form and disperse. We found that the ionising feedback from the newborn massive stars expels the gas from the collapse centre, flipping-up the gravitational potential as a consequence of the mass removal from the inside-out. Since neither the parental clouds nor the formed shells are distributed symmetrically around the HII region, net forces pulling out the stars are present, accelerating them towards the edges of the cavity. We call this mechanism ``gravitational feedback", in which the gravity from the expelled gas appears to be the crucial mechanism producing unbound clusters that expand away from their formation centre in an accelerated way in young stellar clusters. This mechanism naturally explains the "Hubble flow-like" expansion observed in several young clusters.

We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E$_{up}$) were detected. From the detected transitions, we estimate that the measured D/H does not exceed 1$\%$, with a measured CH$_2$DOH/CH$_3$OH and CH$_3$OD/CH$_3$OH of about (32$\pm$8)$\times$10$^{-4}$ and (10$\pm$4)$\times$10$^{-4}$, respectively. This finding is consistent with the hypothesis of a short-time scale formation during the pre-stellar phase. We find a relative abundance ratio CH$_2$DOH/CH$_3$OD of 3.2 $\pm$ 1.5. This result is consistent with a statistical deuteration. We cannot exclude H/D exchanges between water and methanol if water deuteration is of the order 0.1$\%$, as suggested by recent Herschel observations.

Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at {\lambda}=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability.

Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of $\sim$0.05--1.5 \msun, we find 399 binaries, including both RV variables and SB2s. The mass ratio distribution of SB2s is consistent with a uniform for $q<0.95$ with an excess of twins with $q>0.95$. The period distribution is consistent with what has been observed in close binaries ($<10$ AU) in the evolved populations. Three systems are found to have $q\sim$0.1, with a companion located within the brown dwarf desert. There are not any strong trends in the multiplicity fraction (MF) as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at $\Sigma_*\sim30$ stars/pc$^{-2}$, and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of $\sim$2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal mass systems that is less effective in binaries with lower $q$.

Background: low-mass stars are the dominant product of the star formation process, and they trace star formation over the full range of environments, from isolated globules to clusters in the central molecular zone. In the past two decades, our understanding of the spatial distribution and properties of young low-mass stars and protostars has been revolutionized by sensitive space-based observations at X-ray and IR wavelengths. By surveying spatial scales from clusters to molecular clouds, these data provide robust measurements of key star formation properties. Goal: with their large numbers and their presence in diverse environments, censuses of low mass stars and protostars can be used to measure the dependence of star formation on environmental properties, such as the density and temperature of the natal gas, strengths of the magnetic and radiation fields, and the density of stars. Here we summarize how such censuses can answer three basic questions: i.) how is the star formation rate influenced by environment, ii.) does the IMF vary with environment, and iii.) how does the environment shape the formation of bound clusters? Answering these questions is an important step toward understanding star and cluster formation across the extreme range of environments found in the Universe. Requirements: sensitivity and angular resolution improvements will allow us to study the full range of environments found in the Milky Way. High spatial dynamic range (< 1arcsec to > 1degree scales) imaging with space-based telescopes at X-ray, mid-IR, and far-IR and ground-based facilities at near-IR and sub-mm wavelengths are needed to identify and characterize young stars.

We present Gaia-DR2 astrometry of a sample of YSO candidates in Ophiuchus, Serpens Main and Serpens South/W40 in the Aquila Rift, which had been mainly identified by their infrared excess with Spitzer. We compare the Gaia-DR2 parallaxes against published and new parallaxes obtained from our Very Long Baseline Array (VLBA) program GOBELINS. We obtain consistent results between Gaia and the VLBA for the mean parallaxes in each of the regions analyzed here. We see small offsets, when comparing mean values, of a few tens of micro-arcseconds in the parallaxes, which are either introduced by the Gaia zero-point error or due to a selection effect by Gaia toward the brightest, less obscured stars. Gaia-DR2 data alone conclusively places Serpens Main and Serpens South at the same distance, as we first inferred from VLBA data alone in a previous publication. Thus, Serpens Main, Serpens South and W40 are all part of the same complex of molecular clouds, located at a mean distance of 436+/-9 pc. In Ophiuchus, both Gaia and VLBA suggest a small parallax gradient across the cloud, and the distance changes from 144.2+/-1.3 pc to 138.4+/-2.6 pc when going from L1689 to L1688.

We derive the distance and structure of the Perseus molecular cloud by combining trigonometric parallaxes from Very Long Baseline Array (VLBA) observations, taken as part of the GOBELINS survey, and Gaia Data Release 2. Based on our VLBA astrometry, we obtain a distance of 321+/-10 pc for IC 348. This is fully consistent with the mean distance of 320+/-26 measured by Gaia. The VLBA observations toward NGC 1333 are insufficient to claim a successful distance measurement to this cluster. Gaia parallaxes, on the other hand, yield a mean distance of 293+/-22 pc. Hence, the distance along the line of sight between the eastern and western edges of the cloud is ~30 pc, which is significantly smaller than previously inferred. We use Gaia proper motions and published radial velocities to derive the spatial velocities of a selected sample of stars. The average velocity vectors with respect to the LSR are (u,v,w) = (-6.1+/-1.6, 6.8+/-1.1, -0.9+/-1.2) and (-6.4+/-1.0, 2.1+/-1.4, -2.4+/-1.0) km/s for IC 348 and NGC 1333, respectively. Finally, our analysis of the kinematics of the stars has shown that there is no clear evidence of expansion, contraction, or rotational motions within the clusters.

Binary stars make up a significant portion of all stellar systems. Consequently, an understanding of the bulk properties of binary stars is necessary for a full picture of star formation. Binary surveys indicate that both multiplicity fraction and typical orbital separation increase as functions of primary mass. Correlations with higher order architectural parameters such as mass ratio are less well constrained. We seek to identify and characterize double-lined spectroscopic binaries (SB2s) among the 1350 M dwarf ancillary science targets with APOGEE spectra in the SDSS-III Data Release 13. We measure the degree of asymmetry in the APOGEE pipeline cross-correlation functions (CCFs), and use those metrics to identify a sample of 44 high-likelihood candidate SB2s. At least 11 of these SB2s are known, having been previously identified by Deshapnde et al, and/or El Badry et al. We are able to extract radial velocities (RVs) for the components of 36 of these systems from their CCFs. With these RVs, we measure mass ratios for 29 SB2s and 5 SB3s. We use Bayesian techniques to fit maximum likelihood (but still preliminary) orbits for 4 SB2s with 8 or more distinct APOGEE observations. The observed (but incomplete) mass ratio distribution of this sample rises quickly towards unity. Two-sided Kolmogorov-Smirnov tests and probabilities of 18.3% and 18.7%, demonstrating that the mass ratio distribution of our sample is consistent with those measured by Pourbaix et al. and Fernandez et al., respectively.

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.

We present new trigonometric parallaxes and proper motions of young stellar objects in the Taurus molecular cloud complex from observations collected with the Very Long Baseline Array as part of the Gould's Belt Distances Survey (GOBELINS). We detected 26 young stellar objects and derived trigonometric parallaxes for 18 stars with an accuracy of 0.3$\%$ to a few percent. We modeled the orbits of six binaries and determined the dynamical masses of the individual components in four of these systems (V1023 Tau, T Tau S, V807 Tau and V1000 Tau). Our results are consistent with the first trigonometric parallaxes delivered by the Gaia satellite and reveal the existence of significant depth effects. We find that the central portion of the dark cloud Lynds 1495 is located at $d=129.5\pm 0.3$ pc while the B 216 clump in the filamentary structure connected to it is at $d=158.1\pm 1.2$ pc. The closest and remotest stars in our sample are located at $d=126.6\pm 1.7$ pc and $d=162.7\pm 0.8$ pc yielding a distance difference of about 36 pc. We also provide a new distance estimate for HL Tau that was recently imaged. Finally, we compute the spatial velocity of the stars with published radial velocity and investigate the kinematic properties of the various clouds and gas structures in this region.

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 present observations of VLBA 20, a radio source found towards the edge of the Orion Nebula Cluster (ONC). Nonthermal emission dominates the spectral energy distribution of this object from radio to mid-infrared regime, suggesting that VLBA 20 is extragalactic. This source is heavily scattered in the radio regime. Very Long Baseline Array observations resolve it to ~34x19 mas at 5 GHz, and the wavelength dependence of the scattering disk is consistent with $\nu^{-2}$ at other frequencies. The origin of the scattering is most likely the ionized X-ray emitting gas from the winds of the most massive stars of the ONC. The scattering is highly anisotropic, with the axis ratio of 2:1, higher than what is typically observed towards other sources.

We present results from high resolution optical spectra toward 66 young stars in the Orion B molecular cloud to study their kinematics and other properties. Observations of the H$\alpha$ and Li I 6707 {\AA} lines are used to check membership and accretion properties. While the stellar radial velocities of in NGC 2068 and L1622 show good agreement with that of the molecular gas, many of the stars in NGC 2024 show a considerable offset. This could be a signature of either expansion of the cluster, high degree of the ejection of the stars from the cluster through the dynamical interaction, or the acceleration of the gas due to stellar feedback.

NGC 1980 is a young cluster that is located about 0.5 degrees south of the Orion Nebula Cluster (ONC). Recent studies by Bouy et al. and Pillitteri et al. have suggested that NGC 1980 contains an older population of stars compared to a much younger ONC, and that it belongs to a foreground population that may be located in front of the Orion A molecular gas by as much as 40 pc. In this work we present low-resolution spectra towards 148 young stars found towards the NGC 1980 region. We determine the spectral types of these stars, examine accretion signatures and measure the extinction towards them. We determine that based on these observations, the age of the population of NGC 1980 is indistinguishable from L1641, estimated to be ~3 Myr, comparable with the study by Fang et al.

We present the first results of the Gould's Belt Distances Survey (GOBELINS), a project aimed at measuring the proper motion and trigonometric parallax of a large sample of young stars in nearby regions using multi-epoch Very Long Baseline Array (VLBA) radio observations. Enough VLBA detections have now been obtained for 16 stellar systems in Ophiuchus to derive their parallax and proper motion. This leads to distance determinations for individual stars with an accuracy of 0.3 to a few percent. In addition, the orbits of 6 multiple systems were modelled by combining absolute positions with VLBA (and in some cases, near infrared) angular separations. Twelve stellar systems are located in the dark cloud Lynds 1688, the individual distances for this sample are highly consistent with one another, and yield a mean parallax for Lynds 1688 of $\varpi=7.28\pm0.06$ mas, corresponding to a distance $d=137.3\pm1.2$ pc. This represents an accuracy better than 1%. Three systems for which astrometric elements could be measured are located in the eastern streamer (Lynds 1689) and yield an estimate of $\varpi=6.79\pm0.16$ mas, corresponding to a distance $d=147.3\pm3.4$ pc. This suggests that the eastern streamer is located about 10 pc farther than the core, but this conclusion needs to be confirmed by observations (currently being collected) of additional sources in the eastern streamer. From the measured proper motions, we estimate the one-dimensional velocity dispersion in Lynds 1688 to be 2.8$\pm$1.8 and 3.0$\pm$2.0 ${\rm km~s}^{-1}$, in R.A. and DEC., respectively, these are larger than, but still consistent within $1\sigma$, with those found in other studies.

Using multi-epoch VLA observations, covering a time baseline of 29.1 years, we have measured the proper motions of 88 young stars with compact radio emission in the core of the Orion Nebula Cluster (ONC) and the neighboring BN/KL region. Our work increases the number of young stars with measured proper motion at radio frequencies by a factor of 2.5 and enables us to perform a better statistical analysis of the kinematics of the region than was previously possible. Most stars (79 out of 88) have proper motions consistent with a Gaussian distribution centered on $\overline{\mu_{\alpha}\cos{\delta}}=1.07\pm0.09\quad{\rm mas\,yr^{-1}}$, and $\overline{\mu_{\delta}}=-0.84\pm0.16\quad{\rm mas\,yr^{-1}}$, with velocity dispersions of $\sigma_{\alpha}=1.08\pm0.07\quad{\rm mas\,\,yr^{-1}},$ $\sigma_{\delta}=1.27\pm0.15\quad{\rm mas\,\,yr^{-1}}$. We looked for organized movements of these stars but found no clear indication of radial expansion/contraction or rotation. The remaining nine stars in our sample show peculiar proper motions that differ from the mean proper motions of the ONC by more than 3-$\sigma$. One of these stars, V 1326 Ori, could have been expelled from the Orion Trapezium 7,000 years ago. Two could be related to the multi-stellar disintegration in the BN/KL region, in addition to the previously known sources BN, I and n. The others either have high uncertainties (so their anomalous proper motions are not firmly established) or could be foreground objects.

We report on new distances and proper motions to seven stars across the Serpens/Aquila complex. The observations were obtained as part of the Gould's Belt Distances Survey (GOBELINS) project between September 2013 and April 2016 with the Very Long Baseline Array (VLBA). One of our targets is the proto-Herbig AeBe object EC 95, which is a binary system embedded in the Serpens Core. For this system, we combined the GOBELINS observations with previous VLBA data to cover a total period of ~8 years, and derive the orbital elements and an updated source distance. The individual distances to sources in the complex are fully consistent with each other, and the mean value corresponds to a distance of $436.0\pm9.2$~pc for the Serpens/W40 complex. Given this new evidence, we argue that Serpens Main, W40 and Serpens South are physically associated and form a single cloud structure.

We present the results of the Gould's Belt Distances Survey (GOBELINS) of young star forming regions towards the Orion Molecular Cloud Complex. We detected 36 YSOs with the Very Large Baseline Array (VLBA), 27 of which have been observed in at least 3 epochs over the course of 2 years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions towards these stars to study the structure and kinematics of the Complex. We measured a distance of 388$\pm$5 pc towards the Orion Nebula Cluster, 428$\pm$10 pc towards the southern portion L1641, 388$\pm$10 pc towards NGC 2068, and roughly $\sim$420 pc towards NGC 2024. Finally, we observed a strong degree of plasma radio scattering towards $\lambda$ Ori.

We present a near-IR survey for the visual multiples in the Orion molecular clouds region at separations between 100 and 1000 AU. These data were acquired at 1.6~$\mu$m with the NICMOS and WFC3 cameras on the Hubble Space Telescope. Additional photometry was obtained for some of the sources at 2.05~$\mu$m with NICMOS and in the $L'$-band with NSFCAM2 on the IRTF. Towards 129 protostars and 197 pre-main sequence stars with disks observed with WFC3, we detect 21 and 28 candidate companions between the projected separations of 100---1000 AU, of which less than 5 and 8, respectively, are chance line of sight coincidences. The resulting companion fraction ($CF$) after the correction for the line of sight contamination is 14.4$^{+1.1}_{-1.3}$% for protostars and 12.5$^{+1.2}_{-0.8}$% for the pre-main sequence stars. These values are similar to those found for main sequence stars, suggesting that there is little variation in the $CF$ with evolution, although several observational biases may mask a decrease in the $CF$ from protostars to the main sequence stars. After segregating the sample into two populations based on the surrounding surface density of YSOs, we find that the $CF$ in the high stellar density regions ($\Sigma_{YSO} > 45$~pc$^{-2}$) is approximately 50% higher than that found in the low stellar density regions ($\Sigma_{YSO} < 45$~pc$^{-2}$). We interpret this as evidence for the elevated formation of companions at 100 to 1000 AU in the denser environments of Orion. We discuss possible reasons for this elevated formation.

We examine the spectroscopic binary population for two massive nearby regions of clustered star formation, the Orion Nebula Cluster and NGC 2264, supplementing the data presented by Tobin et al. (2009, 2015) with more recent observations and more extensive analysis. The inferred multiplicity fraction up to 10 AU based on these observations is $5.3\pm 1.2$% for NGC 2264 and $5.8\pm 1.1$% for the ONC; they are consistent with the distribution of binaries in the field in the relevant parameter range. Eight of the multiple systems in the sample have enough epochs to make an initial fit for the orbital parameters. Two of these sources are double-lined spectroscopic binaries; for them we determine the mass ratio. Our reanalysis of the distribution of stellar radial velocities towards these clusters presents a significantly better agreement between stellar and gas kinematics than was previously thought.

We present multi-epoch, large-scale ($\sim$ 2000 arcmin${}^2$), fairly deep ($\sim$ 16 $\mu$Jy), high-resolution ($\sim$ 1") radio observations of the Perseus star-forming complex obtained with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. These observations were mainly focused on the clouds NGC 1333 and IC 348, although we also observed several fields in other parts of the Perseus complex. We detect a total of 206 sources, 42 of which are associated with young stellar objects (YSOs). The radio properties of about 60% of the YSOs are compatible with a non-thermal radio emission origin. Based on our sample, we find a fairly clear relation between the prevalence of non-thermal radio emission and evolutionary status of the YSOs. By comparing our results with previously reported X-ray observations, we show that YSOs in Perseus follow a G\"udel-Benz relation with $\kappa$ = 0.03 consistent with other regions of star formation. We argue that most of the sources detected in our observations but not associated with known YSOs are extragalactic, but provide a list of 20 unidentified radio sources whose radio properties are consistent with being YSO candidates. Finally we also detect 5 sources with extended emission features which can clearly be associated with radio galaxies.

We present deep ($\sim 17~\mu$Jy) radio continuum observations of the Serpens molecular cloud, the Serpens south cluster, and the W40 region obtained using the Very Large Array in its A configuration. We detect a total of 146 sources, 29 of which are young stellar objects (YSOs), 2 are BV stars and 5 more are associated with phenomena related to YSOs. Based on their radio variability and spectral index, we propose that about 16 of the remaining 110 unclassified sources are also YSOs. For approximately 65% of the known YSOs detected here as radio sources, the emission is most likely non-thermal, and related to stellar coronal activity. As also recently observed in Ophiuchus, our sample of YSOs with X-ray counterparts lies below the fiducial G\"udel & Benz relation. Finally, we analyze the proper motions of 9 sources in the W40 region. This allows us to better constrain the membership of the radio sources in the region.

We present a multi-epoch radio study of the Taurus-Auriga star-forming complex made with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. We detect a total of 610 sources, 59 of which are related to young stellar objects and 18 to field stars. The properties of 56\% of the young stars are compatible with non-thermal radio emission. We also show that the radio emission of more evolved young stellar objects tends to be more non-thermal in origin and, in general, that their radio properties are compatible with those found in other star forming regions. By comparing our results with previously reported X-ray observations, we notice that young stellar objects in Taurus-Auriga follow a G\"{u}del-Benz relation with $\kappa$=0.03, as we previously suggested for other regions of star formation. In general, young stellar objects in Taurus-Auriga and in all the previous studied regions seem to follow this relation with a dispersion of $\sim1$ dex. Finally, we propose that most of the remaining sources are related with extragalactic objects but provide a list of 46 unidentified radio sources whose radio properties are compatible with a YSO nature.

We present results from a high-sensitivity (60 $\mu$Jy), large-scale (2.26 square degree) survey obtained with the Karl G. Jansky Very Large Array as part of the Gould's Belt Survey program. We detected 374 and 354 sources at 4.5 and 7.5 GHz, respectively. Of these, 148 are associated with previously known Young Stellar Objects (YSOs). Another 86 sources previously unclassified at either optical or infrared wavelengths exhibit radio properties that are consistent with those of young stars. The overall properties of our sources at radio wavelengths such as their variability and radio to X-ray luminosity relation are consistent with previous results from the Gould's Belt Survey. Our detections provide target lists for followup VLBA radio observations to determine their distances as YSOs are located in regions of high nebulosity and extinction, making it difficult to measure optical parallaxes.

Edge-on protostars are valuable for understanding the disk and envelope properties of embedded young stellar objects, since the disk, envelope, and envelope cavities are all distinctly visible in resolved images and well constrained in modeling. Comparing 2MASS, WISE, Spitzer, Herschel, APEX, and IRAM photometry and limits from 1.25 to 1200 microns, Spitzer spectroscopy from 5 to 40 microns, and high-resolution Hubble imaging at 1.60 and 2.05 microns to radiative transfer modeling, we determine envelope and disk properties for the Class I protostar HOPS 136, an edge-on source in Orion's Lynds 1641 region. The source has a bolometric luminosity of 0.8 L_sun, a bolometric temperature of 170 K, and a ratio of submillimeter to bolometric luminosity of 0.8%. Via modeling, we find a total luminosity of 4.7 L_sun (larger than the observed luminosity due to extinction by the disk), an envelope mass of 0.06 M_sun, and a disk radius and mass of 450 AU and 0.002 M_sun. The stellar mass is highly uncertain but is estimated to fall between 0.4 and 0.5 M_sun. To reproduce the flux and wavelength of the near-infrared scattered-light peak in the spectral energy distribution, we require 5.4 x 10^-5 M_sun of gas and dust in each cavity. The disk has a large radius and a mass typical of more evolved T Tauri stars in spite of the significant remaining envelope. HOPS 136 appears to be a key link between the protostellar and optically revealed stages of star formation.

Individual outbursting young stars are important laboratories for studying the physics of episodic accretion and the extent to which this phenomenon can explain the luminosity distribution of protostars. We present new and archival data for V2775 Ori (HOPS 223), a protostar in the L 1641 region of the Orion molecular clouds that was discovered by Caratti o Garatti et al. (2011) to have recently undergone an order-of-magnitude increase in luminosity. Our near-infrared spectra of the source have strong blueshifted He I 10830 absorption, strong H2O and CO absorption, and no H I emission, all typical of FU Orionis sources. With data from IRTF, 2MASS, HST, Spitzer, WISE, Herschel, and APEX that span from 1 to 70 microns pre-outburst and from 1 to 870 microns post-outburst, we estimate that the outburst began between 2005 April and 2007 March. We also model the pre- and post-outburst spectral energy distributions of the source, finding it to be in the late stages of accreting its envelope with a disk-to-star accretion rate that increased from about 2x10^-6 M_sun/yr to about 10^-5 M_sun/yr during the outburst. The post-outburst luminosity at the epoch of the FU Orionis-like near-IR spectra is 28 L_sun, making V2775 Ori the least luminous documented FU Orionis outburster with a protostellar envelope. The existence of low-luminosity outbursts supports the notion that a range of episiodic accretion phenomena can partially explain the observed spread in protostellar luminosities.