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.

Theoretical models and spectroscopic observations of newborn stars suggest that protoplantary disks have an inner "wall" at a distance set by the disk interaction with the star. Around T Tauri stars, the size of this disk hole is expected to be on a 0.1 au scale that is unresolved by current adaptive optics imaging, though some model-dependent constraints have been obtained by near-infrared interferometry. Here we report the first measurement of the inner disk wall around a solar-mass young stellar object, YLW 16B in the rho Ophiuchi star-forming region, by detecting the light-travel time of the variable radiation from the stellar surface to the disk. Consistent time lags were detected on two nights, when the time series in H (1.6 mu m) and K (2.2 mu m) bands were synchronized while the 4.5 mu m emission lagged by 74.5 +/- 3.2 s. Considering the nearly edge-on geometry of the disk, the inner rim should be 0.084 au from the protostar on average, with an error of order 0.01 au. This size is likely larger than the range of magnetospheric truncations and consistent with an optically and geometrically thick disk front at the dust sublimation radius at similar to 1500 K. The widths of the cross-correlation functions between the data in different wavebands place possible new constraints on the geometry of the disk.

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.

Thesis (Ph. D.)--University of California, Riverside, 2008. / Includes abstract. Available via ProQuest Digital Dissertations. Title from first page of PDF file (viewed March 10, 2010). Includes bibliographical references (p. 126-135). Also issued in print.

This thesis is a computational study of circumstellar gas disks, with a special focus on modeling techniques and on numerical methods not only as scientific tools but also as a target of study. In particular, in-depth discussions are included on the main numerical strategy used, namely the moving-mesh method for astrophysical hydrodynamics. In this work, the moving-mesh approach is used to simulate circumstellar disks for the first time. / Astronomy

What processes are responsible for the dispersal of protoplanetary disks? In this dissertation, beginning with a brief Introduction to planet detection, disk dispersal and high-contrast imaging in Chapter 1, I will describe how ground-based adaptive optics (AO) imaging can help to inform these processes. Chapter 2 presents Polarized Differential Imaging (PDI) of the transitional disk SR21 at H-band taken as part of the Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS). These observations were the first to show that transition disk cavities can appear markedly different at different wavelengths. The observation that the sub-mm cavity is absent in NIR scattered light is consistent with grain filtration at a planet-induced gap edge. Chapter 3 presents SEEDS data of the transition disk Oph IRS 48. This highly asymmetrical disk is also most consistent with a planet-induced clearing mechanism. In particular, the images reveal both the disk cavity and a spiral arm/divot that had not been imaged previously. This study demonstrates the power of multiwavelength PDI imaging to verify disk structure and to probe azimuthal variation in grain properties. Chapter 4 presents Magellan visible light adaptive optics imaging of the silhouette disk Orion 218-354. In addition to its technical merits, these observations reveal the surprising fact that this very young disk is optically thin at H-alpha. The simplest explanation for this observation is that significant grain growth has occurred in this disk, which may be responsible for the pre-transitional nature of its SED. Chapter 5 presents brief descriptions of several other works-in-progress that build on my previous work. These include the MagAO Giant Accreting Protoplanet Survey (GAPlanetS), which will probe the inner regions of transition disks at unprecedented resolution in search of young planets in the process of formation. Chapters 6-8 represent my educational research in quantitative literacy, beginning with an introduction to the literature and study motivation in Chapter 6. Chapter 7 describes the development and validation of the Quantitative Reasoning for College Science (QuaRCS) Assessment instrument. Chapter 8 briefly describes the next steps for Phase II of the QuaRCS study.

We present Herschel observations of the debris disk around γ Doradus (HD 27290, HIP 19893) from the Herschel Key Programme DEBRIS (Disc Emission via Bias-free Reconnaissance in the Infrared/Submillimetre). The disk is well-resolved with PACS at 70, 100 and 160 micron and detected with SPIRE at 250 and 350 micron. The 250 micron image is only resolved along the disk's long axis. The SPIRE 500 micron 3 σ detection includes a nearby background source. γ Dor's spectral energy distribution (SED) is sampled in the submillimetre for the first time and modelled with multiple modified-blackbody functions to account for its broad shape. Two approaches are used, both of which reproduce the SED in the same way: a model of two narrow dust rings and a model of an extended, wide dust belt. The former implies the dust rings have temperatures of ~90 and ~40 K, corresponding to blackbody radii of 25 and 135 AU, respectively. The latter model suggests the dust lies in a wide belt extending from 15 to 230 AU. The resolved images, however, show dust extending beyond ~350 AU. This is consistent with other debris disks whose actual radii are observed to be a factor of 2 - 3 times larger than the blackbody radii. Although it is impossible to determine a preferred model from the SED alone, the resolved images suggest that the dust is located in a smooth continuous belt rather than discrete narrow rings. Both models estimate that the dust mass is 6.7 x 10^{-3} Earth masses and that fractional luminosity is 2.5 x 10^{-5}. This amount of dust is within the levels expected from steady state evolution given the age of γ Dor and therefore a transient event is not needed to explain the dust mass. No asymmetries that would hint at a planetary body are evident in the disk at Herschel's resolution. However, the constraints placed on the dust's location suggest that the most likely region to find planets is within 20 AU of the star. / Graduate

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