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.

abstract: Understanding the temperature structure of protoplanetary disks (PPDs) is paramount to modeling disk evolution and future planet formation. PPDs around T Tauri stars have two primary heating sources, protostellar irradiation, which depends on the flaring of the disk, and accretional heating as viscous coupling between annuli dissipate energy. I have written a "1.5-D" radiative transfer code to calculate disk temperatures assuming hydrostatic and radiative equilibrium. The model solves for the temperature at all locations simultaneously using Rybicki's method, converges rapidly at high optical depth, and retains full frequency dependence. The likely cause of accretional heating in PPDs is the magnetorotational instability (MRI), which acts where gas ionization is sufficiently high for gas to couple to the magnetic field. This will occur in surface layers of the disk, leaving the interior portions of the disk inactive ("dead zone"). I calculate temperatures in PPDs undergoing such "layered accretion." Since the accretional heating is concentrated far from the midplane, temperatures in the disk's interior are lower than in PPDs modeled with vertically uniform accretion. The method is used to study for the first time disks evolving via the magnetorotational instability, which operates primarily in surface layers. I find that temperatures in layered accretion disks do not significantly differ from those of "passive disks," where no accretional heating exists. Emergent spectra are insensitive to active layer thickness, making it difficult to observationally identify disks undergoing layered vs. uniform accretion. I also calculate the ionization chemistry in PPDs, using an ionization network including multiple charge states of dust grains. Combined with a criterion for the onset of the MRI, I calculate where the MRI can be initiated and the extent of dead zones in PPDs. After accounting for feedback between temperature and active layer thickness, I find the surface density of the actively accreting layers falls rapidly with distance from the protostar, leading to a net outward flow of mass from ~0.1 to 3 AU. The clearing out of the innermost zones is possibly consistent with the observed behavior of recently discovered "transition disks." / Dissertation/Thesis / Ph.D. Physics 2012

abstract: The lowest-mass stars, known as M-dwarfs, form target samples for upcoming exoplanet searches, and together with lower-mass substellar objects known as brown dwarfs, are among prime targets for detailed study with high-contrast adaptive optics (AO) imaging and sub-millimeter interferometry. In this thesis, I describe results from three studies investigating the companion properties and environments of low-mass systems: (1) The 245-star M-dwarfs in Multiples (MinMs) Survey, a volume-limited survey of field M-dwarf companions within 15 pc, (2) the Taurus Boundary of Stellar/Substellar (TBOSS) Survey, an ongoing study of disk properties for low-mass members within the Taurus star-forming region, and (3) spectroscopy of a brown dwarf companion using the Gemini Planet Imager (GPI). Direct imaging of M-dwarfs is a sensitive technique to identify low-mass companions over a wide range of orbital separation, and the high proper motion of nearby M-dwarfs eases confirmation of new multiple stars. Combining AO and wide-field imaging, the MinMs Survey provides new measurements of the companion star fraction (CSF), separation distribution, and mass ratio distribution for the nearest K7-M6 dwarfs. These results demonstrate the closer orbital separations (~6 AU) and lower frequency (~23% CSF) of M-dwarf binaries relative to higher-mass stars. From the TBOSS project, I report 885µm Atacama Large Millimeter/sub-millimeter Array continuum measurements for 24 Taurus members spanning the stellar/substellar boundary (M4-M7.75). Observations of submillimeter emission from dust grains around the lowest-mass hosts show decreasing disk dust mass for decreasing host star mass, consistent with low frequencies of giant planets around M-dwarfs. Compared to the older stellar association of Upper Scorpius, Taurus disks have a factor of four higher mass in submillimeter-sized grains. From the GPI Exoplanet Survey, I describe near-infrared spectroscopy of an unusually red companion orbiting inside the debris disk of an F5V star. As the second brown dwarf discovered within the innermost region of a debris disk, the properties of this system offer important dynamical constraints for companion-disk interaction and a useful benchmark for brown dwarf and giant planet atmospheric study. / Dissertation/Thesis / Doctoral Dissertation Astrophysics 2017

HD 100453AB is a 10 +/- 2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.'' 05, or similar to 108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, which prevented a proper treatment of the dynamical influence of the companion on the disk. We observed HD. 100453AB between 2015 and 2017, utilizing extreme adaptive optics systems on the Very Large Telescope and the Magellan Clay Telescope. We combined the astrometry from these observations with published data to constrain the parameters of the binary's orbit to a = 1.'' 06 +/- 0.'' 09, e = 0.17 +/- 0.07, and i = 32 degrees.5 +/- 6 degrees.5. We utilized publicly available ALMA (CO)-C-12 data to constrain the inclination of the disk, i(disk) similar to 28 degrees, which is relatively coplanar with the orbit of the companion and consistent with previous estimates from scattered light images. Finally, we input these constraints into hydrodynamic and radiative transfer simulations to model the structural evolution of the disk. We find that the spiral structure and truncation of the circumprimary disk in HD 100453 are consistent with a companion-driven origin. Furthermore, we find that the primary star's rotation, its outer disk, and the companion exhibit roughly the same direction of angular momentum, and thus the system likely formed from the same parent body of material.

We present Magellan adaptive optics Ha imaging of HH 508, which has the highest surface brightness among protostellar jets in the Orion Nebula. We find that HH 508 actually has a shorter component to the west, and a longer and knotty component to the east. The east component has a kink at 0.'' 3 from the jet-driving star theta(1) Ori B-2, so it may have been deflected by the wind/radiation from the nearby theta(1) Ori B1B5. The origin of both components is unclear, but if each of them is a separate jet, then theta(1) Ori B-2 may be a tight binary. Alternatively, HH 508 may be a slow-moving outflow, and each component represents an illuminated cavity wall. The ionization front surrounding theta(1) Ori B2B3 does not directly face theta(1) Ori B1B5, suggesting that the EUV radiation from theta(1) Ori C plays a dominant role in affecting the morphology of proplyds even in the vicinity of theta(1) Ori B1B5. Finally, we report an Ha blob that might be ejected by the binary proplyd LV 1.

In this paper, we investigate the relation between disk mass and mass accretion rate to constrain the mechanism of angular momentum transport in protoplanetary disks. We find a correlation between dust disk mass and mass accretion rate in Chamaeleon I with a slope that is close to linear, similar to the one recently identified in Lupus. We investigate the effect of stellar mass and find that the intrinsic scatter around the best-fit M-dust-M star and M-acc-M star relations is uncorrelated. We simulate synthetic observations of an ensemble of evolving disks using a Monte Carlo approach and find that disks with a constant alpha viscosity can fit the observed relations between dust mass, mass accretion rate, and stellar mass but overpredict the strength of the correlation between disk mass and mass accretion rate when using standard initial conditions. We find two possible solutions. In the first one, the observed scatter in M-dust and M-acc is not primordial, but arises from additional physical processes or uncertainties in estimating the disk gas mass. Most likely grain growth and radial drift affect the observable dust mass, while variability on large timescales affects the mass accretion rates. In the second scenario, the observed scatter is primordial, but disks have not evolved substantially at the age of Lupus and Chamaeleon I owing to a low viscosity or a large initial disk radius. More accurate estimates of the disk mass and gas disk sizes in a large sample of protoplanetary disks, through either direct observations of the gas or spatially resolved multiwavelength observations of the dust with ALMA, are needed to discriminate between both scenarios or to constrain alternative angular momentum transport mechanisms such as MHD disk winds.

The evolution of protoplanetary disks is believed to be driven largely by angular momentum transport resulting from magnetized disk winds and turbulent viscosity. The ionization of the disk that is essential for these processes has been thought to be due to host star coronal X-rays but could also arise from energetic particles produced by coronal flares, or traveling shock waves, and advected by the stellar wind. We have performed test-particle numerical simulations of energetic protons propagating into a realistic T. Tauri stellar wind, including a superposed small-scale magnetostatic turbulence. The isotropic (Kolmogorov power spectrum) turbulent component is synthesized along the individual particle trajectories. We have investigated the energy range [0.1-10] GeV, consistent with expectations from Chandra X-ray observations of large flares on T. Tauri stars and recent indications by the Herschel Space Observatory of a significant contribution of energetic particles to the disk ionization of young stars. In contrast with a previous theoretical study finding a dominance of energetic particles over X-rays in the ionization throughout the disk, we find that the disk ionization is likely dominated by X-rays over much of its area, except within narrow regions where particles are channeled onto the disk by the strongly tangled and turbulent magnetic field. The radial thickness of such regions is 5 stellar radii close to the star and broadens with increasing radial distance. This likely continues out to large distances from the star (10 au or greater), where particles can be copiously advected and diffused by the turbulent wind.

We present multi-wavelength radio observations obtained with the VLA of the protoplanetary disk surrounding the young brown dwarf 2MASS J04442713+2512164 (2M0444) in the Taurus star-forming region. 2M0444 is the brightest known brown dwarf disk at millimeter wavelengths, making this an ideal target to probe radio emission from a young brown dwarf. Thermal emission from dust in the disk is detected at 6.8 and 9.1 mm, whereas the 1.36 cm measured flux is dominated by ionized gas emission. We combine these data with previous observations at shorter sub-mm and mm wavelengths to test the predictions of dust evolution models in gas-rich disks after adapting their parameters to the case of 2M0444. These models show that the radial drift mechanism affecting solids in a gaseous environment has to be either completely made inefficient, or significantly slowed down by very strong gas pressure bumps in order to explain the presence of mm/cm-sized grains in the outer regions of the 2M0444 disk. We also discuss the possible mechanisms for the origin of the ionized gas emission detected at 1.36 cm. The inferred radio luminosity for this emission is in line with the relation between radio and bolometric luminosity valid for for more massive and luminous young stellar objects, and extrapolated down to the very low luminosity of the 2M0444 brown dwarf.

For decades, astronomers have used observations of mature planetary systems to constrain planet formation theories, beginning with our own solar system and now the thousands of known exoplanets. Recent advances in instrumentation have given us a direct view of some steps in the planet formation process, such as large-scale protostar and protoplanetary disk features and evolution. However, understanding the details of how planets accrete and interact with their environment requires direct observations of protoplanets themselves. Transition disks, protoplanetary disks with inner clearings that may be caused by forming planets, are the best targets for these studies. Their large distances, compared to the stars normally targeted for direct imaging of exoplanets, make protoplanet detection difficult and necessitate novel imaging techniques. In this dissertation, I describe the results of using non-redundant masking (NRM) to search for forming planets in transition disk clearings. I first present a data reduction pipeline that I wrote to this end, using example datasets and simulations to demonstrate reduction and imaging optimizations. I discuss two transition disk NRM case studies: T Cha and LkCa 15. In the case of T Cha, while we detect significant asymmetries, the data cannot be explained by orbiting companions. The fluxes and orbital motion of the LkCa 15 companion signals, however, can be naturally explained by protoplanets in the disk clearing. I use these datasets and simulated observations to illustrate the effects of scattered light from transition disk material on NRM protoplanet searches. I then demonstrate the utility of the dual-aperture Large Binocular Telescope Interferometer's NRM mode on the bright B[e] star MWC 349A. I discuss the implications of this work for planet formation studies as well as future prospects for NRM and related techniques on next generation instruments.

Small solids embedded in gaseous protoplanetary disks are subject to strong dust-gas friction. Consequently, tightly coupled dust particles almost follow the gas flow. This near conservation of the dust-to-gas ratio along streamlines is analogous to the near conservation of entropy along flows of (dust-free) gas with weak heating and cooling. We develop this thermodynamic analogy into a framework to study dusty gas dynamics in protoplanetary disks. We show that an isothermal dusty gas behaves like an adiabatic pure gas, and that finite dust-gas coupling may be regarded as effective heating/cooling. We exploit this correspondence to deduce that (1) perfectly coupled, thin dust layers cannot cause axisymmetric instabilities; (2) radial dust edges are unstable if the dust is vertically well-mixed; (3) the streaming instability necessarily involves a gas pressure response that lags behind dust density; and (4) dust-loading introduces buoyancy forces that generally stabilize the vertical shear instability associated with global radial temperature gradients. We also discuss dusty analogs of other hydrodynamic processes (e.g., Rossby wave instability, convective overstability, and zombie vortices) and how to simulate dusty protoplanetary disks with minor tweaks to existing codes for pure gas dynamics.