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.

We present the discovery of seven new proplyds (i.e., sources surrounded by cometary H alpha emission characteristic of offset ionization fronts (IFs)) in NGC 1977, located about 30' north of the Orion Nebula Cluster (ONC) at a distance of similar to 400 pc. Each of these proplyds is situated at projected distances 0.04-0.27 pc from the B1V star 42 Orionis (c Ori), which is the main source of UV photons in the region. In all cases the IFs of the proplyds are clearly pointing toward the common ionizing source, 42 Ori, and six of the seven proplyds clearly show tails pointing away from it. These are the first proplyds to be found around a B star, with previously known examples instead being located around O stars, including those in the ONC around theta(1) Ori C. The radii of the offset IFs in our proplyds are between similar to 200 and 550 au; two objects also contain clearly resolved central sources that we associate with disks of radii 50-70 au. The estimated strength of the FUV radiation field impinging on the proplyds is around 10-30 times less than that incident on the classic proplyds in the ONC. We show that the observed proplyd sizes are however consistent with recent models for FUV photoevaporation in relatively weak FUV radiation fields.

We present Gemini Planet Imager polarized intensity imagery of HD 100453 in Y, J, and K1 bands that reveals an inner gap (9-18 au), an outer disk (18-39 au) with two prominent spiral arms, and two azimuthally localized dark features that are also present in Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) total intensity images. Spectral energy distribution fitting further suggests that the radial gap extends to 1 au. The narrow, wedge-like shape of the dark features appears similar to predictions of shadows cast by an inner disk that is misaligned with respect to the outer disk. Using the Monte Carlo radiative transfer code HOCHUNCK3D, we construct a model of the disk that allows us to determine its physical properties in more detail. From the angular separation of the features, we measure the difference in inclination between the disks (45 degrees) and their major axes, PA = 140 degrees east of north for the outer disk, and 100 degrees for the inner disk. We find an outer-disk inclination of 25 degrees +/- 10 degrees from face-on, in broad agreement with the Wagner et al. measurement of 34 degrees. SPHERE data in J and H bands indicate a reddish disk, which indicates that HD 100453 is evolving into a young debris disk.

We present optical and near-infrared high-contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at Ha. The new GPI H-band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple data sets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 "c" varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD 100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.

L'objectif de cette thèse est d'étudier la formation de tourbillons dans la zone morte des disques protoplanétaires. Un code numérique 3D compressible a été mis au point et utilisé pour cette étude. Deux instabilités hydrodynamiques sont envisagées pour former les tourbillons: l'instabilité de Rossby et l'instabilité barocline.La première entraine la fragmentation d'une sur-densité annulaire en une chaîne de tourbillons qui se rattrapent les uns les autres et finissent par fusionner en un seul tourbillon qui reste stable sur de très longues durées lorsque son rapport d’aspect est suffisamment grand, et possède une structure quasi bidimensionnelle. En revanche, les tourbillons tridimensionnels de petits rapport d'aspect sont affectés par l’instabilité elliptique qui les détruits en quelques rotations. Seuls persistent ceux de grand rapport d'aspect.L'instabilité barocline, fondamentalement non linéaire, produit des tourbillons à partir de perturbations d'amplitude finies ; ces tourbillons sont ensuite amplifiés et fusionnent en tourbillons plus gros si le disque est stratifié de façon instable et s’il permet aussi le transfert de chaleur. Deux types de transfert thermique ont été envisagés pour étudier cette instabilité qui conduit alors à des différences significatives dans la structure des tourbillons formés. Le rapport d'aspect étant lié à la vorticité, l'amplification des tourbillons se traduit par une diminution de leur rapport d'aspect, et les rend donc sujet à l'instabilité elliptique. Cependant, ils ne sont pas détruit et gardent une structure tourbillonnaire grâce à l'amplification barocline. / The objective of this thesis is to study the formation of vortices in the dead-zone of protoplanetary disks. A 3D compressible numerical code has been performed and used for this study. Two hydrodynamical instabilities are considered for vortex formation: the Rossby wave instability and the baroclinic instability.The first one leads tp the fragmentation of an annular bump into a chain of vortices that catch one another and merge in a single vortex; this vortex remains stable on very long durations when its aspect ratio is large enough and has a quasi two-dimensional structure. In contrast, tridimensional small aspect ratios vortices are affected by the elliptical instability and are destroyed in a few rotation periods. Only vortices with large aspect ratios can survive.The baroclinic instability, a basically non-linear one, can produce vortices from small amplitude perturbations; these vortices are then amplified and merge in bigger vortices if the disk is unstably stratified and also permits heat transfer. Two types of heat transfer have been considered leading to significant differences in the structures of the resulting vortices. As aspect ratio and vorticity are strongly related, the baroclinic amplification reduces the aspect ratio and, so, make the vortex sensitive to the elliptical instability. However, such vortices are not destroyed and keep a vertical structure thanks to the baroclinic amplification.

La formation des planètes a lieu dans les disques protoplanétaires constitués de gaz et de poussières. Si ces dernières ne représentent que 1% de la masse totale du disque, elles jouent un rôle fondamental pour l’évolution chimique des disques en agissant comme catalyseurs pour la formation des molécules. Comprendre cette chimie est essentiel pour remonter aux conditions physiques initiales qui ont permis la naissance des planètes.Au cours de ma thèse, j’ai étudié la chimie à la surface des grains de poussières et son impact sur l’évolution chimique du nuage moléculaire, condition initale de la formation du disque, et du disque protoplanétaire. Grâce à des simulations numériques, à l’aide du code de chimie gaz-grain Nautilus, j’ai pu montrer l’importance des réactions de diffusion et des interactions gaz-grain pour les abondances des espèces en phase gazeuse. Les résultats du modèle couplés aux observations ont également mis en évidence les effets de la structure physique (température, densité, AV) sur la distribution des molécules dans les disques. / Planetary formation occurs in the protoplanetary disks of gas and dust. Although dust represents only 1% of the total disk mass, it plays a fundamental role in disk chemical evolution since it acts as a catalyst for the formation of molecules. Understanding this chemistry is therefore essential to determine the initial conditions from which planets form.During my thesis, I studied grain-surface chemistry and its impact on the chemical evolution of molecular cloud, initial condition for disk formation, and protoplanetary disk. Thanks to numerical simulations, using the gas-grain code Nautilus, I showed the importance of diffusion reactions and gas-grain interactions for the abundances of gas-phase species. Model results combined with observations also showed the effects of the physical structure (in temperature, density, AV) on the molecular distribution in disks.

This thesis presents a study of the sub-light year regions of Active Galactic Nuclei (AGN). These environments contain accretion discs that orbit a central super-massive black hole. The luminosity of the AGN inner regions varies over time across all wavelengths with variability at longer wavelengths lagging behind that at shorter wavelengths. Since the AGN themselves are too remote and too compact to resolve directly, I exploit these time lags to infer the physical characteristics of the accretion disc and surrounding gas clouds that emit broad emission lines. These characteristics include the inclination and temperature profile of the accretion disc, and the shape (or light curve) of the luminosity fluctuations that drive the accretion disc variability. This thesis details the work in the first author papers of Starkey et al. (2016, 2017), in which I detail the statistical code, CREAM (Continuum REverberting AGN Markov Chain Monte Carlo), that I developed to analyse AGN accretion disc variability. I apply the code to a set of AGN light curve observations of the Seyfert 1 galaxy NGC 5548 by the AGN STORM collaboration (De Rosa et al., 2015; Edelson et al., 2015; Fausnaugh et al., 2016a; Goad et al., 2016; Starkey et al., 2017). I also present work detailing my variability analysis of the Seyfert galaxies NGC 6814, NGC 2617, MCG 08-11-11 and NGC 4151. This work has contributed to the analysis presented in (Troyer et al. 2016, Fausnaugh et al. submitted). I also investigate the implications of a twin accretion disc structure (Nealon et al., 2015) on the disc time lag measurements across near UV and optical wavelengths. I finish by detailing a modification to CREAM that allows it to merge continuum light curves observed in a common filter, but taken by multiple telescopes with different calibration and instrumental effects to consider.

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

The number of extrasolar planet discoveries has increased dramatically over the last 15 years. Nearly 700 exoplanets have currently been observed through a variety of observation techniques. Most of the currently documented exoplanets differ greatly from the planets in our own Solar System, with various combinations of eccentric orbits, short orbital periods, and masses many times that of Jupiter. More recently, planets belonging to a new class of `distant gas giants' have also been discovered with orbits of 30 to 100 times that of Jupiter. The wide variety of different planet formation outcomes stem from a complex interplay between gravitational interactions, hydrodynamic interactions and competitive accretion among the planets that is not yet fully understood. Simulations performed using a series of modifications to an existing, widely used hydrodynamic code (FARGO) are presented. The main goal is to develop a more rigorous and robust gas accretion scheme that is valid and consistent for the ranges of exolanetary gas giant masses, eccentricities and semimajor axes that have been observed to better understand the mechanisms involved in their formation. The resulting scheme is a more robust and accurate prescription for gas accretion onto planetary cores in a manner that is mostly resolution independent and valid over a large range of masses (less than an Earth mass to multiple Jupiter masses). The modified scheme accounts for multiple, competing, dynamic accretion mechanisms (including atmospheric effects) and their associated time scales between an arbitrary number of protoplanets. This updated accretion scheme provides a means for exploring the entire formation process of gas giants out of a variety of initial conditions in a self-consistent manner. The modifications made to the code as well as simulation results will be discussed and explored.