Aims. We investigate the suitability of alpha-viscosity in self-similar models for self-gravitating disks with a focus on active galactic nuclei (AGN) disks. Methods. We use a self-similar approach to simplify the partial di ff erential equations arising from the evolution equation, which are then solved using numerical standard procedures. Results. We find a self-similar solution for the dynamical evolution of self-gravitating alpha-disks and derive the significant quantities. In the Keplerian part of the disk our model is consistent with standard stationary alpha-disk theory, and self-consistent throughout the self-gravitating regime. Positive accretion rates throughout the disk demand a high degree of self-gravitation. Combined with the temporal decline of the accretion rate and its low amount, the model prohibits the growth of large central masses. Conclusions. alpha-viscosity cannot account for the evolution of the whole mass spectrum of super-massive black holes (SMBH) in AGN. However, considering the involved scales it seems suitable for modelling protoplanetary disks.

We present new Atacama Large Millimeter/sub-millimeter Array (ALMA) 1.3 mm continuum observations of the SR 24S transition disk with an angular resolution less than or similar to 0'.18 (12 au radius). We perform a multi-wavelength investigation by combining new data with previous ALMA data at 0.45 mm. The visibilities and images of the continuum emission at the two wavelengths are well characterized by a ring-like emission. Visibility modeling finds that the ring-like emission is narrower at longer wavelengths, in good agreement with models of dust-trapping in pressure bumps, although there are complex residuals that suggest potentially asymmetric structures. The 0.45 mm emission has a shallower profile inside the central cavity than the 1.3 mm emission. In addition, we find that the (CO)-C-13 and (CO)-O-18 (J = 2-1) emission peaks at the center of the continuum cavity. We do not detect either continuum or gas emission from the northern companion to this system (SR 24N), which is itself a binary system. The upper limit for the dust disk mass of SR 24N is less than or similar to 0.12 M-circle plus, which gives a disk mass ratio in dust between the two components of M-dust,M-SR 24S/M-dust,M-SR 24N greater than or similar to 840. The current ALMA observations may imply that either planets have already formed in the SR 24N disk or that dust growth to millimeter sizes is inhibited there and that only warm gas, as seen by rovibrational CO emission inside the truncation radii of the binary, is present.

Using Keck/HIRES spectra (Delta v similar to 7 km s(-1)) we analyze forbidden lines of [O I] 6300 angstrom, [O I] 5577 angstrom. and [S II] 6731 angstrom. from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high-velocity component (HVC) associated with microjets, we study the properties of the low-velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (magnetohydrodynamic) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material. LVC emission is seen in all 30 stars with detected [O. I] but only in two out of eight with detected [S. II], so our analysis is largely based on the properties of the [O. I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks. The full width at half maximum of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 au and 0.5 to 5 au, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC, however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive.

Variable emission from Sgr A*, the luminous counterpart to the super-massive black hole at the center of our Galaxy, arises from the innermost portions of the accretion flow. Better characterization of the variability is important for constraining models of the low-luminosity accretion mode powering Sgr A*, and could further our ability to use variable emission as a probe of the strong gravitational potential in the vicinity of the 4 x 10(6) M-circle dot black hole. We use the Herschel Spectral and Photometric Imaging Receiver (SPIRE) to monitor Sgr. A* at wavelengths that are difficult or impossible to observe from the ground. We find highly significant variations at 0.25, 0.35, and 0.5 mm, with temporal structure that is highly correlated across these wavelengths. While the variations correspond to < 1% changes in the total intensity in the Herschel beam containing Sgr. A*, comparison to independent, simultaneous observations at 0.85 mm strongly supports the reality of the variations. The lowest point in the light curves, similar to 0.5 Jy below the time-averaged flux density, places a lower bound on the emission of Sgr. A* at 0.25 mm, the first such constraint on the THz portion of the spectral energy distribution. The variability on few hour timescales in the SPIRE light curves is similar to that seen in historical 1.3 mm data, where the longest time series is available, but the distribution of variations in the sub-mm do not show a tail of large-amplitude variations seen at 1.3 mm. Simultaneous X-ray photometry from XMM-Newton shows no significant variation within our observing period, which may explain the lack of very large submillimeter variations in our data if X-ray and submillimeter flares are correlated.

Context. Whether or not a self-gravitating accretion disk fragments is still an open issue. There are many different physical and numerical explanations for fragmentation, but simulations often show a non-convergent behavior for ever better resolution. Aims. We aim to investigate the influence of different numerical limiters in Godunov type schemes on the fragmentation boundary in self-gravitating disks. Methods. We have compared the linear and non-linear outcomes in two-dimensional shearingsheet simulations using the VANLE ER and the SUPERBEE limiter. Results. We show that choosing inappropriate limiting functions to handle shock-capturing in Godunov type schemes can lead to an overestimation of the surface density in regions with shallow density gradients. The effect amplifies itself on timescales comparable to the dynamical timescale even at high resolutions. This is exactly the environment in which clumps are expected to form. The effect is present without, but scaled up by, self-gravity and also does not depend on cooling. Moreover it can be backtracked to a well known effect called oversteepening. If the effect is also observed in the linear case, the fragmentation limit is shifted to larger values of the critical cooling timescale.

Closure phases along different baseline triangles carry a large amount of information regarding the structures of the images of black holes in interferometric observations with the Event Horizon Telescope. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A* to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of the small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the funnel-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the u-v space that show a high level of variability will need to be identified and excised.

In this thesis, we examine the role of accretion disks in the formation of stellar systems, focusing on young massive disks which regulate the flow of material from the parent molecular core down to the star. We study the evolution of disks with high infall rates that develop strong gravitational instabilities. We begin in chapter 1 with a review of the observations and theory which underpin models for the earliest phases of star formation and provide a brief review of basic accretion disk physics, and the numerical methods which we employ. In chapter 2 we outline the current models of binary and multiple star formation, and review their successes and shortcomings from a theoretical and observational perspective. In chapter 3 we begin with a relatively simple analytic model for disks around young, very massive stars, showing that instability in these disks may be responsible for the higher multiplicity fraction of massive stars, and perhaps the upper mass to which they grow. We extend these models in chapter 4 to explore the properties of disks and the formation of binary companions across a broad range of stellar masses. In particular, we model the role of global and local mechanisms for angular momentum transport in regulating the relative masses of disks and stars. We follow the evolution of these disks throughout the main accretion phase of the system, and predict the trajectory of disks through parameter space. We follow up on the predictions made in our analytic models with a series of high resolution, global numerical experiments in chapter 5. Here we propose and test a new parameterization for describing rapidly accreting, gravitationally unstable disks. We find that disk properties and system multiplicity can be mapped out well in this parameter space. Finally, in chapter 6, we address whether our studies of unstable disks are relevant to recently detected massive planets on wide orbits around their central stars.

We perform a spectroscopic survey of the foreground population in Orion. A with MMT/Hectospec. We use these data, along with archival spectroscopic data and photometric data, to derive spectral types, extinction values, and masses for 691 stars. Using the Spitzer Space Telescope data, we characterize the disk properties of these sources. We identify 37 new transition disk (TD) objects, 1 globally depleted disk candidate, and 7 probable young debris disks. We discover an object with a mass of. less than 0.018-0.030 M-circle dot, which harbors a flaring disk. Using the Ha emission line, we characterize the accretion activity of the sources with disks, and confirm that the. fraction of accreting TDs is lower than that of optically thick disks (46% +/- 7% versus 73% +/- 9%, respectively). Using kinematic data from the Sloan Digital Sky Survey and APOGEE INfrared Spectroscopy of the Young Nebulous Clusters program (IN-SYNC), we confirm that the foreground population shows similar kinematics to their local molecular clouds and other young stars in the same regions. Using the isochronal ages, we find that the foreground population has a median age of. around 1-2 Myr, which is similar to that of other young stars in Orion. A. Therefore, our results argue against the presence of a large and old foreground cluster in front of Orion. A.

Significant X-ray variability and flaring has been observed from Sgr A* but is poorly understood from a theoretical standpoint. We perform general relativistic magnetohydrodynamic simulations that take into account a population of non-thermal electrons with energy distributions and injection rates that are motivated by PIC simulations of magnetic reconnection. We explore the effects of including these non-thermal electrons on the predicted broadband variability of Sgr A* and find that X-ray variability is a generic result of localizing non-thermal electrons to highly magnetized regions, where particles are likely to be accelerated via magnetic reconnection. The proximity of these high-field regions to the event horizon forms a natural connection between IR and X-ray variability and accounts for the rapid timescales associated with the X-ray flares. The qualitative nature of this variability is consistent with observations, producing X-ray flares that are always coincident with IR flares, but not vice versa, i.e., there are a number of IR flares without X-ray counterparts.

Low mass, self-gravitating accretion disks admit quasi-steady, "gravito-turbulent" states in which cooling balances turbulent viscous heating. However, numerical simulations show that gravito-turbulence cannot be sustained beyond dynamical timescales when the cooling rate or corresponding turbulent viscosity is too large. The result is disk fragmentation. We motivate and quantify an interpretation of disk fragmentation as the inability to maintain gravito-turbulence due to formal secondary instabilities driven by: (1) cooling, which reduces pressure support; and/or (2) viscosity, which reduces rotational support. We analyze the axisymmetric gravitational stability of viscous, non-adiabatic accretion disks with internal heating, external irradiation, and cooling in the shearing box approximation. We consider parameterized cooling functions in 2D and 3D disks, as well as radiative diffusion in 3D. We show that generally there is no critical cooling rate/viscosity below which the disk is formally stable, although interesting limits appear for unstable modes with lengthscales on the order of the disk thickness. We apply this new linear theory to protoplanetary disks subject to gravito-turbulence modeled as an effective viscosity, and cooling regulated by dust opacity. We find that viscosity renders the disk beyond similar to 60 au dynamically unstable on radial lengthscales a few times the local disk thickness. This is coincident with the empirical condition for disk fragmentation based on a maximum sustainable stress. We suggest turbulent stresses can play an active role in realistic disk fragmentation by removing rotational stabilization against self-gravity, and that the observed transition in behavior from gravito-turbulent to fragmenting may reflect instability of the gravito-turbulent state itself.