Most massive stars form in dense clusters where gravitational interactions with other. stars may be common. The two nearest forming massive stars, the BN object and Source I, located behind the Orion Nebula, were ejected with velocities of similar to 29 and similar to 13 km s(-1) about 500 years ago by such interactions. This event generated an explosion in the gas. New ALMA observations show in unprecedented detail, a roughly spherically symmetric distribution of over a hundred (CO)-C-12 J = 2-1 streamers with velocities extending from V-LSR = -150 to +145 km s(-1) The streamer radial velocities increase (or decrease) linearly with projected distance from the explosion center, forming a '' Hubble Flow '' confined to within 50 ''. of the explosion center. They point toward the high proper-motion, shock-excited H-2 and [Fe II] '' fingertips '' and lower-velocity CO in the H-2 wakes comprising Orion's '' fingers.'' In some directions, the H-2 '' fingers '' extend more than a factor of two farther from the ejection center than the CO streamers. Such deviations from spherical symmetry may be caused by ejecta running into dense gas or the dynamics of the N-body interaction that ejected the stars and produced the explosion. This similar to 10(48) erg event may have been powered by the release of gravitational potential energy associated with the formation of a compact binary or a protostellar merger. Orion may be the prototype for a new class of stellar explosiozn responsible for luminous infrared transients in nearby galaxies.

In this thesis a critical assessment of the ages derived using theoretical pre-main-sequence (pre-MS) stellar evolutionary models is presented by comparing the predictions to the low-mass pre-MS population of 14 young star-forming regions (SFRs) in colour-magnitude diagrams (CMDs). Deriving pre-MS ages requires precise distances and estimates of the reddening. Therefore, the main-sequence (MS) members of the SFRs have been used to derive a self-consistent set of statistically robust ages, distances and reddenings with associated uncertainties using a maximum-likelihood fitting statistic and MS evolutionary models. A photometric method (known as the Q-method) for de-reddening individual stars in regions where the extinction is spatially variable has been updated and is presented. The effects of both the model dependency and the SFR composition on these derived parameters are also discussed. The problem of calibrating photometric observations of red pre-MS stars is examined and it is shown that using observations of MS stars to transform the data into a standard photometric system can introduce significant errors in the position of the pre-MS locus in CMD space. Hence, it is crucial that precise photometric studies (especially of pre- MS objects) be carried out in the natural photometric system of the observations. This therefore requires a robust model of the system responses for the instrument used, and thus the calculated responses for the Wide-Field Camera on the Isaac Newton Telescope are presented. These system responses have been tested using standard star observations and have been shown to be a good representation of the photometric system. A benchmark test for the pre-MS evolutionary models is performed by comparing them to a set of well-calibrated CMDs of the Pleiades in the wavelength regime 0.4−2.5 μm. The masses predicted by these models are also tested against dynamical masses using a sample of MS binaries by calculating the system magnitude in a given photometric band- pass. This analysis shows that for Teff ≤ 4000 K the models systematically overestimate the flux by a factor of 2 at 0.5 μm, though this decreases with wavelength, becoming negligible at 2.2 μm. Thus before the pre-MS models are used to derive ages, a recalibration of the models is performed by incorporating an empirical colour-Teff relation and bolometric corrections based on the Ks-band luminosity of Pleiades members, with theoretical corrections for the dependence on the surface gravity (log g). The recalibrated pre-MS model isochrones are used to derive ages from the pre-MS populations of the SFRs. These ages are then compared with the MS derivations, thus providing a powerful diagnostic tool with which to discriminate between the different pre- MS age scales that arise from a much stronger model dependency in the pre-MS regime. The revised ages assigned to each of the 14 SFRs are up to a factor two older than previous derivations, a result with wide-ranging implications, including that circumstellar discs survive longer and that the average Class II lifetime is greater than currently believed.

Star and stellar cluster formation in spiral galaxies is one of the biggest questions of astrophysics. In this thesis, I study how star formation, and the formation of stellar clusters, proceeds using SPH simulations. These simulations model a region of 400 pc and 107 solar masses. Star formation is modelled through the use of sink particles which represent small groups of stars. Star formation occurs in high density regions, created by galactic spiral arm passage. The spiral shock compresses the gas and generates high density regions. Once these regions attain sufficiently high density, self-gravity becomes dominant and drives collapse and star formation. The regions fragment hierarchically, forming local small groups of stars. These fall together to form clusters, which grow through subsequent mergers and large scale gas infall. As the individual star formation occurs over large distances before forming a stellar cluster, this process can result in significant age spreads of 1-2 Myrs. One protocluster is found to fail to merge due to the large scale tidal forces from the nearby regions, and instead expands forming a dispersed population of young stars such as an OB association.

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.