The transfer of polarized radiation in spectral lines: Formalism and solutions in simple cases

A general treatment of the transfer of polarized radiation in spectral lines assuming a Rayleigh phase function and a general law of frequency redistribution is derived. It is shown how nine families of coupled integral equations for the moments of the radiation field arise which are necessary to fully describe the state of polarization of the emergent radiation from a plane-parallel, semi-infinite atmosphere. The special case of angle independent redistribution functions is derived from the general formalism, and it is shown how the nine families of integral equations reduce to the six linearly independent integral equations derived by Collins (1972). To serve as a test of the formulation, solutions for isothermal atmospheres are given.     

On nonstationary radiative transfer in a spectral line in stellar atmospheres

We consider nonstationary radiative transfer in a line in stellar atmospheres simulated as a stationary semi-infinite plane-parallel medium. We assume complete frequency redistribution in the elementary act of scattering. We assume that the time a photon spends in the medium is determined only by the mean time spent in the absorbed state. We obtain an explicit expression for the resolvent of the nonstationary integral equation of transfer, which is a bilinear expansion with respect to the eigenfunctions found in [12] for the corresponding stationary transfer equation.Translated fromAstrofizika, Vol. 37, No. 3, 1994.I am grateful to the American Astronomical Society for financial support of this work.     

Formation of spectral lines in stellar atmospheres. Linear and nonlinear theory

The present paper represents a short review of some directions of radiative transfer theory in spectral lines (RTSL) concerning investigations of authors in this field.Published in Astrofizika, Vol. 38, No. 4, pp. 565–571, October–December, 1995.     

Polarization of radiation from stellar atmospheres. The grey case

An accurate numerical method is presented for the solution of the transfer equations in a plane-parallel atmosphere in which scattering occurs according to Rayleigh's law. Some results are given for the polarization and limb darkening of both integrated and monochromatic radiation emerging from grey atmospheres with various ratios of scattering to absorption. The method is equally applicable to non-grey atmospheres.     

On the influence of density and temperature fluctuations on the formation of spectral lines in stellar atmospheres

A method taking into account the influence of temperature and density flucuations generated by the velocity field in stellar atmospheres on the formation of spectral lines is presented. The influenced line profile is derived by exchanging the values in a static atmosphere by a mean value and a fluctuating one. The correlations are calculated with the help of the well-know hydrodynamic eqs. It results, that in normal stellar atmospheres the visual lines are only very weakly influenced by such fluctuations due to the small values of the gradients of the pressure and density and of the velocity dispersion.     

Macro-and micro-turbulent filter functions for weak lines in stellar atmospheres

Following a similar discussion given earlier for the solar case (De Jager, 1972) we compute in this paper spectral line profiles for the spatial wavelengths in which a stellar motion field can be decomposed, and thereafter the macro-and micro-turbulent filter functionsf _M(k) andf (k), where is the optical scale height andf ²(k) dk the fraction of the energy of the turbulent motions between wavenumbersk andk+dk of the spectrum of turbulence that contributes to either kind of turbulence. If micro-and macro-turbulent velocity components are known for a certain star, and if the spectrum of turbulence is sharp enough, the ratiof _M/f would enable one to derive the average size of the turbulent elements in the star's atmosphere. The computations apply to weak lines in idealized stellar atmospheres, and refer to two cases: isotropic turbulence, and radial pulsations. These filters can be suitably used in a diagnostic method for the analysis of the motion field in the solar and stellar atmospheres. Some examples of applications to stars of very different kinds are given.     

Diagnostics for regions of formation of atomic lines in stellar atmospheres

We present a procedure to make direct diagnostics of the physical conditions of the regions of formation of the atomic spectral lines in stellar atmospheres using the atomic line widths at half maximum and the number of lines visible of a given atomic series in the observed stellar spectra. This is accomplished using the theoretical widths at half maximum of the atomic lines induced by the broadening produced by thermal energy fluctuations and considering the maximum number of levels that exist in those atoms under the physical conditions of the given system. The procedure is easy to use in any application. As an example we apply the procedure to the observed Lyman lines of hydrogen in the ultraviolet of some stars.     

Non-LTE effects for Na I lines in X-ray illuminated stellar atmospheres

The formation of Na I lines in X-ray illuminated atmospheres is investigated by abandoning the assumption of local thermodynamic equilibrium (LTE). Calculations are performed on the basis of a 21-level Na I model atom for the LTE model atmospheres of irradiated F-G stars obtained with allowance for a reflection effect in the first approximation. The state of extreme “overrecombination” is shown to exist for the populations of all Na I levels in the case of external illumination. Absorption features in the profiles of “cool” and “normal” Na I lines have been found to be enhanced compared to the LTE approximation. Effects of the angle of incidence and intensity of the external radiation on the formation of level populations and line profiles when abandoning LTE are analyzed. The existence of overrecombination for Na I is explained by the small X-ray heating function and the large optical cooling function. Na I level populations are shown to depend weakly on the presence of “overionization” for Na II in the atmospheres of irradiated stars.     

Helium in stellar atmospheres

Helium, which was first discovered on the sun with the help of spectral analysis, plays, together with hydrogen, a principal role in astrophysics. We consider here two fundamental quantities: primordial helium abundance formed during Big Bang nucleosynthesis and the current initial helium abundances in nearby stars. It is shown that stellar atmospheres are enriched in helium during the main-sequence stage. Observational evidence for helium contamination in close OB-binaries is discussed. Stars with strong abundance anomalies are considered, such as chemically peculiar Ap and Bp helium-deficient stars and some types of objects with helium atmospheres.     

Order-of-scattering of partially polarized radiation in inhomogeneous,anisotropically scattering atmospheres

A complete set of transfer equations required for the order-of-scattering analysis of partially polarized radiation in inhomogeneous, anisotropically scattering atmospheres is provided. The equations have been derived for both a local study using the radiative transfer equation and its associated auxiliary equation for the source-matrix, and a global study in terms of the scattering and transmission matrices; they account for the polarity of the scattering medium. Their derivations for the finite order scattering and the finitely cumulative scattering, in particular, have yielded important new equations expressing the invariance principles and the integro-differential recurrences for the scattering and transmission matrices. These novel expressions contain as a special case Bellmanet al's (1972) equations for the simpler case of isotropic scattering of unpolarized light in homogeneous atmospheres.     

Order-of-scattering of partially polarized radiation in inhomogeneous,anisotropically scattering atmospheres

Reciprocity and symmetry relationships, representing local invariants for the scattering phase-matrix, are derived for twelve cases of particle assemblies studied by van de Hulst (1957) including situations of scattering in an arbitrary direction, in the near forward and near backward directions. These relations are used to generate corresponding relations representing global invariants for the scattering and transmission matrices of atmospheres consisting of such assemblies. The latter relations are obtained from the matrix integro-differential equations for scattering and transmission; they apply to single scattering, any finite order of scattering, and after an arbitrary cumulation of scattering orders (finite or infinite). Our results are summarized in Tables I and II for general inhomogeneous atmospheres and for particular inhomogeneous atmospheres that are symmetrical with respect to their central level. The latter case includes homogeneous atmospheres as a special case. The largest set of local relations obtained contains three independent relations (called universal, reversal, exchange) which can further be combined to yield four additional dependent relations. This circumstance happens in three out of the above twelve cases. In the remaining cases fewer relations (both independent and dependent) remain valid. Likewise, a maximal set of three independent global relations is obtained for general inhomogeneous atmospheres; they too can be linearly combined to yield seven other dependent relations. For the symmetrically inhomogeneous atmospheres, three independent and seven dependent additional relations are obtained. On the basis of these tables, it becomes a trivial matter to provide the local and global invariants (both the independent and the dependent relations) for any assembly of particles and atmospheric inhomogeneity. A mixture of Rayleigh-Cabannes scattering by anisotropic molecules or extremely small particles and Mie scattering by large isotropic particles is considered for illustration. Lastly, the group properties of these invariants are studied.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract No. NAS 7-100, sponsored by the U.S. National Aeronautics and Space Administration.     

The Rosseland mean opacity of stellar atmospheres

Rosseland mean opacities $\text{a}\text{?}{}_{\mathrm{R}}\text{(T,}\text{lgP}{}_{\mathrm{e}}\text{)}$ [see Table 2] per unit hydrogen particle are calculated in this paper for the atmosphere of the solar type [see Table 1].The sources of the continuous opacities under consideration in this paper are as follows: (1) H^?, HI, H^?₂, H⁺₂ He^?, HeI, HeII, CI, CII, CIII, NI, NII, NIII, OI, OII, NaI, MgI, MgII, AlI, AlII, SiI, SiII, CII, KI, CaII in the form of absorption. (2) HI, HeI, CI, NI, OI, H₂ in the form of Rayleigh scattering. (3) Free electrons in the form of Thomson scattering.     

Pinchrosphere — The coolest layer in stellar atmospheres

An attempt is made to demonstrate the possibility of the existence of one more layer—pinchrosphere—in the structure of atmospheres of late-type stars—later than G5. Pinchrosphere is the coolest layer with very low electron temperature, 5500K, it is located on the photosphere and near the inner boundary of the chromosphere. Practically all elements in pinchrosphere, particularly hydrogen and magnesium, are in neutral state. The main physical parameters of the pinchrosphere are derived.     

Formation of Zr I and II lines under non-LTE conditions of stellar atmospheres

The formation of Zr I and Zr II lines in stellar atmospheres under non-LTE conditions has been considered for the first time. A model zirconium atom has been composed using 148 Zr I levels, 772 Zr II levels, and the ground Zr III state. Non-LTE calculations have been performed for model atmospheres with T _eff = 5500 and 6000 K, log g = 2.0 and 4.0, [M/H] = −3, −2, −1, 0. In the entire investigated range of parameters, the Zr I levels are shown to be underpopulated relative to their LTE populations in the line formation region. In contrast, the excited Zr II levels are overpopulated, while the ground state and lower excited levels of Zr II retain their LTE populations. Since the non-LTE effects cause the Zr I and Zr II spectral lines being investigated to weaken, the non-LTE corrections to the abundance derived from Zr I and Zr II lines are positive. For Zr II lines, they increase with decreasing metallicity and surface gravity up to 0.34 dex for the model with T _eff = 5500, log g = 2.0, and [M/H] = −2. The non-LTE effects depend weakly on temperature. The non-LTE corrections for Zr I lines reach 0.33 dex for solar-metallicity models. Zr I and Zr II lines in the solar spectrum have been analyzed. The non-LTE zirconium abundances derived from lines in the two ionization stages are shown to agree between themselves within the error limits, while the LTE abundance difference is 0.28 dex. The zirconium abundance in the solar atmosphere (averaged over Zr I and Zr II lines) is log ɛ_Zr,⊙ = 2.63 ± 0.07.     

Convective energy transport in stellar atmospheres

The motion of convective cells in an environment which changes rapidly with depth is examined. In such an environment a cell may move through regions with different levels of ionization and with associated differences in heat capacity. The energy equation is cast in a manner which is independent of the history of these cells. The convective flux at a given level of the atmosphere is written as an average over an ensemble of cells originating at a range of other levels. A procedure for correcting the temperature gradient for these non-local effects is described and results for a model solar atmosphere are given. The principal results are: (1) The rms velocity varies smoothly and is non-zero well into the photosphere (e.g.,v _rsm=1.4 km/sec at =0.2); (2) Convective overshoot reduces the radiative flux to 60% and 90% of the total at =2.5 and 0.2 respectively; and (3) The interior adiabat of the convective envelope is less sensitive to the assumed value of the average cell size than in the usual treatment of convection.Supported in part by the National Science Foundation [GP-9433, GP-9114], the Office of Naval Research [Nonr-220(47)], and Air Force Grant AG-AFOSR-171-67.     

Convective energy transport in stellar atmospheres

We discuss a theoretical method of computing the temperature structure of hot and cool streams in convective stellar atmospheres. The method is based on the model that the streams are due to organized cells whose diameters are greater than the thickness of the photosphere. The excess thermal energy of matter rising from the deeper layers, where the entropy is higher than in the photosphere, is converted to radiation in a steady front. This model, applied to the solar case, exhibits a peak-to-peak contrast of 30–40% between granules and lanes. This contrast agrees with the Stratoscope data reduced by Namba and Diemel (1969). As a necessary part of the theory, we obtain an expression for the perturbation in radiative heat exchange which may be used in a medium with a strongly preferred direction such as a stellar atmosphere.Supported in part by the National Science Foundation [GP-15911 (formerly GP-9433), GP-9114] and the Office of Naval Research [Nonr-220(47)].     

Turbulence influenced line formation in stellar atmospheres

The influence of stochastic velocity fields with finite correlation lengths on the formation of spectral lines is taken into account without restrictions to specific velocity models. To construct a perturbation theory treating the influence of stochastic motion on the radiative transfer we start with the stochastic transfer equation for plane-parallel atmospheres and expand it to an infinite hierarchical system. An appropriate cut-off of the infinite system admits to achieve exactly the micro- and macroturbulent limit at every level of approach. The formalism is derived for n-point correlations of the absorption coefficient and specialised for 2-point correlations.