Structure Of Circumstellar Disks Research Articles

Multi-wavelength VLTI study of the puffed-up inner rim of a circumbinary disc

Context. The presence of stable, compact circumbinary discs of gas and dust around post-asymptotic giant branch (post-AGB) binary systems has been well established. We focus on one such system, IRAS 08544−4431. Aims. We present an interferometric multi-wavelength analysis of the circumstellar environment of IRAS 08544−4431. The aim is to constrain different contributions to the total flux in the H-, K-, L-, and N-bands in the radial direction. Methods. The data obtained with the three current instruments on the Very Large Telescope Interferometer (VLTI), VLTI/PIONIER, VLTI/GRAVITY, and VLTI/MATISSE, range from the near-infrared, where the post-AGB star dominates, to the mid-infrared, where the disc dominates. We fitted the following two geometric models to the visibility data to reproduce the circumbinary disc: a ring with a Gaussian width and a flat disc model with a temperature gradient. The flux contributions from the disc, the primary star (modelled as a point source), and an over-resolved component were recovered along with the radial size of the emission, the temperature of the disc as a function of radius, and the spectral dependencies of the different components. Results. The trends of all visibility data were well reproduced with the geometric models. The near-infrared data were best fitted with a Gaussian ring model, while the mid-infrared data favoured a temperature gradient model. This implies that a vertical structure is present at the disc inner rim, which we attribute to a rounded puffed-up inner rim. The N-to-K size ratio is 2.8, referring to a continuous flat source, analogues to young stellar objects. Conclusions. By combining optical interferometric instruments operating at different wavelengths, we can resolve the complex structure of circumstellar discs and study the wavelength-dependent opacity profile. A detailed radial, vertical, and azimuthal structural analysis awaits a radiative transfer treatment in 3D to capture all non-radial complexity.

Nonideal MHD Simulation of HL Tau Disk: Formation of Rings

Abstract Recent high-resolution observations unveil ring structures in circumstellar disks. The origin of these rings has been widely investigated under various theoretical scenarios. In this work we perform global 3D nonideal MHD simulations including effects from both ohmic resistivity and ambipolar diffusion (AD) to model the HL Tau disk. The nonideal MHD diffusion profiles are calculated based on the global dust evolution calculation including sintering effects. Disk ionization structure changes dramatically across the snow line due to the change of dust size distribution close to the snow line of major volatiles. We find that accretion is mainly driven by disk wind. Gaps and rings can be quickly produced from different accretion rates across the snow line. Furthermore, AD leads to highly preferential accretion at the midplane, followed by magnetic reconnection. This results in a local zone of decretion that drains the mass in the field reconnection area, which leaves a gap and an adjacent ring just outside of it. Overall, under favorable conditions, both snow lines and nonideal MHD effects can lead to gaseous gaps and rings in protoplanetary disks.

Tracing planet-induced structures in circumstellar disks using molecular lines

Circumstellar disks are considered to be the birthplace of planets. Specific structures like spiral arms, gaps, and cavities are characteristic indicators of planet-disk interaction. Investigating these structures can provide insights into the growth of protoplanets and the physical properties of the disk. We investigate the feasibility of using molecular lines to trace planet-induced structures in circumstellar disks. Based on 3D hydrodynamic simulations of planet-disk interactions, we perform self-consistent temperature calculations and produce N-LTE molecular line velocity-channel maps and spectra of these disks using our new N-LTE line radiative transfer code Mol3D. Subsequently, we simulate ALMA observations using the CASA simulator. We consider two nearly face-on inclinations, 5 disk masses, 7 disk radii, and 2 different typical pre-main-sequence host stars (T Tauri, Herbig Ae). We calculate up to 141 individual velocity-channel maps for five molecules/isotopoloques in a total of 32 rotational transitions to investigate the frequency dependence of the structures indicated above. We find that the majority of protoplanetary disks in our parameter space could be detected in the molecular lines considered. However, unlike the continuum case, gap detection is not straightforward in lines. For example, gaps are not seen in symmetric rings but are masked by the pattern caused by the global (Keplerian) velocity field. We identify specific regions in the velocity-channel maps that are characteristic of planet-induced structures. Simulations of high angular resolution molecular line observations demonstrate the potential of ALMA to provide complementary information about the planet-disk interaction as compared to continuum observations. In particular, the detection of planet-induced gaps is possible under certain conditions.(abridged)

An interferometric study of the post-AGB binary 89 Herculis

Binary post-AGB stars are interesting laboratories to study both the evolution of binaries as well as the structure of circumstellar disks. A multiwavelength high angular resolution study of the prototypical object 89 Herculis is performed with the aim of identifying and locating the different emission components seen in the SED. A large interferometric data set, collected over the past decade and covering optical and near-IR wavelengths, is analyzed with simple geometric models. Combining the interferometric constraints with the photometry and the optical spectra, we reassess the energy budget of the post-AGB star and its circumstellar environment. We report the first (direct) detection of a large (35-40%) optical circumstellar flux contribution and spatially resolve its emission region. Given this large amount of reprocessed and/or redistributed optical light, the fitted size of the emission region is rather compact and fits with(in) the inner rim of the circumbinary dust disk. This rim dominates our K band data through thermal emission and is rather compact, emitting significantly already at a radius of twice the orbital separation. We interpret the circumstellar optical flux as due to a scattering process, with the scatterers located in the extremely puffed-up inner rim of the disk and possibly also in a bipolar outflow seen pole-on. A non-LTE gaseous origin in an inner disk cannot be excluded but is considered highly unlikely. This direct detection of a significant amount of circumbinary light at optical wavelengths poses several significant questions regarding our understanding of both post-AGB binaries and the physics in their circumbinary disks. Although the identification of the source of emission/scattering remains inconclusive without further study on this and similar objects, the implications are manifold.

Tracing large-scale structures in circumstellar disks with ALMA

Planets are supposed to form in circumstellar disks. The gravitational potential of a planet perturbs the disk and leads to characteristic structures, i.e. spiral waves and gaps, in the disk's density profile. We perform a large-scale parameter study of the observability of these planet-induced structures in circumstellar disks with ALMA. On the basis of HD and MHD simulations, we calculated the disk temperature structure and (sub)mm images of these systems. These were used to derive simulated ALMA images. Because appropriate objects are frequent in Taurus, we focused on a distance of 140pc and a declination of 20{\deg}. The explored range of star-disk-planet configurations consists of 6 HD simulations (including magnetic fields and different planet masses), 9 disk sizes, 15 total disk masses, 6 different central stars, and two different grain size distributions. On almost all scales and in particular down to a scale of a few AU, ALMA is able to trace disk structures induced by planet-disk interaction or by the influence of magnetic fields on the wavelength range between 0.4 and 2.0mm. In most cases, the optimum angular resolution is limited by the sensitivity. However, within the range of typical masses of protoplanetary disks (0.1-0.001Msun) the disk mass has a minor impact on the observability. It is possible to resolve disks down to 2.67e-6Msun and trace gaps induced by a planet with M_p/M_s = 0.001 in disks with 2.67e-4Msun with a signal-to-noise ratio greater than three. The central star has a major impact on the observability of gaps, as well as the considered maximum grainsize of the dust in the disk. In general, it is more likely to trace planet-induced gaps in our MHD models, because gaps are wider in the presence of magnetic fields. We also find that zonal flows resulting from MRI create gap-like structures in the disk's re-emission radiation, which are observable with ALMA.

EX Lupi in quiescence

EX Lup is the prototype of EXors, a subclass of low-mass pre-main sequence stars whose episodic eruptions are attributed to temporarily increased accretion. In quiescence the optical and near-infrared properties of EX Lup cannot be distinguished from those of normal T Tau stars. Here we investigate whether it is the circumstellar disk structure which makes EX Lup an atypical Class II object. During outburst the disk might undergo structural changes. Our characterization of the quiescent disk is intended to serve as a reference to study the physical changes related to one of EX Lupi's strongest known eruptions in 2008 Jan-Sep. We searched the literature for photometric and spectroscopic observations including ground-based, IRAS, ISO and Spitzer data. After constructing the optical-infrared spectral energy distribution (SED), we compared it with the typical SEDs of other young stellar objects and modeled it using the Monte Carlo radiative transfer code RADMC. A mineralogical decomposition of the 10 micron silicate emission feature and also the description of the optical and near-infrared spectra were performed. The SED is in general similar to that of a typical T Tauri star, though above 7 micron EX Lup emits higher flux. The quiescent phase data suggest low level variability in the optical-mid-infrared domain. Integrating the optical and infrared fluxes we derived a bolometric luminosity of 0.7 L_Sun. The 10 micron silicate profile could be fitted by a mixture consisting of amorphous silicates, no crystalline silicates were found. A modestly flaring disk model with a total mass of 0.025 M_Sun and an outer radius of 150 AU was able to reproduce the observed SED. The derived inner radius of 0.2 AU is larger than the sublimation radius, and this inner gap sets EX Lup aside from typical T Tauri stars.

Gas and Dust Emission at the Outer Edge of Protoplanetary Disks

We investigate the apparent discrepancy between gas and dust outer radii derived from millimeter observations of protoplanetary disks. Using 230 and 345 GHz continuum and CO J = 3–2 data from the Submillimeter Array for four nearby disk systems (HD 163296, TW Hydrae, GM Aurigae, and MWC 480), we examine models of circumstellar disk structure and the effects of their treatment of the outer disk edge. We show that for these disks, models described by power laws in surface density and temperature that are truncated at an outer radius are incapable of reproducing both the gas and dust emission simultaneously: the outer radius derived from the dust continuum emission is always significantly smaller than the extent of the molecular gas disk traced by CO emission. However, a simple model motivated by similarity solutions of the time evolution of accretion disks that includes a tapered exponential edge in the surface density distribution (and the same number of free parameters) does much better at reproducing both the gas and dust emission. While this analysis does not rule out the disparate radii implied by the truncated power law models, a realistic alternative disk model, grounded in the physics of accretion, provides a consistent picture for the extent of both the gas and dust.

Infrared signatures and models of circumstellar dust disks

This contribution describes theoretical models for the temperature and density structure of circumstellar disks and what properties of disks may be determined from the analysis of disk images and spectral energy distributions. I will summarize the contributions from several groups that have developed publicly available radiation transfer codes and grids of models of disk spectra. In addition, tools for fitting spectra from large datasets, such as the growing Spitzer archive, are now available.

Detecting protoplanets with ALMA

Theoretical investigations show that planet-disk interactions cause structures in circumstellar disks, which are usually much larger in size than the planet itself and thus more easily detectable. The specific result of planet-disk interactions depends on the evolutionary stage of the disk. Exemplary signatures of planets embedded in disks are gaps and spiral density waves in the case of young, gas-rich protoplanetary disks and characteristic asymmetric density patterns in debris disks. Numerical simulations convincingly demonstrate that high-resolution imaging performed with observational facilities which are already available or will become available in the near future will allow to trace these “fingerprints” of planets in protoplanetary and debris disks. These observations will provide a deep insight into specific phases of the formation and early evolution of planets in circumstellar disks. In this context, the Atacama Large Millimeter Array (ALMA) will play a crucial role by allowing to trace features in disks which are indicative for various stages of the formation and early evolution of planets in circumstellar disks.

The effects of X-ray photoionization and heating on the structure of circumstellar discs

We present the results of a theoretical study investigating the effects of photoionization and heating by X‐rays on discs around low‐mass stars. In particular, we address the question of whether or not X‐rays can drive a disc wind. First, we construct a one‐dimensional ‘quasi‐hydrostatic’ model, which solves for the vertical structure introduced by X‐ray heating. We consider uniform X‐ray illumination of the disc, but the X‐ray fluxes required to heat the disc significantly are much greater than those seen by recent observations. When the model is extended to consider heating from a central X‐ray source, we find that the one‐dimensional model is only valid very close to the star. We extend our analysis to consider a simple two‐dimensional model, treating the disc as a two‐layered structure and solving for its density profile self‐consistently. For T Tauri stars, we are able to set a crude upper limit on the mass‐loss rate that can be driven by X‐ray photoevaporation, with a value of ≃10−13 g cm‐2 s−1. Our model is designed to maximize this value, and most likely overestimates it significantly. However, we still find a mass‐loss rate which is less than that found in studies of ultraviolet photoevaporation. We conclude that in the presence of a significant ultraviolet field, X‐ray driven disc winds are unlikely to play a significant role in the evolution of discs around low‐mass stars.

NICMOS Coronagraphy

ABSTRACTThe NICMOS coronagraphic capability is a unique resource. It has been used to search for and/or confirm the presence of planetary or brown dwarf companions around nearby stars, to probe the structure of circumstellar disks, and to study the host environment of quasars with damped Lyα systems. These and other programs are expected to continue with NICMOS, using the recently installed NICMOS cooling system. In this paper we present an operational overview of NICMOS coronagraphy from its inception on 1998 February 6 through 1998 December 18. During this time interval, which corresponds to H ubble Space Telescope (HST) cycles 7 and 7N, coronagraphic observations of 329 targets for 22 programs were obtained. The flight software was enhanced to locate the position of the hole as well as the target during onboard acquisitions. Coronagraphic images are affected by changes in the point‐spread function (PSF) due to thermally induced HST focus changes, NICMOS cold‐mask movement, and light scatter around the coronagraphic hole. Proper PSF subtraction requires a well‐planned observing strategy that is appropriate for the target, and a matched PSF target. For example, for point sources, observations with the same filter taken back‐to‐back in the same orbit with a roll of the spacecraft between observations yields good results. For circumstellar disks, contemporaneous observations of a matched PSF target are required.

The Formation and Structure of Circumstellar Disks

AbstractSeveral theories have been proposed to explain the origin of Be star disks. Among them are Wind-Compressed Disks, accretion disks, decretion disks, and “explosive” ejections. In reviewing these mechanisms, I first concentrate on the current status of the Wind-Compressed Disk model. In particular, I discuss how non-radial forces may prevent disk formation and then discuss various physical effects that may restore the disk. Second, I examine the observational evidence and what it tells us about the structure of the disk. Of particular interest is evidence in favor of Keplerian disks. Finally, I discuss theories for Keplerian disk formation and some of the constraints such theories must satisfy.