A Submillimetre Study of Nearby Star Formation using Molecular Line Data

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.University of Exeter studentshi

The JCMT Gould Belt Survey: understanding the influence of outflows on Gould Belt clouds

JOURThis is the final version of the article. It was first published by Oxford Journals for the Royal Astronomical Society via http://dx.doi.org/10.1093/mnrasl/slv202Using James Clerk Maxwell Telescope (JCMT) Gould Belt Survey data from CO J = 3 → 2 isotopologues, we present a meta-analysis of the outflows and energetics of star-forming regions in several Gould Belt clouds. The majority of the regions are strongly gravitationally bound. There is evidence that molecular outflows transport large quantities of momentum and energy. Outflow energies are at least 20 per cent of the total turbulent kinetic energies in all of the regions studied and greater than the turbulent energy in half of the regions. However, we find no evidence that outflows increase levels of turbulence, and there is no correlation between the outflow and turbulent energies. Even though outflows in some regions contribute significantly to maintaining turbulence levels against dissipation, this relies on outflows efficiently coupling to bulk motions. Other mechanisms (e.g. supernovae) must be the main drivers of turbulence in most if not all of these regions

The JCMT Gould Belt Survey: the effect of molecular contamination in SCUBA-2 observations of Orion A

Thermal emission from cold dust grains in giant molecular clouds can be used to probe the physical properties, such as density, temperature and emissivity in star-forming regions. We present the SCUBA-2 shared-risk observations at 450 $\mu$m and 850 $\mu$m of the Orion A molecular cloud complex taken at the James Clerk Maxwell Telescope (JCMT). Previous studies showed that molecular emission lines can contribute significantly to the measured fluxes in those continuum bands. We use the HARP $^{12}$CO J=3-2 integrated intensity map for Orion A in order to evaluate the molecular line contamination and its effects on the SCUBA-2 maps. With the corrected fluxes, we have obtained a new spectral index $\alpha$ map for the thermal emission of dust in the well-known integral-shaped filament. Furthermore, we compare a sample of 33 sources, selected over the Orion A molecular cloud complex for their high $^{12}$CO J=3-2 line contamination, to 27 previously identified clumps in OMC-4. This allows us to quantify the effect of line contamination on the ratio of 850 $\mu$m to 450 $\mu$m flux densities and how it modifies the deduced spectral index of emissivity $\beta$ for the dust grains. We also show that at least one Spitzer-identified protostellar core in OMC-5 has a $^{12}$CO J=3-2 contamination level of 16 %. Furthermore, we find the strongest contamination level (44 %) towards a young star with disk near OMC-2. This work is part of the JCMT Gould Belt Legacy Survey.Comment: 13 pages, 6 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS

The JCMT Gould Belt Survey: the effect of molecular contamination in SCUBA-2 observations of Orion A

JOURThis is the final version of the article. It was first published by Oxford Journals for the Royal Astronomical Society via http://dx.doi.org/ 10.1093/mnras/stv3009Thermal emission from cold dust grains in giant molecular clouds can be used to probe the physical properties, such as density, temperature and emissivity in star-forming regions. We present the Submillimetre Common-User Bolometer Array (SCUBA-2) shared-risk observations at 450 and 850 μm of the Orion A molecular cloud complex taken at the James Clerk Maxwell Telescope (JCMT). Previous studies showed that molecular emission lines can contribute significantly to the measured fluxes in those continuum bands. We use the Heterodyne Array Receiver Programme 12CO J = 3-2 integrated intensity map for Orion A in order to evaluate the molecular line contamination and its effects on the SCUBA-2 maps. With the corrected fluxes, we have obtained a new spectral index α map for the thermal emission of dust in the well-known integral-shaped filament. Furthermore, we compare a sample of 33 sources, selected over the Orion A molecular cloud complex for their high 12CO J = 3-2 line contamination, to 27 previously identified clumps in OMC 4. This allows us to quantify the effect of line contamination on the ratio of 850-450 μm flux densities and how it modifies the deduced spectral index of emissivity β for the dust grains. We also show that at least one Spitzer-identified protostellar core in OMC 5 has a 12CO J = 3-2 contamination level of 16 per cent. Furthermore, we find the strongest contamination level (44 per cent) towards a young star with disc near OMC 2. This work is part of the JCMT Gould Belt Legacy Survey

The James Clerk Maxwell telescope Legacy Survey of the Gould Belt: a molecular line study of the Ophiuchus molecular cloud

CO, 13CO, and C18O J = 3-2 observations are presented of the Ophiuchus molecular cloud. The 13CO and C18O emission is dominated by the Oph A clump, and the Oph B1, B2, C, E, F, and J regions. The optically thin(ner) C18O line is used as a column density tracer, from which the gravitational binding energy is estimated to be 4.5 × 1039 J (2282 M⊙ km2 s-2). The turbulent kinetic energy is 6.3 × 1038 J (320 M⊙ km2 s-2), or seven times less than this, and therefore the Oph cloud as a whole is gravitationally bound. 30 protostars were searched for high-velocity gas, with 8 showing outflows, and 20 more having evidence of high-velocity gas along their lines of sight. The total outflow kinetic energy is 1.3 × 1038 J (67 M⊙ km2 s-2), corresponding to 21 per cent of the cloud's turbulent kinetic energy. Although turbulent injection by outflows is significant, but does not appear to be the dominant source of turbulence in the cloud. 105 dense molecular clumplets were identified, which had radii ˜0.01-0.05 pc, virial masses ˜0.1-12 M⊙, luminosities ˜0.001-0.1 K km s-1 pc-2, and excitation temperatures ˜10-50 K. These are consistent with the standard Giant Molecular Cloud (GMC) based size-linewidth relationships, showing that the scaling laws extend down to size scales of hundredths of a parsec, and to subsolar-mass condensations. There is however no compelling evidence that the majority of clumplets are undergoing free-fall collapse, nor that they are pressure confined

The TOP-SCOPE Survey of Planck Galactic Cold Clumps : Survey Overview and Results of an Exemplar Source, PGCC G26.53+0.17

The low dust temperatures (<14 K) of Planck Galactic cold clumps (PGCCs) make them ideal targets to probe the initial conditions and very early phase of star formation. "TOP-SCOPE" is a joint survey program targeting similar to 2000 PGCCs in J = 1-0 transitions of CO isotopologues and similar to 1000 PGCCs in 850 mu m continuum emission. The objective of the "TOP-SCOPE" survey and the joint surveys (SMT 10 m, KVN 21 m, and NRO 45 m) is to statistically study the initial conditions occurring during star formation and the evolution of molecular clouds, across a wide range of environments. The observations, data analysis, and example science cases for these surveys are introduced with an exemplar source, PGCC G26.53+0.17 (G26), which is a filamentary infrared dark cloud (IRDC). The total mass, length, and mean line mass (M/L) of the G26 filament are similar to 6200 M-circle dot, similar to 12 pc, and similar to 500 M-circle dot pc(-1), respectively. Ten massive clumps, including eight starless ones, are found along the filament. The most massive clump as a whole may still be in global collapse, while its denser part seems to be undergoing expansion owing to outflow feedback. The fragmentation in the G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, nonmagnetized, and turbulent supported cylinder. A bimodal behavior in dust emissivity spectral index (beta) distribution is found in G26, suggesting grain growth along the filament. The G26 filament may be formed owing to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.Peer reviewe

Magnetic Fields toward Ophiuchus-B Derived from SCUBA-2 Polarization Measurements

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope in Hawaii. This work is part of the B-fields in Star-forming Region Observations survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over ~0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically toward the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large-scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630 ± 410 μG in the Oph-B2 sub-clump using a Davis–Chandrasekhar–Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio λ = 1.6 ± 1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical

How do stars gain their mass? A JCMT/SCUBA-2 transient survey of protostars in nearby star-forming regions

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion—prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In A James Clerk Maxwell Telescope/SCUBA-2 Transient Survey of Protostars in Nearby Star-forming Regions, we are monitoring monthly with SCUBA-2 the submillimeter emission in eight fields within nearby (&lt;500 pc) star-forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 deg2 includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125–0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of submillimeter observations to reach a precision of ~2%–3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any submillimeter variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for submillimeter variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars

The JCMT Gould Belt Survey: First results from SCUBA-2 observations of the Cepheus Flare region

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that if modelled as thermally supported Bonnor–Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival 13CO velocity dispersion measurements and find that our cores are typically pressure confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.ISSN:0035-8711ISSN:1365-2966ISSN:1365-871

The JCMT Gould Belt Survey : a first look at SCUBA-2 observations of the Lupus I molecular cloud

This paper presents observations of the Lupus I molecular cloud at 450 and 850 μm with Submillimetre Common User Bolometer Array (SCUBA-2) as part of the James Clerk Maxwell Telescope Gould Belt Survey (JCMT GBS). Nine compact sources, assumed to be the discs of young stellar objects (YSOs), 12 extended protostellar, pre-stellar and starless cores, and one isolated, low-luminosity protostar, are detected in the region. Spectral energy distributions, including submillimetre fluxes, are produced for 15 YSOs, and each is fitted with the models of Robitaille et al. The proportion of Class 0/I protostars is higher than that seen in other Gould Belt regions such as Ophiuchus and Serpens. Circumstellar disc masses are calculated for more evolved sources, while protostellar envelope masses are calculated for protostars. Up to four very low luminosity objects are found; a large fraction when compared to other Spitzer c2d regions. One YSO has a disc mass greater than the minimum mass solar nebula. 12 starless/protostellar cores are detected by SCUBA-2 and their masses are calculated. The stability of these cores is examined using both the thermal Jeans mass and a turbulent virial mass when possible. Two cores in Lupus I are super-Jeans and contain no known YSOs. One of these cores has a virial parameter of 1.1 ± 0.4, and could therefore be pre-stellar. The high ratio of Class 0/I to Class III YSOs (1:1), and the presence of a pre-stellar core candidate, provides support for the hypothesis that a shock recently triggered star formation in Lupus I