Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFastChem: A computer program for efficient complex chemical equilibrium calculations in the neutral/ionized gas phase with applications to stellar and planetary atmospheres
90
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
Interactions between the winds of stars and the magnetospheres and atmospheres of planets involve many processes, including the acceleration of particles, heating of upper atmospheres, and a diverse range of atmospheric loss processes. Winds remove angular momentum from their host stars causing rotational spin-down and a decay in magnetic activity, which protects atmospheres from erosion. While wind interactions are strongly influenced by the X-ray and ultraviolet activity of the star and the chemical composition of the atmosphere, the role of planetary magnetic fields is unclear. In this chapter, I review our knowledge of the properties and evolution of stellar activity and winds and discuss the influences of these processes on the long term evolution of planetary atmospheres. I do not consider the large number of important processes taking place at the surfaces of planets that cause exchanges between the atmosphere and the planets interior.
The evolution of stars and planets is mostly controlled by the properties of their atmosphere. This is particularly true in the case of exoplanets close to their stars, for which one has to account both for an (often intense) irradiation flux, and from an intrinsic flux responsible for the progressive loss of the inner planetary heat. The goals of the present work are to help understanding the coupling between radiative transfer and advection in exoplanetary atmospheres and to provide constraints on the temperatures of the deep atmospheres. This is crucial in assessing whether modifying assumed opacity sources and/or heat transport may explain the inflated sizes of a significant number of giant exoplanets found so far. I use a simple analytical approach inspired by Eddingtons approximation for stellar atmospheres to derive a relation between temperature and optical depth valid for plane-parallel static grey atmospheres which are both transporting an intrinsic heat flux and receiving an outer radiation flux. The model is parameterized as a function of mean visible and thermal opacities, respectively. The model is shown to reproduce relatively well temperature profiles obtained from more sophisticated radiative transfer calculations of exoplanetary atmospheres. It naturally explains why a temperature inversion (stratosphere) appears when the opacity in the optical becomes significant compared to that in the infrared. I further show that the mean equivalent flux (proportional to T^4) is conserved in the presence of horizontal advection on constant optical depth levels. This implies with these hypotheses that the deep atmospheric temperature used as outer boundary for the evolution models should be calculated from models pertaining to the entire planetary atmosphere, not from ones that are relevant to the day side or to the substellar point. In these conditions, present-day models yield deep temperatures that are ~1000K too cold to explain the present size of planet HD 209458b. An tenfold increase in the infrared to visible opacity ratio would be required to slow the planetary cooling and contraction sufficiently to explain its size. However, the mean equivalent flux is not conserved anymore in the presence of opacity variations, or in the case of non-radiative vertical transport of energy: The presence of clouds on the night side or a downward transport of kinetic energy and its dissipation at deep levels would help making the deep atmosphere hotter and may explain the inflated sizes of giant exoplanets.
The Milky Way Galaxy is literally teeming with exoplanets; thousands of planets have been discovered, with thousands more planet candidates identified. Terrestrial-like planets are quite common around other stars, and are expected to be detected in large numbers in the future. Such planets are the primary targets in the search for potentially habitable conditions outside the solar system. Determining the atmospheric composition of exoplanets is mandatory to understand their origin and evolution, as atmospheric processes play crucial roles in many aspects of planetary architecture. In this work we construct and exploit a 1D radiative transfer model based on the discrete-ordinates method in plane-parallel geometry. Radiative results are linked to a convective flux that redistributes energy at any altitude producing atmospheric profiles in radiative-convective equilibrium. The model has been applied to a large number (6250) of closely dry synthetic ce{CO2} atmospheres, and the resulting pressure and thermal profiles have been interpreted in terms of parameter variability. Although less accurate than 3D general circulation models, not properly accounting for e.g., clouds and atmospheric and ocean dynamics, 1D descriptions are computationally inexpensive and retain significant value by allowing multidimensional parameter sweeps with relative ease.
There are many open questions about prebiotic chemistry in both planetary and exoplanetary environments. The increasing number of known exoplanets and other ultra-cool, substellar objects has propelled the desire to detect life and prebiotic chemistry outside the solar system. We present an ion-neutral chemical network constructed from scratch, Stand2015, that treats hydrogen, nitrogen, carbon and oxygen chemistry accurately within a temperature range between 100 K and 30000 K. Formation pathways for glycine and other organic molecules are included. The network is complete up to H6C2N2O3. Stand2015 is successfully tested against atmospheric chemistry models for HD209458b, Jupiter and the present-day Earth using a simple 1D photochemistry/diffusion code. Our results for the early Earth agree with those of Kasting (1993) for CO2, H2, CO and O2, but do not agree for water and atomic oxygen. We use the network to simulate an experiment where varied chemical initial conditions are irradiated by UV light. The result from our simulation is that more glycine is produced when more ammonia and methane is present. Very little glycine is produced in the absence of any molecular nitrogen and oxygen. This suggests that production of glycine is inhibited if a gas is too strongly reducing. Possible applications and limitations of the chemical kinetics network are also discussed.
We describe the incorporation of polarized radiative transfer into the atmospheric radiative transfer modelling code VSTAR (Versatile Software for Transfer of Atmospheric Radiation). Using a vector discrete-ordinate radiative transfer code we are able to generate maps of radiance and polarization across the disc of a planet, and integrate over these to get the full-disc polarization. In this way we are able to obtain disc-resolved, phase-resolved and spectrally-resolved intensity and polarization for any of the wide range of atmopsheres that can be modelled with VSTAR. We have tested the code by reproducing a standard benchmark problem, as well as by comparing with classic calculations of the polarization phase curves of Venus. We apply the code to modelling the polarization phase curves of the hot Jupiter system HD 189733b. We find that the highest polarization amplitudes are produced with optically thick Rayleigh scattering clouds and these would result in a polarization amplitude of 27 ppm for the planetary signal seen in the combined light of the star and planet. A more realistic cloud model consistent with the observed transmission spectrum results is an amplitude of ~20 ppm. Decreasing the optical depth of the cloud, or making the cloud particles more absorbing, both have the effect of increasing the polarization of the reflected light but reducing the amount of reflected light and hence the observed polarization amplitude.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
