This thesis presents my work characterising star formation in Gould Belt molecular clouds using submillimetre observations from SCUBA-2 on the James Clerk Maxwell Telescope (JCMT). I use these observations alongside data from previously published surveys using instruments including the Spitzer Space Telescope. I investigate the effect of including submillimetre data on the numbers, classifications and lifetimes of Young Stellar Objects (YSOs) in Gould Belt molecular clouds, particularly protostars. Following a literature review, I use SCUBA-2 450 and 850 μm observations to characterise star formation in the Lupus I molecular cloud. A total of eleven previously identified YSOs are detected with SCUBA-2, as well as eleven starless cores. Two cores have masses greater than the Jeans mass, and one has a virial parameter of 1.1 0.4, meaning these cores could be unstable against collapse. I use submillimetre emission to calculate disk masses, and find that one YSO has a disk mass greater than the minimum mass solar nebula. I find that Lupus I has a high percentage of both protostars and Very Low Luminosity Objects (VeLLOs). I also fit YSO Spectral Energy Distributions (SEDs) with models, allowing protostellar envelope masses and temperatures to be calculated, and interstellar extinction to be constrained for some YSOs. The signs of recent and future star formation support the hypothesis that a shock has triggered a star forming event in Lupus I. I also use SCUBA-2 data in conjunction with archival Spitzer and Herschel data to produce SEDs for five new candidate First Hydrostatic Cores (FHSCs) in Serpens South. These observations were then fit with models by the first author of this work, Alison Young. This work was able to identify two of the FHSC candidates as probable FHSCs, and constrain the rotation rate and inclination of one of them. I use JCMT Gould Belt Survey (GBS) observations of ten molecular clouds to produce an updated catalogue of protostars in these clouds. I use the FellWalker algorithm to find individual sources in the SCUBA-2 maps, and match them to the Spitzer YSO catalogue of Dunham et al. (2015). I use bolometric temperature to classify 362 out of 592 candidates as Class 0 or Class I protostars - a factor of two increase compared to the Spitzer catalogue due to improved submillimetre coverage. I find that protostellar lifetimes of 0.59 – 0.89 Myr - approximately 25 % longer than previously estimated. I also calculate protostellar luminosities, envelope masses, and envelope temperatures, and examine the distributions. Finally, I newly identify 19 protostars as VeLLOs, and increase the number of known VeLLOs in these clouds by a factor of two.

Le milieu interstellaire (MIS) où se forment les étoiles est constitué de gaz très dilué dominé par l'hydrogène moléculaire, et de grains de poussière de taille submicrométrique. Ces poussières jouent un rôle crucial en atténuant la lumière des étoiles lointaines, protégeant ainsi les molécules du gaz des rayonnements ultra-violets, et en servant de catalyseurs à une chimie hétérogène à très basse température. Outre la synthèse de l'hydrogène moléculaire, la surface des grains permet de former des molécules organiques dites complexes comme le méthanol (CH3OH) à partir de l'hydrogénation (et la deutération) du monoxyde de carbone (CO). Les glaces ainsi formées participent à la complexification moléculaire du MIS et seront à terme intégrées au sein de disques de poussières, berceaux des astéroïdes, comètes et exo-planètes. L'objectif de cette thèse est l'étude des mécanismes d'échanges hydrogène-deuterium sur certains groupements fonctionnels de molécules organiques simples, méthanol par exemple, présentes à la surface ou dans les manteaux des grains interstellaires. La thèse est centrée sur une exploration expérimentale de ces processus en phase condensée, à l'aide d'une expérience de cryogénie synthétisant des glaces à très basse température (15K) couplée à un spectromètre infrarouge. Nous montrons que ces échanges se produisent avant la sublimation du manteau de glace sur des groupes fonctionnels capables d'établir des liaisons hydrogènes avec les molécules d'eau voisines. Le processus catalysant est vraisemblablement la cristallisation de la glace d'eau. Des études cinétiques nous permettent d'évaluer les énergies d'activation du transfert H/D (6745K) et de la transition amorphe-cristalline (8100K), et de déterminer la constante de vitesse d'échange dans le domaine de température 120-140~K. Cette constante de vitesse est, de plus, comparée à des calculs semi-classiques basés sur un traitement ab initio. En marge de ces expériences, des observations millimétriques de la molécule de méthanol en direction de proto-étoiles confirment une variabilité des abondances relatives des isotopologues simplement deutérés de cette molécule en fonction de la masse de la protoétoile. / The interstellar medium where stars are formed consists of a dilute gas which is dominated by molecular hydrogen and dust grains less than a few microm in size. The dust plays a crucial role in the attenuation of light from the stars. They also protect molecules within the gas from UV photons. Furthermore, they serve as heterogeneous catalysts for chemistry at low temperature. The surface of the grains also permit the formation of complex organic molecules such as methanol via the hydrogenation and/or deuteration of carbon monoxide. The ices are formed and subsequently participate in increasing the molecular complexity of the clouds. Finally, they are incorporated into debris disks, asteroids, comets, and exoplanets. The objective of this thesis is to study the mechanism of hydrogen/deuterium exchange within certain functionnal groups of simple organic molecules such as methanol, which are present on the surface of these grain mantles. The thesis is focused on the experimental determination of these processes in the condensed phase. This will be achieved with the aid of a cryogenic synthesis of the ices at very low temperatures coupled with infrared spectrometry. We observe that it is possible for the exchange to proceed before the sublimation of the ice mantles. However, this is only the case when the functional groups within the molecule may form hydrogen bonds with water. From our results we see that this process seems to be catalysed by the crystalization of the water ice. The kinetics study permits us to evalute the activation energy for the H/D exchange (6745 K) and for the transition from amorphous to crystaline ice (8100 K). In addition it also allows us to determine the rate constant for the exchange in the temperature range 120-140 K. In addition we have performed theoretical calulation in an attempt to elucidate the mechanism for the exchange. However, the experimental rate constant for the exchange is much larger in comparison to the one predicted by a semi- classical treatment based on the AB initio potential we have obtained. Further to this observations of methanol towards protostars have been conducted. These observations show that there is a variation in the relative abundance of the CH2DOH and CH3OD. This variation in relative abundance seems to have some dependence upon the mass of the protostar, with high mass stars showing (CH2DOH/CH3OD ≤ 1) and low/intermiediate mass stars showing (CH2DOH/CH3OD >> 3).

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.

Les étoiles se forment par effondrement gravitationnel de condensations pré-stellaires. Le jeune embryon stellaire (phase Classe 0) croît en masse par l'accrétion progressive de l'enveloppe de gaz et de poussières dans lequel il est enfoui. Par conservation du moment cinétique, si le moment du coeur pré-stellaire est totalement transféré à l'embryon pendant la phase d'accrétion, la force gravitationnelle ne peut contrer la force centrifuge et l'embryon se fragmente prématurément. Pour former une étoile comme notre Soleil, l'enveloppe en rotation doit nécessairement réduire son moment cinétique de 5 à 10 ordres de grandeur en l'évacuant ou en le redistribuant. L'un des principaux défis de la formation stellaire est de quantifier l'ampleur de ce "problème du moment cinétique" et d'identifier les mécanismes responsables de la redistribution du moment.L'objectif de cette thèse est d'étudier la cinématique des enveloppes proto-stellaires de Classe 0 afin d'établir leurs distributions de moment cinétique. Pour cela, j'ai utilisé des observations de raies moléculaires à haute résolution angulaire de l'Interféromètre du Plateau de Bure et du télescope de 30m de l'IRAM issues du large programme CALYPSO (Continuum and Lines in Young Protostellar Objects, PI : Ph. André) pour un échantillon de 12 proto-étoiles de Classe 0 à une distance d<400 pc. Cette analyse a permis de mesurer des mouvements de rotation différentielle et d'établir, pour la première fois, des distributions radiales du moment cinétique spécifique sur une grande gamme d'échelles (~50-10000 au) dans 11 des 12 enveloppes proto-stellaires de l'échantillon. Deux régimes distincts ont ainsi été mis en exergue: un profil constant à petites échelles (<1600 au) et une augmentation du moment avec le rayon aux grandes échelles (1600-10000 au).Le profil constant montre que la matière participant directement à la formation de l'étoile possède un moment cinétique spécifique (~5 10^-4 km/s pc, <1600 au) similaire à celui observé dans les petits disques entourant les étoiles T-Tauri.Les gradients de vitesse observés aux grandes échelles (>3000 au), historiquement utilisés pour mesurer la rotation des coeurs et quantifier le problème du moment cinétique, ne sont pas dus à la pure rotation des enveloppe proto-stellaires, mais sont dominés par d'autres mécanismes. Plusieurs scénarios sont donc discutés pour interpréter le changement de régime dans les profils de moment cinétique aux échelles >1600 au: une empreinte des conditions initiales de la phase pré-stellaire, un changement de mécanismes dominants (contre-rotation, transition effondrement-rotation) ou l'influence de la dynamique des filaments interstellaires (turbulence, effondrement, chocs) dans lesquels les proto-étoiles sont enfouies. / Stars form via the gravitational collapse of a pre-stellar condensation. The young stellar embryo (Class 0 phase) mass increases via the progressive accretion of the gaseous and dusty envelope within which it is buried. As a direct consequence of the angular momentum conservation, if the angular momentum of the pre-stellar core is totally transferred to the central embryo during the accretion phase, the gravitational force can not counteract the centrifugal force and the embryo fragments prematurely before reaching the main sequence. To form a star such as our Sun, the rotating envelope needs to reduce its angular momentum by 5 to 10 orders of magnitude by ejecting or redistributing it. One of the main challenges of stellar formation is to quantify the amplitude of this "angular momentum problem" and identify the mechanisms responsible for the angular momentum redistribution.The goal of this PhD thesis is to study the kinematics of Class 0 protostellar envelopes in order to probe the distribution of their angular momentum. To do this, I used high-resolution observations of molecular lines with the Plateau de Bure Interferometer and the 30m telescope at IRAM taken as part of the large programme CALYPSO (Continuum and Lines in Young Protostellar Objects, PI : Ph. André). The sample gathers 12 Class 0 protostars with distances d<400 pc. This analysis allows to measure differential rotation motions and provides, for the first time in a large sample, robust constraints on the radial distributions of specific angular momentum in a large range of scales (~50-10000 au) for 11 of the 12 protostellar envelopes targeted in the sample. Two distinct regimes are revealed: a constant profile at small scales (<1600 au) and an increasing of the angular momentum at larger radii (1600-10000 au).The constant profile shows that the specific angular momentum (~5 10-4 km/s pc, <1600 au) of the material directly involved in the star formation is similar to the value observed in the small disks surrounding the T-Tauri stars.Velocity gradients observed on large scales (>3000 au) - that are historically used to measure the rotation of the core and quantify the angular momentum problem - are not due to pure envelope rotation but can be dominated by other mechanisms. I discuss several scenarios in order to interpret this change of regime in the angular momentum profiles at scales >1600 au: the imprints of the initial conditions of the pre-stellar phase, a change of dominant mechanisms (counter-rotation, transition between infall and rotation) or the influence of the interstellar filament dynamics (turbulence, collapse, shocks) within which protostars are buried.

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.

No description available.