Heating and cooling processes in disks

The research leading to these results has received funding from the European Union Seventh Framework Programme FP7-2011 under grant agreement no 284405. 10th Lecture from Summer School “Protoplanetary Disks: Theory and Modelling Meet Observations”This chapter summarises current theoretical concepts and methods to determine the gas temperature structure in protoplanetary disks by balancing all relevant heating and cooling rates. The processes considered are non-LTE line heating/cooling based on the escape probability method, photo-ionisation heating and recombination cooling, free-free heating/cooling, dust thermal accommodation and high-energy heating processes such as X-ray and cosmic ray heating, dust photoelectric and PAH heating, a number of particular follow-up heating processes starting with the UV excitation of H2, and the release of binding energy in exothermal reactions. The resulting thermal structure of protoplanetary disks is described and discussed.Publisher PDFPeer reviewe

Prospects of using simulations to study the photospheres of brown dwarfs

We discuss prospects of using multi-dimensional time-dependent simulations to study the atmospheres of brown dwarfs and extrasolar giant planets, including the processes of convection, radiation, dust formation, and rotation. We argue that reasonably realistic simulations are feasible, however, separated into two classes of local and global models. Numerical challenges are related to potentially large dynamic ranges, and the treatment of scattering of radiation in multi-D geometries.Comment: 6 pages, 3 figures, to appear in the Proceedings of the IAU Symposium 239 "Convection in Astrophysics", eds. F. Kupka, I.W. Roxburgh, and K.L. Cha

Detectability of dirty dust grains in brown dwarf atmospheres

Dust clouds influence the atmospheric structure of brown dwarfs, and they affect the heat transfer and change the gas-phase chemistry. However, the physics of their formation and evolution is not well understood. In this letter, we predict dust signatures and propose a potential observational test of the physics of dust formation in brown dwarf atmosphere based on the spectral features of the different solid components predicted by dust formation theory. A momentum method for the formation of dirty dust grains (nucleation, growth, evaporation, drift) is used in application to a static brown dwarf atmosphere structure to compute the dust grain properties, in particular the heterogeneous grain composition and the grain size. Effective medium and Mie theory are used to compute the extinction of these spherical grains. Dust formation results in grains whose composition differs from that of grains formed at equilibrium. Our kinetic model predicts that solid amorphous SiO2[s] (silica) is one of the most abundant solid component followed by amorphous MgSiO4[s] and MgSiO3[s], while SiO2[s] is absent in equilibrium models because it is a metastable solid. Solid amorphous SiO2[s] possesses a strong broad absorption feature centered at 8.7mum, while amorphous Mg2SiO4[s]/MgSiO3[s] absorb at 9.7mum beside other absorption features at longer wavelength. Those features at lambda < 15mum are detectable in absorption if grains are small (radius < 0.2mum) in the upper atmosphere as suggested by our model. We suggest that the detection of a feature at 8.7mum in deep infrared spectra could provide evidence for non-equilibrium dust formation that yields grains composed of metastable solids in brown dwarf atmospheres. This feature will shift towards 10mum and broaden if silicates (e.g. fosterite) are much more abundant.Comment: A&A Letter, accepte

Magnetic Flares on Asymptotic Giant Branch Stars

We investigate the consequences of magnetic flares on the surface of asymptotic giant branch (AGB) and similar stars. In contrast to the solar wind, in the winds of AGB stars the gas cooling time is much shorter than the outflow time. As a result, we predict that energetic flaring will not inhibit, and may even enhance, dust formation around AGB stars. If magnetic flares do occur around such stars, we expect some AGB stars to exhibit X-ray emission; indeed certain systems including AGB stars, such as Mira, have been detected in X-rays. However, in these cases, it is difficult to distinguish between potential AGB star X-ray emission and, e.g., X-ray emission from the vicinity of a binary companion. Analysis of an archival ROSAT X-ray spectrum of the Mira system suggests an intrinsic X-ray luminosity 2x10^{29} erg/sec and temperature 10^7 K. These modeling results suggest that magnetic activity, either on the AGB star (Mira A) or on its nearby companion (Mira B), is the source of the X-rays, but do not rule out the possibility that the X-rays are generated by an accretion disk around Mira B.Comment: ApJ, Accepted; revised version of astro-ph/020923

Tracing the Physical Conditions in Active Galactic Nuclei with Time-Dependent Chemistry

We present an extension of the code ProDiMo that allows for a modeling of processes pertinent to active galactic nuclei and to an ambient chemistry that is time dependent. We present a proof-of-concept and focus on a few astrophysically relevant species, e.g., H+, H2+ and H3+; C+ and N+; C and O; CO and H2O; OH+, H2O+ and H3O+; HCN and HCO+. We find that the freeze-out of water is strongly suppressed and that this affects the bulk of the oxygen and carbon chemistry occurring in AGN. The commonly used AGN tracer HCN/HCO+ is strongly time-dependent, with ratios that vary over orders of magnitude for times longer than 10^4 years. Through ALMA observations this ratio can be used to probe how the narrow-line region evolves under large fluctuations in the SMBH accretion rate. Strong evolutionary trends, on time scales of 10^4-10^8 years, are also found in species such as H3O+, CO, and H2O. These reflect, respectively, time dependent effects in the ionization balance, the transient nature of the production of molecular gas, and the freeze-out/sublimation of water.Comment: Accepted for publication at the Journal of Physical Chemistry A "Oka Festschrift: Celebrating 45 Years of Astrochemistry

Preface

Publisher PDFPeer reviewe

The Binarity of Eta Carinae and its Similarity to Related Astrophysical Objects

I examine some aspects of the interaction between the massive star Eta Carinae and its companion, in particular during the eclipse-like event, known as the spectroscopic event or the shell event. The spectroscopic event is thought to occur when near periastron passages the stellar companion induces much higher mass loss rate from the primary star, and/or enters into a much denser environment around the primary star. I find that enhanced mass loss rate during periastron passages, if it occurs, might explain the high eccentricity of the system. However, there is not yet a good model to explain the presumed enhanced mass loss rate during periastron passages. In the region where the winds from the two stars collide, a dense slow flow is formed, such that large dust grains may be formed. Unlike the case during the 19th century Great Eruption, the companion does not accrete mass during most of its orbital motion. However, near periastron passages short accretion episodes may occur, which may lead to pulsed ejection of two jets by the companion. The companion may ionize a non-negligible region in its surrounding, resembling the situation in symbiotic systems. I discuss the relation of some of these processes to other astrophysical objects, by that incorporating Eta Car to a large class of astrophysical bipolar nebulae.Comment: Updated version. ApJ, in pres

Continuum and line modelling of discs around young stars. I. 300000 disc models for Herschel/GASPS

We have combined the thermo-chemical disc code ProDiMo with the Monte Carlo radiative transfer code MCFOST to calculate a grid of ~300000 circumstellar disc models, systematically varying 11 stellar, disc and dust parameters including the total disc mass, several disc shape parameters and the dust-to-gas ratio. For each model, dust continuum and line radiative transfer calculations are carried out for 29 far IR, sub-mm and mm lines of [OI], [CII], 12CO and o/p-H2O under 5 inclinations. The grid allows to study the influence of the input parameters on the observables, to make statistical predictions for different types of circumstellar discs, and to find systematic trends and correlations between the parameters, the continuum fluxes, and the line fluxes. The model grid, comprising the calculated disc temperatures and chemical structures, the computed SEDs, line fluxes and profiles, will be used in particular for the data interpretation of the Herschel open time key programme GASPS. The calculated line fluxes show a strong dependence on the assumed UV excess of the central star, and on the disc flaring. The fraction of models predicting [OI] and [CII] fine-structure lines fluxes above Herschel/PACS and Spica/SAFARI detection limits are calculated as function of disc mass. The possibility of deriving the disc gas mass from line observations is discussed.Comment: accepted by MNRAS. 5 pages, 4 figures, 3 table

Radiation thermo-chemical models of protoplanetary disks II. Line diagnostics

Methods. We use the recently developed disk code ProDiMo to calculate the physico-chemical structure of protoplanetary disks and apply the Monte-Carlo line radiative transfer code RATRAN to predict observable line profiles and fluxes. We consider a series of Herbig Ae type disk models ranging from 10^-6 M_Sun to 2.2 10^-2 M_Sun (between 0.5 and 700 AU) to discuss the dependency of the line fluxes and ratios on disk mass for otherwise fixed disk parameters. Results. We find the [CII] 157.7 mum line to originate in LTE from the surface layers of the disk, where Tg > Td . The total emission is dominated by surface area and hence depends strongly on disk outer radius. The [OI] lines can be very bright (> 10^-16 W/m^2) and form in slightly deeper and closer regions under non-LTE conditions. The high-excitation [OI] 145.5 mum line, which has a larger critical density, decreases more rapidly with disk mass than the 63.2 mum line. Therefore, the [OI] 63.2 mum/145.5 mum ratio is a promising disk mass indicator, especially as it is independent of disk outer radius for Rout > 200 AU. CO is abundant only in deeper layers A_V >~ 0.05. For too low disk masses (M_disk <~10^-4 M_Sun) the dust starts to become transparent, and CO is almost completely photo-dissociated. For masses larger than that the lines are an excellent independent tracer of disk outer radius and can break the outer radius degeneracy in the [OI] 63.2 mum/[CII]157.7 mum line ratio. Conclusions. The far-IR fine-structure lines of [CII] and [OI] observable with Herschel provide a promising tool to measure the disk gas mass, although they are mainly generated in the atomic surface layers. In spatially unresolved observations, none of these lines carry much information about the inner, possibly hot regions < 30 AU.Comment: accepted for publication in A&

Radiation thermo-chemical models of protoplanetary disks I. Hydrostatic disk structure and inner rim

This paper introduces a new disk code, called ProDiMo, to calculate the thermo-chemical structure of protoplanetary disks and to interpret gas emission lines from UV to sub-mm. We combine frequency-dependent 2D dust continuum radiative transfer, kinetic gas-phase and UV photo-chemistry, ice formation, and detailed non-LTE heating & cooling balance with the consistent calculation of the hydrostatic disk structure. We include FeII and CO ro-vibrational line heating/cooling relevant for the high-density gas close to the star, and apply a modified escape probability treatment. The models are characterized by a high degree of consistency between the various physical, chemical and radiative processes, where the mutual feedbacks are solved iteratively. In application to a T Tauri disk extending from 0.5AU to 500AU, the models are featured by a puffed-up inner rim and show that the dense, shielded and cold midplane (z/r<0.1, Tg~Td) is surrounded by a layer of hot (5000K) and thin (10^7 to 10^8 cm^-3) atomic gas which extends radially to about 10AU, and vertically up to z/r~0.5. This layer is predominantly heated by the stellar UV (e.g. PAH-heating) and cools via FeII semi-forbidden and OI 630nm optical line emission. The dust grains in this "halo" scatter the star light back onto the disk which impacts the photo-chemistry. The more distant regions are characterized by a cooler flaring structure. Beyond 100AU, Tgas decouples from Tdust even in the midplane and reaches values of about Tg~2Td. Our models show that the gas energy balance is the key to understand the vertical disk structure. Models calculated with the assumption Tg=Td show a much flatter disk structure.Comment: 24 pages, 14 figures, 120 equations, accepted by A&A, download a high-resolution version from http://www.roe.ac.uk/~ptw/prodimo1_article.pd