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).

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.

Understanding how diffuse molecular clouds at large scales (~10 pc) assemble mass into dense, star-forming cores at small scales (~ 0.1 pc) is crucial to building a holistic theory of star formation. While recent observations suggest that filaments play an important role in the mass assembly of dense cores, detailed gas kinematics studies are still lacking. My dissertation presents three innovative techniques that enable us to study star-forming filaments' complex gas kinematics in unprecedented detail: multi-component spectral fit, multi-dimensional filament identification, and membership assignment of velocity-coherent structures. Through these techniques, I analyzed star-forming filaments in the Perseus Molecular Cloud and unveiled unexpectedly complex velocity structures at scales where filaments are well resolved, to as low as the 0.01 pc scale. Moreover, the correlations I discovered between the various filament properties further suggest a scenario in which thermally supercritical filaments grow continuously via accretion from their surroundings while simultaneously forming cores through fragmentation along their lengths. / Graduate / 2022-01-08

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.

Cette thèse a pour but de faire le lien entre l’évolution des galaxies, leur morphologie et les processus physiques internes, notamment la formation stellaire comme le résultat du milieu interstellaire turbulent et multiphase, en utilisant les simulations cosmologiques zoom-in, les simulations des galaxies isolées et en interaction, et le modèle analytique de la formation stellaire. Dans le chapitre 1, j’explique la motivation pour cette thèse et je passe brièvement en revue le contexte nécessaire lié à la formation des galaxies et la modélisation en utilisant les simulations numériques. Tout d’abord, j’explore l’évolution de la morphologie des galaxies du type de la Voie Lactée dans la série des simulations cosmologiques zoom-in à travers l’analyse des barres. J’analyse l’évolution de la fraction des barres avec le redshift, sa dépendance en fonction de la masse stellaire et l’histoire d’accrétion de galaxies individuelles. Je montre en particulier, que la fraction de barres décroit avec le redshift croissant, en accord avec les observations. Ce travail montre également que les résultats obtenus suggèrent que l’époque de la formation des barres correspond à la transition entre une phase précoce “violente” de la formation de galaxies spirales à z > 1, pendant laquelle elles sont souvent perturbées par les fusions avec les galaxies de masse comparable ou par multiple fusions avec les galaxies de petite masse, mais aussi les instabilités violentes de disque, et une phase "séculaire" tardive à z < 1, quand la morphologie finale est généralement stabilisée vers une structure dominée par le disque. Cette analyse est présentée dans le chapitre 2. Étant donné que ces simulations cosmologiques forment trop d'étoiles trop tôt par rapport aux populations de galaxies observées, je me concentre dans le chapitre 3 sur la formation stellaire dans un échantillon de simulation de galaxies en isolation, à bas redshift, et à résolution du parsec et sous-parsec. J'étudie l'origine physique de leurs relations de formation stellaire avec les cassures, et montre que le seuil de densité surfacique pour une formation stellaire efficace peut être lié à la densité caractéristique d'apparition de turbulence supersonique. Ce résultat s'applique aussi bien aux galaxies qui fusionnent, dans lesquelles l'augmentation de la turbulence compressive déclenchée par les marées compressives les conduit au régime de sursaut de formation d'étoiles. Un modèle analytique idéalisé de formation stellaire liant la densité surfacique de gaz au taux de formation stellaire comme une fonction de la présence de turbulence supersonique et la structure associée du milieu interstellaire est ensuite présenté dans le chapitre 4. Ce modèle prédit une cassure à basse densité de surface qui est suivie par un régime de loi de puissance à haute densité dans différents systèmes en accord avec les relations de formation stellaire des galaxies observées et simulées. La dernière partie de cette thèse est dédiée à la technique alternative de zoom-in cosmologique (Martig et al. 2009) et son implémentation dans le code à raffinement de maillage adaptatif RAMSES. Dans le chapitre 5, je présente les caractéristiques de base de cette technique aussi bien que certains de nos tout premiers résultats dans le contexte de l'accrétion cosmologique diffuse. / This thesis aims at making the link between galaxy evolution, morphology and internal physical processes, namely star formation as the outcome of the turbulent multiphase interstellar medium, using the cosmological zoom-in simulations, simulations of isolated and merging galaxies, and the analytic model of star formation. In Chapter 1, I explain the motivation for this thesis and briefly review the necessary background related to galaxy formation and modeling with the use of numerical simulations. I first explore the evolution of the morphology of Milky-Way-mass galaxies in a suite of zoom-in cosmological simulations through the analysis of bars. I analyze the evolution of the fraction of bars with redshift, its dependence on the stellar mass and accretion history of individual galaxies. I show in particular, that the fraction of bars declines with increasing redshift, in agreement with the observations. This work also shows that the obtained results suggest that the bar formation epoch corresponds to the transition between an early "violent" phase of spiral galaxies formation at z > 1, during which they are often disturbed by major mergers or multiple minor mergers as well as violent disk instabilities, and a late "secular" phase at z < 1, when the final morphology is generally stabilized to a disk-dominated structure. This analysis is presented in Chapter 2. Because such cosmological simulations form too many stars too early compared to observed galaxy populations, I shift the focus in Chapter 3 to star formation in a sample of low-redshift galaxy simulations in isolation at parsec and sub-parsec resolution. I study the physical origin of their star formation relations and breaks and show that the surface density threshold for efficient star formation can be related to the typical density for the onset of supersonic turbulence. This result holds in merging galaxies as well, where increased compressive turbulence triggered by compressive tides during the interaction drives the merger to the regime of starbursts. An idealized analytic model for star formation relating the surface density of gas and star formation rate as a function of the presence of supersonic turbulence and the associated structure of the ISM is then presented in Chapter 4. This model predicts a break at low surface densities that is followed by a power-law regime at high densities in different systems in agreement with star formation relations of observed and simulated galaxies. The last part of this thesis is dedicated to the alternative cosmological zoom-in technique Martig et al. 2009 and its implementation in the Adaptive Mesh Refinement code RAMSES. In Chapter 5, I will present the basic features of this technique as well as some of our very first results in the context of smooth cosmological accretion.