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.

This third paper in the series modeling QSOs and AGN as clusters of accreting black holes studies the accretion flow within an externally supplied cluster. Significant radiation will be emitted by the cluster core, but the black holes in the outer halo, where the flow is considered spherically symmetric, will not contribute much to the overall luminosity of the source because of their large velocities relative to the infalling gas, and therefore their small accretion radii. As a result the scenario discussed in Paper I will refer to the cluster cores, rather than to entire clusters. This will steepen the high frequency region of the spectrum unless inverse Compton scattering is effective. In many cases accretion flow in the central part of the cluster will be optically thick to electron scattering resulting in a spectrum featuring optically thick radiative component in addition to power -law regimes. The fitting of these spectra to QSO and AGN observations is discussed, and application to 3C 273 is worked out as an example.

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.

We develop a simple model to predict the radial distribution of planetesimal formation. The model is based on the observed growth of dust to millimeter-sized particles, which drift radially, pile-up, and form planetesimals where the stopping time and dust-to-gas ratio intersect the allowed region for streaming instability-induced gravitational collapse. Using an approximate analytic treatment, we first show that drifting particles define a track in metallicity-stopping time space whose only substantial dependence is on the disk's angular momentum transport efficiency. Prompt planetesimal formation is feasible for high particle accretion rates (relative to the gas, (M) over dot(p)/(M) over dot greater than or similar to 3 x 10(-2) for alpha = 10(-2)), which could only be sustained for a limited period of time. If it is possible, it would lead to the deposition of a broad and massive belt of planetesimals with a sharp outer edge. Numerically including turbulent diffusion and vapor condensation processes, we find that a modest enhancement of solids near the snow line occurs for centimeter-sized particles, but that this is largely immaterial for planetesimal formation. We note that radial drift couples planetesimal formation across radii in the disk, and suggest that considerations of planetesimal formation favor a model in which the initial deposition of material for giant planet cores occurs well beyond the snow line.

It is currently believed that angular momentum transport in accretion disks is mediated by magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability (MRI). More than 15 years after its discovery, an accretion disk model that incorporates the MRI as the mechanism driving the MHD turbulence is still lacking. This dissertation constitutes the first in a series of steps towards establishing the formalism and methodology needed to move beyond the standard accretion disk model and incorporating the MRI as the mechanism enabling the accretion process. I begin by presenting a local linear stability analysis of a compressible, differentially rotating flow and addressing the evolution of the MRI beyond the weak-field limit when magnetic tension forces due to strong toroidal fields are considered. Then, I derive the first formal analytical proof showing that, during the exponential growth of the instability, the mean total stress produced by correlated MHD fluctuations is positive and leads to a net outward flux of angular momentum. I also show that some characteristics of the MHD stresses that are determined during this initial phase are roughly preserved in the turbulent saturated state observed in local numerical simulations. Motivated by these results, I present the first mean-field MHD model for angular momentum transport driven by the MRI that is able to account for a number of correlations among stresses found in local numerical simulations. I point out the relevance of a new type of correlation that couples the dynamical evolution of the Reynolds and Maxwell stresses and plays a key role in developing and sustaining the MHD turbulence. Finally, I address how the turbulent transport of angular momentum depends on the magnitude of the local shear. I show that turbulent MHD stresses in accretion disks cannot be described in terms of shear-viscosity.

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.

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第22248号 / 理博第4562号 / 新制||理||1655(附属図書館) / 京都大学大学院理学研究科物理学・宇宙物理学専攻 / (主査)准教授 野上 大作, 教授 嶺重 慎, 教授 長田 哲也 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DFAM