We present an implicit numerical method to solve the time-dependent equations of radiation hydrodynamics (RHD) in axial symmetry assuming hydrostatic equilibrium perpendicular to the equatorial plane (1+1D) of a gaseous disk. The equations are formulated in conservative form on an adaptive grid and the corresponding fluxes are calculated by a spacial second order advection scheme. Self-gravity of the disk is included by solving the Possion equation. We test the resulting numerical method through comparison with a simplified analytical solution as well as through the long term viscous evolution of protoplanetary disk when due to viscosity matter is transported towards the central host star and the disk depletes. The importance of the inner boundary conditions on the structural behaviour of disks is demonstrated with several examples.

The collapse of interstellar gaseous clouds towards a protostar leads to the formation of accretion disks around the central star. Such disks can be dynamically stable if they settle in an axisymmetric state. In this letter, we investigate the long-term stability of astrophysical viscous disks around various protostars. We apply an implicit numerical code which solves the equations of radiation hydrodynamics and treats turbulence-induced viscosity according to the $\alpha$-viscosity model. We show how the viscosity is related to the disk mass. A stability criterion to determine the maximum disk mass can be formulated. We analyse such instabilities for a variety of radial points with different orbital distances from the host star and discuss the feedback on the disk in the event of an unstable protoplanetary disk. Additionally, we examine the critical disk-mass for disks with variable outer boundaries and compare them to observations of protostellar disks in the Upper Scorpius OB Association and near the Lupus complex. We derive an easily applicable method to obtain an estimate for maximum disk masses when the outer disk radius ins known.

Context. In tight binary star systems, tidal interactions can significantly influence the rotational and orbital evolution of both stars, and therefore their activity evolution. This can have strong effects on the atmospheric evolution of planets that are orbiting the two stars. Aims. In this paper, we aim to study the evolution of stellar rotation and of X-ray and ultraviolet (XUV) radiation in tight binary systems consisting of two solar mass stars and use our results to study planetary atmosphere evolution in the habitable zones of these systems. Methods. We have applied a rotation model developed for single stars to binary systems, taking into account the effects of tidal interactions on the rotational and orbital evolution of both stars. We used empirical rotation-activity relations to predict XUV evolution tracks for the stars, which we used to model hydrodynamic escape of hydrogen dominated atmospheres. Results. When significant, tidal interactions increase the total amount of XUV energy emitted, and in the most extreme cases by up to factor of $\sim$50. We find that in the systems that we study, habitable zone planets with masses of 1~M$_\oplus$ can lose huge hydrogen atmospheres due to the extended high levels of XUV emission, and the time that is needed to lose these atmospheres depends on the binary orbital separation.For some orbital separations, and when the stars are born as rapid rotators, it is also possible for tidal interactions to protect atmospheres from erosion by quickly spinning down the stars. For very small orbital separations, the loss of orbital angular momentum by stellar winds causes the two stars to merge. We suggest that the merging of the two stars could cause previously frozen planets to become habitable due to the habitable zone boundaries moving outwards.

Planets that form early enough to be embedded in the circumstellar gas disk accumulate thick atmospheres of nebular gas. Models of these atmospheres need to specify the surface luminosity (i.e. energy loss rate) of the planet. This luminosity is usually associated with a continuous inflow of solid bodies, where the gravitational energy released from these bodies is the source of energy. However, if these bodies release energy in the atmosphere instead of at the surface, this assumption might not be justified. Our aim is to explore the interactions of infalling planetesimals with primordial atmospheres at an embedded phase of evolution. We investigate effects of atmospheric interaction on the planetesimals (mass loss) and the atmosphere (heating/cooling). We used atmospheric parameters from a snapshot of time-dependent evolution simulations for embedded atmospheres and simulated purely radial, infall events of siliceous planetesimals in a 1D, explicit code. We implemented energy transfer between friction, radiation transfer by the atmosphere and the body and thermal ablation; this gives us the possibility to examine the effects on the planetesimals and the atmosphere. We find that a significant amount of gravitational energy is indeed dissipated into the atmosphere, especially for larger planetary cores, which consequently cannot contribute to the atmospheric planetary luminosity. Furthermore, we examine that planetesimal infall events for cores, $M_\mathrm{C} > 2$M$_{\oplus}$, which actually result in a local cooling of the atmosphere; this is totally in contradiction with the classical model.

Terrestrial planets formed within gaseous protoplanetary disks can accumulate significant hydrogen envelopes. The evolution of such an atmosphere due to XUV driven evaporation depends on the activity evolution of the host star, which itself depends sensitively on its rotational evolution, and therefore on its initial rotation rate. In this letter, we derive an easily applicable method for calculating planetary atmosphere evaporation that combines models for a hydrostatic lower atmosphere and a hydrodynamic upper atmosphere. We show that the initial rotation rate of the central star is of critical importance for the evolution of planetary atmospheres and can determine if a planet keeps or loses its primordial hydrogen envelope. Our results highlight the need for a detailed treatment of stellar activity evolution when studying the evolution of planetary atmospheres.

We investigate the origin and loss of captured hydrogen envelopes from protoplanets between `sub-Earth'-like bodies of 0.1$M_{\oplus}$ up to `super-Earths' with 5$M_{\oplus}$ in the HZ of a Sun like G star, assuming their rocky cores had formed before the nebula dissipated. We model the gravitational accumulation of nebula gas around a core as a function of protoplanetary luminosity during accretion and calculate the resulting surface temperature by solving the hydrostatic structure equations for the protoplanetary nebula. Depending on nebular properties and resulting luminosities, for planetary bodies of 0.1--1$M_{\oplus}$ we obtain hydrogen envelopes with masses between $\sim 2.5\times 10^{19}$--$1.5\times 10^{26}$ g. For `super-Earths' with masses between 2--5$M_{\oplus}$ hydrogen envelopes within the mass range of $\sim 7.5\times 10^{23}$--$1.5\times 10^{28}$ g can be captured. To study the escape of these hydrogen-dominated protoatmospheres, we apply a hydrodynamic upper atmosphere model and calculate the loss rates due to the heating by the high XUV flux of the young star. Our results indicate that under most nebula conditions `sub-Earth' and Earth-mass planets can lose their envelopes by thermal escape during the first $100$ Myr after the disk dissipated. However, if a nebula has a low dust depletion factor or low accretion rates resulting in low protoplanetary luminosities, it is possible that even protoplanets with Earth-mass cores may keep their hydrogen envelopes during their whole lifetime. In contrast to lower mass protoplanets, `super-Earths' accumulate a huge amount of nebula gas and lose only tiny fractions of their primordial envelopes. Our results agree with the fact that Venus, Earth, and Mars are not surrounded by dense hydrogen envelopes, as well as with the recent discoveries of low density `super-Earths' that most likely could not get rid of their protoatmospheres.

Latest research in planet formation indicate that Mars formed within a few million years (Myr) and remained a planetary embryo that never grew to a more massive planet. It can also be expected from dynamical models, that most of Mars' building blocks consisted of material that formed in orbital locations just beyond the ice line which could have contained ~0.1-0.2 wt. % of H2O. By using these constraints, we estimate the nebula-captured and catastrophically outgassed volatile contents during the solidification of Mars' magma ocean and apply a hydrodynamic upper atmosphere model for the study of the soft X-ray and extreme ultraviolet (XUV) driven thermal escape of the martian protoatmosphere during the early active epoch of the young Sun. The amount of gas that has been captured from the protoplanetary disk into the planetary atmosphere is calculated by solving the hydrostatic structure equations in the protoplanetary nebula. Depending on nebular properties such as the dust grain depletion factor, planetesimal accretion rates and luminosities, hydrogen envelopes with masses >=3x10^{19} g to <=6.5x10^{22} g could have been captured from the nebula around early Mars. Depending of the before mentioned parameters, due to the planets low gravity and a solar XUV flux that was ~100 times stronger compared to the present value, our results indicate that early Mars would have lost its nebular captured hydrogen envelope after the nebula gas evaporated, during a fast period of ~0.1-7.5 Myr. After the solidification of early Mars' magma ocean, catastrophically outgassed volatiles with the amount of ~50-250 bar H2O and ~10-55 bar CO2 could have been lost during ~0.4-12 Myr, if the impact related energy flux of large planetesimals and small embryos to the planet's surface lasted long enough, that the steam atmosphere could have been prevented from condensing. If this was not the case... (continued)