Search NASA⌕ Search

Engineering topics

Badnell, N. R.

Publications and source records attributed to Badnell, N. R..

At least 19 records

Quantitative comparison of opacities calculated using the R -matrix and distorted-wave methods: Fe xvii

ABSTRACT We present here a detailed calculation of opacities for Fe xvii at the physical conditions corresponding to the base of the Solar convection zone. Many ingredients are involved in the calculation of opacities. We review the impact of each ingredient on the final monochromatic and mean opacities (Rosseland and Planck). The necessary atomic data were calculated with the R-matrix and the distorted-wave (DW) methods. We study the effect of broadening, of resolution, of the extent of configuration sets and of configuration interaction to understand the differences between several theoretical predictions as well as the existing large disagreement with measurements. New Dirac R-matrix calculations including all configurations up to the n = 4, 5, and 6 complexes have been performed as well as corresponding Breit–Pauli DW calculations. The DW calculations have been extended to include autoionizing initial levels. A quantitative contrast is made between comparable DW and R-matrix models. We have reached self-convergence with n = 6 R-matrix and DW calculations. Populations in autoionizing initial levels contribute significantly to the opacities and should not be neglected. The R-matrix and DW results are consistent under the similar treatment of resonance broadening. The comparison with the experiment shows a persistent difference in the continuum while the filling of the windows shows some improvement. This study defines our path to the next generation of opacities and opacity tables for stellar modelling.

Delahaye, F.↗

Iso-nuclear tungsten dielectronic recombination rates for use in magnetically-confined fusion plasmas

Under the auspices of the IAEA Atomic and Molecular Data Center and the Korean Atomic Energy Research Institute, our assembled group of authors has reviewed the current state of dielectronic recombination (DR) rate coefficients for various ion stages of tungsten (W). Subsequent recommendations were based upon available experimental data, first-principle calculations carried out in support of this paper and from available recombination data within existing atomic databases. If a recommendation was possible, data were compiled, evaluated and fitted to a functional form with associated uncertainty information retained, where available. This paper also considers the variation of the W fractional abundance due to the underlying atomic data when employing different data sets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bootstrapping dielectronic recombination from second-row elements and the orion nebula

Dielectronic recombination (DR) is the dominant recombination process for most heavy elements in photoionized clouds. Accurate DR rates for a species can be predicted when the positions of autoionizing states are known. Unfortunately such data are not available for most third- and higher-row elements. This introduces an uncertainty that is especially acute for photoionized clouds, where the low temperatures mean that DR occurs energetically through very low-lying autoionizing states. This paper discusses S{sup 2+} → S{sup +} DR, the process that is largely responsible for establishing the [S iii]/[S ii] ratio in nebulae. We derive an empirical rate coefficient using a novel method for second-row ions, which do have accurate data. Photoionization models are used to reproduce the [O iii]/[O ii]/[O i]/[Ne iii] intensity ratios in central regions of the Orion Nebula. O and Ne have accurate atomic data and can be used to derive an empirical S{sup 2+} → S{sup +} DR rate coefficient at ∼10{sup 4} K. We present new calculations of the DR rate coefficient for S{sup 2+} → S{sup +} and quantify how uncertainties in the autoionizing level positions affect it. The empirical and theoretical results are combined and we derive a simple fit to the resulting rate coefficient at all temperatures for incorporation into spectral synthesis codes. This method can be used to derive empirical DR rates for other ions, provided that good observations of several stages of ionization of O and Ne are available.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron-ion recombination of Fe{sup 12+} forming Fe{sup 11+}: Laboratory measurements and theoretical calculations

We have measured dielectronic recombination (DR) for Fe{sup 12+} forming Fe{sup 11+} using the heavy ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. Using our results, we have calculated a plasma rate coefficient from these data that can be used for modeling astrophysical and laboratory plasmas. For the low temperatures characteristic of photoionized plasmas, the experimentally derived rate coefficient is orders of magnitude larger than the previously recommended atomic data. The existing atomic data were also about 40% smaller than our measurements at temperatures relevant for collisionally ionized plasmas. Recent state-of-the-art theory has difficulty reproducing the detailed energy dependence of the DR spectrum. However, for the Maxwellian plasma rate coefficient, recent theoretical results agree with our measurements to within about 30% for both photoionized and collisionally ionized plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron-impact excitation of Fe{sup 2+}: A comparison of intermediate coupling frame transformation, Breit-Pauli and Dirac R-matrix calculations

Modeling the spectral emission of low-charge iron group ions enables the diagnostic determination of the local physical conditions of many cool plasma environments such as those found in H II regions, planetary nebulae, active galactic nuclei, etc. Electron-impact excitation drives the population of the emitting levels and, hence, their emissivities. By carrying-out Breit-Pauli and intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of Fe{sup 2+}, which both use the exact same atomic structure and the same close-coupling expansion, we demonstrate the validity of the application of the powerful ICFT method to low-charge iron group ions. This is in contradiction to the finding of Bautista et al., who carried-out ICFT and Dirac R-matrix calculations for the same ion. We discuss possible reasons.

79 ASTRONOMY AND ASTROPHYSICS↗

OPACITIES AND SPECTRA OF THE r-PROCESS EJECTA FROM NEUTRON STAR MERGERS

Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements through the r-process. The subsequent radioactive decay of the nuclei can power transient electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying such events is an important goal of future optical surveys, offering new perspectives on the origin of r-process nuclei and the astrophysical sources of gravitational waves. Predictions of the transient light curves and spectra, however, have suffered from the uncertain optical properties of heavy ions. Here we argue that the opacity of an expanding r-process material is dominated by bound-bound transitions from those ions with the most complex valence electron structure, namely the lanthanides. For a few representative ions, we run atomic structure models to calculate the radiative transition rates for tens of millions of lines. The resulting r-process opacities are orders of magnitude larger than that of ordinary (e.g., iron-rich) supernova ejecta. Radiative transport calculations using these new opacities suggest that the light curves should be longer, dimmer, and redder than previously thought. The spectra appear to be pseudo-blackbody, with broad absorption features, and peak in the infrared ({approx}1 {mu}m). We discuss uncertainties in the opacities and attempt to quantify their impact on the spectral predictions. The results have important implications for observational strategies to find and study the radioactively powered electromagnetic counterparts to neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

SUPPRESSION OF DIELECTRONIC RECOMBINATION DUE TO FINITE DENSITY EFFECTS

We have developed a general model for determining density-dependent effective dielectronic recombination (DR) rate coefficients in order to explore finite-density effects on the ionization balance of plasmas. Our model consists of multiplying by a suppression factor those highly-accurate total zero-density DR rate coefficients which have been produced from state-of-the-art theoretical calculations and which have been benchmarked by experiment. The suppression factor is based upon earlier detailed collision-radiative calculations which were made for a wide range of ions at various densities and temperatures, but used a simplified treatment of DR. A general suppression formula is then developed as a function of isoelectronic sequence, charge, density, and temperature. These density-dependent effective DR rate coefficients are then used in the plasma simulation code Cloudy to compute ionization balance curves for both collisionally ionized and photoionized plasmas at very low (n{sub e} = 1 cm{sup -3}) and finite (n{sub e} = 10{sup 10} cm{sup -3}) densities. We find that the denser case is significantly more ionized due to suppression of DR, warranting further studies of density effects on DR by detailed collisional-radiative calculations which utilize state-of-the-art partial DR rate coefficients. This is expected to impact the predictions of the ionization balance in denser cosmic gases such as those found in nova and supernova shells, accretion disks, and the broad emission line regions in active galactic nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗

Collisional Ionization Equilibrium for Optically Thin Plasmas

Reliably interpreting spectra from electron-ionized cosmic plasmas requires accurate ionization balance calculations for the plasma in question. However, much of the atomic data needed for these calculations have not been generated using modern theoretical methods and their reliability are often highly suspect. We have utilized state-of-the-art calculations of dielectronic recombination (DR) rate coefficients for the hydrogenic through Na-like ions of all elements from He to Zn. We have also utilized state-of-the-art radiative recombination (RR) rate coefficient calculations for the bare through Na-like ions of all elements from H to Zn. Using our data and the recommended electron impact ionization data of Mazzotta et al. (1998), we have calculated improved collisional ionization equilibrium calculations. We compare our calculated fractional ionic abundances using these data with those presented by Mazzotta et al. (1998) for all elements from H to Ni, and with the fractional abundances derived from the modern DR and RR calculations of Gu (2003a,b, 2004) for Mg, Si, S, Ar, Ca, Fe, and Ni.

Bryans, P.↗

Calculation of Atomic Data for NASA Missions

The interpretation of cosmic spectra relies on a vast sea of atomic data which are not readily obtainable from analytic expressions or simple calculations. Rather, their evaluation typically requires state-of-the-art atomic physics calculations, with the inclusion of weaker effects (spin-orbit and configuration interactions, relaxation, Auger broadening, etc.), to achieve the level of accuracy needed for use by astrophysicists. Our NASA-supported research program is focused on calculating data for three important atomic processes, 1) dielectronic recombination (DR), 2) inner-shell photoabsorption, and 3) fluorescence and Auger decay of inner-shell vacancy states. Some additional details and examples of our recent findings are given.

Gorczyca, T. W.↗

The Strings of Eta Carina: The HST/STIS Spectra and [Ca II]

Long linear, filamentary ejecta, are found to move at very high velocity external to the Homunculus, the circumstellar hourglass-shaped ejecta surrounding Eta Carinae. The origin of the strings is a puzzle. As an example, the Weigelt Blobs have N at 10X solar and C, O at 0.01X solar abundance, along with He/H significantly enhanced. This abundance pattern is evidence for extreme CNO-processing. Similarly, the Strontium Filament has Ti/Ni at 100X solar, presumably due to the lack of oxygen to form Ti-oxide precipitates onto dust grains. We have obtained 2-D spectra with the HST/STIS of the Strontium Filament and a portion of a string. These deep spectral exposures, at moderate dispersion, span much of the near red spectral region from 5000 to 9000A. We have identified twelve emission lines in these spectra with proper velocities and spatial structure of this string and obtained line ratios for [Ca II] (7293/7325A) and [Fe Ill (7157/8619A) which are useful for determining physical conditions in this nebulosity. In an attempt to use the [Ca II] ratio to determine the physical parameters, and ultimately the abundances in the strings, we have constructed a statistical equilibrium model for Ca II , including radiative and collisional rates. These results incorporate our newly calculated atomic data for levels n = 3,4,5 and 6 configurations of Ca II. The aim is to compute the [Ca II] line ratios and use them as a diagnostic of the physical parameters. Using the [Fe II] ratio we find that for Te=10,000 K, the electron density is Ne approx.10(exp 6)/cu cm. We plan to use the [Ca II] ratio to confirm this result. Then, we will extend the use of this multilevel model Ca II atom to study the physical conditions of the Strontium filament where eight lines of Ca II, both allowed and forbidden, had been identified. With the physical conditions determined, we will be able to derive reliable estimates for the gas phase abundances in the strings.

Melendez, M. B.↗