Thesis (Ph.D.) - University of St Andrews, November 2008.

abstract: The composition of planets and their volatile contents are intimately connected to the structure and evolution of their parent protoplanetary disks. The transport of momentum and volatiles is often parameterized by a turbulent viscosity parameter $\alpha$, which is usually assumed to be spatially and temporally uniform across the disk. I show that variable $\alpha$(r,z) (where $r$ is radius, and $z$ is height from the midplane) attributable to angular momentum transport due to MRI can yield disks with significantly different structure, as mass piles up in the 1-10 AU region resulting in steep slopes of p $>$ 2 here (where p is the power law exponent in $\Sigma \propto r^{-p}$). I also show that the transition radius (where bulk mass flow switches from inward to outward) can move as close in as 3 AU; this effect (especially prominent in externally photoevaporated disks) may significantly influence the radial water content available during planet formation. I then investigate the transport of water in disks with different variable α profiles. While radial temperature profile sets the location of the water snowline (i.e., inside of which water is present as vapor; outside of which, as ice on solids), it is the rates of diffusion and drift of small icy solids and diffusion of vapor across the snow line that determine the radial water distribution. All of these processes are highly sensitive to local $\alpha$. I calculate the effect of radially varying α on water transport, by tracking the abundance of vapor in the inner disk, and fraction of ice in particles and larger asteroids beyond the snow line. I find one α profile attributable to winds and hydrodynamical instabilities, and motivated by meteoritic constraints, to show considerable agreement with inferred water contents observed in solar system asteroids. Finally, I calculate the timing of gap formation due to the formation of a planet in disks around different stars. Here, I assume that pebble accretion is the dominant mechanism for planetary growth and that the core of the first protoplanet forms at the water snow line. I discuss the dependence of gap timing to various stellar and disk properties. / Dissertation/Thesis / Doctoral Dissertation Astrophysics 2018

Thesis (Ph. D.)--University of Washington, 2004. / Vita. Includes bibliographical references (p. 103-106).

Class I protostars are thought to represent an early stage in the lifetime of protoplanetary disks, when they are still embedded in their natal envelope. Here we measure the disk masses of 10 Class I protostars in the Taurus Molecular Cloud to constrain the initial mass budget for forming planets in disks. We use radiative transfer modeling to produce synthetic protostar observations and fit the models to a multi-wavelength data set using a Markov Chain Monte Carlo fitting procedure. We fit these models simultaneously to our new Combined Array for Research in Millimeter-wave Astronomy 1.3 mm observations that are sensitive to the wide range of spatial scales that are expected from protostellar disks and envelopes so as to be able to distinguish each component, as well as broadband spectral energy distributions compiled from the literature. We find a median disk mass of 0.018 M-circle dot on average, more massive than the Taurus Class II disks, which have median disk mass of similar to 0.0025 M-circle dot. This decrease in disk mass can be explained if dust grains have grown by a factor of 75 in grain size, indicating that by the Class II stage, at a few Myr, a significant amount of dust grain processing has occurred. However, there is evidence that significant dust processing has occurred even during the Class I stage, so it is likely that the initial mass budget is higher than the value quoted here.

We present Karl G. Jansky Very Large Array 1.3, 3.6, and 6 cm continuum maps of compact radio sources in the Orion Nebular Cluster (ONC). We mosaicked 34 arcmin(2) at 1.3 cm, 70 arcmin(2) at 3.6 cm and 109 arcmin(2) at 6 cm, containing 778 near-infrared detected young stellar objects and 190 Hubble Space Telescope-identified proplyds (with significant overlap between those characterizations). We detected radio emission from 175 compact radio sources in the ONC, including 26 sources that were detected for the first time at these wavelengths. For each detected source, we fitted a simple free-free and dust emission model to characterize the radio emission. We extrapolate the free-free emission spectrum model for each source to ALMA bands to illustrate how these measurements could be used to correctly measure protoplanetary disk dust masses from submillimeter flux measurements. Finally, we compare the fluxes measured in this survey with previously measured fluxes for our targets, as well as four separate epochs of 1.3 cm data, to search for and quantify the variability of our sources.

We study the formation of planetesimals in protoplanetary disks from the gravitational collapse of solid over-densities generated via the streaming instability. To carry out these studies, we implement and test a particle-mesh self-gravity module for the ATHENA code that enables the simulation of aerodynamically coupled systems of gas and collisionless self-gravitating solid particles. Upon employment of our algorithm to planetesimal formation simulations, we find that (when a direct comparison is possible) the ATHENA simulations yield predicted planetesimal properties that agree well with those found in prior work using different numerical techniques. In particular, the gravitational collapse of streaming-initiated clumps leads to an initial planetesimal mass function that is well-represented by a power law, dN / dM(p) proportional to M-p(-p), with p similar or equal to 1.6 +/- 0.1, which equates to a differential size distribution of dN / dR(p) proportional to R-p(-q), with q similar or equal to 2.8 +/- 0.1. We find no significant trends with resolution from a convergence study of up to 512(3) grid zones and N-par approximate to 1.5 x 10(8) particles. Likewise, the power-law slope appears indifferent to changes in the relative strength of self-gravity and tidal shear, and to the time when (for reasons of numerical economy) self-gravity is turned on, though the strength of these claims is limited by small number statistics. For a typically assumed radial distribution of minimum mass solar nebula solids (assumed here to have dimensionless stopping time tau = 0.3), our results support the hypothesis that bodies on the scale of large asteroids or Kuiper Belt Objects could have formed as the high-mass tail of a primordial planetesimal population.

Context. Understanding the diversity of planets requires studying the morphology and physical conditions in the protoplanetary disks in which they form. Aims. We aim to study the structure of the similar to 10 Myr old protoplanetary disk HD 100453, to detect features that can trace disk evolution and to understand the mechanisms that drive these features. Methods. We observed HD100453 in polarized scattered light with VLT/SPHERE at optical (0.6 mu m, 0.8 mu m) and near-infrared (1.2 mu m) wavelengths, reaching an angular resolution of similar to 0.02 '', and an inner working angle of similar to 0.09 ''. Results. We spatially resolve the disk around HD 100453, and detect polarized scattered light up to similar to 0.42 '' (similar to 48 au). We detect a cavity, a rim with azimuthal brightness variations at an inclination of similar to 38 degrees with respect to our line of sight, two shadows and two symmetric spiral arms. The spiral arms originate near the location of the shadows, close to the semi major axis. We detect a faint feature in the SW that can be interpreted as the scattering surface of the bottom side of the disk, if the disk is tidally truncated by the M-dwarf companion currently seen at a projected distance of similar to 119 au. We construct a radiative transfer model that accounts for the main characteristics of the features with an inner and outer disk misaligned by similar to 72 degrees. The azimuthal brightness variations along the rim are well reproduced with the scattering phase function of the model. While spirals can be triggered by the tidal interaction with the companion, the close proximity of the spirals to the shadows suggests that the shadows could also play a role. The change in stellar illumination along the rim induces an azimuthal variation of the scale height that can contribute to the brightness variations. Conclusions. Dark regions in polarized images of transition disks are now detected in a handful of disks and often interpreted as shadows due to a misaligned inner disk. However, the origin of such a misalignment in HD100453, and of the spirals, is still unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit. Observations over a few years will allow us to measure the spiral pattern speed, and determine if the shadows are fixed or moving, which may constrain their origin.

We present the highest spatial resolution ALMA observations to date of the Class I protostar WL 17 in the rho Ophiuchus L1688 molecular cloud complex, which show that it has a 12 au hole in the center of its disk. We consider whether WL 17 is actually a Class II disk being extincted by foreground material, but find that such models do not provide a good fit to the broadband spectral energy distribution (SED) and also require such high extinction that it would presumably arise from dense material close to the source, such as a remnant envelope. Self-consistent models of a disk embedded in a rotating collapsing envelope can nicely reproduce both the ALMA 3 mm observations and the broadband SED of WL 17. This suggests that WL 17 is a disk in the early stages of its formation, and yet even at this young age the inner disk has been depleted. Although there are multiple pathways for such a hole to be created in a disk, if this hole was produced by the formation of planets it could place constraints on the timescale for the growth of planets in protoplanetary disks.

The properties of disks around brown dwarfs and very low mass stars (hereafter VLMOs) provide important boundary conditions on the process of planet formation and inform us about the numbers and masses of planets than can form in this regime. We use the Herschel Space Observatory PACS spectrometer to measure the continuum and [O I] 63 mu m line emission toward 11 VLMOs with known disks in the Taurus and Chamaeleon I star-forming regions. We fit radiative transfer models to the spectral energy distributions of these sources. Additionally, we carry out a grid of radiative transfer models run in a regime that connects the luminosity of our sources with brighter T Tauri stars. We find that VLMO disks with sizes 1.3-78 au, smaller than typical T Tauri disks, fit well the spectral energy distributions assuming that disk geometry and dust properties are stellar mass independent. Reducing the disk size increases the disk temperature, and we show that VLMOs do not follow previously derived disk temperature-stellar luminosity relationships if the disk outer radius scales with stellar mass. Only 2 out of 11 sources are detected in [O I] despite a better sensitivity than was achieved for T Tauri stars, suggesting that VLMO disks are underluminous. Using thermochemical models, we show that smaller disks can lead to the unexpected [O I] 63 mu m nondetections in our sample. The disk outer radius is an important factor in determining the gas and dust observables. Hence, spatially resolved observations with ALMA-to establish if and how disk radii scale with stellar mass-should be pursued further.

Thesis (Ph. D.)--University of Texas at Austin, 2001. / Vita. Includes bibliographical references (leaves 84-91). Available also in a digital version from Dissertation Abstracts.