User:Henrikjonsson/Sandbox

The photosphere denotes those solar or stellar surface layers from which optical radiation escapes. These stellar outer layers can be modeled by different computer programs. Often, calculated models are used, together with other programs, to calculate synthetic spectra for stars. For example, in varying the assumed abundance of a chemical element, and comparing the synthetic spectra to observed ones, the abundance of that element in that particular star can be determined. As computers have evolved, the complexity of the models has deepened, becoming more realistic in including more physical data and excluding more of the simplifying assumptions. This evolution of the models has also made them applicable to different kinds of stars.

Local thermodynamic equilibrium
This assumption (LTE) means that within any local computational volume, the state of thermodynamical equilibrium is assumed:
 * The inflow of radiation is determined by a blackbody spectrum set by the local temperature only. This radiation then interacts with the matter inside the volume.
 * The number of atoms or molecules occupying different excited energy states is determined by the Maxwell-Boltzmann distribution. This distribution is determined by the atomic excitation energies, and the local temperature.
 * The number of atoms in different ionization states is determined by the Saha equation. This distribution is determined by the atomic ionization energy, and the local temperature.

Plane parallel and spherical atmospheres
A common simplifying assumption is that the atmosphere is plane parallel, meaning that physical variables depend on one space coordinate only: the vertical depth (i.e., one assumes that we see the stellar atmosphere “head-on”, ignoring the curved portions towards the limbs). In stars where the photosphere is relatively thick compared to the stellar diameter, this is not a good approximation and an assumption of a spherical atmosphere is more appropriate.

Hydrostatic equilibrium
This means that the star is currently not undergoing any radical changes in structure involving large scale pulsations, flows or mass loss.

Mixing length and microturbulence
This assumption means that the convective motions in the atmosphere are described by the mixing-length theory, modeled as parcels of gas rising and disintegrating. To account for some of the small-scale effects in convective motions, a parameter called microturbulence is often used. The microturbulence corresponds to the motions of atoms or molecules on scales smaller than the photon mean free path.

Different methods of treating opacity
To fully model the photosphere one would need to include every absorption line of every element present. This is not feasible because it would be computationally extremely demanding, and also all spectra are not fully known. Therefore one needs to simplify the treatment of opacity. Methods used in photospheric models include: Opacity sampling means that the radiative transfer is evaluated for a number of optical wavelengths spread across the interesting parts of the spectrum. Although the model would improve with more frequencies included, opacity sampling uses as few as practical, to still get a realistic model, thereby minimizing calculation time. In using opacity distribution functions, the spectra are divided into subsections, within which the absorption probabilities are rearranged and simplified to one smooth function. Similar to the opacity sampling method, this is improved by adding more intervals but at the cost of prolonging the computation time.
 * Opacity sampling (OS)
 * Opacity distribution functions (ODF)

Different models
There are several different computer codes available modeling stellar photospheres. Some of them are described here and some of them are linked under “External links” below.

ATLAS
The ATLAS code was originally presented in 1970 by Robert Kurucz using the assumption of LTE and hydrostatic and plane parallel atmospheres. Since the source code is publicly available on the web, it has been amended by different persons numerous times over the years and nowadays exists in many versions. There are both plane parallel and spherical versions as well as those using opacity sampling or opacity distribution functions.

MARCS
The MARCS (Model Atmospheres in Radiative and Convective Scheme) code was originally presented in 1975 by Bengt Gustafsson and others. The original code simulated stellar spectra assuming the atmosphere to be in hydrostatic equilibrium, plane parallel, with convection described by mixing-length theory. The evolution of the code has since involved better modeling of the line opacity (opacity sampling instead of opacity distribution functions), spherical modeling and including an increasing number of physical data. Nowadays a large grid of different models is available on the web.

PHOENIX
The PHOENIX code is a model “risen from the ashes” from an earlier model called SNIRIS and mainly developed by Peter Hauschildt from 1992 onwards; it is regularly updated and made available on the web. Currently it uses spherical atmospheres and allows deviations from LTE for several elements but mostly assumes LTE.

3D hydrodynamical models
There are efforts to construct models not assuming LTE, and/or computing the detailed hydrodynamic motions instead of hydrostatic assumptions. These models are physically more realistic but also require more physical data such as cross sections and probabilities for different atomic processes. Such models are computationally rather demanding, and have not yet reached a stage of broader distribution.