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.

Despite all advances in multi-dimensional hydrodynamics, investigations of stellar evolution and stellar pulsations still depend on one-dimensional computations. The present work devises an alternative to the mixing length theory or turbulence models usually adopted for the modelling of convective transport in such studies. Assuming that the largest convective patterns generate the majority of convective transport, the convective velocity field is described using two parallel radial columns to represent up- and downstream flows. Horizontal exchange in the form of fluid flow and radiation over their connecting interface couples the two columns and allows a simple circulating motion. The main parameters of this convective description have a straightforward geometrical meaning, namely the diameter of the columns (corresponding to the size of the convective cells) and the ratio of cross section between up- and downdrafts. For this geometrical setup, the time-dependent solution of the equations of radiation hydrodynamics is computed from an implicit scheme which has the advantage of being not affected by the Courant-Friedrichs-Lewy time step limit. In order to demonstrate the approach, results for the example of convection zones in Cepheids are presented.

Galactic Cepheids in the vicinity of the Sun have a residual line-of-sight velocity, or gamma-velocity, which shows a systematic blueshift of about 2 km/s compared to an axisymmetric rotation model of the Milky Way. This term is either related to the space motion of the star and, consequently, to the kinematic structure of our Galaxy, or it is the result of the dynamical structure of the Cepheids' atmosphere. We aim to show that these residual gamma-velocities are an intrinsic property of Cepheids. We observed nine galactic Cepheids with the HARPS spectroscope, focusing specifically on 17 spectral lines. For each spectral line of each star, we computed the gamma-velocity (resp. gamma-asymmetry) as an average value of the interpolated radial velocity (resp. line asymmetry) curve. For each Cepheid in our sample, a linear relation is found between the gamma-velocities of the various spectral lines and their corresponding gamma-asymmetries, showing that residual gamma-velocities stem from the intrinsic properties of Cepheids. We also provide a physical reference to the stellar gamma-velocity: it should be zero when the gamma-asymmetry is zero. Following this definition, we provide very precise and physically calibrated estimates of the gamma-velocities for all stars of our sample. To understand this very subtle gamma-asymmetry effect, further numerical studies are needed. Cepheids' atmosphere are strongly affected by pulsational dynamics, convective flows, nonlinear physics, and complex radiative transport. Hence, all of these effects have to be incorporated simultaneously and consistently into the numerical models to reproduce the observed line profiles in detail.