For the calculation of complex neutral/ionized gas phase chemical equilibria, we present a semi-analytical versatile and efficient computer program, called FastChem. The applied method is based on the solution of a system of coupled nonlinear (and linear) algebraic equations, namely the law of mass action and the element conservation equations including charge balance, in many variables. Specifically, the system of equations is decomposed into a set of coupled nonlinear equations in one variable each, which are solved analytically whenever feasible to reduce computation time. Notably, the electron density is determined by using the method of Nelder and Mead at low temperatures. The program is written in object-oriented C++ which makes it easy to couple the code with other programs, although a stand-alone version is provided. FastChem can be used in parallel or sequentially and is available under the GNU General Public License version 3 at https://github.com/exoclime/FastChem together with several sample applications. The code has been successfully validated against previous studies and its convergence behavior has been tested even for extreme physical parameter ranges down to 100 K and up to 1000 bar. FastChem converges stable and robust in even most demanding chemical situations, which posed sometimes extreme challenges for previous algorithms.

During solar flares fast electron beams generated in the solar corona are non-thermal radio sources in terms of type III bursts. Sometimes they can enter into the interplanetary space, where they can be observed by in-situ measurements as it is done e.g. by the WIND spacecraft. On the other hand, they can be the source of non-thermal X-ray radiation as e.g. observed by RHESSI, if they precipitate toward the dense chromosphere due to bremsstrahlung. Since these energetic electrons are generated in the corona and observed at another site, the study of transport of such electrons in the corona and interplanetary space is of special interest. The transport of electrons is influenced by the global magnetic and electric field as well as local Coulomb collisions with the particles in the background plasma.

The dust formed in extended circumstellar envelopes of long-period variables and Miras has a strong influence on the envelope dynamics. A radiatively driven instability caused by the formation of dust leads to the development of an autonomous dynamics characterised by a set of distinct frequencies. We study the interplay between the envelope's internal dynamics and an external excitation by a pulsating star.

Long-period variables (LPVs) and Miras are highly evolved stars at the upper region of the asymptotic giant branch. They exhibit a pronounced variability of their luminosity with a more or less well defined period and suffer large mass loss in the form of stellar winds. Due to this extensive mass loss, LPVs and Miras are surrounded by extended circumstellar dust shells (CDSs). The dynamics of these envelopes is the result of a complex interplay via an external excitation by the pulsating central star, dust formation and radiative transfer. Our study is aimed at an understanding of the dynamics of CDSs around carbon-rich, standard luminous LPVs and Miras. These shells often show multiperiodicity with secondary periods as high as a few 10^4 d superimposed on a main period which is in the range of approximately 10^2-10^3 d. Such secondary periods may be caused at least in part by the presence of dust. We consider an excitation of the CDSs by either a harmonic force, provided by the oscillation of the central star, or by a stochastic force with a continuous power spectrum. The resulting numerically computed dynamical behaviour of the shell is analysed with the help of Fourier analysis and stroboscopic maps. CDSs may be described as multioscillatory systems which are driven by the pulsating stars. A set of normal modes can be identified. The obtained periods of these modes are of the order of some 10^3 d, which is a characteristic timescale for dust nucleation, growth and elemental enrichment in the dust formation zone. Depending on the oscillation period and strength of the central star, the envelope reacts periodically, multi- periodically or irregularly.

We investigate the mass loss of highly evolved, low- and intermediate mass stars and stellar samples with subsolar metallicity. We give a qualitative as well as quantitative description which can be applied to LMC/SMC-type stellar populations. For that purpose we apply the same approach as we did for solar metallicity stars and calculate hydrodynamical wind models including dust formation with LMC and SMC abundances under consideration of an adapted model assumption. In particular, we improved the treatment of the radiative transfer problem in order to accommodate larger non-local contributions occurring with smaller opacities. For each wind model we determine an averaged mass-loss rate. The resulting, approximate mass-loss formulae are then applied to well-tested and calibrated stellar evolution calculations in order to quantify the stellar mass loss. The dynamical models for LMC and SMC metallicity result in mass-loss rates of the same order of magnitude as the solar metallicity models which is in this basic approach in agreement with observations. The hydrodynamical properties like e.g. the outflow velocity differ (for fixed C/O abundance ratio) noticeably, though. While critical luminosities of LMC and solar metallicity models fairly coincide, the SMC models need higher luminosities to develop dust-driven winds.

Substellar atmospheres are observed to be irregularly variable for which the formation of dust clouds is the most promising candidate explanation. The atmospheric gas is convectively unstable and, last but not least, colliding convective cells are seen as cause for a turbulent fluid field. Since dust formation depends on the local properties of the fluid, turbulence influences the dust formation process and may even allow the dust formation in an initially dust-hostile gas. A regime-wise investigation of dust forming substellar atmospheric situations reveals that the largest scales are determined by the interplay between gravitational settling and convective replenishment which results in a dust-stratified atmosphere. The regime of small scales is determined by the interaction of turbulent fluctuations. Resulting lane-like and curled dust distributions combine to larger and larger structures. We compile necessary criteria for a subgrid model in the frame of large scale simulations as result of our study on small scale turbulence in dust forming gases.

Dust formation in brown dwarf atmospheres is studied by utilising a model for driven turbulence in the mesoscopic scale regime. We apply a pseudo-spectral method where waves are created and superimposed within a limited wavenumber interval. The turbulent kinetic energy distribution follows the Kolmogoroff spectrum which is assumed to be the most likely value. Such superimposed, stochastic waves may occur in a convectively active environment. They cause nucleation fronts and nucleation events and thereby initiate the dust formation process which continues until all condensible material is consumed. Small disturbances are found to have a large impact on the dust forming system. An initially dust-hostile region, which may originally be optically thin, becomes optically thick in a patchy way showing considerable variations in the dust properties during the formation process. The dust appears in lanes and curls as a result of the interaction with waves, i.e. turbulence, which form larger and larger structures with time. Aiming on a physical understanding of the variability of brown dwarfs, related to structure formation in substellar atmospheres, we work out first necessary criteria for small-scale closure models to be applied in macroscopic simulations of dust forming astrophysical systems.

In this paper, we present mean gas and dust opacities relevant to the physical conditions typical of protoplanetary discs. As the principal absorber for temperatures below ~1,500 K, we consider spherical and aggregate dust particles of various sizes, chemical structure, and porosity, consisting of ice, organics, troilite, silicates, and iron. For higher temperatures, ions, atoms, molecules, and electrons are included as the main opacity sources. Rosseland and Planck mean opacities are calculated for temperatures between 5 K and 10,000 K and gas densities ranging from 10^{-18} g/ccm to 10^{-7} g/ccm. The dependence on the adopted model of dust grains is investigated. We compare our results with recent opacity tables and show how different opacity models affect the calculated hydrodynamical structure of accretion discs.