We generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength domain of the absorption. Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the greenhouse effect saturates and any Earth-like planet has a maximal temperature which depends on the type of and distance to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the visible, could have at most a surface temperature of 1200-1300K irrespective of the nature of the greenhouse gas. We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the value of the wavelength above which the atmosphere becomes optically thick. Unless this wavelength is small and resides in the optical region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. This has the observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, saturation takes place by radiating through windows in the FIR.

We present a new implicit numerical algorithm for the calculation of the time dependent non-Local Thermodynamic Equilibrium of a gas in an external radiation field that is accurate, fast and unconditionally stable for all spatial and temporal increments. The method is presented as a backward difference scheme in 1-D but can be readily generalised to 3-D. We apply the method for calculating the evolution of ionisation domains in a hydrogen plasma with plane-parallel Gaussian density enhancements illuminated by sources of UV radiation. We calculate the speed of propagation of ionising fronts through different ambient densities and the interaction of such ionising fronts with density enhancements. We show that for a typical UV source that may be present in the early universe, the introduction of a density enhancement of a factor ~10 above an ambient density 10^{-4} atoms/cm^3 could delay the outward propagation of an ionisation front by millions of years. Our calculations show that within the lifetime of a single source (~a few million years), and for ambient intergalactic densities appropriate to redshifts z~6-20, neutral fractions of of ~10^{-3}-10^{-5} can be achieved within its zone of influence. We also present calculations which demonstrate that once started, ionisation will proceed very efficiently as multiple sources are subsequently introduced, even if the time between the appearence of such sources may be much longer than their lifetimes.

A method for the fast and accurate solution of the radiative transfer equation in plane-parallel media with coherent isotropic scattering is presented. This largely analytical approach uses the formalism of meromorphic functions in order to obtain the total solution, i.e. the angle and depth variation of the radiation field. A discretization of space and angle coordinates is not required. For the further application in the accretion disks a particular case of a finite slab whose the midplane is the symmetry plane is considered.