Search

Wilhelm et al. have recently called attention to the unidentified broad emission features near 117 nm in the solar spectrum. They discuss the observed properties of these features in detail but do not identify the source of this emission. We show that the broad autoionizing transitions of neutral sulfur are responsible for these emission features. Autoionizing lines of \ion{S}{i} occur throughout the spectrum between Lyman alpha and the Lyman limit. Sulfur is a normal contributor to stellar spectra. We use non-LTE chromospheric model calculations with line data from the Kurucz 2004 \ion{S}{i} line list to simulate the solar spectrum in the range 116 to 118 nm. We compare the results with SUMER disk-center observations from Curdt et al. and limb observations from Wilhelm et al. Our calculations generally agree with the SUMER observations of the broad autoionizing \ion{S}{i} emission features, the narrow \ion{S}{i} emission lines, and the continuum in this wavelength region, and agree with basic characteristics of the center-to-limb observations. In addition to modeling the average spectrum, we show that a change of $\pm$ 200 K in the temperature distribution causes the intensity to change by a factor of 4. This exceeds the observed intensity variations 1) with time in quiet regions at these wavelengths, and 2) with position from cell centers to bright network. These results do not seem compatible with current dynamical models that have temporal variations of 1000 K or more in the low chromosphere., Comment: 5 pages, 4 figures, 1 table (Accepted for Astronomy & Astrophysics)

This grant supported the research and publication of a major 26-page paper in The Astrophysical Journal, by Fontenla, Avrett, & Loeser (2002): 'Energy Balance in the Solar Transition Region. IV. Hydrogen and Helium Mass Flows with Diffusion.' This paper extended our previous modeling of the chromosphere-corona transition region to include cases with particle and mass flows. Inflows and outflows were shown to produce striking changes in the profiles of hydrogen and helium lines. An important conclusion is that line shifts are much less significant than the changes in line intensity and central reversal due to the influence of flows on the excitation and ionization of atoms in the solar atmosphere. This modeling effort at SAO is the only current one being undertaken anywhere to simulate in detail the full range of non-LTE absorption, emission, and scattering processes in the solar atmosphere to account for the entire solar spectrum from radio waves to X-rays. This effort is being continued with internal SAO funding at a relatively slow pace. Further NASA support in the future would yield results of great value for the interpretation of solar observations from NASA spacecraft.

During the last three years we have continued the development of extensive computer programs for constructing realistic models of the solar atmosphere and for calculating detailed spectra to use in the interpretation of solar observations. This research involves two major interrelated efforts: work by Avrett and Loeser on the Pandora computer program for optically thick non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Kurucz on the detailed high-resolution synthesis of the solar spectrum using data for over 58 million atomic and molecular lines. Our objective is to construct atmospheric models from which the calculated spectra agree as well as possible with high-and low-resolution observations over a wide wavelength range. Such modeling leads to an improved understanding of the physical processes responsible for the structure and behavior of the atmosphere.

This grant funded an observational and theoretical program to study the structure and dynamics of the solar photosphere and low chromosphere, and the spectral signatures that result. The overall goal is to learn about mechanisms that cause heating of the overlying atmosphere, and produce variability of solar emission in spectral regions important for astrophysics and space physics. The program exploited two new ground-based observational capabilities: one using the Swedish Solar Telescope on La Palma for very high angular resolution observations of the photospheric intensity field (granulation) and proxies of the magnetic field (G-band images); and the other using the Near Infrared Magnetograph at the McMath-Pierce Solar Facility to map the spatial variation and dynamic behavior of the solar temperature minimum region using infrared CO lines. We have interpreted these data using a variety of theoretical and modelling approaches, some developed especially for this project. Previous annual reports cover the work done up to 31 May 1997. This final report summarizes our work for the entire period, including the period of no-cost extension from 1 June 1997 through September 30 1997. In Section 2 we discuss observations and modelling of the photospheric flowfields and their consequences for heating of the overlying atmosphere, and in Section 3 we discuss imaging spectroscopy of the CO lines at 4.67 mu.

Since the early 1970s we have been developing the extensive computer programs needed to construct models of the solar atmosphere and to calculate detailed spectra for use in the interpretation of solar observations. This research involves two major related efforts: work by Avrett and Loeser on the Pandora computer program for non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Rurucz on the detailed synthesis of the solar spectrum based on opacity data or over 58 million atomic and molecular lines. our goals are: to determine models of the various features observed on the Sun (sunspots, different components of quiet and active regions, and flares) by means of physically realistic models, and to calculate detailed spectra at all wavelengths that match observations of those features. These two goals are interrelated: discrepancies between calculated and observed spectra are used to determine improvements in the structure of the models, and in the detailed physical processes used in both the model calculations and the spectrum calculations. The atmospheric models obtained in this way provide not only the depth variation of various atmospheric parameters, but also a description of the internal physical processes that are responsible for non-radiative heating, and for solar activity in general.

One-dimensional hydrostatic models of quiet and active solar regions can be constructed that generally account for the observed intensities of lines and continue throughout the spectrum, except for the infrared CO lines. There is an apparent conflict between: (1) observations of the strongest infrared CO lines formed in LTE at low-chromospheric heights but at temperatures much cooler than the average chromospheric values; and (2) observations of Ca II, UV (ultraviolet), and microwave intensities that originate from the same chromospheric heights but at the much higher temperatures characteristic of the average chromosphere. A model M(sub CO) has been constructed which gives a good fit to the full range of mean CO line profiles (averaged over the central area of the solar disk and over time) but this model conflicts with other observations of average quiet regions. A model L(sub CO) which is approximately 100 K cooler than M(sub CO) combined with a very bright network model F in the proportions 0.6 L(sub CO) + 0.4 F is found to be generally consistent with the CO, Ca II, UV, and microwave observations. Ayres, Testerman, and Brault found that models COOLC and FLUXT in the proportions 0.925 and 0.075 account for the CO and Ca II lines, but these combined models give an average UV intensity at 140 nm about 20 times larger than observed. The 0.6 L(sub CO) + 0.4 F result may give a better description of the cool and hot components that produce the space- and time-averaged spectra. Recent observations carried out by Uitenbroek, Noyes, and Rabine with high spatial and temporal resolution indicate that the faintest intensities in the strong CO lines measured at given locations usually become much brighter within 1 to 3 minutes. The cool regions thus seem to be mostly the low-temperature portions of oscillatory waves rather than cool structures that are stationary.

Since the early 1970s we have been developing the extensive computer programs needed to construct models of the solar atmosphere and to calculate detailed spectra for use in the interpretation of solar observations. This research involves two major related efforts: work by Avrett and Loeser on the Pandora computer program for non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Kurucz on the detailed synthesis of the solar spectrum based on opacity data for over 58 million atomic and molecular lines. Our goals are to determine models of the various features observed on the sun (sunspots, different components of quiet and active regions, and flares) by means of physically realistic models, and to calculate detailed spectra at all wavelengths that match observations of those features. These two goals are interrelated: discrepancies between calculated and observed spectra are used to determine improvements in the structure of the models, and in the detailed physical processes used in both the model calculations and the spectrum calculations. The atmospheric models obtained in this way provide not only the depth variation of various atmospheric parameters, but also a description of the internal physical processes that are responsible for nonradiative heating, and for solar activity in general.

The resonance lines of singly ionized magnesium, the MgII h&k lines at about 280 nm, are two of the small number of lines in the solar spectrum that are optically thick in the chromospheric part of the solar atmosphere. Potentially these lines contain information on the initial temperature rise that occurs at the top of the photosphere. Unfortunately, few good observations of the lines exist due to their wavelength near 280 nm the ultraviolet. However, a fair number of observations (on the order of 200) are available from the data base of the UltraViolet Polarimeter and Spectrometer (UVSP) instrument that flew on board of NASA's Solar Maximum Mission (SMM) satellite. In addition, this data base contains a number of spectra that include the Mg I resonance line at (lambda)285.2nm, just longward of the h&k lines. The neutral magnesium line is not as strong as its ionic counterparts and samples slightly lower parts of the atmosphere. Its width is a sensitive diagnostic of the ionization balance between neutral and singly ionized magnesium, which determines the opacity scale (and formation height) of other diagnostically important MgI lines like the 457.1 nm intercombination line, the magnesium b lines and the infrared MgI emission lines near 12 microns. Analysis of the observed line profiles shows that it is necessary to include the effects of partial frequency redistribution (PRD) in the formation of the line as in the case of the h&k lines. This implies that the core of the line is very sensitive to the way scattering is treated in the modeling of the line, and in turn this allows us to separate the uncertain effects in the atomic data (viz. the Van der Waals broadening) from the uncertainties in the underlying atmospheric model. The main objective of this research was to compare observed spectra of the magnesium resonance lines against theoretical line profiles calculated from recent models of the solar atmosphere by Fontenla et al., hereafter called FAL. These models extend earlier solar models by Avrett et al. and provide yet a better match between observed and calculated spectra, in particular of the hydrogen Lyman line, by including the effects of ambipolar diffusion of hydrogen and helium in the mass and energy transport equations for the chromosphere-corona transition region. Although the transition region models can be built self-consistently (without additional observational input) once the proper boundary conditions are given at the top of the chromosphere, there are uncertainties in the models precisely because these boundary conditions are not well known. With the analysis of UVSP magnesium observations we strive to constrain these uncertainties.

Since the early 1970s we have been developing the extensive computer programs needed to construct models of the solar atmosphere and to calculate detailed spectra for use in the interpretation of solar observations. This research involves two major related efforts: work by Avrett and Loeser on the Pandora computer program for non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Kurucz on the detailed synthesis of the solar spectrum based on opacity data for over 58 million atomic and molecular lines. Our goals are to determine models of the various features observed on the Sun (sunspots, different components of quiet and active regions, and flares) by means of physically realistic models, and to calculate detailed spectra at all wavelengths that match observations of those features. These two goals are interrelated: discrepancies between calculated and observed spectra are used to determine improvements in the structure of the models, and in the detailed physical processes used in both the model calculations and the spectrum calculations. The atmospheric models obtained in this way provide not only the depth variation of various atmospheric parameters, but also a description of the internal physical processes that are responsible for non-radiative heating, and for solar activity in general.

Our energy-balance models of the lower transit ion region were presented in a previous paper (New Models of the Chromosphere and Transition Region). Here we show the influence of particle diffusion on the calculated hydrogen and helium number densities for a given temperature-density model (model C in the preceding paper). We have also solved the statistical equilibrium and radiative transfer equations for a 13-level He I atom (22 radiative transitions) and a 6-level He II ion (15 radiative transitions) together with He III. The resulting He I and He II level-1 number densities and He III density are shown as functions of temperature. Diffusion substantially increases eta(sub HeI) for T greater than 35,000 K and decreases eta(sub HeI) (while increasing eta(sub(HeII)) between 9,000 and 25,000 K. Including the effects of diffusion also increases eta(sub HeIII) for T less than 60,000 K.

The present characterization of the solar photosphere as a radiative boundary for the solar interior gives attention to both semiempirical and a priori models for the structure of emergent radiation, and emphasizes that any discussion of photospheric structures smaller than 300 km is easily compromised by the radiant integration over such scale lengths. Plane-parallel modes have indicated that Boltzmann-Saha statistics are largely valid over the photosphere save within 100 km of the temperature minimum. An evaluation is made of observational evidence for horizontal inhomogeneities and such characteristics as granular and temperature minimum region fluctuations, as well as the recent demonstration of relatively cold CO-absorbing regions.

Models are presented for the two continuum-emitting kernels observed in the white-light flare (WLF) of June 15, 1982. They are the first semiempirical models of a WLF which are consistent not only with observations of the continuum emission level but also with a set of spectral lines having heights of formation which span the chromosphere and upper photosphere. It is shown that the models are not compatible with the hypothesis that the continuum emission is caused by enhanced Balmer and Paschen hydrogen continua, and they present strong evidence instead that the emission is of photospheric origin and that its source is due to H(-). The observation and the models derived from them show that white-light emission can occur in areas of the active region where there is no chromospheric emission and in particular no H-alpha emission. This fact seems to rule out the viability of downward transport mechanisms as the source of the energy required for the WLF.

The CO cooling rate for models of the solar atmosphere using the detailed line-by-line CO opacity in the fundamental band, and carrying out a full radiative transfer calculation for each line is computed. The importance of the different assumptions that have been made to obtain the CO cooling rate and find that when detailed optical depth effects are taken into account, the calculated CO cooling rate at line optical depths near unity can be smaller than optically thin estimates by more than an order of magnitude is studied. It is found that CO cooling does not account for the missing source of radiative cooling in the temperature minimum region of the quiet sun.

Various methods of inverting solar Mg I 4571 and 5173 spectral line observations are examined to find the best method of using these lines to calculate the vertical temperature and electron density structure around the temperature minimum region. Following a perturbation analysis by Mein (1971), a Fredholm integral equation of the first kind is obtained which can be inverted to yield these temperature and density structures as a function of time. Several inversion methods are tested and compared. The methods are used to test data as well as to a subset of observations of these absorption lines taken on February 3, 1986 before and during a solar flare. A small but significant increase is found in the temperature and a relatively large increase in the electron density during this flare. The observations are inconsistent with heating and ionization by an intense beam of electrons and with ionization by UV photoionization of Si I.

Line fluxes and profiles are computed for a sequence of spherically symmetric T Tauri wind models. The calculations indicate that the H-alpha emission of T Tauri stars arises in an extended and probably turbulent circumstellar envelope at temperatures above about 8000 K. The models predict that Mg II resonance line emission should be strongly correlated with H-alpha fluxes; observed Mg II/H-alpha ratios are inconsistent with the models unless extinction corrections have been underestimated. The models predict that most of the Ca II resonance line and IR triplet emission arises in dense layers close to the star rather than in the wind. H-alpha emission levels suggest mass loss rates of about 10 to the -8th solar mass/yr for most T Tauri stars, in reasonable agreement with independent analysis of forbidden emission lines. These results should be useful for interpreting observed line profiles in terms of wind densities, temperatures, and velocity fields.

A semiempirical chromospheric model is proposed for TX Psc which is a prototype for the N-type carbon stars. Observational data imply that the chromospheric temperature rise must begin at a low density, that the temperature gradient in the lower chromosphere must be steep, that partial redistribution must be employed in the Mg II calculation, and that the lower chromosphere is expanding away from the photosphere with a velocity of close to 50 km/s. The present model also shows that the microturbulent velocity is about 7 km/s at the temperature minimum region, dropping to 5 km/s in the chromosphere, and that the Lyman lines are optically thick in the chromosphere.

A nine-level atomic model is presented for C I lines synthesis in the sun. All the atomic parameters are based on recent experimental and theoretical data. Profiles of the multiplets at 1561 and 1657 A are computed from a model of the solar chromosphere. The sensitivity of the results to changes in the atomic parameters is examined. Partial frequency redistribution must be included in these lines to obtain agreement between calculated and observed profiles.

This paper presents a complete atomic model for Mg I line synthesis, where all the atomic parameters are based on recent experimental and theoretical data. It is shown how the computed profiles at 4571 A and 5173 A are influenced by the choice of these parameters and the number of levels included in the model atom. In addition, observed profiles of the 5173 A b2 line and theoretical profiles for comparison (based on a recent atmospheric model for the average quiet sun) are presented.