Search

Observational studies of star formation reveal spatial distributions of Young Stellar Objects (YSOs) that are `snapshots' of an ongoing star formation process. Using methods from spatial statistics it is possible to test the likelihood that a given distribution process could produce the observed patterns of YSOs. I determine the sensitivity of the spatial statistical tests Diggle's G function (G), the `free-space' function (F), Ripley's K and O-ring for application to astrophysical data. To do this I applied each test to simulated data containing 2D Gaussian clusters projected on a random distribution of background stars. By varying the number of stars within the Gaussian cluster and the number of background stars I determined the ability of the tests to reject complete spatial randomness (CSR) with changing signal-to-noise. Ripley's K and O-ring were shown to be much more sensitive to Gaussian clusters than G and F. I then apply the O-ring test to determine if column density alone is sufficient to explain the locations of Class 0/I YSOs within Serpens South, Serpens Core, Ophiuchus, NGC1333 and IC348. Star formation is known to occur more readily where more raw materials are available, a relationship that is often expressed in the form of a 'Kennicutt--Schmidt' relation where the surface density of Young Stellar Objects (YSOs) is proportional to column density to some power, μ. Using the O-ring test as a summary statistic, confidence envelopes were produced for different values of μ from probability models made using the Herschel column density maps. The YSOs were tested against four distribution models: the best-estimate of μ for the region, μ = 0 (i.e. random) above a column density threshold and zero probability elsewhere, μ = 1, and the power-law that best represents the five regions as a collective, μ = 2.05 ± 0.20. Serpens South and NGC1333 rejected the μ = 2.05 model on small scales of ~ 0.15 pc which implies that small-scale interactions may be influencing their distribution. On scales above 0.15 pc, the positions of YSOs in all five regions can be well described using column density alone.

Radio continuum emission has been demonstrated to be an excellent tracer of recent star formation, yet current understanding remains empirical, and is primarily based upon integrated multi-frequency or resolved mono-chromatic radio continuum studies. In this thesis, we gain further insight into the relationship between the observed radio continuum emission and recent star formation through the use of spatially resolved, multi-frequency radio continuum observations, focusing our analysis on the relatively simple star-forming dwarf irregular galaxies, IC10 and NGC1569. High-resolution (∼1 pc) 1.5GHz and 5.0GHz e-MERLIN observations of IC10, combined with a multi-wavelength classification scheme, reveal that the majority of the compact star-formation products found within low mass, star forming galaxies are HII regions, with radio supernova and supernova remnants being observedmore rarely. By taking a census of these star-formation products, we derive useful lower limits on the instantaneous star-formation rate in this simple stellar system. We further demonstrate that high-resolution observations can be used to remove the contribution of contaminating background galaxies to integrated galaxy radio spectra, which can be a significant correction in very low mass systems. On larger spatial scales (∼100 pc), we developed a Bayesian Markov-Chain Monte- Carlo approach to separate the thermal and non-thermal radio continuum emission components on a resolved basis. We applied this procedure to VLA observations of NGC1569 that span a wide frequency range (1-34GHz). During our fitting procedure, we use Ha maps found in the literature to constrain the thermal radio continuumemission component. We find that NGC1569 exhibits a high average thermal fraction at 1GHz (∼25%) which is in line with recent integrated studies, and that the average non-thermal spectral index for the main disk is a = −0.53, which is consistent with the injection spectrum found in galactic supernova remnants. Taken together, these results indicate that NGC1569 has recently finished a star formation burst, which is in close agreement with the literature where this galaxy is classified as being in a post- starburst phase, and that supernova remnants are the predominant source of cosmic ray electrons in this system. We further demonstrate that the resulting maps from our separation procedure can be used to further constrain additional quantities, such as the line of sight reddening at arcsecond resolution, providing for example a measure of internal extinction by dust in NGC1569's HII regions as well as spatially resolved equipartition magnetic field strengths. With our low-resolution VLA observations, we found that the supernova remnant, NGC1569-38, exhibits a break in its radio spectral energy distribution at ∼ 7GHz. We followed this observation up with e-MERLIN observations of NGC1569 and found that the surrounding InterstellarMediumsignificantly contributes to the observed spectral energy distribution in the lower resolution maps. Based on NGC1569-38's observed surface brightness and minimum energy magnetic field strength compared to a sample of galactic and extra-galactic supernova remnants,we determine thatNGC1569- 38 is young and is just entering the Sedov-Taylor phase of expansion. We demonstrate that the observed break in a spatially resolved radio continuum spectral energy distribution of a young supernova Remnant, even at resolutions of about 100 pc, is due to contamination by emission from the surrounding Interstellar medium. To reveal effects such as cosmic ray electron ageing, higher resolution, matched to the size of the supernova remnant, multi-frequency radio spectra are required.

Some of the most fundamental measurements we can make of the Universe are where and when stars formed in galaxies. In recent years, astronomers have converged on a picture in which the star formation rate density of the Universe peaks at approximately redshift (z) 2, when the Universe is around a quarter of its present age. There, star-forming galaxies harbour large reservoirs of molecular gas, assemble stellar mass rapidly, and typically display disordered morphologies. In this thesis, I study the evolution of galaxies on large and small scales, with a particular focus on the epoch around the peak of cosmic star formation. My overarching aim is to understand the physical processes that drive and quench star formation in galaxies over cosmic time. In the first half of this thesis, I focus on global measurements of star formation, using the High-z Emission line survey (HiZELS), a deep, near-infrared narrow-band survey, which identifies star-forming galaxies at z=0.8-2.2. I characterise the dark matter halo environments of these galaxies via a clustering analysis, along with a Halo Occupation Distribution model fitting procedure, then study the relationships between host dark matter halo environment and galaxy properties. I show that the clustering strength and the host dark matter halo masses of the HiZELS galaxies increase linearly with H-alpha luminosity (and, by implication, star formation rate) at all three redshifts. The typical galaxies in our samples are star-forming centrals, residing in dark matter haloes of mass ~10^12M_sol. I find a remarkably tight redshift-independent relation between the H-alpha luminosity scaled by the characteristic luminosity, L_H-alpha/L_H-alpha*(z), and the host dark matter halo mass of central galaxies. Simple analytic modelling suggests that this is consistent with a model in which the dark matter halo environment is a strong driver of galaxy star formation rate and therefore of the evolution of the star formation rate density in the Universe. I investigate this further by distinguishing the stellar mass and star formation rate dependencies of the clustering of HiZELS galaxies. I compare my observational results to the predictions of a pioneering cosmological hydrodynamical simulation, the Virgo Consortium's Evolution and Assembly of GaLaxies and their Environments project, known as EAGLE. In the subsequent chapters of this thesis, I focus more heavily on simulations of galaxy formation, which are important tools for constraining and understanding the physics at play in galaxies. I use EAGLE to investigate the quenching of star formation in simulated galaxies via a novel application of Principal Component Analysis. I show that the key relations between halo mass, stellar mass and star formation rate are in good agreement with observed low-redshift galaxies. Having studied the global properties of star-forming galaxies, I then turn to smaller scales, investigating what we can learn from spatially-resolved imaging. I present a detailed study of the spatially-resolved dust continuum emission of realistic simulated high-redshift galaxies. These galaxies, drawn from the FIRE-2 simulations, reach Milky Way masses by z~2. Post-processing them using radiative transfer techniques, I obtain predictions for the full rest-frame far-ultraviolet to far-infrared Spectral Energy Distributions of these simulated galaxies, as well as maps of their emission across the wavelength spectrum. As has been observed in distant galaxies, the rest-frame far-infrared emission of the simulated galaxies is compact, spanning half-light radii of ~0.5-4kpc. The derived morphologies of simulated galaxies are notably different in different wavebands; a galaxy can appear clumpy and extended in the far-ultraviolet yet compact at far-infrared wavelengths. Finally, I perform a multi-wavelength study of a single observed galaxy, SHiZELS-14 (z=2.24), drawn from the HiZELS survey and subsequently imaged at 0.15'' resolution at multiple wavelengths. The data comprise kpc-resolution imaging in three different widely used tracers of star formation: the H-alpha emission line (from SINFONI/VLT), rest-frame far-ultraviolet continuum (from HST F606W imaging), and the rest-frame far-infrared (from ALMA), as well as the rest-frame optical (from HST F140W imaging). SHiZELS-14 displays a compact, dusty centre, as well as extended emission in both H-alpha and the rest-frame FIR. The ultraviolet emission is spatially offset from the extended dust emission, and appears to trace holes in the dust distribution. I find that the dust attenuation varies across the spatial extent of the galaxy, reaching up to ~5 magnitudes of extinction at H-alpha wavelengths in the most dusty regions. Global star formation rates inferred using standard calibrations to the different tracers vary from ~10-1000M_sol, and are particularly discrepant in the galaxy's dusty centre. This galaxy highlights the biased view of galaxy evolution provided by short-wavelength data in the absence of long-wavelength data, and is in line with my simulations.

This dissertation explores cosmological simulations of galaxy formation, focusing on the asymptotic efficiency of conversion of baryons into stars and asking how this quantity varies in different cosmological models. By exploring the behaviour into the future and in counter-factual universes with altered cosmological parameters, I hope to gain insight into the robustness of galaxy formation codes for their behaviour in different parameter regimes, addressing the concern that the models may be fine-tuned to fit observations, rather than being genuinely predictive. Through the calibration of multiple star formation and feedback parameters within zoom simulations using the Enzo code, I obtain a set of parameter values able to reproduce the observational constraint of baryon makeup of haloes of masses between 1010 and 1012M⊙. Comparing two different star formation setups, I show that feedback is self-consistent with a higher feedback energy injection associated with a low star formation efficiency, and vice versa. I also explore the reproducibility of the simulation results and conclude that operational differences in the implementation of the numerical code can create approximately 10% - 25% deviation in the stellar and baryon mass respectively. Using the best star formation and feedback prescription from the zoom simulations, I attempt a cosmological box simulation in a standard cold dark matter cosmology beyond z = 0. Comparing with previous work, I find a similar elongation in the general evolution of the distribution of temperature and density of the gas into the future. I then extend the analysis to the evolution of the halo mass function, the equation of state of the intergalactic medium and star formation rate into the future to determine the level of convergence with various resolutions. Interestingly, I identify a turn-around in the cosmic star formation, deviating from the extrapolation obtained from observations. There is also a cross over in the fraction of baryon in the form of stars and in gas around the period of 'freeze out' in both the zoom and box simulations. Lastly, I present results from simulations of counter-factual simulations, exploring the effects of different values of Λ on structure formation and evolution. These simulations allow the study of anthropic explanations for the small observed value of Λ, considering the relative frequency of observers within a multiverse. It appears that a distinct peak in the star formation rate may be due to the presence of Λ. For a higher value of Λ, the peak in star formation is of a lower value, and it occurs earlier in the evolution of the universe. Given these differences in the star formation rate, the UV background must be calculated self consistently in these models. This corrected UV background will provide a more realistic representation of the counter-factual universe. I show that almost 54% of possible observers reside in universes with a value of Λ similar to ours, implying that anthropic principle can explain why Λ is of such a small value.

The following dissertation details the work undertaken to compile the most complete compact and ultracompact H II region catalogue to date via multi-wavelength inspection of survey data. We utilise data from the recently available SASSy 850 μm survey and explain the data reduction process leading to the final Outer Galaxy Survey. Combined with the earlier SASSy component that focuses on the Perseus arm, we were able to use the catalogue to identify massive star forming clumps in the Outer Galaxy (RGC > 8.5 kpc) and cross-match with infrared and radio data of known UC H II regions from the RMS database. We sought to compare this sample with previous existing catalogues in the Inner Galaxy to examine massive star formation as a function of galactocentric radius and particularly focusing on two distinct environmental regions (Inner and Outer Galaxy) with overall different average metallicities and gas-to-dust ratios. For this Inner Galaxy sample (RGC < 8.5 kpc), we adopted the compact H II regions from works that used similar methods to cross-match ATLASGAL with either CORNISH or RMS, depending on the location within the Galactic plane. We present a new UC H II region catalogue that more than doubles the original sample size of previous work, totalling 536 embedded H II regions and 445 host clumps. We examine the distance independent values of NLy/M and Lbol/M as proxies for massive star formation efficiency and overall star formation efficiency, respectively. We find a significant trend showing that Lbol/M decreases with increasing RGC, suggesting that the overall star formation per unit mass of the molecular clump is less in the Outer Galaxy. Further investigations with the VLA into the Outer Galaxy for 200 unmatched SASSy molecular clumps revealed seven additional 5 GHz continuum emission sources embedded within the clumps. We conclude that massive star formation is not dependent on Galactic position or environmental factors such as metallicity and the gas-to-dust ratio that are known to vary with galactocentric radius. These results will be relevant to future models attempting to use the Milky Way as a template for evolutionary and formation models of nearby galaxies where similar conditions are found but not observable to the degree of detail as seen in our own Galaxy.

Remarkable progress has been made over the last few decades in furthering our understanding of the growth of cosmic structure. Nonetheless, there remains a great deal of uncertainty regarding the precise details of the complex baryonic physics that regulate galaxy formation. Any theory of star formation in galaxies must encompass the radiative cooling of gas into dark matter haloes, the formation of a turbulent, multiphase interstellar medium (ISM), the efficiency with which molecular gas is able to collapse into cores and ultimately stars, and the subsequent interaction of those stars with the gas through ionizing radiation, winds and supernova (SN) explosions. Given the highly non-linear nature of the problem, numerical simulations provide an invaluable tool with which to study galaxy formation. Yet, even with contemporary computational resources, the inherently large dynamical range of spatial scales that must be tackled makes the development of such models extremely challenging, inevitably leading to the adoption of `subgrid' approximations at some scale. In this thesis, I explore new methods of incorporating the physics of star formation and stellar feedback into high resolution hydrodynamic simulations of galaxies. I first describe a new implementation of star formation and SN feedback that I have developed for the state-of-the-art moving mesh code Arepo. I carry out a detailed study into various classes of subgrid SN feedback schemes commonly adopted in the literature, including injections of thermal and/or kinetic energy, two parametrizations of delayed cooling feedback and a 'mechanical' feedback scheme that injects the appropriate amount of momentum depending on the relevant scale of the SN remnant (SNR) resolved. All schemes make use of individually time-resolved SN events. Adopting isolated disk galaxy setups at different resolutions, with the highest resolution runs reasonably resolving the Sedov-Taylor phase of the SNR, I demonstrate that the mechanical scheme is the only physically well-posed method of those examined, is efficient at suppressing star formation, agrees well with the Kennicutt-Schmidt relation and leads to converged star formation rates and galaxy morphologies with increasing resolution without fine tuning any parameters. However, I find that it is difficult to produce outflows with high enough mass loading factors at all but the highest resolution. I discuss the various possible solutions to this effect, including improved modelling of star formation. Moving on to a more self-consistent setup, I carry out a suite of cosmological zoom-in simulations of low mass haloes at very high resolution, performed to z = 4, to investigate the ability of SN feedback models to produce realistic galaxies. The haloes are selected in a variety of environments, ranging from voids to crowded locations. In the majority of cases, SN feedback alone has little impact at early times even in low mass haloes ($\sim10^{10}\,\mathrm{M_\odot}$ at z = 0). This appears to be due largely to the build up of very dense gas prior to SN events, suggesting that other mechanisms (such as other stellar feedback processes) are required to regulate ISM properties before SNe occur. The effectiveness of the feedback also appears to be strongly dependent on the merger history of the halo. Finally, I describe a new scheme to drive turbulence in isolated galaxy setups. The turbulent structure of the ISM very likely regulates star formation efficiencies on small scales, as well as affecting the clustering of SNe. The large range of potential drivers of ISM turbulence are not fully understood and are, in any case, unlikely to arise ab initio in a whole galaxy simulation. I therefore neglect these details and adopt a highly idealised approach, artificially driving turbulence to produce an ISM structure of my choice. This enables me to study the effects of a given level of ISM turbulence on global galaxy properties, such as the fragmentation scale of the disk and the impact on SN feedback efficiencies. I demonstrate this technique in the context of simulations of isolated dwarfs, finding that moderate levels of turbulent driving in combination with SN feedback can produce a steady-state of star formation rates and global galaxy properties, rather than the extremely violent SN feedback that is produced by a rapidly fragmenting disk.

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.

Previous studies have shown that thermal feedback from protostars plays a key role in the process of low-mass star formation. In this thesis, we explore the effects of protostellar heating on the formation of stellar clusters. We describe new methods for modelling protostellar accretion luminosities and protostellar evolution in calculations of star formation. We then present results of a series of numerical simulations of stellar cluster formation which include these effects, and examine their impact. We begin by investigating the dependence of stellar properties on the initial density of molecular clouds. We find that the dependence of the median stellar mass on the initial density of the cloud is weaker than the dependence of the thermal Jeans mass when radiative effects are included. We suggest that including protostellar accretion luminosities and protostellar evolution may weaken this dependence further, and may account for the observed invariance of the median stellar mass in Galactic star-forming regions. Next, we investigate the effects of including accretion feedback from sink particles on the formation of small stellar groups. We find that including accretion feedback in calculations suppresses fragmentation even further than calculations that only include radiative transfer within the gas. Including feedback also produces a higher median stellar mass, which is insensitive to the sink particle accretion radius used. Finally, we compare calculations of small stellar clusters which model the evolution of protostars using a live stellar model with those which use a fixed stellar structure. We find that the dynamics of the clusters are primarily determined by the accretion luminosities of protostars, but that the relative effects of protostellar evolution depend on the accretion rate and advection of energy into the protostar. We also demonstrate how such calculations may be used to study the properties of young stellar populations.

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.

Stars form inside clouds of molecular gas and dust. In the early stages of stellar evolution the remainders of the initial cloud form a circumstellar disk. For the next few million years the disk will slowly dissipate via accretion, outflows, photoevaporation and planet growth while the star makes its way onto the Main Sequence. This stage of a star's life is referred to as the T Tauri phase and is characterised by high-level spectrophotometric variability. This thesis aims to study and map out the environments of T Tauri stars down to the very low mass regime by the means of time-domain monitoring. Different physical processes in the system manifest themselves as variability on different time- scales as well as produce characteristic spectroscopic and photometric features at various wave- lengths. In order to study young stellar objects in depth, the observing campaigns presented in this work were designed to cover a large range of time-scales - minutes, hours, days and months. Combining all the data, this thesis establishes a baseline of over a decade for some objects. The observations also cover a wide range of wavelengths from the optical to the mid-infrared part of the spectrum. The star RW Aur experienced two long-lasting dimming events in 2010 and 2014. This thesis presents a large collection of spectral and photometric measurements carried out just before and during the 2014 event. Spectral accretion signatures indicate no change in the accretion activity of the system. Photometry indicates that parallel to the dimming in the optical the star becomes brighter in the mid-infrared. The observations in this work combined with literature data suggest that the origin of the 2014 event is most likely obscuration of the star by hot dust from the disk being lifted into the disk wind. Very low mass stars (< 0.4 M⊙) are the most common type of star in the Galaxy. In order to understand the early stages of stellar evolution we must study young very low mass stars. This work investigates the photometric and spectroscopic variability of seven brown dwarfs in star forming regions near σ Ori and ε Ori. All targets exhibit optical photometric variability between from 0.1 to over 1.0 magnitude that persists on a time-scale of at least one decade. Despite the photometric variability no change in the spectral type is measured. In the cases where the stars are accreting, modelling of the spectral changes suggest the accretion flow is more homogeneous and less funnelled compared to Sun-like T Tauri stars. The non-accreting variables are more plausibly explained by obscuration by circumstellar material, possibly a ring made out of multiple clouds of dust grains and pebbles with varying optical depths. The star-disk systems studied in this thesis have some broader implications for star and planet formation theory. The case-study of RW Aur has unambiguously demonstrated that the planet- forming environment is very dynamic and can change dramatically on short time-scales, which in turn would have implications for the diversity of planetary systems found in the Galaxy. The Orion stars have shown that the current theory for the T Tauri stage of stellar evolution is valid down to the very low mass regime. The seven dwarfs are a good example for the evolutionary path of circumstellar disks, showing the transition from gas-high, flared accretion disks (σ Ori) to dust-dominated, depleted, structured debris disks (ε Ori).

I study mid-infrared and star formation properties of AGN samples using infrared observations, and star-forming galaxies using radio observations in order to investigate the link between star formation, AGN activity and radio luminosity. I present the results of these investigations in this thesis. I carried out an analysis of four complete samples of radio-loud AGN (3CRR, 2Jy, 6CE and 7CE) using near- and mid-IR data taken by the Wide-Field Infrared Survey Explorer (WISE). The combined sample consists of quasars and radio galaxies, and covers a redshift range 0:003 < z < 3:395. The dichotomy in the mid-IR properties of low- and high-excitation radio galaxies (LERGs - HERGs) is analysed using large complete samples. The results show that a division in the accretion modes of powerful LERGs and HERGs clearly stands out in the mid-IR radio plane. Evaluation of the positions of the sample objects in WISE colour-colour diagrams shows that widely used WISE colour cuts are not completely reliable in selecting AGN. I examined the link between AGN activity and star formation by constructing matched samples of local (0 < z < 0:6) radio-loud and radio-quiet AGN in the Herschel-ATLAS fields. AGN accretion and jet powers in these active galaxies are traced by [OIII] emission-line and radio luminosity, respectively. Star formation properties were derived using Herschel 250-_m and stellar mass measurements are taken from the SDSS-MPA-JHU catalogue. The stacking analyses show that star formation rates (SFRs) and specific star formation rate (SSFRs) of both radio-loud and radio-quiet AGN increase with increasing AGN power but that radio-loud AGN tend to have lower SFR. Additionally, radio-quiet AGN are found to have approximately an order of magnitude higher SSFRs than radio-loud AGN for a given level of AGN power. The difference between the star formation properties of radio-loud and -quiet AGN is also seen in samples matched in stellar mass. I also investigated the relationship between SFR and low-frequency radio luminosity observed in star-forming galaxies. I used a sample of star-forming galaxies in the 19 local Universe selected from the SDSS-MPA-JHU catalogue. LOFAR observations of the Herschel-ATLAS North Galactic Pole field (NGP) were carried out as part of the LOFAR surveys Key Science Project at an effective frequency of 150 MHz, which provided low-frequency radio luminosity of sample galaxies. SFRs of galaxies in the sample were derived using MAGPHYS spectral energy distribution (SED) fitting. The results of this study show that the slope of L150/SFR is less than unity and not universal for all star-forming galaxies (SFGs) in the local Universe (0 < z < 0:3). The slope of the L150/SFR relation is also found to be steeper than the L1:4/SFR relation, probably due to the contribution from thermal radio emission at 1.4 GHz. If the L150=SFR relation for strongly star-forming objects is explained naively by electron calorimetry, I conclude that low luminosity sources are not ideal calorimeters and differ from the main locus of SFGs at low redshifts. The different gradients we obtain for the far- IR/radio correlation using samples selected at different frequencies reveal the selection effects on relations derived in this thesis.

The young moving groups are collections of nearby (<200 pc), young (5-150 Myr) pre-main sequence stars; these stars offer us one of the best opportunities to characterise stellar multiplicity, sub-stellar phenomena, disc evolution and planet formation. Here we present results from a series of multiplicity studies aimed at producing comprehensive multiplicity statistics of the young moving groups. The aim was to compare the derived statistics of the young moving groups to other populations in order to investigate whether the abundance and properties of multiple systems are environment independent. We have combined high-resolution spectroscopy, AO-imaging and direct imaging to identify and characterise multiple systems across a huge range of orbital periods (1- 10e10 day). The observational techniques also allow us to constrain the abundance of multiple systems in these populations by calculating detection limits. We found many similarities (frequency of spectroscopic binaries; frequency, mass-ratio and physical separation of visual binaries) between the young moving groups and both younger and older regions, for multiple systems with physical separations smaller than 1000 au. We did, however, identify a significant number of new wide (>1000 au) companions. We reconciled the apparent excess of wide binary systems, when compared to the field population, by arguing that the wide systems are weakly bound and most likely decaying. By comparing the multiplicity statistics in one particular moving group we showed that the dynamical evolution of non-hierarchical protostars could lead to the population of wide binaries we can observe today. Our results indicate that the majority of low-mass stars form in small groups with 3 or 4 components that undergo significant dynamical evolution. The multiplicity properties of the young nearby moving groups are statistically similar to many other populations, supporting the environment-independent formation of multiple systems.

Active galaxies have been in the forefront of astronomic research since their first discovery, at least 50 years ago (e.g. Schmidt, 1963; Matthews & Sandage, 1963). The putative supermassive black hole (SMBH) at their center characterizes their properties and regulates the evolution of these objects. In this thesis, I study the 'demographics' and 'ecology' of active galactic nuclei (AGN) in the context of their evolution and the interaction with their environments (mainly their host galaxy). The number density of AGN has been found to peak at 1 < z < 3 (e.g. Ueda et al., 2003; Hasinger et al., 2005; Richards et al., 2005; Aird et al., 2010), similar to the star formation history (e.g. Silverman et al., 2008a; Aird et al., 2010). However, when taking into account obscuration, faint AGN are found to peak at lower redshift (z ≤ 2) than that of bright AGN (z ≈ 2 - 3; e.g. Hasinger et al., 2005; Hopkins et al., 2007; Xue et al., 2011). This qualitative behaviour is also broadly seen in star-forming galaxies (e.g. Cowie et al., 1996) and is often referred to as 'cosmic downsizing', although this term has developed a number of usages with respect to galaxies (e.g. Bundy et al., 2006; Cimatti et al., 2006; Faber et al., 2007; Fontanot et al., 2009). Though this behaviour is well established up to z ≈ 3, the nature of how and when the initial seed of these AGNs were formed remains an open question. For this study, I use Chandra surveys to study some of the most distant AGN in the Universe (z > 3). The combination of two different size and depth Chandra surveys (Chandra-COSMOS and ChaMP) provides me with the largest to-date z > 3 AGN sample, over a wide range of rest-frame 2-10 keV luminosities [log (Lₓ/erg s⁻¹) = 43.3-46.0] and obscuration (NH = 10²⁰ - 10²³ cm⁻²). I find strong evidence about a strong decline in number density of X-ray AGN above z ≈ 3, and also the association of this decline with a luminosity-dependent density evolution (LDDE; e.g. Gilli et al., 2007). Especially at high redshifts, the different evolution models predict quite different numbers of AGNs. The large size and the wide X-ray luminosity range of this sample reduces the uncertainties of previous studies at similar redshifts making it possible to distinguish between the different models and suggest that observations appear to favour the LDDE model. The observed AGN downsizing behaviour seen via the measured X-ray luminosity function (XLF) could arise due to changes in the mass of the typical active SMBH and/or changes in the typical accretion rate. But how does the growth of SMBHs over cosmic time influence its environment? A powerful way to address this question is to compare the host galaxy properties over a wide range of AGN and accretion rate types. Radio-jets are one of the most prominent constituents of AGN as they can interact directly with the host galaxy. Although AGN with radio jets are rare (they make up to 10 per cent of the total AGN population) radio galaxies make up over 30 per cent of the massive galaxy population and it is likely that all massive galaxies go through a radio-loud phase, as the activity is expected to be cyclical (e.g Best et al., 2005). It is therefore, important to investigate the impact of radio jets on the host galaxy and particularly the star formation. The method I follow focuses on the comparison of the host galaxy properties between optically selected quasar samples, with and without strong radio emission associated with powerful radio-jets, matched in AGN luminosity. Herschel far-infrared observations are used to trace the star formation in the host galaxy, providing minimal AGN contamination. In my first approach, I have constructed a sample of radio-loud and radio-quiet quasars from the Faint Images Radio Sky at Twenty-one centimetres (FIRST) and the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), over the H-ATLAS Phase 1 Area (9h, 12h and 14.5h). The main result of this work is that RLQs at lower AGN luminosities tend to have on average higher FIR and 250-μm luminosity with respect to RQQs matched in AGN luminosity and redshift. However, evolution effects could be strong as the quasars in this sample cover a wide range of redshifts (0.4 < z < 5). Therefore, I follow a second approach with the advantage of a QSO sample selection at a single redshift epoch, decomposing the evolution effects from the AGN/star-formation study. The results indicate that radio-jets in powerful QSOs can both suppress and enhance the star formation in their host galaxies. These fundings are consistent with a galaxy mass and jet-power dependence model. Then we expect more massive galaxies to have more star-formation for a given jet-power because their star-formation is more enhanced by the jet. Although radio-jets are the best candidates for a direct AGN impact to the host galaxy, many models refer to an AGN feedback associated with energetic AGN winds and outflows which are expected to suppress the star formation in powerful AGN when compared to the overall galaxy population. My results do not suggest star formation is suppressed in the hosts of optically selected QSOs at z ≈ 1, with more than 30 per cent of them being associated with strong star formation rates (SFR ≈ 350 M⊙ yr⁻¹). Although different interpretations are possible, this result can be explained through periods of enhanced AGN activity and star-forming bursts, possibly through major mergers. However, optical QSOs comprise only a small fraction of the total AGN population. Even if the 'unified model' predicts that the host galaxy properties should not be affected by the viewing angle (type-1 vs. type-2 AGN), several studies have shown results supporting a scenario departing from the basic model. Investigating star formation in the hosts of 24 μm selected type-1&2 AGN, I found that the type-2 AGNs display on average higher star-formation rate than type-1 AGNs. This result is in agreement with previous studies suggesting an undergoing transition between a hidden growth phase and an unobscured AGN phase.

Stars form via the gravitational collapse of molecular clouds during which time the protostellar object contracts by over seven orders of magnitude. If all the angular momentum present in the natal cloud was conserved during collapse, stars would approach rotational velocities rapid enough to tear themselves apart within just a few Myr. In contrast to this, observations of pre-main sequence rotation rates are relatively slow (∼ 1 − 15 days) indicating that significant quantities of angular momentum must be removed from the star. I use observations of fully convective pre-main sequence stars in two well-studied, nearby regions of star formation (namely the Orion Nebula Cluster and Taurus-Auriga) to determine the removal rate of stellar angular momentum. I find the accretion disc-hosting stars to be rotating at a slower rate and contain less specific angular momentum than the disc-less stars. I interpret this as indicating a period of accretion disc-regulated angular momentum evolution followed by near-constant rotational evolution following disc dispersal. Furthermore, assuming that the age spread inferred from the Hertzsprung-Russell diagram constructed for the star forming region is real, I find that the removal rate of angular momentum during the accretion-disc hosting phase to be more rapid than that expected from simple disc-locking theory whereby contraction occurs at a fixed rotation period. This indicates a more efficient process of angular momentum removal must operate, most likely in the form of an accretion-driven stellar wind or outflow emanating from the star-disc interaction. The initial circumstellar envelope that surrounds a protostellar object during the earliest stages of star formation is rotationally flattened into a disc as the star contracts. An effective viscosity, present within the disc, enables the disc to evolve: mass accretes inwards through the disc and onto the star while momentum migrates outwards, forcing the outer regions of the disc to expand. I used spatially resolved submillimetre detections of the dust and gas components of protoplanetary discs, gathered from the literature, to measure the radial extent of discs around low-mass pre-main sequence stars of ∼ 1−10 Myr and probe their viscous evolution. I find no clear observational evidence for the radial expansion of the dust component. However, I find tentative evidence for the expansion ofthe gas component. This suggests that the evolution of the gas and dust components of protoplanetary discs are likely governed by different astrophysical processes. Observations of jets and outflows emanating from protostars and pre-main sequence stars highlight that it may also be possible to remove angular momentum from the circumstellar material. Using the sample of spatially resolved protoplanetary discs, I find no evidence for angular momentum removal during disc evolution. I also use the spatially resolved debris discs from the Submillimetre Common-User Bolometer Array-2 Observations of Nearby Stars survey to constrain the amount of angular momentum retained within planetary systems. This sample is compared to the protoplanetary disc angular momenta and to the angular momentum contained within pre-stellar cores. I find that significant quantities of angular momentum must be removed during disc formation and disc dispersal. This likely occurs via magnetic braking during the formation of the disc, via the launching of a disc or photo-evaporative wind, and/or via ejection of planetary material following dynamical interactions.

In this thesis I present work done on the subject of star formation through the study of infrared dark clouds. We studied the velocity fields in several IRDCs using spectral line mapping. We also performed observations of a high density tracer and its deuterated counterpart. These observations allow me to assess the kinematics of these clouds and the evolutionary state of the observed targets. The sample observed is an important starting point for the search of early and quiescent high-mass regions. We mapped several IRDCs using the 22m ATNF Mopra Telescope in high-density molecular tracers at 3 mm, HNC (1-0) and N2H+ (1-0). We present integrated intensity emission and velocity field maps of these IRDCs. The molecular emission in the maps matches well with the extinction seen in the mid-IR. For an IRDC-complex we see connecting emission in the whole filament. We calculate kinematic distances and masses of the IRDCs. The IRDCs typically display an ordered velocity field within the clouds. The mean velocity gradient of the sample was 0.4 km/s/pc. We show how this velocity gradient can mean gas flows within the cloud into the central regions in order to feed the central cores. We observed 54 cores in IRDCs using N2H+ (1-0) and (3-2) to determine the kinematics of the densest material, where stars will form. We also observed N2D+ (3-2) towards 29 of the brightest peaks to analyze the level of deuteration which is an excellent probe of the quiescent of the early stages of star formation. There were 13 detections of N2D+ (3-2). This is one of the largest samples of IRDCs yet observed in these species. The deuteration ratio in the sources with detected N2D+ (3-2) has a mean of 0.024 and reaches a maximum value of 0.14. For most of the sources the material traced by N2D+ and N2H+ (3-2) still has significant turbulent motions, however three objects show subthermal N2D+ velocity dispersion. Surprisingly the presence or absence of an embedded 70 micron source shows no correlation with the detection of N2D+ (3-2), nor does it correlate with any change in velocity dispersion or excitation temperature. Comparison with recent models of deuteration suggest evolutionary timescales of these regions of several freefall times or less.

This thesis presents a detailed study of the impact of radio-AGN on star formation and the interstellar medium (ISM) of galaxies across cosmic time. To do this, this thesis uses far-IR/sub-mm data from the Herschel Space Observatory. I create a well-selected sample of 1599 radio sources using the NRAO VLA Sky Survey (NVSS) and Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) data in combination with the UKIRT Infrared Deep Sky Survey - Large Area Survey (UKIDSS - LAS) data. I find that the far-IR luminosities and dust temperatures of radio galaxies are lower in comparison to those of non-radio-detected galaxies. This luminosity deficit grows with increasing stellar mass. I argue that the reasons for these differences is probably due to indirect radio-AGN feedback, i.e. radio jets mechanically heat the halo-environment, preventing external sources of cold gas from entering the host and forming stars. The far-IR luminosity and dust temperature is found to decrease as a function of radio source size. I find the most likely explanation for this is jet-induced star formation while the jets are confined to the ISM. Finally, a method for identifying reliable high-z, star-bursting radio sources in the H-ATLAS is described with which statistically significant studies of radio-jet induced star formation may be undertaken.

I have used blue near-infrared colours to select a group of UKIDSS T dwarfs with spectral types later than T4. From amongst these I identify two kinematic halo T-dwarf candi- dates. Blue near-infrared colours have been attributed to collisionally-induced hydrogen absorption, which is enhanced by either high surface gravity or low metallicity. Proper motions are measured and distances estimated, allowing the determination of tangential velocities. U and V components are estimated for our objects by assuming Vrad = 0. From this, ULAS J0926+0835 is found to have U = 62 kms−1 and V = −140 kms−1 and ULAS J1319+1209 is found to have U = 192 kms−1 and V = −92 kms−1. These values are consistent with potential halo membership. However, surprisingly, these are not the bluest objects in the selection. The bluest is ULAS J1233+1219, with J −K = −1.16±0.07, and surprisingly this object is found to have thin disc-like U and V . Our sample also contains Hip 73786B, which I find to be a companion to the star Hip 73786. Hip 73786 is a metal- poor star, with [Fe/H]= −0.3 ± 0.1 and is located at a distance of 19±0.7 pc. U, V,W space velocity components are calculated for Hip 73786A and B, finding that U = −48±7 kms−1, V = −75 ± 4 kms−1 and W = −44 ± 8 kms−1. From the properties of the pri- mary, Hip 73786B is found to be at least 1.6Gyr old. As a metal poor object, Hip 73786B represents an important addition to the sample of known T dwarf benchmarks. Using mid-infrared data from WISE, I also identify T dwarfs with abnormally-red H − W2 and consider possible causes for their extreme colours. In particular I exam- ine three prominent examples of this phenomenon, ULAS J1416+1348B, 2MASS J0939- 2448 and BD+01o 2920B. A plot of spectral type against MW2-magnitude suggests that ULAS J1416+1348B is potentially an unresolved binary, similar to 2MASS J0939-2448. However, the plot also indicates that BD+01o 2920B is not an unresolved binary. I also present new FIRE spectroscopy for ULAS J1416+1348B and 2MASS J0939-2448. These data show that ULAS J1416+1348B has a similar shape to the Y -band spectrum to that of BD+01o 2920B, thus suggesting that the two objects have a similar metallicity, whereas 2MASS J0939-2448 appears to be a more metal-rich object. Using a new parallactic dis- tance, I derive a luminosity of (6.9±0.7)×1020W for ULAS J1416+1348B. I also find a radial velocity of −39 ± 1 kms−1 for this object. The agreement between this and that of the L dwarf SDSS J1416+1348A confirms that these two objects are physically-associated. I also present a set of simulated unresolved binaries; the colours of these systems do not appear to redden significantly with the addition of cooler companions. From this, I suggest that the colours of ULAS J1416+1348B and BD+01o 2920B cannot be solely attributed to any possible unresolved companions; for these two objects, composition and/or surface gravity must be playing a substantial role. Consideration of model predictions provides extra evidence for this argument, showing as it does that high log g and low metallicity can redden H − W2 colours by as much as »0.5mag as compared to a high-metallicity and low log g object of the same effective temperature. I also present kinematics and photometry for several more new candidate low-metallicity T dwarfs. Spectra are also presented, where available. In addition I provide new follow-up JHK spectroscopy for ULAS J0926+0835, ULAS J1233+1219 and ULAS J1319+1209. These new spectra allow full JHK-based spectral typing for these objects.

This thesis details the development of the radiation transport code torus for radiation hydrodynamic applications and its subsequent use in investigating problems regarding radiative feedback. The code couples Monte Carlo photoionization with grid-based hydrodynamics and has the advantage that all of the features available to a dedicated radiation transport code are at its disposal in RHD applications. I discuss the development of the code, including the hydrodynamics scheme, the adaptive mesh refinement (AMR) framework and the coupling of radiation transport with hydrodynamics. Extensive testing of the resulting code is also presented. The main application involves the study of radiatively driven implosion (RDI), a mechanism where the expanding ionized region about a massive star impacts nearby clumps, potentially triggering star formation. Firstly I investigate the way in which the radiation field is treated, isolating the relative impacts of polychromatic and diffuse field radiation on the evolution of radiation hydrodynamic RDI models. I also produce synthetic SEDs, radio, Hα and forbidden line images of the bright rimmed clouds (BRCs) resulting from the RDI models, on which I perform standard diagnostics that are used by observers to obtain the cloud conditions. I test the accuracy of the diagnostics and show that considering the pressure difference between the neutral cloud and surrounding ionized layer can be used to infer whether or not RDI is occurring. Finally I use more synthetic observations to investigate the accuracy of molecular line diagnostics and the nature of line profiles of BRCs. I show that the previously unexplained lack of dominant blue-asymmetry (a blue-asymmetry is the expected signature of a collapsing cloud) in the line profiles of BRCs can be explained by the shell of material, swept up by the expanding ionized region, that drives into the cloud. The work in this thesis combines to help resolve the difficulties in understanding radiative feedback, which is a non–linear process that happens on small astrophysical timescales, by improving numerical models and the way in which they are compared with observations.

Within this thesis, we present an extended multiwavelength analysis of the rich massive Galactic star-forming complex G305. We have focused our attention on studying the both the embedded massive star-forming population within G305, while also identifying the intermediate-, to lowmass content of the region also. Though massive stars play an important role in the shaping and evolution of their host galaxies, the physics of their formation still remains unclear. We have therefore set out to studying the nature of star formation within this complex, and also identify the impact that such a population has on the evolution of G305. We firstly present a Herschel far-infrared study towards G305, utilising PACS 70, 160 μm and SPIRE 250, 350, and 500 μm observations from the Hi-GAL survey of the Galactic plane. The focus of this study is to identify the embedded massive star-forming population within G305, by combining far-infrared data with radio continuum, H2O maser, methanolmaser,MIPS, and Red MSX Source survey data available from previous studies. From this sample we identify some 16 candidate associations are identified as embedded massive star-forming regions, and derive a two-selection colour criterion from this sample of log(F70/F500)! 1 and log(F160/F350)! 1.6 to identify an additional 31 embedded massive star candidates with no associated starformation tracers. Using this result, we are able to derive a star formation rate (SFR) of 0.01 - 0.02 M! yr−1. Comparing this resolved star formation rate, to extragalactic star formation rate tracers (based on the Kennicutt-Schmidt relation), we find the star formation activity is underestimated by a factor of !2 in comparison to the SFR derived from the YSO population. By next combining data available from 2MASS and VVV, Spitzer GLIMPSE and MIPSGAL, MSX, and Herschel Hi-GAL, we are able to identify the low-, to intermediate-mass YSOs present within the complex. Employing a series of stringent colour selection criteria and fitting reddened stellar atmosphere models, we are able remove a significant amount of contaminating sources from our sample, leaving us with a highly reliable sample of some 599 candidate YSOs. From this sample, we derive a present-day SFR of 0.005±0.001M! yr−1, and find the YSOmass function (YMF) of G305 to be significantly steeper than the standard Salpeter-Kroupa IMF. We find evidence of mass segregation towards G305, with a significant variation of the YMF both with the active star-forming region, and the outer region. The spatial distribution, and age gradient, of our 601 candidate YSOs also seem to rule out the scenario of propagating star formation within G305, with a more likely scenario of punctuated star formation over the lifetime of the complex.

This thesis presents molecular line and radio continuum observations of the giant molecular cloud (GMC) complex known as G305. The energy input from high-mass stars in the form of powerful winds and ionising radiation is one of the primary feedback mechanisms in GMCs. This feedback is thought to play a dual role both dispersing and destroying the natal environment but also sweeping up and compressing molecular gas and potentially triggering new episodes of star formation. Despite their importance to the evolution of GMCs and galaxies as a whole, the physical processes behind the formation and evolution of high-mass stars remains poorly understood. We therefore set out to obtain wide-field observations of the ionised and molecular environment to study the impact of high-mass stars on the evolution of G305. Observations conducted with the Mopra telescope of the molecular gas traced by NH3 in the (1,1), (2,2) and (3,3) transition and CO (12CO, 13CO and C18O J = 1–0) reveals the reservoir for future star formation in G305 and allows the physical properties and kinematics of the region to be studied. We identify 15 large molecular clouds and 57 smaller molecular clumps towards G305. The physical properties of the molecular gas are consistent with G305 being amongst the most massive a vigorous star forming regions in the Galaxy. We find a total molecular gas mass of 2:5–6:5 105M indicating that there is a large reservoir for future star formation. By considering virial equilibrium within the molecular clumps we discover that only 14% of the molecular clumps in G305 are gravitationally unstable, however these clumps contain > 30% of the molecular mass in G305 suggesting there is scope for considerable future star formation. To study the ionised environment towards G305 we have obtained some of the largest and most detailed wide-area mosaics with the Australia Telescope Compact Array to date. These radio continuum observations were performed simultaneously at 5.5 and 8.8 GHz and by applying two imaging techniques we are able to resolve HII regions from the ultra-compact to classical evolutionary phase. This has allowed high-mass star formation within G305 to be traced over the extent and lifetime of the complex. We discover that more than half of the observable total ionising flux in G305 is associated with embedded high-mass star formation around the periphery of a central cavity that has been driven into the molecular gas by a cluster of optically visible massive stars. By considering the contribution of embedded and visible massive stars to the observed radio continuum we suggest that more than 45 massive stars exist within G305. Combination of these two studies and recent and ongoing star formation provides the most in depth view of G305 to date and allows the star formation history and impact of high-mass stars to be investigated. We find compelling morphological evidence that suggests triggering is responsible for at least some of the observed high-mass star formation and construct a star formation history for the region.

Despite its importance the fundamental question of how massive stars form remains unanswered, with improvements to both models and observations having crucial roles to play. To quote Bate et al. (2003) computational models of star formation are limited because “conditions in molecular clouds are not sufficiently well understood to be able to select a representative sample of cloud cores for the initial conditions”. It is this notion that motivates the study of the environments within Giant Molecular Clouds (GMCs) and Infrared Dark Clouds (IRDCs), known sites of massive star formation, at the clump and core level. By studying large populations of these objects, it is possible to make conclusions based on global properties. With this in mind I study the dense molecular clumps within one of the most massive GMCs in the Galaxy: the W51 GMC. New observations of the W51 GMC in the 12 CO, 13 CO and C18 O (3 – 2) transitions using the HARP instrument on the JCMT are presented. With the help of the clump finding algorithm CLUMPFIND a total of 1575 dense clumps are identified of which 1130 are associated with the W51 GMC, yielding a dense mass reservoir of 1.5 × 105 M contained within these clumps. Of these clumps only 1% by number are found to be super-critical, yielding a super-critical clump formation efficiency of 0.5%, below current SFE estimates of the region. This indicates star formation within the W51 GMC will diminish over time although evidence from the first search for molecular outflows presents the W51 GMC in an active light with a lower limit of 14 outflows. The distribution of the outflows within the region searched found them concentrated towards the W51A region. Having much smaller sizes and masses, obtaining global properties of clumps and cores within IRDCs required studying a large sample of these objects. To do this pre-existing data from the SCUBA Legacy Catalogue was utilised to study IRDCs within a catalogues based on 8 µm data. This data identifies 154 IRDC cores that are detected at 850 µm and 51 cores that were not. This work suggests that cores not detected at 850 µm are low mass, low column density and low temperature cores that are below the sensitivity limit of SCUBA at 850 µm. Utilising observations at 24 µm from the Spitzer space telescope, allows for an investigation of current star formation by looking for warm embedded objects within the cores. This work reveals 69% of the IRDC cores have 24 µm embedded objects. IRDC cores without associated 24 µm emission (“starless” IRDC cores) may have yet to form stars, or may contain low mass YSOs below the detection limit. If it is assumed that cores without 24 µm embedded sources are at an earlier evolutionary stage to cores with embedded objects a statistical lifetime for the quiescent phase of a few 103 – 104 years is derived.

The primary focus of my PhD is to quantify the spatial distribution of star-forming environments from optical to radio wavelengths using data from the Hubble Space Telescope, the Very Large Telescope, the Spitzer Space Telescope, the Herschel Space Observatory, and the Caltech Submillimeter Observatory. Towards the end of my PhD study I have developed theoretical models. With these observational and theoretical avenues I have led a series of research projects to (1) quantify the initial spatial structure of pre-stellar cores and proto-stars, (2) test whether massive stars can form in isolation or not, (3) and develop a theoretical model on how young massive clusters form. These research projects have been fruitful as my collaborators and I have shown that pre-stellar cores and stars form in a smooth continuum of surface densities from a few to thousands of stars per pc^2. These two works have important implications on our understanding of what a young stellar cluster is and how star forming environments can evolve to form field star populations or gravitationally bound clusters. In my second study my collaborators and I found evidence for isolated massive star formation in the active star forming region 30 Doradus, in the Large Magellanic Cloud. The result impacts the field of the initial mass function and star formation models. Massive stars forming in isolation is consistent with a stochastically sampled initial mass function. Additionally, the result would put constraints on theoretical models on massive star formation. Continuing my work on massive star forming environments my collaborators and I have developed a theoretical model on how young massive clusters form. From the models we argue that feedback energies can be contained by the gravitational potential well of the massive progenitors. Furthermore, we predict the physical properties the massive cluster progenitors in terms of initial gas mass, radii and flux brightness to enable a search for these objects in Galactic plane surveys and upcoming telescopes. Using the common thread of spatial distribution analysis of star formation I describe my future research plans, which entails studies on extragalactic scales in the conclusion.

In this thesis, I explore both obscured and unobscured star formation over a large fraction of cosmic time. I use the HAWK-I Y -band science verification data over GOODS-South, in conjunction with optical and infrared data to search for Lyman-break galaxies at z >∼ 6.5 (i.e. within the first billion years of the Universe). I find four possible (two robust) z′-drop candidates (z >∼ 6.5) and four possible (but no robust) Y -drop candidates (z >∼ 7). I use my results to place constraints on the luminosity function at z ∼ 6.5 and find significant evolution in the population of Lyman-break galaxies between 3 < z <∼ 6.5. I also explore obscured star formation with a population of 70μm selected galaxies over the COSMOS field. I use AAT spectroscopy in conjunction with other available spectroscopic redshifts for my sample, and photometric redshifts otherwise, to calculate the total infrared luminosity of each galaxy. Two libraries of spectral energy distributions are considered; Siebenmorgen & Krügel (2007) templates and Chary & Elbaz (2001) models. We have supplemented our data with that of Huynh et al. (2007) collected over the GOODS-North field and adapted it to directly compare with the results of this work. The far-infrared luminosity function is then determined using the 1/Vmax technique. A double power law parameterisation is found to provide the best fit to the data. The far-infrared luminosity function was fitted for all parameters and the evolution was measured out to z ∼ 1. Three different types of evolution were allowed, pure luminosity, pure density and luminosity dependent density evolution. In all cases strong positive evolution was evident with the best-fit case being pure luminosity evolution where p = 2.4+0.6 −0.7. Due to the larger volume surveyed compared to previous studies, this work provides better constraints on the bright end of the far-infrared luminosity function displaying a shallower bright end slope (α2 ∼ −1.6) than previously determined, implying a higher number density of the most luminous objects and thereby a greater contribution from these objects to the total infrared energy density. However the shallower slope determined here can be reconciled with other work if the Chary & Elbaz (2001) models are used instead of the Siebenmorgen & Krügel templates; demonstrating that spectral energy distribution model selection is a key component in determining luminosity functions at far-infrared wavelengths. The far-infrared–radio correlation (FIRC; qIR) was determined for the sample of 70μm selected star-forming galaxies using 1.4GHz radio data over the COSMOS field, and no evolution was found out to z ∼ 2. The 70μm monochromatic evolution in the FIRC was also examined (q70) and no evolution was found in this parameter with redshift.

Recent observations at the James Clerk Maxwell Telescope (JCMT) and elsewhere have identified a class of very deeply embedded, possibly protostellar, sources which are not associated with any of the traditional indicators of star formation, such as HII regions and near-infrared emission, but which do lie close to otherwise isolated H2O masers. This thesis describes a search, based on catalogues of known water maser positions, for new deeply embedded cores similar to those found in S106 and M17. In addition to millimetre molecular line and submillimetre continuum observations, 22 GHz and 8 GHz radio observations have been made of a number of the sources in order to obtain more accurate maser positions and to search for any associated compact HII regions. Observing sources such as these in less active star forming regions provides a cleaner environment in which to examine the maser excitation and the ongoing process of star formation. A sample of 44 water maser sources was observed from which submillimetre continuum emission was detected from 40 (91 per cent). The most striking feature of the data is the close association of the masers with the submillimetre cores: the data are consistent with masers occurring within 6000 AU of the embedded core. The results can be summarised as follows: o High temperature gas has been detected, and most of the submillimetre cores have mean densities greater than 10^6 /cm^3. o The masers have low velocities with respect to the molecular cloud and are uniformly distributed within 6600 AU of the submillimetre core. o The isotropic maser luminosity is proportional to the far-infrared luminosity over more than 5 orders of magnitude. o There is no obvious correlation between the near-infrared spectral class and the maser emission. o Where a radio spectral index is known the majority of sources are optically thin HII regions. The bulk of the remainder are undetected and have a flux density less than 1 mJy. o Masers not associated directly with a submillimetre core show explicit evidence for shocks. Embedded cores \emph{have} been detected with this survey and the maser emission is consistent with collisionally excited pump models.