Aims. We numerically studied close encounters between a young stellar system hosting a massive, gravitationally fragmenting disk and an intruder diskless star with the aim of determining the evolution of fragments that have formed in the disk prior to the encounter. Methods. Numerical hydrodynamics simulations in the non-inertial frame of reference of the host star were employed to simulate the prograde and retrograde co-planar encounters. The initial configuration of the target system (star plus disk) was obtained via a separate numerical simulation featuring the gravitational collapse of a solar-mass pre-stellar core. Results. We found that close encounters can lead to the ejection of fragments that have formed in the disk of the target prior to collision. In particular, prograde encounters are more efficient in ejecting the fragments than the retrograde encounters. The masses of ejected fragments are in the brown-dwarf mass regime. They also carry away an appreciable amount of gas in their gravitational radius of influence, implying that these objects may possess extended disks or envelopes, as also previously suggested. Close encounters can also lead to the ejection of entire spiral arms, followed by fragmentation and formation of freely-floating objects straddling the planetary mass limit. However, numerical simulations with a higher resolution are needed to confirm this finding.

We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II], and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass. outflow rates, (M) over dot(w). Thereby we observe a strong correlation of (M) over dot(w) with bolometric luminosity, and with the inferred mass accretion rates of the central objects, (M) over dot(a), which continues through the Class 0 range the trend observed in Class II young stellar objects. Along this trend from large to small mass. flow rates, the different classes of young stellar objects lie in the sequence Class 0-Class I/flat-spectrum-Class II, indicating that the trend is an evolutionary sequence in which (M) over dot(a) and (M) over dot(w) decrease together with increasing age, while maintaining rough proportionality. The survey results include two that. are key tests of magnetocentrifugal outflow-acceleration mechanisms: the distribution of the outflow/accretion branching ratio b = (M) over dot(w)/(M) over dot(a), and limits on the distribution of outflow speeds. Neither rules out any of the three leading outflow-acceleration, angular-momentum-ejection mechanisms, but they provide some evidence that disk winds and accretion-powered stellar winds (APSWs) operate in many protostars. An upper edge observed in the branching-ratio distribution is consistent with the upper bound of b = 0.6 found in models of APSWs, and a large fraction (31%) of the sample have a. branching ratio sufficiently small that only disk winds, launched on scales as large as several au, have been demonstrated to account for them.

We have imaged GM Aurigae with the Hubble Space Telescope, detected its disk in scattered light at 1400 and 1650 angstrom, and compared these with observations at 3300 angstrom, 5550 angstrom, 1.1 mu m, and 1.6 mu m. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300 angstrom shows no evidence of the 24 au radius cavity that has been previously observed in submillimeter observations. Comparison with dust grain opacity models indicates that. the surface of the entire disk is populated with submicron grains. We have compiled a. spectral energy distribution from 0.1 mu m to 1 mm. and used it to constrain a model of the star + disk system that includes the submillimeter cavity using the Monte Carlo radiative transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of submicron grains interior to the submillimeter cavity wall. We suggest one explanation for this that. could be due to a planet of mass <9 M-J interior to 24 au. A unique cylindrical structure is detected in the far-UV data from the Advanced Camera for Surveys/ Solar Blind Channel. It is aligned along the system semiminor axis, but does not resemble an accretion-driven jet. The structure is limb. brightened and extends 190 +/- 35 au above the disk midplane. The inner radius of the limb. brightening is 40 +/- 10 au, just beyond the submillimeter cavity wall.

The notion that low- to intermediate-mass young stellar objects (YSOs) gain mass at a constant rate during the early stages of their evolution appears to be challenged by observations of YSOs suffering sudden increases of the rate at which they gain mass from their circumstellar discs. Also, this idea that stars spend most of their lifetime with a low accretion rate and gain most of their final mass during short-lived episodes of high accretion bursts, helps to solve some long-standing problems in stellar evolution. The original classification of eruptive variables divides them in two separate subclasses known as FU Orionis stars (FUors) and EX Lupi stars (EXors). In this classical view FUors are at an early evolutionary stage and are still gaining mass from their parent envelopes, whilst EXors are thought to be older objects only surrounded by an accretion disc. The problem with this classical view is that it excludes younger protostars which have higher accretion rates but are too deeply embedded in circumstellar matter to be observed at optical wavelengths. Optically invisible protostars have been observed to display large variability in the near-infrared. These and some recent discoveries of new eruptive variables, show characteristics that can be attributed to both of the optically-defined subclasses of eruptive variables. The new objects have been proposed to be part of a new class of eruptive variables. However, a more accepted scenario is that in fact the original classes only represent two extremes of the same phenomena. In this sense eruptive variability could be explained as arising from one physical mechanism, i.e. unsteady accretion, where a variation in the parameters of such mechanism can cause the different characteristics observed in the members of this class. With the aim of studying the incidence of episodic accretion among young stellar objects, and to characterize the nature of these eruptive variables we searched for high amplitude variability in two multi-epoch infrared surveys: the UKIDSS Galactic Plane Survey (GPS) and the Vista Variables in the Via Lactea (VVV). In order to further investigate the nature of the selected variable stars, we use photometric information arising from public surveys at near- to farinfrared wavelengths. In addition we have performed spectroscopic and photometric follow-up for a large subset of the samples arising from GPS and VVV. We analyse the widely separated two-epoch K-band photometry in the 5th, 7th and 8th data releases of the UKIDSS Galactic Plane Survey. We find 71 stars with _K > 1 mag, including 2 previously known OH/IR stars and a Nova. Even though the mid-plane is mostly excluded from the dataset, we find the majority (66%) of our sample to be within known star forming regions (SFRs), with two large concentrations in the Serpens OB2 association (11 stars) and the Cygnus-X complex (27 stars). The analysis of the multi-epoch K-band photometry of 2010-2012 data from VVV covering the Galactic disc at |b| < 1◦ yields 816 high amplitude variables, which include known variables of different classes such as high mass X-ray binaries, Novae and eclipsing binaries among others. Remarkably, 65% of the sample are found concentrated towards areas of star formation, similar to the results from GPS. In both surveys, sources in SFRs show spectral energy distributions (SEDs) that support classification as YSOs. This indicates that YSOs dominate the Galactic population of high amplitude infrared variable stars at low luminosities and therefore likely dominate the total high amplitude population. Spectroscopic follow-up allows us to confirm the pre-main sequence nature of several GPS and VVV Objects. Most objects in both samples show spectroscopic signatures that can be attributed to YSOs undergoing high states of accretion, such as veiling of photospheric features and CO emission, or show FUor-like spectra. We also find a large fraction of objects with 2.12 μm H2 emission that can be explained as arising from shock-excited emission caused by molecular outflows. Whether these molecular outflows are related to outbursts events cannot be confirmed from our data. Adding the GPS and VVV spectroscopic results, we find that between 6 and 14 objects are new additions to the FUor class from their close resemblance to the near-infrared spectra of FUors, and at least 23 more objects are new additions to the eruptive variable class. For most of these we are unable to classify them into any of the original definitions for this variable class. In any case, we are adding up to 37 new stars to the eruptive variable class which would double the current number of known objects. We note that most objects are found to be deeply embedded optically invisible stars, thus increasing the number of objects belonging to this subclass by a much larger factor. In general, objects in our samples which are found to be likely eruptive variable stars show a mixture of characteristics that can be attributed to both of the optically-defined classes. This agrees well with the recent discoveries in the literature. Finally, we are able to derive a first rough estimate on the incidence of episodic accretion among class I YSOs in the star-forming complex G305. We find that _ 9% of such objects are in a state of high accretion. This number is in agreement with previous theoretical and observational estimates among class I YSOs.

This thesis primarily uses submillimetre molecular line data from HARP, a heterodyne array on the James Clerk Maxwell Telescope (JCMT), to further investigate star formation in the Ophiuchus L1688 cloud. HARP was used to observe CO J = 3-2 isotopologues: 12CO, 13CO and C18O; and the dense gas tracer HCO+ J = 4-3. A method for calculating molecular line contamination in the SCUBA-2 450 and 850 μm dust continuum data was developed, which can be used to convert 12CO J =6-5and J =3-2 maps of integrated intensity (K km s−1) to molecular line flux (mJy beam−1) contaminating the continuum emission. Using HARP maps of 12CO J = 3-2, I quantified the amount of molecular line contamination found in the SCUBA-2 850 μm maps of three different regions, including NGC 1333 of Perseus and NGC 2071 and NGC 2024 of Orion B. Regions with ‘significant’ (i.e. > 20%) molecular line contamination correspond to molecular outflows. This method is now being used to remove molecular line contamination from regions with both SCUBA-2 dust continuum and HARP 12CO map coverage in the Gould Belt Legacy Survey (GBS). The Ophiuchus L1688 cloud was observed in all three CO J = 3-2 isotopologues. I carried out a molecular outflow analysis in the region on a list of 30 sources from the Spitzer ‘c2d’ survey [Evans et al., 2009]. Out of the 30 sources, 8 had confirmed bipolar outflows, 20 sources had ‘confused’ outflow detections and 2 sources did not have outflow detections. The Ophiuchus cloud was found to be gravitationally bound with the turbulent kinetic energy a factor of 7 lower than the gravitational binding energy. The high-velocity outflowing gas was found to be only 21% of the turbulence in the cloud, suggesting outflows are significant but not the dominant source of turbulence in the region. Other factors were found to influence the global high-velocity outflowing gas in addition to molecular outflows, including hot dust from nearby B-type stars, outflow remnants from less embedded sources and stellar winds from the Upper Scorpius OB association. To trace high density gas in the Ophiuchus L1688 cloud, HCO+ J = 4-3 was observed to further investigate the relationship between high column density and high density in the molecular cloud. Non-LTE codes RADEX and TORUS were used to develop density models corresponding to the HCO+ emission. The models involved both constant density and peaked density profiles. RADEX [van der Tak et al., 2007] models used a constant density model along the line-of-sight and indicated the HCO+ traced densities that were predominantly subthermally excited with den- sities ranging from 10^3–10^5 cm^−3. Line-of-sight estimates ranged from several parsecs to 90 pc, which was unrealistic for the Ophiuchus cloud. This lead to the implementation of peaked density profiles using the TORUS non-LTE radiative transfer code. Initial models used a ‘triangle’ density profile and a more complicated log-normal density probability density function (PDF) profile was subsequently implemented. Peaked density models were relatively successful at fitting the HCO+ data. Triangle models had density fits ranging from 0.2–2.0×10^6 cm^−3 and 0.1–0.3×10^6 cm^−3 for the 0.2 and 0.3 pc cloud length models re- spectively. Log-normal density models with constant-σ had peak density ranges from 0.2–1.0 ×10^5 cm^−3 and 0.6–2.0×10^5 cm^−3 for 0.2 and 0.3 pc models respectively. Similarly, log-normal models with varying-σ had lower and upper density limits corresponding to the range of FWHM velocities. Densities (lower and upper limits) ranged from 0.1–1.0 ×10^6 and 0.5–3.0 ×10^5 cm^-3 for the 0.2 and 0.3 pc models respectively. The result of the HCO+ density modelling indicated the distributions of starless, prestellar and protostellar cores do not have a preference for higher densities with respect to the rest of the cloud. This is contrary to past research suggesting the probability of finding a submillimetre core steeply rises as a function of column density (i.e. density; Belloche et al. 2011; Hatchell et al. 2005). Since the majority of sources are less embedded (i.e Class II/III), it is possible the evolutionary state of Ophiuchus is the main reason the small sample of Class 0/I protostars do not appear to have a preference for higher densities in the cloud.