Hydrogen line blanketed model stellar atmospheres
Hydrogen line blanketed stellar model atmospheres
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Hydrogen line blanketed stellar model atmospheres
Several grids of theoretical H II region models are computed by photoionization modeling in order to determine the extent to which the choice of the ionizing stellar atmosphere model affects the calibration of emission-line diagnostic diagrams of Evans and Dopita (1985) and the semiempirical H II region abundance sequence calibration of Evans and Dopita. Emission-line diagnostic diagrams are presented and compared for model nebulae ionized by Hummer and Mihalas (1970) unblanketed LTE atmospheres, Kurucz (1979) line-blanketed LTE atmospheres, Mihalas (1972) unblanketed non-LTE, and a truncated blackbody spectrum. The models demonstrate that for solar nebular and atmospheric abundances, there are only minor differences between H II models ionized by the Hummer and Mihalas atmospheres. The unblanketed non-LTE stellar atmosphere models of Mihalas and truncated blackbody spectra are shown to be unsuitable for general H II region modelling.
The luminosity, mass, and elemental abundances, as well as other properties of each star are studied in order to locate them in an evolutionary pattern. A method for determining the flux, gravity, and abundances at the stellar surface is the construction of theoretical stellar atmospheric models that predict the observed energy distribution and detailed stellar spectrum.
I discuss errors in theory and in interpreting observations that are produced by the failure to consider resolution in space, time, and energy. I discuss convection in stellar model atmospheres and in stars. Large errors in abundances are possible such as the factor of ten error in the Li abundance for extreme Population II stars. Finally I discuss the variation of microturbulent velocity with depth, effective temperature, gravity, and abundance. These variations must be dealt with in computing models and grids and in any type of photometric calibration. I have also developed a new opacity-sampling version of my model atmosphere program called ATLAS12. It recognizes more than 1000 atomic and molecular species, each in up to 10 isotopic forms. It can treat all ions of the elements up through Zn and the first 5 ions of heavier elements up through Es. The elemental and isotopic abundances are treated as variables with depth. The fluxes predicted by ATLAS12 are not accurate in intermediate or narrow bandpass intervals because the sample size is too small. A special stripped version of the spectrum synthesis program SYNTHE is used to generate the surface flux for the converged model using the line data on CD-ROMs 1 and 15. ATLAS12 can be used to produce improved models for Am and Ap stars. It should be very useful for investigating diffusion effects in atmospheres. It can be used to model exciting stars for H II regions with abundances consistent with those of the H II region. These programs and line files will be distributed on CD-ROMs.
Two approaches to accelerating the method of complete linearization for calculating NLTE model stellar atmospheres are suggested. The first one, the so-called Kantorovich variant of the Newton-Raphson method, consists of keeping the Jacobi matrix of the system fixed, which allows us to calculate the costly matrix inversions only a few times and then keep them fixed during the subsequent computations. The second method is an application of the Ng acceleration. Both methods are extremely easy to implement with any model atmosphere code based on complete linearization. It is demonstrated that both methods, and especially their combination, yield a rapidly and globally convergent algorithm, which takes 2 to 5 times less computer time, depending on the model at hand and the required accuracy, than the ordinary complete linearization. Generally, the time gain is more significant for more complicated models. The methods were tested for a broad range of atmospheric parameters, and in all cases they exhibited similar behavior. Ng acceleration applied on the Kantorovich variant thus offers a significant improvement of the standard complete-linearization method, and may now be used for calculating relatively involved NLTE model stellar atmospheres.
Monte Carlo relaxation method for physically self consistent model stellar atmospheres, discussing stability, accuracy and convergence of solution
The primary goal of this analysis is to determine whether the effects of atomic bound-bound transitions on stellar atmospheric structure can be represented well in models. The investigation is based on an approach which is called the method of artificial absorption edges. The method is described, developed, tested, and applied to the problem of fitting a model stellar atmosphere to Sirius. It is shown that the main features of the entire observed spectrum of Sirius can be reproduced to within the observational uncertainty by a blanketed flux-constant model with T sub eff = 9700 K and Log g = 4.26. The profile of H sub gamma is reproduced completely within the standard deviations of the measurements except near line center, where non-LTE effects are expected to be significant. The equivalent width of H sub gamma, the Paschen slope, the Balmer jump, and the absolute flux at 5550 A all agree with the observed values.
The method of artificial absorption edges is applied to a grid of model atmospheres to include the effects of hydrogen and metal line blanketing on the atmospheric structure. A best-fit model of Sirius is obtained which is free of the discrepancies inherent in unblanketed or partially blanketed models. All major spectral features calculated from the model agree with the corresponding observed features to within the observational uncertainty. Some additional properties of the models are discussed.
Line profiles of 36 multiplets of silicon (Si) II, III, and IV were computed for a grid of model atmospheres covering the range from 15,000 to 35,000 K in effective temperature and 2.5 to 4.5 in log (gravity). The computations involved simultaneous solution of the steady-state statistical equilibrium equations for the populations and of the equation of radiative transfer in the lines. The variables were linearized, and successive corrections were computed until a minimal accuracy of 1/1000 in the line intensities was reached. The common assumption of local thermodynamic equilibrium (LTE) was dropped. The model atmospheres used also were computed by non-LTE methods. Some effects that were incorporated into the calculations were the depression of the continuum by free electrons, hydrogen and ionized helium line blocking, and auto-ionization and dielectronic recombination, which later were found to be insignificant. Use of radiation damping and detailed electron (quadratic Stark) damping constants had small but significant effects on the strong resonance lines of Si III and IV. For weak and intermediate-strength lines, large differences with respect to LTE computations, the results of which are also presented, were found in line shapes and strengths. For the strong lines the differences are generally small, except for the models at the hot, low-gravity extreme of our range. These computations should be useful in the interpretation of the spectra of stars in the spectral range B0-B5, luminosity classes III, IV, and V.
In this paper we study the influence of line-merging regions at the intermediate long-wavelength side of a continuum threshold on the computed model atmosphere structure and predicted spectrum. In order to model these regions sufficiently accurately, we have developed two concepts. First, we have extended the occupation probability formalism of Hummer and Mihalas to non-local thermodynamic equilibrium (non-LTE) plasmas. Second, in order to treat the very complicated opacity in the line merging region, we have generalized the concept of opacity distribution functions to treat non-LTE situations. All Rydberg states are consistently included within this framework, so that no arbitrary cutoff of high (LTE) levels is made. We have calculated several pure hydrogen models atmospheres for two effective temperatures, T(sub eff) = 20000 and 35000 K, and discussed the differences between models calculated with various treatments of the line merging. In particular, we have shown that the error in the predicted profiles of Balmer lines resulting from the neglect of line merging is typically of the order of 3-4%, while the errors in the far-UV portion of the Balmer continuum reaches 15-35%. The errors generally decrease with increasing effective temperature. At the same time, the internal accuracy of the models is shown to be about or below 0.5% for all predicted spectral features. We conclude that for interpreting current high-accuracy spectrophotometric observations models including the line merging are necessary, and that the formalism developed in this paper is capable of providing a sufficiently accurate and robust modeling technique.
Stellar atmosphere models that assume interior energy sources, thermodynamic and mechanical equilibrium, and absorption of radiant energy only in continuous spectrum
Convection role on integrated absorption lines in stellar atmosphere models
Active galactic nuclei (AGN) accretion disk spectra were calculated using non-LTE stellar atmosphere models for Kerr and Schwarzschild geometries. It is found that the Lyman limit absorption edge, probably the most conclusive observational evidence for the accretion disk, would be drastically distorted and displaced by the relativistic effects from the large gravitational field of the central black hole and strong Doppler motion of emitting material on the disk surface. These effects are especially pronounced in the Kerr geometry. The strength of the Lyman limit absorption is very sensitive to the surface gravity in the stellar atmosphere models used. For models at the same temperature but different surface gravities, the strength of the Lyman edge exhibits an almost exponential decrease as the surface gravity approach the Eddington limit, which should approximate the thin disk atmosphere. The relativistic effects as well as the vanishing of the Lyman edge at the Eddington gravity may be the reasons that not many Lyman edges in the rest frames of AGNs and quasars are found.
Using the results of the most recent stellar atmosphere models applied to a sample of hot stars, we construct calibrations of effective temperature (T(sub eff)), and gravity (log(sub g)) with a spectral type and luminosity class for Galactic 0-type and early B-type stars. From the model results we also derive an empirical relation between the bolometric correction and T(sub eff) and log g. Using a sample of stars with known distances located in OB associations in the Galaxy and the Large Magellanic Cloud, we derive a new calibration of M(sub v) with spectral class. With these new calibrations and the stellar atmosphere models of Kurucz, we calculate the physical parameters and ionizing photon luminosities in the H(0) and He(0) continua for O and early B-type stars. We find substantial differences between our values of the Lyman- continuum luminosity and those reported in the literature. We also discuss the systematic discrepancy between O-type stellar masses derived from spectroscopic models and those derived from evolutionary tracks. Most likely, the cause of this 'mass discrepancy' lies primarily in the atmospheric models, which are plane parallel and hydrostatic and therefore do not account for an extended atmosphere and the velocity fields in a stellar wind. Finally, we present a new computation of the Lyman-continuum luminosity from 429 known O stars located within 2.5 kpc of the Sun. We find the total ionizing luminosity from this population ((Q(sub 0)(sup T(sub ot))) = 7.0 x 10(exp 51) photons/s) to be 47% larger than that determined using the Lyman continuum values tabulated by Panagia.
An intercomparison of model stellar atmospheres and observations of real stars can lead to a better understanding of the relationship between the physical properties of stars and their observed radiative flux. In this spirit we have determined wide-band and narrow-band magnitudes and colors for a subset of models of K and M giant and supergiant stars selected from the grid of 40 models by Johnson, Bernat and Krupp (1980) (hereafter referred to as JBK). The 24 models selected have effective temperatures of 4000, 3800, 3600, 3400, 3200, 3000, 2750 and 2500 K and log g = 0, 1 or 2. Emergent energy fluxes (erg/ sq cm s A) were calculated at 9140 wavelengths for each model. These computed flux curves were folded through the transmission functions of Wing's 8-color system (Wing, 1971; White and Wing, 1978) and through Johnson's (1965) wide-band (BVRIJKLM) system. The calibration of the resultant magnitudes was made by using the absolute calibration of the flux curve of Vega by Schild, et al. (1971).
The main data need for stellar atmosphere and spectrum modeling remains atomic and molecular transition data, particularly energy levels and transition cross-sections. We emphasize that data is needed for bound-free (b - f) as well as bound-bound (b - b), and collisional as well as radiative transitions. Data is now needed for polyatomic molecules as well as atoms, ions, and diatomic molecules. In addition, data for the formation of, and extinction due to, liquid and solid phase dust grains is needed. A prioritization of species and data types is presented, and gives emphasis to Fe group elements, and elements important for the investigation of nucleosynthesis and Galactic chemical evolution, such as the -elements and n-capture elements. Special data needs for topical problems in the modeling of cool stars and brown dwarfs are described.
Ultraviolet (UV) fluxes and other data for a sample of nine novalike cataclysmic variables are assembled from the literature. The UV fluxes and colors are compared to spectra of steady state model accretion disks, constructed using either Planck functions or model stellar atmosphere spectra as elementary radiators. Deficiencies are found with both sets of models. The restriction to steady state radial temperature profiles is relaxed for the 'stellar atmosphere' disks, but these generalized models still fail to account simultaneously for the absolute UV flux and the UV color of the observed spectra. The conclusion is that stellar atmosphere models do not reflect the physics of disks. More appropriate models probably can and definitely should be made.