Model photosphere

Last updated
Some essential steps in model-atmosphere analysis for determining stellar abundances (Figure by Bengt Gustafsson, Astronomical Observatory, Uppsala). Model atmospheres Bengt Gustafsson.jpg
Some essential steps in model-atmosphere analysis for determining stellar abundances (Figure by Bengt Gustafsson, Astronomical Observatory, Uppsala).

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.

Contents

Common assumptions and computational methods

Local Thermodynamic Equilibrium (LTE)

This assumption (LTE) means that within any local computational volume, the state of thermodynamical equilibrium is assumed:

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.

Expanding atmospheres

Many stars lose mass in the form of a stellar wind. Especially for stars which are very hot (photospheric temperatures > 10,000 Kelvin) and very luminous, these winds can be so dense that major parts of the emergent spectrum are formed in an "expanding atmosphere", i.e. in layers that are moving outward with a high speed that can reach a few 1000 km/s.

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.

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, Roger Bell 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 "risen from the ashes" of an earlier code called SNIRIS and mainly developed by Peter Hauschildt (Hamburger Sternwarte) from 1992 onwards; it is regularly updated and made available on the web. It runs in two different spatial configuration modes: the "classic" one-dimensional mode, assuming spherical symmetry, and the three-dimensional mode. It allows for calculations for many different astrophysical objects, i.e. supernovae, novae, stars and planets. It considers scattering and dust and allows for non-LTE computations over many atomic species, plus LTE over atoms and molecules.

PoWR

The PoWR (Potsdam Wolf-Rayet) code is designed for expanding stellar atmospheres, i.e. for stars with a stellar wind. It has been developed since the 1990s by Wolf-Rainer Hamann and collaborators at the Universität Potsdam (Germany) especially for simulating Wolf-Rayet stars, which are hot stars with very strong mass loss. Adopting spherical symmetry and stationarity, the program computes the occupation numbers of the atomic energy states, including the ionization balance, in non-LTE, and consistently solves the radiative transfer problem in the comoving frame. The stellar wind parameters (mass-loss rate, wind speed) can be specified as free parameter, or, alternatively, calculated from the hydrodynamic equation consistently. As the PoWR code treats consistently the static and expanding layers of the stellar atmosphere, it is applicable for any types of hot stars. The code as such is not yet public, but large sets of models for Wolf-Rayet stars are available on the web.

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.

Applications of Model Photospheres

Model Atmospheres, while interesting in their own right, are frequently used as part of input recipes and tools for studying other astrophysical problems.

Stellar Evolution

As a result of stellar evolution, changes in the internal structure of stars manifest themselves in the photosphere.

Synthetic Spectra

Spectral synthesis programs (e.g. Moog (code)) often use previously generated model photospheres to describe the physical conditions (temperature, pressure, etc...) through which photons must travel to escape the stellar atmosphere. Together with a list of absorption lines and an elemental abundance table, spectral synthesis programs generate synthetic spectra. By comparing these synthetic spectra to observed spectra of distant stars, astronomers can determine the properties (temperature, age, chemical composition, etc...) of these stars.

See also

Stellar structure

Related Research Articles

Main sequence Continuous band of stars that appears on plots of stellar color versus brightness

In astronomy, the main sequence is a continuous and distinctive band of stars that appears on plots of stellar color versus brightness. These color-magnitude plots are known as Hertzsprung–Russell diagrams after their co-developers, Ejnar Hertzsprung and Henry Norris Russell. Stars on this band are known as main-sequence stars or dwarf stars. These are the most numerous true stars in the universe, and include the Earth's Sun.

Star Astronomical spheroid of plasma

A star is an astronomical object consisting of a luminous spheroid of plasma held together by its own gravity. The nearest star to Earth is the Sun. Many other stars are visible to the naked eye at night, but due to their immense distance from Earth they appear as fixed points of light in the sky. The most prominent stars are grouped into constellations and asterisms, and many of the brightest stars have proper names. Astronomers have assembled star catalogues that identify the known stars and provide standardized stellar designations. The observable universe contains an estimated 1022 to 1024 stars, but most are invisible to the naked eye from Earth, including all individual stars outside our galaxy, the Milky Way.

Spectroscopy Study involving matter and electromagnetic radiation

Spectroscopy is the study of the interaction between matter and electromagnetic radiation as a function of the wavelength or frequency of the radiation. In simpler terms, spectroscopy is the precise study of color as generalized from visible light to all bands of the electromagnetic spectrum; indeed, historically, spectroscopy originated as the study of the wavelength dependence of the absorption by gas phase matter of visible light dispersed by a prism. Matter waves and acoustic waves can also be considered forms of radiative energy, and recently gravitational waves have been associated with a spectral signature in the context of the Laser Interferometer Gravitational-Wave Observatory (LIGO).

Stellar classification Classification of stars based on their spectral characteristics

In astronomy, stellar classification is the classification of stars based on their spectral characteristics. Electromagnetic radiation from the star is analyzed by splitting it with a prism or diffraction grating into a spectrum exhibiting the rainbow of colors interspersed with spectral lines. Each line indicates a particular chemical element or molecule, with the line strength indicating the abundance of that element. The strengths of the different spectral lines vary mainly due to the temperature of the photosphere, although in some cases there are true abundance differences. The spectral class of a star is a short code primarily summarizing the ionization state, giving an objective measure of the photosphere's temperature.

Wolf–Rayet star Stars with unusual spectra showing prominent broad emission lines of highly ionised helium and nitrogen or carbon

Wolf–Rayet stars, often abbreviated as WR stars, are a rare heterogeneous set of stars with unusual spectra showing prominent broad emission lines of ionised helium and highly ionised nitrogen or carbon. The spectra indicate very high surface enhancement of heavy elements, depletion of hydrogen, and strong stellar winds. The surface temperatures of known Wolf-Rayet stars range from 20,000 K to around 210,000 K, hotter than almost all other kinds of stars. They were previously called W-type stars referring to their spectral classification.

Asteroseismology Study of oscillations in stars

Asteroseismology is the study of oscillations in stars. Because a star's different oscillation modes are sensitive to different parts of the star, they inform astronomers about the internal structure of the star, which is otherwise not directly possible from overall properties like brightness and surface temperature. Asteroseismology is closely related to helioseismology, the study of stellar oscillations specifically in the Sun. Though both are based on the same underlying physics, more and qualitatively different information is available for the Sun because its surface can be resolved.

Stellar structure Structure of stars

Stellar structure models describe the internal structure of a star in detail and make predictions about the luminosity, the color and the future evolution of the star. Different classes and ages of stars have different internal structures, reflecting their elemental makeup and energy transport mechanisms.

B-type main-sequence star Stellar classification distinguished by bright blue luminosity

A B-type main-sequence star is a main-sequence (hydrogen-burning) star of spectral type B and luminosity class V. These stars have from 2 to 16 times the mass of the Sun and surface temperatures between 10,000 and 30,000 K. B-type stars are extremely luminous and blue. Their spectra have neutral helium, which are most prominent at the B2 subclass, and moderate hydrogen lines. Examples include Regulus and Algol A.

Radiative transfer is the physical phenomenon of energy transfer in the form of electromagnetic radiation. The propagation of radiation through a medium is affected by absorption, emission, and scattering processes. The equation of radiative transfer describes these interactions mathematically. Equations of radiative transfer have application in a wide variety of subjects including optics, astrophysics, atmospheric science, and remote sensing. Analytic solutions to the radiative transfer equation (RTE) exist for simple cases but for more realistic media, with complex multiple scattering effects, numerical methods are required. The present article is largely focused on the condition of radiative equilibrium.

Hayashi track

The Hayashi track is a luminosity–temperature relationship obeyed by infant stars of less than 3 M in the pre-main-sequence phase of stellar evolution. It is named after Japanese astrophysicist Chushiro Hayashi. On the Hertzsprung–Russell diagram, which plots luminosity against temperature, the track is a nearly vertical curve. After a protostar ends its phase of rapid contraction and becomes a T Tauri star, it is extremely luminous. The star continues to contract, but much more slowly. While slowly contracting, the star follows the Hayashi track downwards, becoming several times less luminous but staying at roughly the same surface temperature, until either a radiative zone develops, at which point the star starts following the Henyey track, or nuclear fusion begins, marking its entry onto the main sequence.

The standard solar model (SSM) is a mathematical treatment of the Sun as a spherical ball of gas. This model, technically the spherically symmetric quasi-static model of a star, has stellar structure described by several differential equations derived from basic physical principles. The model is constrained by boundary conditions, namely the luminosity, radius, age and composition of the Sun, which are well determined. The age of the Sun cannot be measured directly; one way to estimate it is from the age of the oldest meteorites, and models of the evolution of the Solar System. The composition in the photosphere of the modern-day Sun, by mass, is 74.9% hydrogen and 23.8% helium. All heavier elements, called metals in astronomy, account for less than 2 percent of the mass. The SSM is used to test the validity of stellar evolution theory. In fact, the only way to determine the two free parameters of the stellar evolution model, the helium abundance and the mixing length parameter, are to adjust the SSM to "fit" the observed Sun.

Hypergiant Rare star with tremendous luminosity and high rates of mass loss by stellar winds

A hypergiant (luminosity class 0 or Ia+) is a very rare type of star that has an extremely high luminosity, mass, size and mass loss because of its extreme stellar winds. The term hypergiant is defined as luminosity class 0 (zero) in the MKK system. However, this is rarely seen in the literature or in published spectral classifications, except for specific well-defined groups such as the yellow hypergiants, RSG (red supergiants), or blue B(e) supergiants with emission spectra. More commonly, hypergiants are classed as Ia-0 or Ia+, but red supergiants are rarely assigned these spectral classifications. Astronomers are interested in these stars because they relate to understanding stellar evolution, especially with star formation, stability, and their expected demise as supernovae.

Red giant Type of large cool star that has exhausted its core hydrogen

A red giant is a luminous giant star of low or intermediate mass in a late phase of stellar evolution. The outer atmosphere is inflated and tenuous, making the radius large and the surface temperature around 5,000 K or lower. The appearance of the red giant is from yellow-orange to red, including the spectral types K and M, but also class S stars and most carbon stars.

R136a1 Wolf–Rayet star with one of the highest mass and luminosity of any known star

RMC 136a1 is one of the most massive and luminous stars known, at 215 M and nearly 6.2 million L, and is also one of the hottest, at around 46,000 K. It is a Wolf–Rayet star at the center of R136, the central concentration of stars of the large NGC 2070 open cluster in the Tarantula Nebula in the Large Magellanic Cloud. The cluster can be seen in the far southern celestial hemisphere with binoculars or a small telescope, at magnitude 7.25. R136a1 itself is 10,000 times fainter and can only be resolved using speckle interferometry.

AB7 Binary star in the Small Magellanic Cloud in the constellation Tucana

AB7, also known as SMC WR7, is a binary star in the Small Magellanic Cloud. A Wolf–Rayet star and a supergiant companion of spectral type O orbit in a period of 19.56 days. The system is surrounded by a ring-shaped nebula known as a bubble nebula.

In spectroscopy, collision-induced absorption and emission refers to spectral features generated by inelastic collisions of molecules in a gas. Such inelastic collisions may induce quantum transitions in the molecules, or the molecules may form transient supramolecular complexes with spectral features different from the underlying molecules. Collision-induced absorption and emission is particularly important in dense gases, such as hydrogen and helium clouds found in astronomical systems.

Stellar molecules are molecules that exist or form in stars. Such formations can take place when the temperature is low enough for molecules to form – typically around 6000 K or cooler. Otherwise the stellar matter is restricted to atoms in the forms of gas or – at very high temperatures – plasma.

Melnick 34 Binary wolf rayet star in the constellation Dorado

BAT99-116 is a binary Wolf–Rayet star near R136 in the 30 Doradus complex in the Large Magellanic Cloud. Both components are amongst the most massive and most luminous stars known, and the system is the most massive-known binary system.

AB8 (star) Star in the constrellation Hydrus

AB8, also known as SMC WR8, is a binary star in the Small Magellanic Cloud (SMC). A Wolf-Rayet star and a main sequence companion of spectral type O orbit in a period of 16.638 days. It is one of only nine known WO stars, the only Wolf-Rayet star in the SMC not on the nitrogen sequence, and the only Wolf-Rayet star in the SMC outside the main bar.

Schwarzschild's equation is used to calculate radiative transfer – energy transfer – through a medium in local thermodynamic equilibrium that both absorbs and emits electromagnetic radiation.

References