We present the results of a J, H, and K photometric variability survey of the central 0.78 square degrees of the young OB association Cygnus OB2. We used data observed with the Wide-Field CAMera at the United Kingdom Infrared Telescope in 2007 (spanning 217 days) to investigate the light curves of 5083 low mass candidate members in the association and explore the occurrence and main characteristics of their near-infrared variability. We identified 2529 stars ($\sim$50$\%$ of the sample) with significant variability with time-scales ranging from days to months. We classified the variable stars into the following three groups according to their light curve morphology: periodic variability (1697 stars), occultation variability (124 stars), and other types of variability (726 stars). We verified that the disk-bearing stars in our sample are significantly more variable in the near-infrared than diskless stars, with a steep increase in the disk-fraction among stars with higher variability amplitude. We investigated the trajectories described by variable stars in the color-space and measured slopes for 335 stars describing linear trajectories. Based on the trajectories in the color-space, we inferred that the sample analyzed is composed of a mix of young stars presenting variability due to hot and cold spots, extinction by circumstellar material, and changes in the disk emission in the near-infrared. We contemplated using the use of near-infrared variability to identify disk-bearing stars and verified that 53.4$\%$ of the known disk-bearing stars in our sample could have been identified as such based solely on their variability. We present 18 newly identified disk-bearing stars and 14 eclipsing binary candidates among CygOB2 lower-mass members.

We present high-resolution $H$- and $K$-band spectroscopic observations of the archetypal T Tauri star T Tau North. Synthetic spectral modeling is used to derive the $K$-band temperature, surface gravity, magnetic field strength, and rotational velocity for this star. The $K$-band spectroscopic temperature measured is $\rm T_{K\text{-band}} = 3976 \pm 90 K$, which is $\sim 1000 \, \rm K$ cooler than the temperature measured from optical observations. Our $K$-band temperature measurement for T Tau N is confirmed using equivalent width line ratio vs. temperature relations in the $H$-band, from which a $\rm T_{H\text{-band}} = 4085 \pm 155 K$ is derived. This optical vs. IR temperature difference is interpreted as cool or hot spots, or both, covering a significant part of the surface of T Tau N. The gravity derived for T Tau N $\log{g}=3.45 \pm 0.14$ is lower than the gravity of nearly every other star in a sample of 24 classical T Tauri stars in Taurus. Combining these temperature and gravity results with magnetic stellar evolutionary models, we find the age of T Tau N to be less than 1 Myr old. These results suggest that T Tau N is in an earlier evolutionary stage than most classical T Tauri stars in Taurus, arguing that it is a protostar ejected from the embedded southern binary system shortly after its formation.

The bipolar jet HH 212, among the finest collimated jets known, has so far been detected only in near-infrared H$_2$ emission. Here we present deep optical images that show two of the major bow shocks weakly detected in optical [SII] emission, as expected for a bona fide Herbig-Haro jet. We present widefield H$_2$ images which reveal two more bow shocks located symmetrically around the source and along the main jet axis. Additionally, examination of Spitzer 4.5 $\mu$m images reveals yet another bright bow shock further to the north along the jet axis; no corresponding bow shock is seen to the south. In total, the HH 212 flow has an extent of 1050 arcsec, corresponding to a projected dimension of 2.0 pc. HH 212 thus joins the growing group of parsec-scale Herbig-Haro jets. Proper motion measurements indicate a velocity of about 170 km/sec, highly symmetric around the source, with an uncertainty of $\sim$30 km/sec, suggesting a probable age of the giant HH 212 flow of about 7000 yr. The jet is driven by a deeply embedded source, known as IRAS 05413-0104. We draw attention to a Spitzer near- and mid-infrared source, which we call IRS-B, located only 7 arcsec from the driving source, towards the outskirts of the dense cloud core. Infrared photometry and spectroscopy suggests that IRS-B is a K-type star with a substantial infrared excess, except that for an extinction of A$_V$ = 44 the star would have only a weakinfrared excess, and so in principle it could be a K-giant at a distance of about 2 kpc.

We report on the source Gaia 19ajj, identifying it as a young star associated with a little-studied star-forming region seen along a complex line-of-sight through the Gum Nebula. The optical lightcurve recently recorded by Gaia exhibits a slow and unsteady 5.5 mag rise over about 3 years, while the mid-infrared lightcurve from NEOWISE over the same time period shows a 1.2 mag rise having similar structure. Available color information is inconsistent with pure extinction reduction as the cause for the photometric brightening. Optical spectroscopic characteristics in the current bright phase include: little in the way of absorption except for the hallmark Li I 6707 A signature of youth plus weak e.g. Ca I and notably Ba II; strong wind/outflow in Ca II, Mg I b, Na I D, Halpha, K I, O I; jet signatures in [O I], [S II], [Ca II], [Fe II], and [Ni II]; and narrow rest-velocity emission in neutral species such as Fe I, Ni I, and Mg I. The infrared spectrum is also characterized by outflow and emission, including: a hot He I wind, jet lines such as [Fe II] and H2; and weak narrow rest-velocity atomic line emission. The 12CO bandheads are weakly in emission, but there is also broad H2O absorption. Gaia 19ajj exhibited a previous bright state in the 2010-2012 time frame. The body of photometric and spectroscopic evidence suggests that the source bears resemblance to V2492 Cyg (PTF 10nvg) and PV Cep, both of which similarly experience bright phases that recur on long timescales, with large-amplitude photometric variations and emission-dominated spectra. We interpret the behavior of Gaia 19ajj as caused by cycles of enhanced disk accretion accompanied by reduced extinction.

We report photometry and spectroscopy of the outburst of the young stellar object ESO-Halpha 99. The outburst was first noticed in Gaia alert Gaia18dvc and later by ATLAS. We have established the outburst light curve with archival ATLAS ``Orange'' filter photometry, Gaia data, new V-band photometry, and J, H, and K_s photometry from IRIS and UKIRT. The brightness has fluctuated several times near the light curve maximum. The TESS satellite observed ESO-Halpha 99 with high cadence during one of these minor minima and found brightness fluctuations on timescales of days and hours. Imaging with UKIRT shows the outline of an outflow cavity, and we find one knot of H_2 1-0 S(1) emission, now named MHO 1520, on the symmetry axis of this nebula, indicating recent collimated outflow activity from ESO-Halpha 99. Its pre-outburst SED shows a flat FIR spectrum, confirming its early evolutionary state and its similarity to other deeply embedded objects in the broader EXor class. The pre-outburst luminosity is 34 L_sun, a much higher luminosity than typical EXors, indicating that ESO-Halpha 99 may be a star of intermediate mass. Infrared and optical spectroscopy show a rich emission line spectrum, including HI lines, strong red CaII emission, as well as infrared CO bandhead emission, all characteristic EXors in the broadest sense. Comparison of the present spectra with an optical spectrum obtained in 1993, presumably in the quiescent state of the object, shows that during the present outburst the continuum component of the spectrum has increased notably more than the emission lines. The Halpha equivalent width during the outburst is down to one half of its 1993 level, and shock-excited emission lines are much less prominent.

A deep objective-prism survey for H-alpha emission stars towards the Canis Major star-forming clouds was performed. A total of 398 Halpha emitters were detected, 353 of which are new detections. There is a strong concentration of these H-alpha emitters towards the molecular clouds surrounding the CMa~OB1 association, and it is likely that these stars are young stellar objects recently born in the clouds. An additional population of H-alpha emitters is scattered all across the region, and probably includes unrelated foreground dMe stars and background Be stars. About 90% of the H-alpha emitters are detected by WISE, of which 75% was detected with usable photometry. When plotted in a WISE colour-colour diagram it appears that the majority are Class II YSOs. Coordinates and finding charts are provided for all the new stars, and coordinates for all the detections. We searched the Gaia-DR2 catalogue and from 334 Halpha emission stars with useful parallaxes, we selected a subset of 98 stars that have parallax errors of less than 20% and nominal distances in the interval 1050 to 1350 pc that surrounds a strong peak at 1185 pc in the distance distribution. Similarly, Gaia distances were obtained for 51 OB-stars located towards Canis Major and selected with the same parallax errors as the H-alpha stars. We find a median distance for the OB stars of 1182 pc, in excellent correspondence with the distance from the H-alpha stars. Two known runaway stars are confirmed as members of the association. Finally, two new Herbig-Haro objects are identified.

While it has been suggested that there is a connection between the magnetic properties and the internal structure of young stars, there have not been enough magnetic measurements to firmly establish such a correlation at the earliest ages. Here, we contribute to this endeavor by presenting stellar parameters and magnetic field strength measurements of BP Tau and V347 Aur, both stars observed with the near-infrared spectrograph iSHELL. We first test the accuracy of our method by fitting synthetic stellar spectra to a sample of nine main and post-main-sequence stars. We report uncertainties of $\sigma_{\rm Teff}$ = 91 K in temperature and $\sigma_{\rm log(g)}$ = 0.14 in gravity. We then apply the modeling technique to BP Tau and measure a surface magnetic field strength of $\langle \rm B \rangle $ = 2.5$^{+0.15}_{-0.16}$ kG, confirming literature results. For this star, however, we obtain a much lower temperature value than previous optical studies ($\Delta \rm T \sim 400$ K) and interpret this significant temperature difference as due to the relatively higher impact of starspots at near-infrared wavelengths than at optical wavelengths. We further apply this technique to the class I protostellar source V347 Aur and measure for the first time its magnetic field strength $\langle \rm B \rangle = $ 1.36$^{+0.06}_{-0.05}$ kG and its surface gravity log(g) = 3.25$^{+0.14}_{-0.14}$. Lastly, we combine our measurements with pre-main-sequence stellar evolutionary models and illustrate the effects produced by starspots on the retrieved masses and ages of young stars.

We investigate the silicate feature of the two Class I components of HH250-IRS, a resolved binary system with a separation of $0''53$ driving a Herbig-Haro flow. Each component has its own circumstellar envelope, and the system is surrounded by a circumbinary disk. We have carried out low resolution spectroscopy in the 8-13$\mu$m range using VISIR at ESO's Very Large Telescope. The silicate features of both sources are clearly different. The NW component has a broad, smooth absorption profile lacking structure. The SE component shows the silicate feature in emission, with structure longwards of 9.5$\mu$m indicating the presence of crystalline dust in the dominant form of forsterite. The apparent lack of an absorption feature caused by foreground dust is probably due to the filling of the band with emission by amorphous silicates in the envelope of the object. Despite their virtually certain coevality, the differences in the components of the HH250-IRS binary are most likely due to markedly different circumstellar environments. The NW component displays an unevolved envelope, whereas dust growth and crystallization has taken place in the SE component. The weak or absent signatures of enstatite in the latter are fairly unusual among envelopes with crystalline dust, and we tentatively relate it to a possible wide gap or an inner truncation of the disk already hinted in previous observations by a drop in the $L'$-band flux, which might indicate that the SE component could actually be a very close binary. We speculate that the clear differences between the silicate feature spectra of both components of HH250-IRS may be due either to disk evolution sped up by multiplicity, or by accretion variability leading to episodes of crystal formation.

We present new results on the young star population of the Lynds~1188 molecular cloud, associated with the Cepheus Bubble, a giant interstellar shell around the association Cep~OB\,2. In order to reveal the star-forming scenario of the molecular cloud located on the supershell, and understand the history of star formation in the region, we identified and characterized young star candidates based on an H$\alpha$ emission survey and various published photometric datasets. Using Gaia DR2 astrometry we studied the spatial distribution of the young star candidates and isolated three groups based on their distances. We constructed spectral energy distributions of our target stars, based on Pan-STARRS, 2MASS, {\it Spitzer} and {\it WISE} photometric data, estimating their spectral types, extinctions, and luminosities. We estimated masses by means of pre-main-sequence evolutionary models, and derived accretion rates from the equivalent width of the H$\alpha$ line. We studied the structure of the cloud by constructing a new extinction map, based on Pan-STARRS data. Our results show that the distribution of low-mass young stars in L1188 is well correlated with that of the dust and mole\-cular gas. We identified two small, compact clusters and a loose aggregate of young stars. We found that star formation in L1188 started about 5 million years ago. The apparent age gradient of young stars across the cloud and the ammonia cores located to the east of the optically visible young stellar groups support the scenario of star formation propagating away from the centre of the Cepheus Bubble.

We have conducted a homogenous near-IR spectroscopic survey of 33 objects with varying degrees of similarity to FU~Orionis. Common spectroscopic features that are characteristic of the three classical FUors FU Ori, V1057 Cyg, and V1515 Cyg are: strong CO absorption, weak metal absorption, strong water bands, low gravity, strong blue shifted He I absorption, and few (if any) emission lines. Based on these criteria, we classify the 33 objects as either bona fide FUors (eruption observed), FUor-like objects (eruption not observed), or peculiar objects with some FUor-like characteristics, and present a spectral atlas of 14 bona-fide FUors, 10 FUor-like objects, and 9 peculiar objects. All objects that we classify as FUors or FUor-like have very similar near-IR spectra. We use this spectral similarity to determine the extinction to each source, and correlate the extinction to the depth of the 3 micron ice band. All bona fide FUors still today maintain the spectrum of a FUor, despite the eruption occurring up to 80 years ago. Most FUors and FUor-like objects occupy a unique space on a plot of Na+Ca vs. CO equivalent widths, whereas the peculiar objects tend to be found mostly elsewhere. Since most FUors show a reflection nebula, we also present an atlas of K-band images of each target. We found that the near-IR spectra of FUors and young brown dwarfs can be extremely similar, a distinguishing feature being the Paschen beta absorption in the spectra of FUors. Although V1647 Ori, AR 6a, and V346 Normae had been previously classified as candidate FUors, we classify them as peculiar objects with some FUor-like properties since their spectra now differ significantly from bona fide FUors. We confirm two new FUor-like objects that were initially identified as candidates based on their near-IR morphology.

Herbig-Haro flows are signposts of recent major accretion and outflow episodes. We aim to determine the nature and properties of the little-known outflow source HH 250-IRS, which is embedded in the Aquila clouds. We have obtained adaptive optics-assisted L-band images with the NACO instrument on the Very Large Telescope (VLT), together with N- and Q-band imaging with VISIR also on the VLT. Using the SINFONI instrument on the VLT we carried out H- and K-band integral field spectroscopy of HH 250-IRS, complemented with spectra obtained with the SpeX instrument at the InfraRed Telescope Facility (IRTF) in the JHKL bands. Finally, the SubMillimeter Array (SMA) interferometer was used to study the circumstellar environment of HH 250-IRS at 225 and 351 GHz with CO (2-1) and CO (3-2) maps and 0.9 mm and 1.3 mm continuum images. The HH 250-IRS source is resolved into a binary with 0''53 separation, corresponding to 120 AU at the adopted distance of 225 pc. The individual components show heavily veiled spectra with weak CO absorption indicative of late-type stars. Both are Class I sources, but their spectral energy distributions between 1.5 $\mu$m and 19 $\mu$m differ markedly and suggest the existence of a large cavity around one of the components. The millimeter interferometric observations indicate that the gas mainly traces a circumbinary envelope or disk, while the dust emission is dominated by one of the circumstellar envelopes. HH 250-IRS is a new addition to the handful of multiple systems where the individual stellar components, the circumstellar disks and a circumbinary disk can be studied in detail, and a rare case among those systems in which a Herbig-Haro flow is present.

We describe a new method for determining proper motions of extended objects, and a pipeline developed for the application of this method. We then apply this method to an analysis of four epochs of [S~II] HST images of the HH~1 jet (covering a period of $\sim 20$~yr). We determine the proper motions of the knots along the jet, and make a reconstruction of the past ejection velocity time-variability (assuming ballistic knot motions). This reconstruction shows an "acceleration" of the ejection velocities of the jet knots, with higher velocities at more recent times. This acceleration will result in an eventual merging of the knots in $\sim 450$~yr and at a distance of $\sim 80"$ from the outflow source, close to the present-day position of HH~1.

We have studied the environment of the FU Ori type star V582 Aur. Our aim is to explore the star-forming region associated with this young eruptive star. Using slitless spectroscopy we searched for H alpha emission stars within a field of 11.5arcmin \times 11.5arcmin, centred on V582 Aur. Based on UKIDSS and Spitzer Space Telescope data we further selected infrared-excess young stellar object candidates. In all, we identified 68 candidate low-mass young stars, 16 of which exhibited H alpha emission in the slitless spectroscopic images. The colour-magnitude diagram of the selected objects, based on IPHAS data, suggests that they are low-mass pre-main-sequence stars associated with the Aur OB 1 association, located at a distance of 1.3 kpc from the Sun. The bright-rimmed globules in the local environment of V582 Aur probably belong to the dark cloud LDN~1516. Our results suggest that star formation in these globules might have been triggered by the radiation field of a few hot members of Aur OB 1. The bolometric luminosity of V582 Aur, based on archival photometric data and on the adopted distance, is 150-320 Lsun.

We have analyzed four epochs of H$\alpha$ and [S~II] HST images of the HH~1/2 outflow (covering a time interval from 1994 to 2014) to determine proper motions and emission line fluxes of the knots of HH~2. We find that our new proper motions agree surprisingly well with the motions measured by Herbig \& Jones (1981), although there is partial evidence for a slight deceleration of the motion of the HH~2 knots from 1945 to 2014. We also measure the time-variability of the H$\alpha$ intensities and the [S~II]/H$\alpha$ line ratios, and find that knots H and A have the largest intensity variabilities (in $1994\to 2014$). Knot H (which now dominates the HH~2 emission) has strengthened substantially, while keeping an approximately constant [S~II]/H$\alpha$ ratio. Knot A has dramatically faded, and at the same time has had a substantial increase in its [S~II]/H$\alpha$ ratio. Possible interpretations of these results are discussed.

BHR 160 is a virtually unstudied cometary globule within the Sco OB4 association in Scorpius at a distance of 1600pc. It is part of a system of cometary clouds which face the luminous O star HD155806. BHR 160 is special because it has an intense bright rim. We attempt to derive physical parameters for BHR 160 and to understand its structure and the origin of its peculiar bright rim. BHR 160 was mapped in the $^{12}$CO, $^{13}$CO and C$^{18}$O (2-1) and (1-0) and CS (3-2) and (2-1) lines. These data, augmented with stellar photometry derived from the ESO VVV survey, were used to derive the mass and distribution of molecular material in BHR 160 and its surroundings. Archival mid-infrared data from the WISE satellite was used to find IR excess stars in the globule and its neighbourhood. An elongated 1' by 0.6' core lies adjacent to the globule bright rim. $^{12}$CO emission covers the whole globule, but the $^{13}$CO, C$^{18}$O and CS emission is more concentrated to the core. The $^{12}$CO line profiles indicate the presence of outflowing material near the core, but the spatial resolution of the mm data is not sufficient for a detailed spatial analysis. The BHR 160 mass estimated from the C$^{18}$O mapping is 100$\pm$50Msun(d/1.6kpc)$^2$ where d is the distance to the globule. Approximately 70 percent of the mass lies in the dense core. The total mass of molecular gas in the direction of BHR 160 is 210$\pm$(d/1.6kpc)$^2$ Msun when estimated from the more extended VVV NIR photometry. We argue that the bright rim of BHR 160 is produced by a close-by early B-type star, HD 319648, that was likely recently born in the globule. This star is likely to have triggered the formation of a source, IRS 1, that is embedded within the core of the globule and detected only in Ks and by WISE and IRAS.

We present an analysis of four epochs of H$\alpha$ and [S II] $\lambda\lambda$ 6716/6731 HST images of HH 1. For determining proper motions we explore a new method based on analysis of spatially degraded images obtained convolving the images with wavelet functions of chosen widths. With this procedure we are able to generate maps of proper motion velocities along and across the outflow axis, as well as (angularly integrated) proper motion velocity distributions. From the four available epochs, we find the time evolution of the velocities, intensities and spatial distribution of the line emission. We find that over the last two decades HH 1 shows a clear acceleration. Also, the H$\alpha$ and [S II] intensities have first dropped, and then recovered in the more recent (2014) images. Finally, we show a comparison between the two available HST epochs of [O III] $\lambda$ 5007 (1994 and 2014), in which we see a clear drop in the value of the [O III]/H$\alpha$ ratio.

We have performed an optical spectroscopic and photometric search for young stellar objects associated with the molecular cloud Lynds 1340, and examined the structure of the cloud by constructing an extinction map, based on SDSS data. The new extinction map suggests a shallow, strongly fragmented cloud, having a mass of some 3700~Msun. Longslit spectroscopic observations of the brightest stars over the area of L1340 revealed that the most massive star associated with L1340 is a B4 type, about 5 solar mass star. The new spectroscopic and photometric data of the intermediate mass members led to a revised distance of 825 (+110 /-80) pc, and revealed seven members of the young stellar population with M > 2 solar masses. Our search for H alpha emission line stars, conducted with the Wide Field Grism Spectrograph 2 on the 2.2-meter telescope of the University of Hawaii and covering a 30 arcmin x 40 arcmin area, resulted in the detection of 75 candidate low-mass pre-main sequence stars, 58 of which are new. We constructed spectral energy distributions of our target stars, based on SDSS, 2MASS, Spitzer, and WISE photometric data, derived their spectral types, extinctions, and luminosities from BVRIJ fluxes, estimated masses by means of pre-main sequence evolutionary models, and examined the disk properties utilizing the 2-24 micron interval of the spectral energy distribution. We measured the equivalent width of the H alpha lines and derived accretion rates. The optically selected sample of pre-main sequence stars has a median effective temperature of 3970 K, stellar mass 0.7 Msun, and accretion rate of 7.6 10^{-9} Msun/yr.

Using the United Kingdom Infrared Telescope on Mauna Kea, we have carried out a new near-infrared J, H, K monitoring survey of almost a square degree of the star-forming Orion Nebula Cluster with observations on 120 nights over three observing seasons, spanning a total of 894 days. We monitored ~15,000 stars down to J=20 using the WFCAM instrument, and have extracted 1203 significantly variable stars from our data. By studying variability in young stellar objects (YSOs) in the H-K, K color-magnitude diagram, we are able to distinguish between physical mechanisms of variability. Many variables show color behavior indicating either dust-extinction or disk/accretion activity, but we find that when monitored for longer periods of time, a number of stars shift between these two variability mechanisms. Further, we show that the intrinsic timescale of disk/accretion variability in young stars is longer than that of dust-extinction variability. We confirm that variability amplitude is statistically correlated with evolutionary class in all bands and colors. Our investigations of these 1203 variables have revealed 73 periodic AA Tau type variables, many large-amplitude and long-period (P > 15 day) YSOs, including three stars showing widely-spaced periodic brightening events consistent with circumbinary disk activity, and four new eclipsing binaries. These phenomena and others indicate the activity of long-term disk/accretion variability processes taking place in young stars. We have made the light curves and associated data for these 1203 variables available online.

The luminous Class I protostar HBC 494, embedded in the Orion A cloud, is associated with a pair of reflection nebulae, Re50 and Re50N, which appeared sometime between 1955 and 1979. We have found that a dramatic brightening of Re50N has taken place sometime between 2006 and 2014. This could result if the embedded source is undergoing a FUor eruption. However, the near-infrared spectrum shows a featureless very red continuum, in contrast to the strong CO bandhead absorption displayed by FUors. Such heavy veiling, and the high luminosity of the protostar, is indicative of strong accretion but seemingly not in the manner of typical FUors. We favor the alternative explanation that the major brightening of Re50N and the simultaneous fading of Re50 is caused by curtains of obscuring material that cast patterns of illumination and shadows across the surface of the molecular cloud. This is likely occurring as an outflow cavity surrounding the embedded protostar breaks through to the surface of the molecular cloud. Several Herbig-Haro objects are found in the region.

We have carried out 200,000 N-body simulations of three identical stellar embryos with masses from a Chabrier IMF and embedded in a molecular core. The bodies are initially non-hierarchical and undergo chaotic motions, while accreting using Bondi-Hoyle accretion. The coupling of dynamics and accretion often leads to one or two dominant bodies controlling the center of the cloud core, while banishing the other(s) to the lower-density outskirts, leading to stunted growth. Eventually each system transforms either to a bound hierarchical configuration or breaks apart into separate single and binary components. The orbital motion is followed for 100 Myr. To illustrate the simulations we introduce the 'triple diagnostic diagram', which plots two dimensionless numbers against each other, representing the binary mass ratio and the mass ratio of the third body to the total system mass. Numerous freefloating BD binaries are formed in these simulations. The separation distribution function is in good correspondence with observations, showing a steep rise at close separations, peaking around 13 AU and then declining more gently. Unresolved BD triple systems may appear as wider BD binaries. Mass ratios are strongly peaked towards unity, as observed, but this is partially due to the initial assumptions. Eccentricities gradually increase towards higher values, due to the lack of viscous interactions in the simulations, which would both shrink the orbits and decrease their eccentricities. The main threat to newly born triple systems is internal instabilities, not external perturbations. Dynamical interactions in newborn triple systems of stellar embryos embedded in and accreting from a cloud core naturally form a population of freefloating BD binaries, and this mechanism may constitute a significant pathway for the formation of BD binaries. (Abstract abbreviated).

We present new H$\alpha$ and H$\beta$ images of the HH~1/2 system, and we find that the H$\alpha$/H$\beta$ ratio has high values in ridges along the leading edges of the HH~1 bow shock and of the brighter condensations of HH~2. These ridges have H$\alpha$/H$\beta=4\to 6$, which is consistent with collisional excitation from the $n=1$ to the $n=3$ and 4 levels of hydrogen in a gas of temperatures $T=1.5\to 10\times 10^4$~K. This is therefore the first direct proof that the collisional excitation/ionization region of hydrogen right behind Herbig-Haro shock fronts is detected.

We present a deep survey of Halpha emission-line stars in the M42 region using wide-field objective prism films. A total of 1699 Halpha emission-line stars were identified, of which 1025 were previously unknown, within an area of 5.5 x 5.5 degrees centred on the Trapezium Cluster. We present Halpha strength estimates, positions, and JHKs photometry extracted from 2MASS, and comparisons to previous surveys. The spatial distribution of the bulk of the stars follows the molecular cloud as seen in CO and these stars are likely to belong to the very young population of stars associated with the Orion Nebula Cluster. Additionally, there is a scattered population of Halpha emission-line stars distributed all over the region surveyed, which may consist partly of foreground stars associated with the young NGC 1980 cluster, as well as some foreground and background dMe or Be stars. The present catalogue adds a large number of candidate low-mass young stars belonging to the Orion population, selected independently of their infrared excess or X-ray emission.

We present observations of near-infrared 2.12 micro-meter molecular hydrogen outflows emerging from 1.1 mm dust continuum clumps in the North America and Pelican Nebula (NAP) complex selected from the Bolocam Galactic Plane Survey (BGPS). Hundreds of individual shocks powered by over 50 outflows from young stars are identified, indicating that the dusty molecular clumps surrounding the NGC 7000 / IC 5070 / W80 HII region are among the most active sites of on-going star formation in the Solar vicinity. A spectacular X-shaped outflow, MHO 3400, emerges from a young star system embedded in a dense clump more than a parsec from the ionization front associated with the Pelican Nebula (IC 5070). Suspected to be a binary, the source drives a pair of outflows with orientations differing by 80 degrees. Each flow exhibits S-shaped symmetry and multiple shocks indicating a pulsed and precessing jet. The `Gulf of Mexico' located south of the North America Nebula (NGC 7000), contains a dense cluster of molecular hydrogen objects (MHOs), Herbig-Haro (HH) objects, and over 300 YSOs, indicating a recent burst of star formation. The largest outflow detected thus far in the North America and Pelican Nebula complex, the 1.6 parsec long MHO 3417 flow, emerges from a 500 Solar mass BGPS clump and may be powered by a forming massive star. Several prominent outflows such as MHO 3427 appear to be powered by highly embedded YSOs only visible at a wavelength > 70 micro-meters. An `activity index' formed by dividing the number of shocks by the mass of the cloud containing their source stars is used to estimate the relative evolutionary states of Bolocam clumps. Outflows can be used as indicators of the evolutionary state of clumps detected in mm and sub-mm dust continuum surveys.

We present high angular resolution, high sensitivity 8.46 GHz (3.6 cm) radio continuum observations made toward the core of the HH~92 outflow with the Very Large Array in 2002-2003 and with the Expanded Very Large Array in 2011. We detect a group of three compact sources distributed in a region 2$''$ in extension and discuss their nature. We conclude that one of the objects (VLA 1) is the exciting source of the giant outflow associated with HH~92. In the case of HH~34 we present new 43.3 GHz (7 mm) observations that reveal the presence of a structure associated with the exciting source and elongated perpendicular to the highly collimated optical jet in the region. We propose that this 7 mm source is a circumstellar disk with radius of $\sim$80 AU and mass of $\sim$0.21 $M_\odot$.

Observations from optical to centimeter wavelengths have demonstrated that multiple systems of two or more bodies is the norm at all stellar evolutionary stages. Multiple systems are widely agreed to result from the collapse and fragmentation of cloud cores, despite the inhibiting influence of magnetic fields. Surveys of Class 0 protostars with mm interferometers have revealed a very high multiplicity frequency of about 2/3, even though there are observational difficulties in resolving close protobinaries, thus supporting the possibility that all stars could be born in multiple systems. Near-infrared adaptive optics observations of Class I protostars show a lower binary frequency relative to the Class 0 phase, a declining trend that continues through the Class II/III stages to the field population. This loss of companions is a natural consequence of dynamical interplay in small multiple systems, leading to ejection of members. We discuss observational consequences of this dynamical evolution, and its influence on circumstellar disks, and we review the evolution of circumbinary disks and their role in defining binary mass ratios. Special attention is paid to eclipsing PMS binaries, which allow for observational tests of evolutionary models of early stellar evolution. Many stars are born in clusters and small groups, and we discuss how interactions in dense stellar environments can significantly alter the distribution of binary separations through dissolution of wider binaries. The binaries and multiples we find in the field are the survivors of these internal and external destructive processes, and we provide a detailed overview of the multiplicity statistics of the field, which form a boundary condition for all models of binary evolution. Finally we discuss various formation mechanisms for massive binaries, and the properties of massive trapezia.

We report photometric observations for comet C/2012 S1 (ISON) obtained during the time period immediately after discovery (r=6.28 AU) until it moved into solar conjunction in mid-2013 June using the UH2.2m, and Gemini North 8-m telescopes on Mauna Kea, the Lowell 1.8m in Flagstaff, the Calar Alto 1.2m telescope in Spain, the VYSOS-5 telescopes on Mauna Loa Hawaii and data from the CARA network. Additional pre-discovery data from the Pan STARRS1 survey extends the light curve back to 2011 September 30 (r=9.4 AU). The images showed a similar tail morphology due to small micron sized particles throughout 2013. Observations at sub-mm wavelengths using the JCMT on 15 nights between 2013 March 9 (r=4.52 AU) and June 16 (r=3.35 AU) were used to search for CO and HCN rotation lines. No gas was detected, with upper limits for CO ranging between (3.5-4.5)E27 molec/s. Combined with published water production rate estimates we have generated ice sublimation models consistent with the photometric light curve. The inbound light curve is likely controlled by sublimation of CO2. At these distances water is not a strong contributor to the outgassing. We also infer that there was a long slow outburst of activity beginning in late 2011 peaking in mid-2013 January (r~5 AU) at which point the activity decreased again through 2013 June. We suggest that this outburst was driven by CO injecting large water ice grains into the coma. Observations as the comet came out of solar conjunction seem to confirm our models.

Our general understanding of multiple star and planet formation is primarily based on observations of young multiple systems in low density regions like Tau-Aur and Oph. Since many, if not most, of the stars are born in clusters, observational constraints from young binaries in those environments are fundamental for understanding both the formation of multiple systems and planets in multiple systems throughout the Galaxy. We build upon the largest survey for young binaries in the Orion Nebula Cluster (ONC) which is based on Hubble Space Telescope observations to derive both stellar and circumstellar properties of newborn binary systems in this cluster environment. We present Adaptive Optics spatially-resolved JHKL'-band photometry and K-band R$\sim$\,5000 spectra for a sample of 8 ONC binary systems from this database. We characterize the stellar properties of binary components and obtain a census of protoplanetary disks through K-L' color excess. For a combined sample of ONC binaries including 7 additional systems with NIR spectroscopy from the literature, we derive mass ratio and relative age distributions. We compare the stellar and circumstellar properties of binaries in ONC with those in Tau-Aur and Oph from samples of binaries with stellar properties derived for each component from spectra and/or visual photometry and with a disk census obtained through K-L color excess. The mass ratio distribution of ONC binaries is found to be indistinguishable from that of Tau-Aur and, to some extent, to that of Oph in the separation range 85-560\,AU and for primary mass in the range 0.15 to 0.8\,M$_{\sun}$.A trend toward a lower mass ratio with larger separation is suggested in ONC binaries which is not seen in Tau-Aur binaries.The components of ONC binaries are found to be significantly more coeval than the overall ONC population and as coeval as components of binaries in Tau-Aur and Oph[...]

We have searched the Lynds 988 dark cloud complex for optical (\ha{} and \sii{}) and near-IR (\mh{} 2.12 $\mu$m) shocks from protostellar outflows. We find 20 new Herbig-Haro objects and 6 new \mh{} shocks (MHO objects), 3 of which are cross detections. Using the morphology in the optical and near-IR, we connect several of these shocks into at least 5 distinct outflow systems and identify their source protostars from catalogs of infrared sources. Two outflows in the cloud, from IRAS 21014+5001 and IRAS 21007+4951, are in excess of 1 pc in length. The IRAS 21007+4951 outflow has carved a large cavity in the cloud through which background stars can be seen. Also, we have found an optical shock which is the counterflow to the previously discovered "northwest outflow" from LkHa 324SE.

The formation of very wide binaries, such as the alpha Cen system with Proxima (also known as alpha Centauri C) separated from alpha Centauri (which itself is a close binary A/B) by 15000 AU, challenges current theories of star formation, because their separation can exceed the typical size of a collapsing cloud core. Various hypotheses have been proposed to overcome this problem, including the suggestion that ultra-wide binaries result from the dissolution of a star cluster - when a cluster star gravitationally captures another, distant, cluster star. Recent observations have shown that very wide binaries are frequently members of triple systems and that close binaries often have a distant third companion. Here we report Nbody simulations of the dynamical evolution of newborn triple systems still embedded in their nascent cloud cores that match observations of very wide systems. We find that although the triple systems are born very compact - and therefore initially are more protected against disruption by passing stars - they can develop extreme hierarchical architectures on timescales of millions of years as one component is dynamically scattered into a very distant orbit. The energy of ejection comes from shrinking the orbits of the other two stars, often making them look from a distance like a single star. Such loosely bound triple systems will therefore appear to be very wide binaries.

We present an optical/infrared study of the dense molecular cloud, L935, dubbed "The Gulf of Mexico", which separates the North America and the Pelican nebulae, and we demonstrate that this area is a very active star forming region. A wide-field imaging study with interference filters has revealed 35 new Herbig-Haro objects in the Gulf of Mexico. A grism survey has identified 41 Halpha emission-line stars, 30 of them new. A small cluster of partly embedded pre-main sequence stars is located around the known LkHalpha 185-189 group of stars, which includes the recently erupting FUor HBC 722.

The young stellar cluster illuminating the W40 H II region, one of the nearest massive star forming regions, has been observed with the ACIS detector on board the Chandra X-ray Observatory. Due to its high obscuration, this is a poorly-studied stellar cluster with only a handful of bright stars visible in the optical band, including three OB stars identified as primary excitation sources. We detect 225 X-ray sources, of which 85% are confidently identified as young stellar members of the region. Two potential distances of the cluster, 260 pc and 600 pc, are used in the paper. Supposing the X-ray luminosity function to be universal, it supports a 600 pc distance as a lower limit for W40 and a total population of at least 600 stars down to 0.1 Mo under the assumption of a coeval population with a uniform obscuration. In fact, there is strong spatial variation in Ks-band-excess disk fraction and non-uniform obscuration due to a dust lane that is identified in absorption in optical, infrared and X-ray. The dust lane is likely part of a ring of material which includes the molecular core within W40. In contrast to the likely ongoing star formation in the dust lane, the molecular core is inactive. The star cluster has a spherical morphology, an isothermal sphere density profile, and mass segregation down to 1.5 Mo. However, other cluster properties, including a \leq{1} Myr age estimate and ongoing star formation, indicate that the cluster is not dynamically relaxed. X-ray diffuse emission and a powerful flare from a young stellar object are also reported.

We explore the origin of a population of distant companions (~1000 - 5000 AU) to Class I protostellar sources recently found by Connelley and co-workers, who noted that the companion fraction diminished as the sources evolved. Here we present N-body simulations of unstable triple systems embedded in dense cloud cores. Many companions are ejected into unbound orbits and quickly escape, but others are ejected with insufficient momentum to climb out of the potential well of the cloud core and associated binary. These loosely bound companions reach distances of many thousands of AU before falling back and eventually being ejected into escapes as the cloud cores gradually disappear. We use the term orphans to denote protostellar objects that are dynamically ejected from their placental cloud cores, either escaping or for a time being tenuously bound at large separations. Half of all triple systems are found to disintegrate during the protostellar stage, so if multiple systems are a frequent outcome of the collapse of a cloud core, then orphans should be common. Bound orphans are associated with embedded close protostellar binaries, but escaping orphans can travel as far as ~0.2 pc during the protostellar phase. The steep climb out of a potential well ensures that orphans are not kinematically distinct from young stars born with a less violent pre-history. The identification of orphans outside their heavily extincted cloud cores will allow the detailed study of protostars high up on their Hayashi tracks at near-infrared and in some cases even at optical wavelengths.

We present deep 3.6 cm radio continuum observations of the HII region W40 obtained using the Very Large Array in its A and B configurations. We detect a total of 20 compact radio sources in a region of 4 x 4 arcmin, with 11 of them concentrated in a band with 30 arcsec of extent. We also present JHK photometry of the W40 cluster taken with the QUIRC instrument on the University of Hawaii 2.2 meter telescope. These data reveal that 15 of the 20 VLA sources have infrared counterparts, and 10 show radio variability with periods less than 20 days. Based on these combined radio and IR data, we propose that 8 of the radio sources are candidate ultracompact HII regions, 7 are likely to be young stellar objects, and 2 may be shocked interstellar gas.

In early 2008, the young low-mass star EX Lupi, the prototype of the EXor class of eruptive variables, optically brightened by over five magnitudes for a period of 7 months. The previous time a change of such amplitude had been observed in EX Lup was over 50 years ago. In this Letter we present new optical and near-IR high resolution spectroscopy of EX~Lup during the 2008 outburst. We investigate the physical characteristics of the outburst both soon after it began and some four months later, and consider the energetics and kinematics observed. Emission line strengths, widths, and profiles changed significantly between the two observations. Also, modeling of the 2.2935 um CO overtone bandhead emission suggests that an inner gap in the circumstellar gas disk around the star may be present and it is from the inner edge of the gas disk that the CO overtone emission probably arises. We derive a mass accretion luminosity and rate during the extreme outburst of ~2+-0.5~Lsun and ~2+-0.5x10^-7 Msun yr^-1, respectively, which suggests that this outburst was indeed one of the strongest witnessed in EX Lup, yet not as strong as those observed in FU Orionis stars.

In this work we study Herbig-Haro objects located in the region around the head of the cometary globule CG 30. Two sets of optical images are presented. The first set was obtained with the 3.5 m New Technology Telescope in 1995 in three emission lines: Halpha, [SII]6731,6716 A and [OII]3729 A. The second set is an Halpha image of the CG 30/31/38 complex obtained in 2006 with the 8 m Subaru telescope. A proper motion study of the HH objects in the region was performed using the Halpha images from both epochs. Due to the high resolution of our images we were able to, for the first time, resolve the HH 120 object into ten knots and measure proper motions for some of them. We discover several new HH objects and a large bipolar jet, HH 950, emerging from the head of CG 30. We suggest that two previously known submillimeter sources are the driving sources for the HH 120 and HH 950 flows.

Cometary globule 12 is a relatively little investigated medium- and low mass star forming region 210 pc above the Galactic plane. NIR J, H, and Ks imaging and stellar photometry is used to analyse the stellar content and the structure of CG 12. Several new members and member candidates of the CG 12 stellar cluster were found. The new members include in particular a highly embedded source with a circumstellar disk or shell and a variable star with a circumstellar disk which forms a binary with a previously known A spectral type cluster member. The central source of the known collimated molecular outflow in CG 12 and an associated "hourglass"-shaped object due to reflected light from the source were also detected. HIRES-enhanced IRAS images are used together with SOFI J,H,Ks imaging to study the two associated IRAS point sources, 13546-3941 and 13547-3944. Two new 12 micrometer sources coinciding with NIR excess stars were detected in the direction of IRAS 13546-3941. The IRAS 13547-3944 emission at 12 and 25 micrometers originates in the Herbig AeBe star h4636n and the 60 and 100 micrometer emission from an adjacent cold source.

We present results from an optical imaging and spectroscopic monitoring campaign on the young, low-mass eruptive variable star V1647 Orionis. The star and associated nebulosity (McNeil's Nebula) were observed over the period February 2004 to February 2006 with observations commencing a few months after the original outburst event occurred. Using the Gemini North telescope, we obtained multi-band optical imaging photometry and medium-resolution long-slit spectroscopy of V1647 Ori on an approximately monthly interval. During this period, V1647 Ori remained at, or close to, peak brightness and then faded by 5 magnitudes to close to its pre-outburst brightness. This implies an outburst timescale of around 27 months. Spectral features seen in both emission and absorption varied considerably during the monitoring period. For example, the Halpha line changed significantly in both intensity and profile. We present and discuss the observed photometric and spectroscopic changes and consider how this eruptive event relates to the early formative stages of low-mass stars.

In order to test the hypothesis that Class I protostellar binary stars are a product of ejections during the dynamical decay of non-hierarchical multiple systems, we combined the results of new adaptive optics (AO) observations of Class I protostars with our previously published AO data to investigate whether Class I protostars with a widely separated companion (r>200 AU) are more likely to also have a close companion (r<200 AU). In total, we observed 47 embedded young stellar objects (YSOs) with either the Subaru natural guide star AO system or the Keck laser guide star AO system. We found that targets with a widely separated companion within 5,000 AU are not more likely to have a close companion. However, targets with another YSO within a projected separation of 25,000 AU are much more likely to have a close companion. Most importantly, every target with a close companion has another YSO within a projected separation of 25,000 AU. We came to the same conclusions after considering a restricted sample of targets within 500 pc and close companions wider than 50 AU to minimize incompleteness effects. The Orion star forming region was found to have an excess of both close binaries and YSOs within 25,000 AU compared to other star forming regions. We interpret these observations as strong evidence that many close Class I binary stars form via ejections and that many of the ejected stars become unbound during the Class I phase.

In late 2003, the young eruptive variable star V1647 Orionis optically brightened by over 5 magnitudes, stayed bright for around 26 months, and then decline to its pre-outburst level. In August 2008 the star was reported to have unexpectedly brightened yet again and we herein present the first detailed observations of this new outburst. Photometrically, the star is now as bright as it ever was following the 2003 eruption. Spectroscopically, a pronounced P Cygni profile is again seen in Halpha with an absorption trough extending to -700 km/s. In the near-infrared, the spectrum now possesses very weak CO overtone bandhead absorption in contrast to the strong bandhead emission seen soon after the 2003 event. Water vapor absorption is also much stronger than previously seen. We discuss the current outburst below and relate it to the earlier event.

We have mapped the central region of the HH 111 protostellar system in 1.33 mm continuum, C18O(J=2-1), 13CO (J=2-1), and SO (N_J=5_6-4_5) emission at ~3" resolution with the Submillimeter Array. There are two sources, VLA 1 (=IRAS 05491+0247) and VLA 2, with the VLA 1 source driving the HH 111 jet. Thermal emission is seen in 1.33 mm continuum tracing the dust in the envelope and the putative disks around the sources. A flattened, torus-like envelope is seen in C18O and 13CO around the VLA 1 source surrounding the dust lane perpendicular to the jet axis, with an inner radius of ~ 400 AU (1"), an outer radius of ~ 3200 AU (8"), and a thickness of ~ 1000 AU (2.5"). It seems to be infalling toward the center with conservation of specific angular momentum rather than with a Keplerian rotation as assumed by Yang et al. 1997. An inner envelope is seen in SO, with a radius of ~ 500 AU (1.3"). The inner part of this inner envelope, which is spatially coincident with the dust lane, seems to have a differential rotation and thus may have formed a rotationally supported disk. The outer part of this inner envelope, however, may have a rotation velocity decreasing toward the center and thus represent a region where an infalling envelope is in transition to a rotationally supported disk. A brief comparison with a collapsing model suggests that the flattened, torus-like envelope seen in C18O and 13CO could result from a collapse of a magnetized rotating toroid.

We present new Keck II NIRSPEC high-spectral resolution 2 um echelle observations of the young eruptive variable star V1647 Orionis. This star went into outburst in late 2003 and faded to its pre-outburst brightness after approximately 26 months. V1647 Orionis is the illuminating star of McNeil's Nebula and is located near M 78 in the Lynds 1630 dark cloud. Our spectra have a resolving power of approximately 18,000 and allow us to study in detail the weak absorption features present on the strong near-IR veiled continuum. An analysis of the echelle orders containing Mg I (2.1066 um) and Al I (2.1099 um), Br-gamma (2.1661 um), the Na I doublet (2.206 and 2.209 um), and the CO overtone bandhead (2.2935 um) gives us considerable information on the physical and geometric characteristics of the regions producing these spectral features. We find that, at high-spectral resolution, V1647 Orionis in quiescence resembles a significant number of FU Orionis type eruptive variables and does not appear similar to the quiescent EX Lupi variables observed. This correspondence is discussed and implications for the evolutionary state of the star are considered.

We present sensitive, high angular resolution ($0\rlap.{''}05$) VLA continuum observations made at 7 mm of the core of the HH 111/121 quadrupolar outflow. We estimate that at this wavelength the continuum emission is dominated by dust, although a significant free-free contribution ($\sim$30%) is still present. The observed structure is formed by two overlapping, elongated sources approximately perpendicular to each other as viewed from Earth. We interpret this structure as either tracing two circumstellar disks that exist around each of the protostars of the close binary source at the core of this quadrupolar outflow or a disk and a jet perpendicular to it. Both interpretations have advantages and disadvantages, and future high angular resolution spectroscopic millimeter observations are required to favor one of them in a more conclusive way.

The W40 complex is a nearby site of recent massive star formation composed of a dense molecular cloud adjacent to an HII region that contains an embedded OB star cluster. The HII region is beginning to blister out and break free from its envelope of molecular gas, but our line of sight to the central stars is largely obscured by intervening dust. Several bright OB stars in W40 - visible at optical, infrared, or cm wavelengths - are providing the ionizing flux that heats the HII region. The known stellar component of W40 is dominated by a small number of partly or fully embedded OB stars which have been studied at various wavelengths, but the lower mass stellar population remains largely unexamined. Despite its modest optical appearance, at 600pc W40 is one of the nearest massive star forming regions, and with a UV flux of about 1/10th of the Orion Nebula Cluster, this neglected region deserves detailed investigation.

We present the observational results of a near-infrared survey of a large sample of Class I protostars designed to determine the Class I binary separation distribution from ~100 AU to ~5000 AU. We have selected targets from a new sample of 267 nearby candidate Class I objects. This sample is well understood, consists of mostly Class I young stellar objects (YSOs) within 1 kpc, has targets selected from the whole sky, and is not biased by previous studies of star formation. We have observed 189 Class I YSOs north of -40 degrees declination at H, K and L'-bands, with a median angular resolution of 0.33" at L'. We determine our detection limit for close binary companions by observing artificial binaries. We choose a contrast limit and an outer detection limit to minimize contamination and to ensure that a candidate companion is gravitationally bound. Our survey uses observations at L' rather than K-band for the detection of binary companions since there is less scattered light and better seeing at L'. This paper presents the positions of our targets, the near-IR photometry of sources detected in our fields at L', as well as the observed properties of the 89 detected companions (73 of which are newly discovered). Although we have chosen contrast and separation limits to minimize contamination, we expect that there are ~6 stars identified as binary companions that are due to contamination. Finder charts at L' for each field are shown to facilitate future studies of these objects.

We present the Class I protostellar binary separation distribution based on the data tabulated in the companion paper. We verify the excess of Class I binary stars over solar-type main-sequence stars, especially at separations beyond 500 AU. Although our sources are in nearby star forming regions distributed across the entire sky (including Orion), none of our objects are in a high stellar density environment. The binary separation distribution changes significantly during the Class I phase, and the binary frequency at separations greater than 1000 AU declines steadily with respect to spectral index. Despite these changes, the binary frequency remains constant until the end of the Class I phase, when it drops sharply. We propose a scenario to account for the changes in the Class I binary separation distribution. This scenario postulates that companions with a separation greater than ~1000 AU were ejected during the Class 0 phase, but remain gravitationally bound due to the mass of the envelope. As the envelope dissipates, these companions become unbound and the binary frequency at wide separations declines. This scenario predicts that a large number of Class 0 objects should be non-hierarchical multiple systems, and that many Class I YSOs with a widely separated companion should also have a very close companion. We also find that Class I protostars are not dynamically pristine, and have experienced dynamical evolution before they are visible as Class I objects. For the first time, evidence is presented showing that the Class I binary frequency and the binary separation distribution strongly depend on the star forming environment. The reason for this dependence remains unclear.

We present new high resolution (R=18,000) near-infrared spectroscopic observations of a sample of classical FU Orionis stars (FUors) and other young stars with FUor characteristics that are sources of Herbig-Haro flows. Spectra are presented for the region 2.203 - 2.236 microns which is rich in absorption lines sensitive to both effective temperatures and surface gravities of stars. Both FUors and FUor-like stars show numerous broad and weak unidentified spectral features in this region. Spectra of the 2.280 - 2.300 micron region are also presented, with the 2.2935 micron v=2-0 CO absorption bandhead being clearly the strongest feature seen in the spectra all FUors and Fuor-like stars. A cross-correlation analysis shows that FUor and FUor-like spectra in the 2.203 - 2.236 micron region are not consistent with late-type dwarfs, giants, nor embedded protostars. The cross-correlations also show that the observed FUor-like Herbig-Haro energy sources have spectra that are substantively similar to those of FUors. Both object groups also have similar near-infrared colors. The large line widths and double-peaked nature of the spectra of the FUor-like stars are consistent with the established accretion disk model for FUors, also consistent with their near-infrared colors. It appears that young stars with FUor-like characteristics may be more common than projected from the relatively few known classical FUors.

We present new optical, near-IR, and mid-IR observations of the young eruptive variable star V1647 Orionis that went into outburst in late 2004 for approximately two years. Our observations were taken one year after the star had faded to its pre-outburst optical brightness and show that V1647Ori is still actively accreting circumstellar material. We compare and contrast these data with existing observations of the source from both pre-outburst and outburst phases. From near-IR spectroscopy we identify photospheric absorption features for the first time that allow us to constrain the classification of the young star itself. Our best fit spectral type is M0+-2 sub-classes with a visual extinction of 19+-2 magnitudes and a K-band veiling of rK~1.5+-0.2. We estimate that V1647Ori has a quiescent bolometric luminosity of ~9.5Lsun and a mass accretion rate of ~1.10^-6Msun yr^-1. Our derived mass and age, from comparison with evolutionary models, are 0.8+-0.2 Msun and ~0.5Myrs, respectively. The presence towards the star of shock excited optical [S II] and [Fe II] emission as well as near-IR H2 and [Fe II] emission perhaps suggests that a new Herbig-Haro flow is becoming visible close to the star.

We have carried out a major survey for visual binaries towards the Orion Nebula Cluster using HST images obtained with an H-alpha filter. Among 781 likely ONC members more than 60" from theta-1 Ori C, we find 78 multiple systems (75 binaries and 3 triples), of which 55 are new discoveries, in the range from 0.1" to 1.5". About 9 binaries are likely line-of-sight associations. We find a binary fraction of 8.8%+-1.1% within the limited separation range from 67.5 to 675 AU. The field binary fraction in the same range is a factor 1.5 higher. Within the range 150 AU to 675 AU we find that T Tauri associations have a factor 2.2 more binaries than the ONC. The binary separation distribution function of the ONC shows unusual structure, with a sudden steep decrease in the number of binaries as the separation increases beyond 0.5", corresponding to 225 AU. We have measured the ratio of binaries wider than 0.5" to binaries closer than 0.5" as a function of distance from the Trapezium, and find that this ratio is significantly depressed in the inner region of the ONC. The deficit of wide binaries in the central part of the cluster is likely due to dissolution or orbital change during their passage through the potential well of the inner cluster region. Many of the companions are likely to be brown dwarfs.

We present high spatial resolution optical integral field spectroscopy of a collimated Herbig-Haro jet viewed nearly edge-on. Maps of the line emission, velocity centroid, and velocity dispersion were generated for the H$\alpha$ and [S II] emission features from the inner collimated jet and exciting source region of the HH 34 outflow. The kinematic structure of the jet shows several maxima and minima in both velocity centroid value and velocity dispersion along the jet axis. Perpendicular to the flow direction the velocity decreases outward from the axis to the limb of the jet, but the velocity dispersion increases. Maps of the electron density structure were derived from the line ratio of [S II] 6731/6716 emission. We have found that the jet exhibits a pronounced ``striped'' pattern in electron density; the high $n_e$ regions are at the leading side of each of the emission knots in the collimated jet, and low $n_e$ regions in the down-flow direction. On average, the measured electron density decreases outward from the inner regions of the jet, but the highest $n_e$ found in the outflow is spatially offset from the nominal position of the exciting star. The results of our high spatial resolution optical integral field spectroscopy show very good agreement with the kinematics and electron density structure predicted by the existing internal working surface models of the HH~34 outflow.

We present a K-band atlas of 106 reflection nebulae, 41 of which are new discoveries. We observed these nebulae with the UH 2.2 m telescope in the course of an imaging survey of 197 objects that were selected to be nearby young Class I sources. K-band images and flux calibrated surface brightness contour plots of each nebula are presented. We found that the near-IR luminosities and physical sizes of the nebulae increase with the bolometric luminosity of the illuminating sources. Only 22 nebulae, about 10% of these candidate Class I sources, have indications of shocked H_2 emission. The great variety of nebulae that we observed prevented us from classifying them based on morphology. However, we note that as the spectral index decreases, the central star is more frequently visible at K-band and the flux from the central star tends to be dominant over the flux from the nebula. For objects that have a higher spectral index, most of the K-band flux is from the reflection nebula, and the central star is less frequently visible. The nebula around IRAS 05450+0019 has a unique morphology, and we speculate that it may be an example of a disk shadow being projected into the surrounding cloud. We present J, H, and K-band images of this object with surface brightness contours, as well as its SED from 1.2 microns to 100 microns.

We present millimeter line observations of the HH 111 outflow and its driving source. The molecular gas emission observed with IRAM 30m and the CSO reveals a small condensation of cold and dense gas. The low-velocity outflow has been mapped with the IRAM PdBI interferometer. The cold gas is distributed in a hollow cylinder surrounding the optical jet. The formation of this cavity and its kinematics are well accounted for in the frame of outflow gas entrainment by jet bow shocks. Evidence of gas acceleration is found along the cavity walls, correlated with the presence of optical bow shocks. The cavity has been expanding with a mean velocity of 4 km/s on a timescale of 8700 yr, similar to the dynamical age of the optical jet. The separation of the inner walls reaches 8"-10", which matches the transverse size of the wings in the bow shock. CSO observations of the J=7-6 line show evidence of a high-velocity and hot gas component (T=300-1000 K) with a low filling factor, associated with shocked molecular gas in the jet. [CI] observations are consistent with C-type non-dissociative shocks. Mapping of the high-velocity molecular bullets B1-B3 located beyond the optical jet, with the PdBI, reveals small structures of 3" by 7" flattened perpendicular to the flow direction. They are made of cold gas of moderate density(a few 10^3 cm-3). The bullets appear to expand into the low-density surrounding medium. We conclude that they are probably shocked gas knots resulting from past time-variable ejections in the jet.

We review five mechanisms for forming brown dwarfs: (i) turbulent fragmentation of molecular clouds, producing very low-mass prestellar cores by shock compression; (ii) collapse and fragmentation of more massive prestellar cores; (iii) disc fragmentation; (iv) premature ejection of protostellar embryos from their natal cores; and (v) photo-erosion of pre-existing cores overrun by HII regions. These mechanisms are not mutually exclusive. Their relative importance probably depends on environment, and should be judged by their ability to reproduce the brown-dwarf IMF, the distribution and kinematics of newly formed brown dwarfs, the binary statistics of brown dwarfs, the ability of brown dwarfs to retain discs, and hence their ability to sustain accretion and outflows. This will require more sophisticated numerical modelling than is presently possible, in particular more realistic initial conditions and more realistic treatments of radiation transport, angular momentum transport and magnetic fields. We discuss the minimum mass for brown dwarfs, and how brown dwarfs should be distinguished from planets.

We present a proper motion study of the shock waves within the classic stellar jet HH 47 based on Hubble Space Telescope H-alpha and [S II] images of the region taken over two epochs. Individual knots within the jet and in the bow shock/Mach disk working surface of HH 47A move significantly in the five years that separate the images, and the excellent spatial resolution of HST makes it possible to measure the proper motions with enough precision to easily observe differential motions throughout the flow. The bright portion of the jet emerges at 37.5 +/- 2.5 degrees from the plane of the sky with an average velocity of 300 km/s. Dynamical ages of the shock waves in the jet range from a few decades for knots recently ejected by the source to ~ 1300 years for the faint extended bow shock HH 47D. The jet curves, but motions of knots in the jet are directed radially away from the exciting source, and velocity variability in the flow drives the shock waves that heat the jet internally. The jet orientation angle varies with time by about 15 degrees, and currently points to the northwestern portion of a cavity outlined by a reflection nebula, where a quasi-stationary shock deflects the jet. The major working surface HH 47A is more complex than a simple bow shock/Mach disk, and contains numerous clumps that move relative to one another with velocities of ~ +/- 40 km/s. Small clumps or instabilities affect the Mach disk, and dense clumps may move all the way through the working surface to cause the bumpy morphology seen at the bow shock. A localized area between the bow shock and Mach disk varies significantly between the two sets of images.

We present new results on the kinematics of the jet HH 110. New proper motion measurements have been calculated from [SII] CCD images obtained with a time baseline of nearly fifteen years. HH 110 proper motions show a strong asymmetry with respect to the outflow axis, with a general trend of pointing towards the west of the axis direction. Spatial velocities have been obtained by combining the proper motions and radial velocities from Fabry-Perot data. Velocities decrease by a factor ~3 over a distance of ~10$^{18}$ cm, much shorter than the distances expected for the braking caused by the jet/environment interaction. Our results show evidence of an anomalously strong interaction between the outflow and the surrounding environment, and are compatible with the scenario in which HH 110 emerges from a deflection in a jet/cloud collision.

We present a study of McNeil's Nebula, a newly appeared reflection nebula in the L1630 cloud, together with photometry and spectroscopy of its source. New IR photometry compared to earlier 2MASS data shows that the star has brightened by about 3 magnitudes in the near-infrared, changing its location in a J-H/H-K diagram precisely along a reddening vector. A Gemini NIRI K-band spectrum shows strong CO-bandhead emission and Br-gamma is in emission, indicative of strong accretion. A Gemini GMOS optical spectrum shows only a red, heavily veiled continuum, with H-alpha strongly in emission and displaying a pronounced P Cygni profile, with an absorption trough reaching velocities up to 600 km s-1. This implies significant mass loss in a powerful wind. However, no evidence is found for any shocks, as commonly seen in collimated outflows from young stars. Apparently the eruption has dispersed a layer of extinction and this, together with the intrinsic brightening of the IRAS source, has allowed an earlier outflow cavity to be flooded with light, thus creating McNeil's Nebula.

We present images of the young stellar object ASR 41 in the NGC 1333 star forming region at the wavelengths of H_alpha and [SII] and in the I, J, H, and K-bands. ASR 41 has the near-infrared morphology of an edge-on disk object, but appears an order of magnitude larger than typical systems of this kind. We also present detailed models of the scattering and radiative transfer in systems consisting of a young star surrounded by a proto-planetary disk, and the whole system being embedded in either an infalling envelope or a uniform molecular cloud. The best fit to the observed morphology can be achieved with a disk of approx. 200 AU diameter, immersed in a low density cloud. The low cloud density is necessary to stay below the sub-mm flux upper limits and to preserve the shadow cast by the disk via single scattering. The results demonstrate that ASR 41 is probably not inherently different from typical edge-on disk objects, and that its large apparent size is due to the shadow of a much smaller disk being projected into the surrounding dusty molecular material.

Adaptive optics imaging of the bright visual T Tauri binary AS 353 with the Subaru Telescope shows that it is a hierarchical triple system. The secondary component, located 5.6" south of AS 353A, is resolved into a subarcsecond binary, AS 353Ba and Bb, separated by 0.24". Resolved spectroscopy of the two close components shows that both have nearly identical spectral types of about M1.5. Whereas AS 353A and Ba show clear evidence for an infrared excess, AS 353Bb does not. We discuss the possible role of multiplicity in launching the large Herbig-Haro flow associated with AS 353A.

We present deep 3.6 cm radio continuum observations of the H II region NGC 2024 in Orion B obtained using the Very Large Array in its A-configuration, with $0\rlap.{''}2$ angular resolution. We detect a total of 25 compact radio sources in a region of $4' \times 4'$. We discuss the nature of these sources and its relation with the infrared and X-ray objects in the region. At least two of the radio sources are obscured proplyds whose morphology can be used to restrict the location of the main ionizing source of the region. This cluster of radio sources is compared with others that have been found in regions of recent star formation.

We have obtained a Halpha position-velocity cube from Fabry-Perot interferometric observations of the HH 110 flow. We analyze the results in terms of anisotropic wavelet transforms, from which we derive the spatial distribution of the knots as well as their characteristic sizes (along and across the outflow axis). We then study the spatial behaviour of the line width and the central radial velocity. The results are interpreted in terms of a simple ``mean flow+turbulent eddy'' jet/wake model. We find that most of the observed kinematics appear to be a direct result of the mean flow, on which are superposed low amplitude (35 km/s) turbulent velocities.

We discuss a scenario in which brown dwarfs are formed like stars, except that their full collapse phases are interrupted through dynamical interactions in small multiple systems, leading to the ejection of the lightest member. This disintegration is a stochastic process, often resulting in the expulsion of newborn low mass stars, but when it occurs early enough the ejected stellar embryo will be a substellar object. This process may be so common at early ages that a large fraction of the ubiquitous brown dwarfs could have formed in this manner. Detailed gas dynamical simulations are required in order to better understand the details of the decay of small newborn multiple systems. We discuss the observational consequences of the ejection hypothesis, noting especially the importance of binaries with brown dwarf components as an observational test. Finally, we note that brown dwarfs that have recently been ejected may be so disturbed, by infall from the collapsing core and also by heavy accretion from perturbed circumstellar disks, that traditional spectral and luminosity criteria may fail to identify their substellar nature.

NGC 1333, a highly active star formation region within the Perseus molecular cloud complex, has been observed with the ACIS-I detector on board the Chandra X-ray Observatory. In our image with a sensitivity limit of ~ 1e28 erg/s, we detect 127 X-ray sources, of which, 95 are known cluster members. The X-ray luminosity function of the discovered YSO population spans a range of 1e28-3.2e31 erg/s in the 0.5-8 keV band. Comparison with K-band source counts indicates that we detect all of the known cluster members with K < 12, and about half of members with K > 12. We detect seven X-ray emitting YSOs which drive optically visible jets as well as one deeply embedded object that has not been detected in near-IR observations. The presence or absence of an outflows does not appear to produce any difference in X-ray properties of YSOs. We also find no systematic differences in X-ray luminosity distributions between two complete subsamples of CTTs and WTTs. That suggests that there is no difference in the astrophysical mechanism of WTT and CTT X-ray emission production. Additional results include: the X-ray emission from two late-B stars which illuminate the reflection nebula originates from unresolved late-type companions; two T Tauri stars are discovered as previously unknown components of visual binaries; and the X-ray counterpart of SVS 16 has the column density much lower than that expected from near-IR photometry and thus its X-ray luminosity is not anomalously high, as has been previously suggested.

We present a long-slit spectrum of the Herbig-Haro object HH 111, obtained with SITS on the Hubble Space Telescope. This spectrum has a spectral resolution of approx. 50 km/s, so it gives a good picture of the kinematical properties of the observed object at very high angular resolution. We find that some of the knots along the jet are associated with sudden drops in the radial velocity (modulus), confirming that the emission from the knots is formed in shocks. We interprt the observed position-velocity diagrams in terms of a model of a jet from a variable source, and we attempt to carry out a reconstruction of the ejection variability necessary for reproducing the observed kinematical structure.

Brown dwarfs may have such low masses because they are prematurely ejected from small unstable multiple systems, while the members are still actively building up their masses. We demonstrate that this scenario is consistent with all currently existing observations of brown dwarfs, and propose further observational tests. We review the status of the latest realistic numerical simulations of disintegrating small N clusters, which show that many of the ejected members end up with masses that are substellar, drifting away from their birth region with velocities rarely exceeding 2 km/s.

We have observed the Orion Molecular Clouds 2 and 3 (OMC-2 and OMC-3) with the Chandra X-ray Observatory (CXO). The northern part of OMC-3 is found to be particularly rich in new X-ray features; four hard X-ray sources are located in and along the filament of cloud cores. Two sources coincide positionally with the sub$mm$-$mm$ dust condensations of MMS 2 and 3 or an outflow radio source VLA 1, which are in a very early phase of star formation. The X-ray spectra of these sources show an absorption column of (1-3) x 10^23 H cm-2. Assuming a moderate temperature plasma, the X-ray luminosity in the 0.5-10 keV band is estimated to be ~10^30 erg s^-1 at a distance of 450 pc. From the large absorption, positional coincidence and moderate luminosity, we infer that the hard X-rays are coming from very young stellar objects embedded in the molecular cloud cores. We found another hard X-ray source near the edge of the dust filament. The extremely high absorption of 3 x 10^23 H cm^-2 indicates that the source must be surrounded by dense gas, suggesting that it is either a YSO in an early accretion phase or a Type II AGN (e.g. a Seyfert 2), although no counterpart is found at any other wavelength. In contrast to the hard X-ray sources, soft X-ray sources are found spread around the dust filaments, most of which are identified with IR sources in the T Tauri phase.

We conjecture that brown dwarfs are substellar objects because they have been ejected from small newborn multiple systems which have decayed in dynamical interactions. In this view, brown dwarfs are stellar embryos for which the star formation process was aborted before the hydrostatic cores could build up enough mass to eventually start hydrogen burning. The disintegration of a small multiple system is a stochastic process, which can be described only in terms of the half-life of the decay. A stellar embryo competes with its siblings in order to accrete infalling matter, and the one that grows slowest is most likely to be ejected. With better luck, a brown dwarf would therefore have become a normal star. This interpretation of brown dwarfs readily explains the rarity of brown dwarfs as companions to normal stars (aka the ``brown dwarf desert''), the absence of wide brown dwarf binaries, and the flattening of the low mass end of the initial mass function. Possible observational tests of this scenario include statistics of brown dwarfs near Class 0 sources, and the kinematics of brown dwarfs in star forming regions while they still retain a kinematic signature of their expulsion. Because the ejection process limits the amount of gas brought along in a disk, it is predicted that substellar equivalents to the classical T Tauri stars should be very rare.

We report the discovery of a dozen Herbig-Haro jets illuminated by the Lyman continuum and/or softer far-ultraviolet radiation fields of nearby high mass stars. Five irradiated outflows lie in the outer parts of the Orion Nebula and seven lie near the reflection nebula NGC 1333 in the Perseus molecular cloud. We propose that the UV radiation field has eroded residual material left over from the formation of these young source stars. Many of the irradiated jets exhibit unusual C-shaped bending. In the outskirts of the Orion Nebula, most irradiated jets bend away from the core of the nebula. On the other hand, in NGC 1333, the C-shaped jets bend back towards the cluster center. Jet bending in the Orion Nebula may be dominated by either the outflow of material from the nebular core or by the rocket effect pushing on the irradiated portion of a mostly neutral jet beam. But in NGC 1333, jet bending may indicate that the source stars have been ejected from the cluster core. Many irradiated jets are asymmetric with one beam much brighter than the other. When fully photo-ionized, irradiated jets may provide unique insights into the physical conditions within outflows powered by young stars, permitting the determination of the density and location of stellar ejecta even in the absence of shocks. We present a model for the photo-ionization of these outflows by external radiation fields and discuss possible mechanisms for producing the observed asymmetries. We demonstrate that the UV radiation field may alter the amount of cloud material entrained by the jet. We also report the discovery of some large scale bow shocks which face the core of the Orion Nebula and which surround visible young stars. These wind-wind collision fronts provide further evidence for a large scale mass flow originating near the nebular core.

We present a millimeter wave study of a cluster of bipolar CO outflows embedded in the western end of the Circinus molecular cloud complex, G317-4, that is traced by very high optical extinction. For an assumed distance of 700 pc, the entire Circinus cloud is estimated to have a mass of about 5E4 solar masses. The opaque western portion that was mapped in this study has a mass of about 10E3 solar masses, contains a number of embedded infrared sources and various compact 1.3 mm continuum sources, and has a remarkable filamentary structure with numerous cavities which appears to be the fossil remnants of past star formation activity. The profusion of outflows in this region are disentangled and linked to driving IRAS sources. Thus, the mapped portion of Circinus contains at least 10 CO emitting molecular outflows. Assuming that star formation has continued at a steady rate for the last several hundred thousand years, the Circinus cloud is expected to have produced dozens of young stars. Their outflows have severely altered the structure and kinematics of this cloud as evidenced by the multitude of prominent cavities and dust filaments that outline their boundaries. This level of star formation activity is consistent with the numerous post-outflow phase H-alpha emission line stars that have been found in this region. The Circinus cloud complex is an archetypical case where star formation activity may have profoundly affected the structure of a molecular cloud, producing its strikingly filamentary and cavitated appearance and providing further evidence that star formation may be a self regulated process.

We have observed an outburst of the T Tauri star EX Lup in March 1994. We present both photometric (BVR) and spectroscopic (low and medium resolution) observations carried out during the decline after outburst. The star appears much bluer during outburst due to an increased emission of a hot continuum. This is accompanied by a strong increase of the veiling of photospheric lines. We observe inverse P Cygni profiles of many emission lines over a large brightness range of EX Lup. We briefly discuss these features towards the model of magnetospherically supported accretion of disk material.

We present an overview of what can be done in the field of star formation using the VLT and its advanced instrumentation, and especially with the very high angular resolution provided by the VLT in the VLTI mode. Herbig-Haro objects and early stellar evolution, pre-main sequence binaries, circumstellar disks and surface properties of young stars are investigated. We show that the VLT and VISA will be not only very useful to better understand star formation processes, but indeed necessary in order to advance in certain subjects.