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.

Cosmic rays (CRs) are transported out of the galaxy by diffusion and advection due to streaming along magnetic field lines and resonant scattering off self-excited Magneto-Hydro-Dynamic (MHD) waves. Thus momentum is transferred to the plasma via the frozen-in waves as a mediator assisting the thermal pressure in driving a galactic wind. Galactic CRs (GCRs) are accelerated by shock waves generated in supernova remnants (SNRs), and they propagate from the disc into the halo. Therefore CR acceleration in the halo strongly depends on the inner disc boundary conditions. We performed hydrodynamical simulations of galactic winds in flux tube geometry appropriate for disc galaxies, describing the CR diffusive-advective transport in a hydrodynamical fashion (by taking appropriate moments of the Fokker-Planck equation) along with the energy exchange with self-generated MHD waves. Our time-dependent CR hydrodynamic simulations confirm that the evolution of galactic winds with feedback depends on the structure of the galactic halo. In case of a wind-structured halo, the wind breaks down after the last super nova (SN) has exploded. The mechanism described here offers a natural and elegant solution to explain the power-law distribution of CRs between the `knee' and the `ankle'. The transition will be naturally smooth, because the Galactic CRs accelerated at SN shocks will be `post-accelerated' by shocks generated at the inner boundary and travelling through the halo.

We modeled the transit signatures in the Lya line of a putative Earth-sized planet orbiting in the HZ of the M dwarf GJ436. We estimated the transit depth in the Lya line for an exo-Earth with three types of atmospheres: a hydrogen-dominated atmosphere, a nitrogen-dominated atmosphere, and a nitrogen-dominated atmosphere with an amount of hydrogen equal to that of the Earth. We calculated the in-transit absorption they would produce in the Lya line. We applied it to the out-of-transit Lya observations of GJ 436 obtained by the HST and compared the calculated in-transit absorption with observational uncertainties to determine if it would be detectable. To validate the model, we also used our method to simulate the deep absorption signature observed during the transit of GJ 436b and showed that our model is capable of reproducing the observations. We used a DSMC code to model the planetary exospheres. The code includes several species and traces neutral particles and ions. At the lower boundary of the DSMC model we assumed an atmosphere density, temperature, and velocity obtained with a hydrodynamic model for the lower atmosphere. We showed that for a small rocky Earth-like planet orbiting in the HZ of GJ436 only the hydrogen-dominated atmosphere is marginally detectable with the STIS/HST. Neither a pure nitrogen atmosphere nor a nitrogen-dominated atmosphere with an Earth-like hydrogen concentration in the upper atmosphere are detectable. We also showed that the Lya observations of GJ436b can be reproduced reasonably well assuming a hydrogen-dominated atmosphere, both in the blue and red wings of the Lya line, which indicates that warm Neptune-like planets are a suitable target for Lya observations. Terrestrial planets can be observed in the Lya line if they orbit very nearby stars, or if several observational visits are available.

Aims. The carbon-rich asymptotic giant branch (AGB) star TX Piscium (TX Psc) has been observed multiple times during multiple epochs and at different wavelengths and resolutions, showing a complex molecular CO line profile and a ring-like structure in thermal dust emission. We investigate the molecular counterpart in high resolution, aiming to resolve the ring-like structure and identify its origin. Methods. Atacama Large Millimeter/submillimeter Array (ALMA) observations have been carried out to map the circumstellar envelope (CSE) of TX Psc in CO(2-1) emission and investigate the counterpart to the ring-like dust structure. Results. We report the detection of a thin, irregular, and elliptical detached molecular shell around TX Psc, which coincides with the dust emission. This is the first discovery of a non-spherically symmetric detached shell, raising questions about the shaping of detached shells. Conclusions. We investigate possible shaping mechanisms for elliptical detached shells and find that in the case of TX Psc, stellar rotation of 2 km/s can lead to a non-uniform mass-loss rate and velocity distribution from stellar pole to equator, recreating the elliptical CSE. We discuss the possible scenarios for increased stellar momentum, enabling the rotation rates needed to reproduce the ellipticity of our observations, and come to the conclusion that momentum transfer of an orbiting object with the mass of a brown dwarf would be sufficient.

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.

Young terrestrial planets, when they are still embedded in a circumstellar disk, accumulate an atmosphere of nebula gas. The evolution and eventual evaporation of the protoplanetary disk affect the structure and dynamics of the planetary atmosphere. These processes, combined with other mass loss mechanisms, such as thermal escape driven by extreme ultraviolet and soft X-ray radiation (XUV) from the young host star, determine how much of the primary atmosphere, if anything at all, survives into later stages of planetary evolution. Our aim is to explore the structure and the dynamic outflow processes of nebula-accreted atmospheres in dependency on changes in the planetary environment. We integrate stationary hydrostatic models and perform time-dependent dynamical simulations to investigate the effect of a changing nebula environment on the atmospheric structure and the timescales on which the protoatmosphere reacts to these changes. We find that the behavior of the atmospheres strongly depends on the mass of the planetary core. For planets of about Mars-mass the atmospheric structure, and in particular the atmospheric mass, changes drastically and on very short timescales whereas atmospheres around higher mass planets are much more robust and inert.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

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)

Context. High-resolution numerical methods have been developed for nonlinear, discontinuous problems as they appear in simulations of astrophysical objects. One of the strategies applied is the concept of artificial viscosity. Aims. Grid-based numerical simulations ideally utilize problem-oriented grids in order to minimize the necessary number of cells at a given (desired) spatial resolution. We want to propose a modified tensor of artificial viscosity which is employable for generally comoving, curvilinear grids. Methods. We study a differential geometrically consistent artificial viscosity analytically and visualize a comparison of our result to previous implementations by applying it to a simple self-similar velocity field. We give a general introduction to artificial viscosity first and motivate its application in numerical analysis. Then we present how a tensor of artificial viscosity has to be designed when going beyond common static Eulerian or Lagrangian comoving rectangular grids. Results. We find that in comoving, curvilinear coordinates the isotropic (pressure) part of the tensor of artificial viscosity has to be modified metrically in order for it to fulfill all its desired properties.

The formation of chondrules in the protoplanetary nebulae causes many questions concerning the formation process, the source of energy for melting the rims, and the composition of the origin material. The aim of this work is to explore the heating of the chondrule in a single precursor as is typical for radiation hydrodynamical shock waves. We take into account the gas-particle friction for the duration of the shock transition and calculate the heat conduction into the chondrules. These processes are located in the protoplanetary nebulae at a region around 2.5 AU, which is considered to be the most likely place of chondrule formation. The present models are a first step towards computing radiative shock waves occurring in a particle-rich environment. We calculated the shock waves using one-dimensional, time-independent equations of radiation hydrodynamics involving realistic gas and dust opacities and gas-particle friction. The evolution of spherical chondrules was followed by solving the heat conduction equation on an adaptive grid. The results for the shock-heating event are consistent with the cosmochemical constraints of chondrule properties. The calculations yield a relative narrow range for density or temperature to meet the requested heating rates of $R > 10^4\,K\, h^{-1}$ as extracted from cosmochemical constraints. Molecular gas, opacities with dust, and a protoplanetary nebula with accretion are necessary requirements for a fast heating process. The thermal structure in the far post-shock region is not fully consistent with experimental constraints on chondrule formation since the models do not include additional molecular cooling processes.

Cosmic rays are transported out of the galaxy by diffusion and advection due to streaming along magnetic field lines and resonant scattering off self-excited MHD waves. Thus momentum is transferred to the plasma via the frozen-in waves as a mediator assisting the thermal pressure in driving a galactic wind. The bulk of the Galactic CRs are accelerated by shock waves generated in SNRs, a significant fraction of which occur in OB associations on a timescale of several $10^7$ years. We examine the effect of changing boundary conditions at the base of the galactic wind due to sequential SN explosions on the outflow. Thus pressure waves will steepen into shock waves leading to in situ post-acceleration of GCRs. We performed simulations of galactic winds in flux tube geometry appropriate for disk galaxies, describing the CR diffusive-advective transport in a hydrodynamical fashion along with the energy exchange with self-generated MHD waves. Our time-dependent CR hydrodynamic simulations confirm the existence of time asymptotic outflow solutions (for constant boundary conditions). It is also found that high-energy particles escaping from the Galaxy and having a power-law distribution in energy ($\propto E^{-2.7}$) similar to the Milky Way with an upper energy cut-off at $\sim 10^{15}$ eV are subjected to efficient and rapid post-SNR acceleration in the lower galactic halo up to energies of $10^{17} - 10^{18}$ eV by multiple shock waves propagating through the halo. The particles can gain energy within less than $3\,$kpc from the galactic plane corresponding to flow times less than $5\cdot 10^6\,$years. The mechanism described here offers a natural solution to explain the power-law distribution of CRs between the "knee" and the "ankle". The mechanism described here offers a natural and elegant solution to explain the power-law distribution of CRs between the "knee" and the "ankle".