-
Satellite lines from autoionizing states of Fe XVI and the problems with the X-ray Fe XVII lines
Authors:
G. Del Zanna,
N. R. Badnell,
P. J. Storey
Abstract:
We present new calculations of atomic data needed to model autoionizing states of Fe XVI. We compare the state energies, radiative and excitation data with a sample of results from previous literature. We find a large scatter of results, the most significant ones in the autoionization rates, which are very sensitive to the configuration interaction and state mixing. We find relatively good agreeme…
▽ More
We present new calculations of atomic data needed to model autoionizing states of Fe XVI. We compare the state energies, radiative and excitation data with a sample of results from previous literature. We find a large scatter of results, the most significant ones in the autoionization rates, which are very sensitive to the configuration interaction and state mixing. We find relatively good agreement between the autoionization rates and the collisional excitation rates calculated with the R-matrix suite of programs and autostructure. The largest model, which includes J-resolved states up to n=10, produces ab-initio wavelengths and intensities of the satellite lines which agree well with solar high-resolution spectra of active regions, with few minor wavelength adjustements. We review previous literature, finding many incorrect identifications, most notably those in the NIST database. We provide several new tentative identifications in the 15-15.7 A range, and several new ones at shorter wavelengths, where previous lines were unidentified. Compared to the previous CHIANTI model, the present one has an increased flux in the 15--15.7 A range at 2 MK of a factor of 1.9, resolving the discrepancies found in the analysis of the Marshall Grazing Incidence X-Ray Spectrometer (MaGIXS) observation. It appears that the satellite lines also resolve the long-standing discrepancy in the intensity of the important Fe XVII 3D line at 15.26 A.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
R-matrix electron-impact excitation data for the H- and He-like ions with Z=6-30
Authors:
Junjie Mao,
G. Del Zanna,
Liyi Gu,
Chunyu Zhang,
N. R. Badnell
Abstract:
Plasma models built on extensive atomic data are essential to interpreting the observed cosmic spectra. H-like Lyman series and He-like triplets observable in the X-ray band are powerful diagnostic lines to measure the physical properties of various types of astrophysical plasmas. Electron-impact excitation is a fundamental atomic process for the formation of H-like and He-like key diagnostic line…
▽ More
Plasma models built on extensive atomic data are essential to interpreting the observed cosmic spectra. H-like Lyman series and He-like triplets observable in the X-ray band are powerful diagnostic lines to measure the physical properties of various types of astrophysical plasmas. Electron-impact excitation is a fundamental atomic process for the formation of H-like and He-like key diagnostic lines. Electron-impact excitation data adopted by the widely used plasma codes (AtomDB, CHIANTI, and SPEX) do not necessarily agree with each other. Here we present a systematic calculation of electron-impact excitation data of H-like and He-like ions with the atomic number Z=6-30 (i.e., C to Zn). Radiation damped R-matrix intermediate coupling frame transformation calculation was performed for each ion with configurations up to $n=6$. We compare the present work with the above three plasma codes and literature to assess the quality of the new data, which are relevant for current and future high-resolution X-ray spectrometers.
△ Less
Submitted 19 October, 2022;
originally announced October 2022.
-
Spectroscopic r-Process Abundance Retrieval for Kilonovae I: The Inferred Abundance Pattern of Early Emission from GW170817
Authors:
Nicholas Vieira,
John J. Ruan,
Daryl Haggard,
Nicole Ford,
Maria R. Drout,
Rodrigo Fernández,
N. R. Badnell
Abstract:
Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computation…
▽ More
Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computationally challenging. We introduce Spectroscopic r-Process Abundance Retrieval for Kilonovae (SPARK), a modular framework to perform Bayesian inference on kilonova spectra with the goals of inferring elemental abundance patterns and identifying absorption features at early times. SPARK inputs an atomic line list and abundance patterns from reaction network calculations into the TARDIS radiative transfer code. It then performs fast Bayesian inference on observed kilonova spectra by training a Gaussian process surrogate for the approximate posteriors of kilonova ejecta parameters, via active learning. We use the spectrum of GW170817 at 1.4 days to perform the first inference on a kilonova spectrum, and recover a complete abundance pattern. Our inference shows that this ejecta was generated by an r-process with either (1) high electron fraction Y_e ~ 0.35 and high entropy s/k_B ~ 25, or, (2) a more moderate Y_e ~ 0.30 and s/k_B ~ 14. These parameters are consistent with a shocked, polar dynamical component, and a viscously-driven outflow from a remnant accretion disk, respectively. We also recover previous identifications of strontium absorption at ~8000 AA, and tentatively identify yttrium and/or zirconium at < 4500 AA. Our approach will enable computationally-tractable inference on the spectra of future kilonovae discovered through multi-messenger observations.
△ Less
Submitted 22 December, 2022; v1 submitted 14 September, 2022;
originally announced September 2022.
-
H, He-like recombination spectra IV; H, He-like recombination spectra IV: clarification and refinement of methodology for $l$-changing collisions
Authors:
N. R. Badnell,
F. Guzmán,
S. Brodie,
R. J. R. Williams,
P. A. M. van Hoof,
M. Chatzikos,
G. J. Ferland
Abstract:
Precise spectral diagnostic modelling of H~{\sc i} and He~{\sc ii} recombination spectra can constrain theoretical models which describe many astrophysical environments. Simple analytic expressions are of interest for collisional $l$-changing rate coefficients that are used by large-scale population modelling codes. We review, clarify and improve-upon the modified Pengelly \& Seaton formulae of Gu…
▽ More
Precise spectral diagnostic modelling of H~{\sc i} and He~{\sc ii} recombination spectra can constrain theoretical models which describe many astrophysical environments. Simple analytic expressions are of interest for collisional $l$-changing rate coefficients that are used by large-scale population modelling codes. We review, clarify and improve-upon the modified Pengelly \& Seaton formulae of Guzmán \etal We show that the recent poor results for it shown by Vrinceanu \etal are due to their misinterpretation of its usage. We also detail efficient numerical algorithms which should enable the full quantum mechanical expression for such rate coefficients to be used much more routinely by modelling codes. We illustrate with some collisional-radiative population modelling for hydrogen.
△ Less
Submitted 20 August, 2021;
originally announced August 2021.
-
R-matrix electron-impact excitation data for the O-like iso-electronic sequence
Authors:
Junjie Mao,
N. R. Badnell,
G. Del Zanna
Abstract:
Astrophysical plasma codes are built on atomic databases. In the current atomic databases, R-matrix electron-impact excitation data of O-like ions are limited. The accuracy of plasma diagnostics with O-like ions depends on the availability and accuracy of the atomic data. This is particularly relevant in the context of future observatories equipped with the next generation of high-resolution spect…
▽ More
Astrophysical plasma codes are built on atomic databases. In the current atomic databases, R-matrix electron-impact excitation data of O-like ions are limited. The accuracy of plasma diagnostics with O-like ions depends on the availability and accuracy of the atomic data. This is particularly relevant in the context of future observatories equipped with the next generation of high-resolution spectrometers. To obtain level-resolved effective collision strengths of O-like ions from \ion{Ne}{III} to \ion{Zn}{XXIII} (i.e., Ne$^{2+}$ to Zn$^{22+}$) over a wide range of temperatures. This includes transitions up to $nl=5d$ for each ion. We also aim to assess the accuracy of the new data, as well as their impact on solar atmosphere plasma diagnostics, compared to those available within the CHIANTI database. A large-scale R-matrix intermediate coupling frame transformation calculations were performed systematically for the O-like iso-electronic sequence. For each ion, 630 fine-structure levels were included in both the configuration interaction target and close-coupling collision expansions. The present results (energy levels, oscillator strengths, and effective collision strengths) of selected ions across the iso-electronic sequence are compared with those in archival databases and the literature. For selected ions across the iso-electronic sequence. We find general agreement with the few previous R-matrix calculations of collision strengths. We illustrate the improvements for a few solar plasma diagnostics over existing CHIANTI atomic models based on distorted wave data. The electron-impact excitation data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf04 and will be available in OPEN-ADAS.
△ Less
Submitted 14 July, 2021;
originally announced July 2021.
-
The influence of photo-induced processes and charge transfer on carbon and oxygen in the lower solar atmosphere
Authors:
R. P. Dufresne,
G. Del Zanna,
N. R. Badnell
Abstract:
To predict line emission in the solar atmosphere requires models which are fundamentally different depending on whether the emission is from the chromosphere or the corona. At some point between the two regions, there must be a change between the two modelling regimes. Recent extensions to the coronal modelling for carbon and oxygen lines in the solar transition region have shown improvements in t…
▽ More
To predict line emission in the solar atmosphere requires models which are fundamentally different depending on whether the emission is from the chromosphere or the corona. At some point between the two regions, there must be a change between the two modelling regimes. Recent extensions to the coronal modelling for carbon and oxygen lines in the solar transition region have shown improvements in the emission of singly- and doubly-charged ions, along with Li-like ions. However, discrepancies still remain, particularly for singly-charged ions and intercombination lines. The aim of this work is to explore additional atomic processes that could further alter the charge state distribution and the level populations within ions, in order to resolve some of the discrepancies. To this end, excitation and ionisation caused by both the radiation field and by atom-ion collisions have been included, along with recombination through charge transfer. The modelling is carried out using conditions which would be present in the quiet Sun, which allows an assessment of the part atomic processes play in changing coronal modelling, separately from dynamic and transient events taking place in the plasma. The effect the processes have on the fractional ion populations are presented, as well as the change in level populations brought about by the new excitation mechanisms. Contribution functions of selected lines from low charge states are also shown, to demonstrate the extent to which line emission in the lower atmosphere could be affected by the new modelling.
△ Less
Submitted 18 February, 2021;
originally announced February 2021.
-
Quantitative comparison of opacities calculated using the $R$-matrix and Distorted-Wave methods: Fe XVII
Authors:
F. Delahaye,
C. P. Connor,
R. T. Smith,
N. R. Badnell
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-…
▽ More
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. The present study defines our path to the next generation of opacities and opacity tables for stellar modeling.
△ Less
Submitted 5 January, 2021;
originally announced January 2021.
-
R-matrix electron-impact excitation data for the N-like iso-electronic sequence
Authors:
Junjie Mao,
N. R. Badnell,
G. Del Zanna
Abstract:
Spectral lines from N-like ions can be used to measure the temperature and density of various types of astrophysical plasmas. The atomic databases of astrophysical plasma modelling codes still have room for improvement in their electron-impact excitation data sets for N-like ions, especially $R$-matrix data. This is particularly relevant for future observatories (e.g. Arcus) which will host high-r…
▽ More
Spectral lines from N-like ions can be used to measure the temperature and density of various types of astrophysical plasmas. The atomic databases of astrophysical plasma modelling codes still have room for improvement in their electron-impact excitation data sets for N-like ions, especially $R$-matrix data. This is particularly relevant for future observatories (e.g. Arcus) which will host high-resolution spectrometers. We aim to obtain level-resolved effective collision strengths for all transitions up to $nl=5d$ over a wide range of temperatures for N-like ions from O II to Zn XXIV (i.e., O$^{+}$ to Zn$^{23+}$) and to assess the accuracy of the present work. We also examine the impact of our new data on plasma diagnostics by modelling solar observations with CHIANTI. We have carried-out systematic $R$-matrix calculations for N-like ions which included 725 fine-structure target levels in both the configuration interaction target and close-coupling collision expansions. The $R$-matrix intermediate coupling frame transformation method was used to calculate the collision strengths, while the AUTOSTRUCTURE code was used for the atomic structures. We compare the present results for selected ions with those in archival databases and the literature. The comparison covers energy levels, oscillator strengths, and effective collision strengths. We show examples of improved plasma diagnostics when compared to CHIANTI models which use only distorted wave data as well as some which use previous $R$-matrix data. The electron-impact excitation data are archived according to the Atomic Data and Analysis Structure (ADAS) data class it adf04 and will be available in OPEN-ADAS. The data can be used to improve the atomic databases for astrophysical plasma diagnostics.
△ Less
Submitted 20 September, 2020;
originally announced September 2020.
-
Effects of density on the oxygen ionisation equilibrium in collisional plasmas
Authors:
R. P. Dufresne,
G. Del Zanna,
N. R. Badnell
Abstract:
The ion populations most frequently adopted for diagnostics in collisional plasmas are derived from the density independent, coronal approximation. In higher density, lower temperature conditions, ionisation rates are enhanced once metastable levels become populated, and recombination rates are suppressed if ions recombine into Rydberg levels. As a result, the formation temperatures of ions shift,…
▽ More
The ion populations most frequently adopted for diagnostics in collisional plasmas are derived from the density independent, coronal approximation. In higher density, lower temperature conditions, ionisation rates are enhanced once metastable levels become populated, and recombination rates are suppressed if ions recombine into Rydberg levels. As a result, the formation temperatures of ions shift, altering the diagnostics of the plasma. To accurately model the effect of ionisation from metastable levels, new electron impact, ionisation cross sections have been calculated for oxygen, both for direct ionisation and excitation--auto-ionisation of the ground and metastable levels. The results have been incorporated into collisional radiative modelling to show how the ionisation equilibrium of oxygen changes once metastable levels become populated. Suppression of dielectronic recombination has been estimated and also included in the modelling, demonstrating the shifts with density in comparison to the coronal approximation. The final results for the ionisation equilibrium are used in differential emission measure modelling to predict line intensities for many lines emitted by O II-VI in the solar transition region. The predictions show improved agreement by 15-40% for O II, O VI and the inter-combination lines of O III-V, when compared to results from coronal approximation modelling. While there are still discrepancies with observations of these lines, this could, to a large part, be explained by variability in the observations.
△ Less
Submitted 1 July, 2020;
originally announced July 2020.
-
Helium line emissivities in the solar corona
Authors:
G. Del Zanna,
P. J. Storey,
N. R. Badnell,
V. Andretta
Abstract:
We present new collisional-radiative models (CRMs) for helium in the quiescent solar corona, and predict the emissivities of the He and He$^+$ lines to be observed by DKIST, Solar Orbiter, and Proba-3. We discuss in detail the rates we selected for these models, highlighting several shortcomings we have found in previous work. As no previous complete and self-consistent coronal CRM for helium exis…
▽ More
We present new collisional-radiative models (CRMs) for helium in the quiescent solar corona, and predict the emissivities of the He and He$^+$ lines to be observed by DKIST, Solar Orbiter, and Proba-3. We discuss in detail the rates we selected for these models, highlighting several shortcomings we have found in previous work. As no previous complete and self-consistent coronal CRM for helium existed, we have benchmarked our largest model at a density of 10$^{6}$ cm$^{-3}$ and temperature of 20,000 K against recent CRMs developed for photoionised nebulae. We then present results for the outer solar corona, using new dielectronic recombination rates we have calculated, which increase the abundance of neutral helium by about a factor of two. We also find that all the optical triplet He I lines, and in particular the well known He I 10830 and 5876 A lines are strongly affected by both photo-excitation and photo-ionisation from the disk radiation, and that extensive CRM models are required to obtain correct estimates. Close to the Sun, at an electron density of 10$^{8}$ cm$^{-3}$ and temperature of 1 MK, we predict the emissivity of the He I 10830 A to be comparable to that of the strong Fe XIII coronal line at 10798 A. However, we expect the He I emissivity to sharply fall in the outer corona, with respect to Fe XIII. We confirm that the He$^+$ Lyman $α$ at 304 A is also significantly affected by photo-excitation and is expected to be detectable as a strong coronal line up to several solar radii.
△ Less
Submitted 16 June, 2020;
originally announced June 2020.
-
Fe III emission in quasars: evidence for a dense turbulent medium
Authors:
Matthew J. Temple,
Gary J. Ferland,
Amy L. Rankine,
Paul C. Hewett,
N. R. Badnell,
C. P. Ballance,
G. Del Zanna,
R. P. Dufresne
Abstract:
Recent improvements to atomic energy-level data allow, for the first time, accurate predictions to be made for the Fe III line emission strengths in the spectra of luminous, $L_\text{bol}=10^{46}-10^{48}$ erg/s, Active Galactic Nuclei. The Fe III emitting gas must be primarily photoionized, consistent with observations of line reverberation. We use CLOUDY models exploring a wide range of parameter…
▽ More
Recent improvements to atomic energy-level data allow, for the first time, accurate predictions to be made for the Fe III line emission strengths in the spectra of luminous, $L_\text{bol}=10^{46}-10^{48}$ erg/s, Active Galactic Nuclei. The Fe III emitting gas must be primarily photoionized, consistent with observations of line reverberation. We use CLOUDY models exploring a wide range of parameter space, together with 26,500 rest-frame ultraviolet spectra from the Sloan Digital Sky Survey, to constrain the physical conditions of the line emitting gas. The observed Fe III emission is best accounted for by dense ($n_H=10^{14}$ cm$^{-3}$) gas which is microturbulent, leading to smaller line optical depths and fluorescent excitation. Such high density gas appears to be present in the central regions of the majority of luminous quasars. Using our favoured model, we present theoretical predictions for the relative strengths of the Fe III UV34 $λλ$1895,1914,1926 multiplet. This multiplet is blended with the Si III] $λ$1892 and C III] $λ$1909 emission lines and an accurate subtraction of UV34 is essential when using these lines to infer information about the physics of the broad line region in quasars.
△ Less
Submitted 30 June, 2020; v1 submitted 15 June, 2020;
originally announced June 2020.
-
R-matrix electron-impact excitation data for the C-like iso-electronic sequence
Authors:
Junjie Mao,
N. R. Badnell,
G. Del Zanna
Abstract:
Emission and absorption features from C-like ions serve as temperature and density diagnostics of astrophysical plasmas. $R$-matrix electron-impact excitation data sets for C-like ions in the literature merely cover a few ions, and often only for the ground configuration. Our goal is to obtain level-resolved effective collision strength over a wide temperature range for C-like ions from \ion{N}{II…
▽ More
Emission and absorption features from C-like ions serve as temperature and density diagnostics of astrophysical plasmas. $R$-matrix electron-impact excitation data sets for C-like ions in the literature merely cover a few ions, and often only for the ground configuration. Our goal is to obtain level-resolved effective collision strength over a wide temperature range for C-like ions from \ion{N}{II} to \ion{Kr}{XXXI} (i.e., N$^{+}$ to Kr$^{30+}$) with a systematic set of $R$-matrix calculations. We also aim to assess their accuracy. For each ion, we included a total of 590 fine-structure levels in both the configuration interaction target and close-coupling collision expansion. These levels arise from 24 configurations $2l^3 nl^{\prime}$ with $n=2-4$, $l=0-1$, and $l^{\prime}=0-3$ plus the three configurations $2s^22p5l$ with $l=0-2$. The AUTOSTRUCTURE code was used to calculate the target structure. Additionally, the $R$-matrix intermediate coupling frame transformation method was used to calculate the collision strengths. We compare the present results of selected ions with archival databases and results in the literature. The comparison covers energy levels, transition rates, and effective collision strengths. We illustrate the impact of using the present results on an \ion{Ar}{xiii} density diagnostic for the solar corona. The electron-impact excitation data is archived according to the Atomic Data and Analysis Structure (ADAS) data class adf04 and will be available in OPEN-ADAS. The data will be incorporated into spectral codes, such as CHIANTI and SPEX, for plasma diagnostics.
△ Less
Submitted 18 December, 2019;
originally announced December 2019.
-
H-, He-like recombination spectra III: $n$-changing collisions in highly-excited Rydberg states and their impact on the radio, IR and optical recombination lines
Authors:
F. Guzmán,
M. Chatzikos,
P. A. M. van Hoof,
Dana S. Balser,
M. Dehghanian,
N. R. Badnell,
G. J. Ferland
Abstract:
At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact para…
▽ More
At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact parameter calculations or semi-empirical fits and the rate coefficients agree within a factor of two. We show in this paper that these uncertainties cause errors lower than 5% in the emission of radio recombination lines (RRL) of most ionized plasmas of typical nebulae. However, in special circumstances where the transitions between Rydberg levels are amplified by maser effects, the errors can increase up to 20%. We present simulations of the optical depth and H$nα$ line emission of Active Galactic Nuclei (AGN) Broad Line Regions (BLRs) and the Orion Nebula Blister to showcase our findings.
△ Less
Submitted 18 March, 2019; v1 submitted 13 March, 2019;
originally announced March 2019.
-
Uncertainties on atomic data. A case study: N IV
Authors:
G. Del Zanna,
L. Fernandez-Menchero,
N. R. Badnell
Abstract:
We consider three recent large-scale calculations for the radiative and electron-impact excitation data of N IV, carried out with different methods and codes. The scattering calculations employed the relativistic Dirac $R$-matrix (DARC) method, the intermediate coupling frame transformation (ICFT) $R$-matrix method, and the B-spline $R$-matrix (BSR) method. These are all large-scale scattering cal…
▽ More
We consider three recent large-scale calculations for the radiative and electron-impact excitation data of N IV, carried out with different methods and codes. The scattering calculations employed the relativistic Dirac $R$-matrix (DARC) method, the intermediate coupling frame transformation (ICFT) $R$-matrix method, and the B-spline $R$-matrix (BSR) method. These are all large-scale scattering calculations with well-tested and sophisticated codes, which use the same set of target states. One concern raised in previous literature is related to the increasingly large discrepancies in the effective collision strengths between the three sets of calculations for increasingly weak and/or high-lying transitions. We have built three model ions and calculated the intensities of all the main spectral lines in this ion. We have found that, despite such large differences, excellent agreement (to within $\pm$~20\%) exists between all the spectroscopically-relevant line intensities. This provides confidence in the reliability of the calculations for plasma diagnostics. We have used the differences in the radiative and excitation rates amongst the three sets of calculations to obtain a measure of the uncertainty in each rate. Using a Monte Carlo approach, we have shown how these uncertainties affect the main theoretical ratios which are used to measure electron densities and temperatures.
△ Less
Submitted 24 January, 2019;
originally announced January 2019.
-
A Quantitative Comparison of Opacities Calculated Using the Distorted- Wave and $\boldsymbol{R}$-Matrix Methods
Authors:
F. Delahaye,
N. R. Badnell,
C. P. Ballance,
P. Palmeri,
S. Preval,
P. Quinet,
C. Ramsbottom,
R. T. Smyth,
M. Turkington,
C. J. Zeippen
Abstract:
The present debate on the reliability of astrophysical opacities has reached a new climax with the recent measurements of Fe opacities on the Z-machine at the Sandia National Laboratory \citep{Bailey2015}. To understand the differences between theoretical results, on the one hand, and experiments on the other, as well as the differences among the various theoretical results, detailed comparisons a…
▽ More
The present debate on the reliability of astrophysical opacities has reached a new climax with the recent measurements of Fe opacities on the Z-machine at the Sandia National Laboratory \citep{Bailey2015}. To understand the differences between theoretical results, on the one hand, and experiments on the other, as well as the differences among the various theoretical results, detailed comparisons are needed. Many ingredients are involved in the calculation of opacities; deconstructing the whole process and comparing the differences at each step are necessary to quantify their importance and impact on the final results. We present here such a comparison using the two main approaches to calculate the required atomic data, the $R$-Matrix and distorted-wave methods, as well as sets of configurations and coupling schemes to quantify the effects on the opacities for the $Fe\ XVII$ and $Ni\ XIV$ ions.
△ Less
Submitted 10 April, 2018;
originally announced April 2018.
-
The 2017 Release of Cloudy
Authors:
G. J. Ferland,
M. Chatzikos,
F. Guzmán,
M. L. Lykins,
P. A. M. van Hoof,
R. J. R. Williams,
N. P. Abel,
N. R. Badnell,
F. P. Keenan,
R. L. Porter,
P. C. Stancil
Abstract:
We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully…
▽ More
We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds.
△ Less
Submitted 5 October, 2023; v1 submitted 30 May, 2017;
originally announced May 2017.
-
Thermodynamically-consistent semi-classical $\ell$-changing rates
Authors:
R. J. R. Williams,
F. Guzmán,
N. R. Badnell,
P. A. M. van Hoof,
M. Chatzikos,
G. J. Ferland
Abstract:
We compare the results of the semi-classical (SC) and quantum-mechanical (QM) formalisms for angular-momentum changing transitions in Rydberg atom collisions given by Vrinceanu & Flannery, J. Phys. B 34, L1 (2001), and Vrinceanu, Onofrio & Sadeghpour, ApJ 747, 56 (2012), with those of the SC formalism using a modified Monte Carlo realization. We find that this revised SC formalism agrees well with…
▽ More
We compare the results of the semi-classical (SC) and quantum-mechanical (QM) formalisms for angular-momentum changing transitions in Rydberg atom collisions given by Vrinceanu & Flannery, J. Phys. B 34, L1 (2001), and Vrinceanu, Onofrio & Sadeghpour, ApJ 747, 56 (2012), with those of the SC formalism using a modified Monte Carlo realization. We find that this revised SC formalism agrees well with the QM results. This provides further evidence that the rates derived from the QM treatment are appropriate to be used when modelling recombination through Rydberg cascades, an important process in understanding the state of material in the early universe. The rates for $Δ\ell=\pm1$ derived from the QM formalism diverge when integrated to sufficiently large impact parameter, $b$. Further to the empirical limits to the $b$ integration suggested by Pengelly & Seaton, MNRAS 127, 165 (1964), we suggest that the fundamental issue causing this divergence in the theory is that it does not fully cater for the finite time taken for such distant collisions to complete.
△ Less
Submitted 27 April, 2017;
originally announced April 2017.
-
Testing atomic collision theory with the two-photon continuum of astrophysical nebulae
Authors:
F. Guzmán,
N. R. Badnell,
M. Chatzikos,
P. A. M. van Hoof,
R. J. R. Williams,
G. J. Ferland
Abstract:
Accurate rates for energy-degenerate l-changing collisions are needed to determine cosmological abundances and recombination. There are now several competing theories for the treatment of this process, and it is not possible to test these experimentally. We show that the H I two-photon continuum produced by astrophysical nebulae is strongly affected by l-changing collisions. We perform an analysis…
▽ More
Accurate rates for energy-degenerate l-changing collisions are needed to determine cosmological abundances and recombination. There are now several competing theories for the treatment of this process, and it is not possible to test these experimentally. We show that the H I two-photon continuum produced by astrophysical nebulae is strongly affected by l-changing collisions. We perform an analysis of the different underlying atomic processes and simulate the recombination and two-photon spectrum of a nebula containing H and He. We provide an extended set of effective recombination coefficients and updated l-changing 2s-2p transition rates using several competing theories. In principle, accurate astronomical observations could determine which theory is correct.
△ Less
Submitted 26 January, 2017;
originally announced January 2017.
-
The electron energy-loss rate due to radiative recombination
Authors:
Junjie Mao,
Jelle Kaastra,
N. R. Badnell
Abstract:
For photoionized plasmas, electron energy-loss rates due to radiative recombination (RR) are required for thermal equilib- rium calculations, which assume a local balance between the energy gain and loss. While many calculations of total and/or partial RR rates are available from literature, specific calculations of associated RR electron energy-loss rates are lacking. Here we focus on electron en…
▽ More
For photoionized plasmas, electron energy-loss rates due to radiative recombination (RR) are required for thermal equilib- rium calculations, which assume a local balance between the energy gain and loss. While many calculations of total and/or partial RR rates are available from literature, specific calculations of associated RR electron energy-loss rates are lacking. Here we focus on electron energy-loss rates due to radiative recombination of H-like to Ne-like ions for all the elements up to and including zinc (Z = 30), over a wide temperature range. We use the AUTOSTRUCTURE code to calculate the level-resolved photoionization cross section and modify the ADASRR code so that we can simultaneously obtain level-resolved RR rate coefficients and associated RR electron energy-loss rate coefficients. The total RR rates and electron energy-loss rates of H i and He i are compared with those found in literature. Furthermore, we utilize and parameterize the weighted electron energy-loss factors (dimensionless) to characterize total electron energy-loss rates due to RR. The RR electron energy-loss data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf48. The RR electron energy-loss data are also incorporated into the SPEX code for detailed modelling of photoionized plamsas.
△ Less
Submitted 3 December, 2016;
originally announced December 2016.
-
Hot DA white dwarf model atmosphere calculations: Including improved Ni PI cross sections
Authors:
S. P. Preval,
M. A. Barstow,
N. R. Badnell,
I. Hubeny,
J. B. Holberg
Abstract:
To calculate realistic models of objects with Ni in their atmospheres, accurate atomic data for the relevant ionization stages needs to be included in model atmosphere calculations. In the context of white dwarf stars, we investigate the effect of changing the Ni {\sc iv}-{\sc vi} bound-bound and bound-free atomic data has on model atmosphere calculations. Models including PICS calculated with {\s…
▽ More
To calculate realistic models of objects with Ni in their atmospheres, accurate atomic data for the relevant ionization stages needs to be included in model atmosphere calculations. In the context of white dwarf stars, we investigate the effect of changing the Ni {\sc iv}-{\sc vi} bound-bound and bound-free atomic data has on model atmosphere calculations. Models including PICS calculated with {\sc autostructure} show significant flux attenuation of up to $\sim 80$\% shortward of 180Å\, in the EUV region compared to a model using hydrogenic PICS. Comparatively, models including a larger set of Ni transitions left the EUV, UV, and optical continua unaffected. We use models calculated with permutations of this atomic data to test for potential changes to measured metal abundances of the hot DA white dwarf G191-B2B. Models including {\sc autostructure} PICS were found to change the abundances of N and O by as much as $\sim 22$\% compared to models using hydrogenic PICS, but heavier species were relatively unaffected. Models including {\sc autostructure} PICS caused the abundances of N/O {\sc iv} and {\sc v} to diverge. This is because the increased opacity in the {\sc autostructure} PICS model causes these charge states to form higher in the atmosphere, moreso for N/O {\sc v}. Models using an extended line list caused significant changes to the Ni {\sc iv}-{\sc v} abundances. While both PICS and an extended line list cause changes in both synthetic spectra and measured abundances, the biggest changes are caused by using {\sc autostructure} PICS for Ni.
△ Less
Submitted 30 October, 2016;
originally announced October 2016.
-
H, He-like recombination spectra II: $l$-changing collisions for He Rydberg states
Authors:
F. Guzmán,
N. R. Badnell,
R. J. R. Williams,
P. A. M. van Hoof,
M. Chatzikos,
G. J. Ferland
Abstract:
Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which foc…
▽ More
Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $\sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy.
△ Less
Submitted 31 August, 2016;
originally announced August 2016.
-
Scaling of collision strengths for highly-excited states of ions of the H- and He-like sequences
Authors:
L. Fernandez-Menchero,
G. Del Zanna,
N. R. Badnell
Abstract:
Emission lines from highly-excited states (n >= 5) of H- and He-like ions have been detected in astrophysical sources and fusion plasmas. For such excited states, R-matrix or distorted wave calculations for electron-impact excitation are very limited, due to the large size of the atomic basis set needed to describe them. Calculations for n >= 6 are also not generally available. We study the behavi…
▽ More
Emission lines from highly-excited states (n >= 5) of H- and He-like ions have been detected in astrophysical sources and fusion plasmas. For such excited states, R-matrix or distorted wave calculations for electron-impact excitation are very limited, due to the large size of the atomic basis set needed to describe them. Calculations for n >= 6 are also not generally available. We study the behaviour of the electron-impact excitation collision strengths and effective collision strengths for the most important transitions used to model electron collision dominated astrophysical plasmas, solar, for example. We investigate the dependence on the relevant parameters: the principal quantum number n or the nuclear charge Z. We also estimate the importance of coupling to highly-excited states and the continuum by comparing the results of different sized calculations. We provide analytic formulae to calculate the electron-impact excitation collision strengths and effective collision strengths to highly-excited states (n >= 8) of H- and He-like ions. These extrapolated effective collision strengths can be used to interpret astrophysical and fusion plasma via collisional-radiative modelling.
△ Less
Submitted 8 July, 2016;
originally announced July 2016.
-
H, He-like recombination spectra I: $l$-changing collisions for hydrogen
Authors:
F. Guzmán,
N. R. Badnell,
R. J. R. Williams,
P. A. M. van Hoof,
M. Chatzikos,
G. J. Ferland
Abstract:
Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe di…
▽ More
Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe dipole $l$-changing collisions ($nl\rightarrow nl'=l\pm 1$) between energy-degenerate states within an $n$-shell. The work of Pengelly & Seaton (1964) has, for half-a-century, been considered the definitive study which "solved" the problem. Recent work by Vrinceanu et al.(2012) recommended the use of rate coefficients from a semi-classical approximation which are nearly an order of magnitude smaller than those of Pengelly & Seaton (1964), with the result that significantly higher densities are needed for the $nl$ populations to come into local thermodynamic equilibrium.
Here, we compare predicted H~I emissivities from the two works and find widespread differences, of up to $\approx 10$%. This far exceeds the 1% precision required to obtain the primordial He/H abundance ratio from observations so as to constrain Big Bang cosmologies. We recommend using the rate coefficients of Pengelly & Seaton (1964) for $l$-changing collisions, to describe the H recombination spectrum, based-on their quantum mechanical representation of the long-range dipole interaction.
△ Less
Submitted 13 April, 2016;
originally announced April 2016.
-
Importance of the completeness of the configuration interaction and close coupling expansions in R-matrix calculations for highly-charged ions: electron-impact excitation of Fe $^{20+}$
Authors:
L Fernandez-Menchero,
A S Giunta,
G Del Zanna,
N R Badnell
Abstract:
We have carried-out two intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of C-like Fe $^{20+}$, both of which use the same expansions for their configuration interaction (CI) and close-coupling (CC) representations. The first expansion arises from the configurations 2s$^2$ 2p$^2$, 2s 2p$^3$, 2p$^4$, {2s$^2$ 2p, 2s 2p$^2$, 2p$^3$} nl, with n…
▽ More
We have carried-out two intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of C-like Fe $^{20+}$, both of which use the same expansions for their configuration interaction (CI) and close-coupling (CC) representations. The first expansion arises from the configurations 2s$^2$ 2p$^2$, 2s 2p$^3$, 2p$^4$, {2s$^2$ 2p, 2s 2p$^2$, 2p$^3$} nl, with n=3,4, for l=0-3, which give rise to 564 CI/CC levels. The second adds configurations 2s$^2$ 2p 5l, for l=0-2, which give rise to 590 CI/CC levels in total. Comparison of oscillator strengths and effective collision strengths from these two calculations demonstrates the lack of convergence in data for n=4 from the smaller one. Comparison of results for the 564 CI/CC level calculation with an earlier ICFT R-matrix calculation which used the exact same CI expansion but truncated the CC expansion to only 200 levels demonstrates the lack of convergence of the earlier data, particularly for n=3 levels. Also, we find that the results of our 590 CC R-matrix calculation are significantly and systematically larger than those of an earlier comparable Distorted Wave-plus-resonances calculation. Thus, it is important still to take note of the (lack of) convergence in both atomic structural and collisional data, even in such a highly-charged ion as Fe $^{20+}$, and to treat resonances non-perturbatively. This is of particular importance for Fe ions given their importance in the spectroscopic diagnostic modelling of astrophysical plasmas.
△ Less
Submitted 11 April, 2016; v1 submitted 18 November, 2015;
originally announced November 2015.
-
Bootstrapping dielectronic recombination from second-row elements and the Orion Nebula
Authors:
N. R. Badnell,
G. J. Ferland,
T. W. Gorczyca,
D. Nikolic,
G. A. Wagle
Abstract:
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 tem…
▽ More
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$^{2+} \rightarrow$ S$^+$ 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$^{2+} \rightarrow$ S$^+$ DR rate coefficient at $\sim 10^{4}$~K. We present new calculations of the DR rate coefficient for S$^{2+} \rightarrow$ S$^+$ 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.
△ Less
Submitted 19 March, 2015;
originally announced March 2015.
-
Validity of the ICFT R-matrix method: Be-like Al$^{9+}$ a case study
Authors:
L. Fernández-Menchero,
G. Del Zanna,
N. R. Badnell
Abstract:
We have carried-out 98-level configuration-interaction / close-coupling (CI/CC) intermediate coupling frame transformation (ICFT) and Breit-Pauli R-matrix calculations for the electron-impact excitation of Be-like Al$^{9+}$. The close agreement that we find between the two sets of effective collision strengths demonstrates the continued robustness of the ICFT method. On the other hand, a compari…
▽ More
We have carried-out 98-level configuration-interaction / close-coupling (CI/CC) intermediate coupling frame transformation (ICFT) and Breit-Pauli R-matrix calculations for the electron-impact excitation of Be-like Al$^{9+}$. The close agreement that we find between the two sets of effective collision strengths demonstrates the continued robustness of the ICFT method. On the other hand, a comparison of this data with previous 238-level CI/CC ICFT effective collision strengths shows that the results for excitation up to n=4 levels are systematically and increasingly underestimated over a wide range of temperatures by R-matrix calculations whose close-coupling expansion extends only to n=4 (98-levels). Thus, we find to be false a recent conjecture that the ICFT approach may not be completely robust. The conjecture was based upon a comparison of 98-level CI/CC Dirac R-matrix effective collision strengths for Al$^{9+}$ with those from the 238-level CI/CC ICFT R-matrix calculations. The disagreement found recently is due to a lack of convergence of the close-coupling expansion in the 98-level CI/CC Dirac work. The earlier 238-level CI/CC ICFT work has a superior target to the 98-level CI/CC Dirac one and provides more accurate atomic data. Similar considerations need to be made for other Be-like ions and for other sequences.
△ Less
Submitted 20 April, 2015; v1 submitted 17 March, 2015;
originally announced March 2015.
-
A possible solution to the Lyman/Balmer line problem in hot DA white dwarfs
Authors:
S. P. Preval,
M. A. Barstow,
N. R. Badnell,
J. B. Holberg,
I. Hubeny
Abstract:
Arguably, the best method for determining the effective temperature ($T_{\mathrm{eff}}$) and surface gravity (log $g$) of a DA white dwarf is by fitting the Hydrogen Lyman and Balmer absorption features. However, as has been shown for white dwarfs with $T_{\mathrm{eff}}$>50,000K, the calculated value from the Lyman and Balmer lines are discrepant, which worsens with increasing temperature. Many di…
▽ More
Arguably, the best method for determining the effective temperature ($T_{\mathrm{eff}}$) and surface gravity (log $g$) of a DA white dwarf is by fitting the Hydrogen Lyman and Balmer absorption features. However, as has been shown for white dwarfs with $T_{\mathrm{eff}}$>50,000K, the calculated value from the Lyman and Balmer lines are discrepant, which worsens with increasing temperature. Many different solutions have been suggested, ranging from the input physics used to calculate the models, to interstellar reddening. We will focus on the former, and consider three variables. The first is the atomic data used, namely the number of transitions included in line blanketing treatments and the photoionization cross sections. The second is the stark broadening treatment used to synthesise the Lyman and Balmer line profiles, namely the calculations performed by Lemke (1997) and Tremblay & Bergeron (2009). Finally, the third is the atmospheric content. The model grids are calculated with a pure H composition, and a metal polluted composition using the abundances of Preval et al. (2013). We present the preliminary results of our analysis, whereby we have determined the $T_{\mathrm{eff}}$ for a small selection of white dwarfs. We plan to extend our analysis by allowing metallicity to vary in future model grids.
△ Less
Submitted 3 October, 2014;
originally announced October 2014.
-
Electron-impact excitation of Fe$^{2+}$: a comparison of intermediate coupling frame transformation, Breit-Pauli and Dirac R-matrix calculations
Authors:
N. R. Badnell,
C. P. Ballance
Abstract:
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 cou…
▽ More
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$^{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. [Ap.J.Lett, 718, L189, (2010)] who carried-out ICFT and Dirac R-matrix calculations for the same ion. We discuss possible reasons.
△ Less
Submitted 19 February, 2014;
originally announced February 2014.
-
Collision Strengths for [O III] Optical and Infrared Lines
Authors:
P. J. Storey,
Taha Sochi,
N. R. Badnell
Abstract:
We present electron collision strengths and their thermally averaged values for the forbidden lines of the astronomically abundant doubly-ionized oxygen ion, O^{2+}, in an intermediate coupling scheme using the Breit-Pauli relativistic terms as implemented in an R-matrix atomic scattering code. We use several atomic targets for the R-matrix scattering calculations including one with 72 atomic term…
▽ More
We present electron collision strengths and their thermally averaged values for the forbidden lines of the astronomically abundant doubly-ionized oxygen ion, O^{2+}, in an intermediate coupling scheme using the Breit-Pauli relativistic terms as implemented in an R-matrix atomic scattering code. We use several atomic targets for the R-matrix scattering calculations including one with 72 atomic terms. We also compare with new results obtained using the intermediate coupling frame transformation method. We find spectroscopically significant differences against a recent Breit-Pauli calculation for the excitation of the [O III] lambda 4363 transition but confirm the results of earlier calculations.
△ Less
Submitted 25 November, 2013;
originally announced November 2013.
-
Opacities and Spectra of the r-process Ejecta from Neutron Star Mergers
Authors:
Daniel Kasen,
N. R. Badnell,
Jennifer Barnes
Abstract:
Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements by the r-process. The subsequent radioactive decay of the nuclei can power electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying such events is an important goal of future transient surveys, offering new perspectives on the origin of r…
▽ More
Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements by the r-process. The subsequent radioactive decay of the nuclei can power electromagnetic emission similar to, but significantly dimmer than, an ordinary supernova. Identifying such events is an important goal of future transient 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 consider the opacity of expanding r-process material and argue that it is dominated by line 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 radiative data for tens of millions of lines. We find that 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 indicate that the transient emission should be dimmer and redder than previously thought. The spectra appear pseudo-blackbody, with broad absorption features, and peak in the infrared (~1 micron). 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 compact object mergers.
△ Less
Submitted 22 March, 2013;
originally announced March 2013.
-
Suppression of Dielectronic Recombination Due to Finite Density Effects
Authors:
D. Nikolić,
T. W. Gorczyca,
K. T. Korista,
G. J. Ferland,
N. R. Badnell
Abstract:
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…
▽ More
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 (ne = 1 cm^-3) and finite (ne=10^10 cm^-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.
△ Less
Submitted 10 March, 2013;
originally announced March 2013.
-
Electron-ion Recombination of Fe XII forming Fe XI: Laboratory Measurements and Theoretical Calculations
Authors:
O. Novotný,
N. R. Badnell,
D. Bernhardt,
M. Grieser,
M. Hahn,
C. Krantz,
M. Lestinsky,
A. Müller,
R. Repnow,
S. Schippers,
A. Wolf,
D. W. Savin
Abstract:
We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC…
▽ More
We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC based on a comparison of the ion beam decay rate with and without the electron beam on. For energies below 75 eV, the spectrum is dominated by dielectronic recombination (DR) resonances associated with 3s-3p and 3p-3d core excitations. At higher energies we observe contributions from 3-N' and 2-N' core excitations DR. We compare our experimental results to state-of-the-art multi-configuration Breit-Pauli (MCBP) calculations and find significant differences, both in resonance energies and strengths. We have extracted the DR contributions from the measured MBRRC data and transformed them into a plasma recombination rate coefficient (PRRC) for temperatures in the range of 10^3 to 10^7 K. We show that the previously recommended DR data for Fe XII significantly underestimate the PRRC at temperatures relevant for both photoionized plasmas (PPs) and collisionaly ionized plasmas (CPs). This is to be contrasted with our MCBP PRRC results which agree with the experiment to within 30% at PP temperatures and even better at CP temperatures. We find this agreement despite the disagreement shown by the detailed comparison between our MCBP and experimental MBRRC results. Lastly, we present a simple parameterized form of the experimentally derived PRRC for easy use in astrophysical modelling codes.
△ Less
Submitted 27 April, 2012;
originally announced April 2012.
-
R-matrix electron-impact excitation data for astrophysically abundant sulphur ions
Authors:
G. Y. Liang,
N. R. Badnell,
G. Zhao,
J. Y. Zhong,
F. L. Wang
Abstract:
We present results for the electron-impact excitation of highly-charged sulphur ions (S8+ - S11+) obtained using the intermediate-coupling frame transformation R-matrix approach. A detailed comparison of the target structure has been made for the four ions to assess the uncertainty on collision strengths from the target structure. Effective collision strengths (Υs) are presented at temperatures ra…
▽ More
We present results for the electron-impact excitation of highly-charged sulphur ions (S8+ - S11+) obtained using the intermediate-coupling frame transformation R-matrix approach. A detailed comparison of the target structure has been made for the four ions to assess the uncertainty on collision strengths from the target structure. Effective collision strengths (Υs) are presented at temperatures ranging from 2\times10^2(z+1)^2 K to 2\times10^6(z+1)^2 K (where z is the residual charge of ions. Detailed comparisons for the Υare made with the results of previous calculations for these ions, which will pose insight on the uncertainty in their usage by astrophysical and fusion modelling codes.
△ Less
Submitted 3 August, 2011; v1 submitted 30 July, 2011;
originally announced August 2011.
-
Revised Results for Non-thermal Recombination Flare Hard X-Ray Emission
Authors:
J. C. Brown,
P. C. V. Mallik,
N. R. Badnell
Abstract:
Brown and Mallik (BM) recently showed that, for hot sources, recombination of non-thermal electrons (NTR) onto highly ionised heavy ions is not negligible compared to non-thermal bremsstrahlung (NTB) as a source of flare hard X-rays (HXRs) and so should be included in modelling non-thermal HXR flare emission. In view of major discrepancies between BM results for the THERMAL continua and those of…
▽ More
Brown and Mallik (BM) recently showed that, for hot sources, recombination of non-thermal electrons (NTR) onto highly ionised heavy ions is not negligible compared to non-thermal bremsstrahlung (NTB) as a source of flare hard X-rays (HXRs) and so should be included in modelling non-thermal HXR flare emission. In view of major discrepancies between BM results for the THERMAL continua and those of the Chianti code and of RHESSI solar data, we critically re-examine and correct the BM analysis and modify the conclusions concerning the importance of NTR. Although the analytic Kramers expression used by BM is correct for the purely hydrogenic recombination cross section, the heuristic expressions used by BM to extend the Kramers expression beyond the `bare nucleus' case to which it applies had serious errors. BM results have therefore been recalculated using corrected expressions, which have been validated against the results of detailed calculations. At T ~ 10-30 MK the dominant ions are Fe 22+, 23+, 24+ for which BM erroneously overestimated NTR emission by around an order of magnitude. Contrary to the BM claim, NTR in hot flare plasmas does NOT dominate over NTB, although in some cases it can be comparable and so still very important in inversions of photon spectra to derive electron spectra, especially as NTR includes sharp edge features. The BM claim of dominance of NTR over NTB in deka-keV emission is incorrect due to a serious error in their analysis. However, the NTR contribution can still be large enough to demand inclusion in spectral fitting, the spectral edges having potentially serious effects on inversion of HXR spectra to infer fast electron spectra.
△ Less
Submitted 24 February, 2010; v1 submitted 17 December, 2009;
originally announced December 2009.
-
Radiative transition rates and collision strengths for Si II
Authors:
M. A. Bautista,
P. Quinet,
P. Palmeri,
N. R. Badnell,
J. Dunn,
N. Arav
Abstract:
Aims. This work reports radiative transition rates and electron impact excitation collision strengths for levels of the 3s23p, 3s3p2, 3s24s, and 3s23d configurations of Siii. Methods. The radiative data were computed using the Thomas-Fermi-Dirac-Amaldi central potential, but with the modifications introduced by Bautista (2008) that account for the effects of electron-electron interactions. We al…
▽ More
Aims. This work reports radiative transition rates and electron impact excitation collision strengths for levels of the 3s23p, 3s3p2, 3s24s, and 3s23d configurations of Siii. Methods. The radiative data were computed using the Thomas-Fermi-Dirac-Amaldi central potential, but with the modifications introduced by Bautista (2008) that account for the effects of electron-electron interactions. We also introduce new schemes for the optimization of the variational parameters of the potential. Additional calculations were carried out with the Relativistic Hartree-Fock and the multiconfiguration Dirac-Fock methods. Collision strengths in LS-coupling were calculated in the close coupling approximation with the R-matrix method. Then, fine structure collision strengths were obtained by means of the intermediate-coupling frame transformation (ICFT) method which accounts for spin-orbit coupling effects. Results. We present extensive comparisons between the results of different approximations and with the most recent calculations and experiment available in the literature. From these comparisons we derive a recommended set of gf- values and radiative transition rates with their corresponding estimated uncertainties. We also study the effects of different approximations in the representation of the target ion on the electron-impact collision strengths. Our most accurate set of collision strengths were integrated over a Maxwellian distribution of electron energies and the resulting effective collision strengths are given for a wide range of temperatures. Our results present significant differences from recent calculations with the B-spline non-orthogonal R-matrix method. We discuss the sources of the differences.
△ Less
Submitted 28 October, 2009;
originally announced October 2009.
-
Electron-ion Recombination of Fe X forming Fe IX and of Fe XI forming Fe X: Laboratory Measurements and Theoretical Calculations
Authors:
M. Lestinsky,
N. R. Badnell,
D. Bernhardt,
M. Grieser,
J. Hoffmann,
D. Lukić,
A. Müller,
D. A. Orlov,
R. Repnow,
D. W. Savin,
E. W. Schmidt,
M. Schnell,
S. Schippers,
A. Wolf,
D. Yu
Abstract:
We have measured electron-ion recombination for Fe$^{9+}$ forming Fe$^{8+}$ and for Fe$^{10+}$ forming Fe$^{9+}$ using merged beams at the TSR heavy-ion storage-ring in Heidelberg. The measured merged beams recombination rate coefficients (MBRRC) for relative energies from 0 to 75 eV are presented, covering all dielectronic recombination (DR) resonances associated with 3s->3p and 3p->3d core tra…
▽ More
We have measured electron-ion recombination for Fe$^{9+}$ forming Fe$^{8+}$ and for Fe$^{10+}$ forming Fe$^{9+}$ using merged beams at the TSR heavy-ion storage-ring in Heidelberg. The measured merged beams recombination rate coefficients (MBRRC) for relative energies from 0 to 75 eV are presented, covering all dielectronic recombination (DR) resonances associated with 3s->3p and 3p->3d core transitions in the spectroscopic species Fe X and Fe XI, respectively. We compare our experimental results to multi-configuration Breit-Pauli (MCBP) calculations and find significant differences. From the measured MBRRC we have extracted the DR contributions and transform them into plasma recombination rate coefficients (PRRC) for astrophysical plasmas with temperatures from 10^2 to 10^7 K. This spans across the regimes where each ion forms in photoionized or in collisionally ionized plasmas. For both temperature regimes the experimental uncertainties are 25% at a 90% confidence level. The formerly recommended DR data severely underestimated the rate coefficient at temperatures relevant for photoionized gas. At the temperatures relevant for photoionized gas, we find agreement between our experimental results and MCBP theory. At the higher temperatures relevant for collisionally ionized gas, the MCBP calculations yield a Fe XI DR rate coefficent which is significantly larger than the experimentally derived one. We present parameterized fits to our experimentally derived DR PRRC.
△ Less
Submitted 1 April, 2009;
originally announced April 2009.
-
R-matrix inner-shell electron-impact excitation of Fe$^{15+}$ including Auger-plus-radiation damping
Authors:
G Y Liang,
A D Whiteford,
N R Badnell
Abstract:
We present results for the inner-shell electron-impact excitation of Fe$^{15+}$ using the intermediate-coupling frame transformation {\it R}-matrix approach in which Auger-plus-radiation damping has been included. The target and close-coupling expansions are both taken to be the 134 levels belonging to the configurations ${\rm 2s^22p^63}l$, ${\rm 2s^22p^53s3}l$, ${\rm 2s^22p^53p^2}$ and…
▽ More
We present results for the inner-shell electron-impact excitation of Fe$^{15+}$ using the intermediate-coupling frame transformation {\it R}-matrix approach in which Auger-plus-radiation damping has been included. The target and close-coupling expansions are both taken to be the 134 levels belonging to the configurations ${\rm 2s^22p^63}l$, ${\rm 2s^22p^53s3}l$, ${\rm 2s^22p^53p^2}$ and ${\rm 2s^22p^53p3d}$. The comparison of Maxwell-averaged effective collision strengths with and without damping shows that the damping reduction is about 30-40% for many transitions at low temperatures, but up to 80% for a few transitions. As a consequence, the results of previous Dirac $R$-matrix calculations (Aggarwal and Keenan, 2008) overestimate the effective collision strengths due to their omission of Auger-plus-radiation damping.
△ Less
Submitted 5 November, 2008;
originally announced November 2008.
-
Dielectronic recombination and stability of warm gas in AGN
Authors:
Susmita Chakravorty,
Ajit K. Kembhavi,
Martin Elvis,
Gary Ferland,
N. R. Badnell
Abstract:
Understanding the thermal equilibrium (stability) curve may offer insights into the nature of the warm absorbers often found in active galactic nuclei. Its shape is determined by factors like the spectrum of the ionizing continuum and the chemical composition of the gas. We find that the stability curves obtained under the same set of the above mentioned physical factors, but using recently deri…
▽ More
Understanding the thermal equilibrium (stability) curve may offer insights into the nature of the warm absorbers often found in active galactic nuclei. Its shape is determined by factors like the spectrum of the ionizing continuum and the chemical composition of the gas. We find that the stability curves obtained under the same set of the above mentioned physical factors, but using recently derived dielectronic recombination rates, give significantly different results, especially in the regions corresponding to warm absorbers, leading to different physical predictions. Using the current rates we find a larger probability of having thermally stable warm absorber at $10^5 \kel$ than previous predictions and also a greater possibility for its multiphase nature. the results obtained with the current dielectronic recombination rate coefficients are more reliable because the warm absorber models along the stability curve have computed coefficient values, whereas previous calculations relied on guessed averages for the same due to lack of available data.
△ Less
Submitted 13 November, 2008; v1 submitted 7 November, 2007;
originally announced November 2007.
-
Atomic data from the Iron Project. LXIV. Radiative transition rates and collision strengths for Ca II
Authors:
M. Melendez,
M. A. Bautista,
N. R. Badnell
Abstract:
This work reports radiative transition rates and electron impact excitation rate coefficients for levels of the n= 3, 4, 5, 6, 7, 8 configurations of Ca II. The radiative data were computed using the Thomas-Fermi-Dirac central potential method in the frozen core approximation and includes the polarization interaction between the valence electron and the core using a model potential. This method…
▽ More
This work reports radiative transition rates and electron impact excitation rate coefficients for levels of the n= 3, 4, 5, 6, 7, 8 configurations of Ca II. The radiative data were computed using the Thomas-Fermi-Dirac central potential method in the frozen core approximation and includes the polarization interaction between the valence electron and the core using a model potential. This method allows for configuration interactions (CI) and relativistic effects in the Breit-Pauli formalism. Collision strengths in LS-coupling were calculated in the close coupling approximation with the R-matrix method. Then, fine structure collision strengths were obtained by means of the intermediate-coupling frame transformation (ICFT) method which accounts for spin-orbit coupling effects. We present extensive comparisons with the most recent calculations and measurements for Ca II as well as a comparison between the core polarization results and the "unpolarized" values. We find that core polarization affects the computed lifetimes by up to 20%. Our results are in very close agreement with recent measurements for the lifetimes of metastable levels. The present collision strengths were integrated over a Maxwellian distribution of electron energies and the resulting effective collision strengths are given for a wide range of temperatures. Our effective collision strengths for the resonance transitions are within ~11% from previous values derived from experimental measurements, but disagree with latter computations using the distorted wave approximation.
△ Less
Submitted 28 April, 2007;
originally announced April 2007.
-
Dielectronic Recombination of Fe XV forming Fe XIV: Laboratory Measurements and Theoretical Calculations
Authors:
D. V. Lukić,
M. Schnell,
D. W. Savin,
C. Brandau,
E. W. Schmidt,
S. Böhm,
A. Müller,
S. Schippers,
M. Lestinsky,
F. Sprenger,
A. Wolf,
Z. Altun,
N. R. Badnell
Abstract:
We have measured resonance strengths and energies for dielectronic recombination (DR) of Mg-like Fe XV forming Al-like Fe XIV via N=3 -> N' = 3 core excitations in the electron-ion collision energy range 0-45 eV. All measurements were carried out using the heavy-ion Test Storage Ring at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. We have also carried out new multiconfigu…
▽ More
We have measured resonance strengths and energies for dielectronic recombination (DR) of Mg-like Fe XV forming Al-like Fe XIV via N=3 -> N' = 3 core excitations in the electron-ion collision energy range 0-45 eV. All measurements were carried out using the heavy-ion Test Storage Ring at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. We have also carried out new multiconfiguration Breit-Pauli (MCBP) calculations using the AUTOSTRUCTURE code. For electron-ion collision energies < 25 eV we find poor agreement between our experimental and theoretical resonance energies and strengths. From 25 to 42 eV we find good agreement between the two for resonance energies. But in this energy range the theoretical resonance strengths are ~ 31% larger than the experimental results. This is larger than our estimated total experimental uncertainty in this energy range of +/- 26% (at a 90% confidence level). Above 42 eV the difference in the shape between the calculated and measured 3s3p(^1P_1)nl DR series limit we attribute partly to the nl dependence of the detection probabilities of high Rydberg states in the experiment. We have used our measurements, supplemented by our AUTOSTRUCTURE calculations, to produce a Maxwellian-averaged 3 -> 3 DR rate coefficient for Fe XV forming Fe XIV. The resulting rate coefficient is estimated to be accurate to better than +/- 29% (at a 90% confidence level) for k_BT_e > 1 eV. At temperatures of k_BT_e ~ 2.5-15 eV, where Fe XV is predicted to form in photoionized plasmas, significant discrepancies are found between our experimentally-derived rate coefficient and previously published theoretical results. Our new MCBP plasma rate coefficient is 19-28% smaller than our experimental results over this temperature range.
△ Less
Submitted 6 April, 2007;
originally announced April 2007.
-
Dielectronic recombination of Fe 3p^q ions: a key ingredient for describing X-ray absorption in active galactic nuclei
Authors:
N. R. Badnell
Abstract:
We have carried-out multi-configuration Breit-Pauli AUTOSTRUCTURE calculations for the dielectronic recombination (DR) of Fe^{8+} - Fe^{12+} ions. We obtain total DR rate coefficients for the initial ground-level which are an order of magnitude larger than those corresponding to radiative recombination (RR), at temperatures where Fe 3p^q (q=2-6) ions are abundant in photoionized plasmas. The res…
▽ More
We have carried-out multi-configuration Breit-Pauli AUTOSTRUCTURE calculations for the dielectronic recombination (DR) of Fe^{8+} - Fe^{12+} ions. We obtain total DR rate coefficients for the initial ground-level which are an order of magnitude larger than those corresponding to radiative recombination (RR), at temperatures where Fe 3p^q (q=2-6) ions are abundant in photoionized plasmas. The resultant total (DR+RR) rate coefficients are then an order of magnitude larger than those currently in use by photoionized plasma modeling codes such as CLOUDY, ION and XSTAR. These rate coefficients, together with our previous results for q=0 and 1, are critical for determining the ionization balance of the M-shell Fe ions which give rise to the prominent unresolved-transition-array X-ray absorption feature found in the spectrum of many active galactic nuclei. This feature is poorly described by CLOUDY and ION, necessitating an ad hoc modification to the low-temperature DR rate coefficients. Such modifications are no longer necessary and a rigorous approach to such modeling can now take place using these data.
△ Less
Submitted 16 August, 2006;
originally announced August 2006.
-
Dielectronic recombination of Fe^{13+}: benchmarking the M-shell
Authors:
N. R. Badnell
Abstract:
We have carried-out a series of multi-configuration Breit-Pauli AUTOSTRUCTURE calculations for the dielectronic recombination of Fe^{13+}. We present a detailed comparison of the results with the high-energy resolution measurements reported recently from the Heidelberg Test Storage Ring by Schmidt et al. Many Rydberg series contribute significantly from this initial 3s^2 3p M-shell ion, resultin…
▽ More
We have carried-out a series of multi-configuration Breit-Pauli AUTOSTRUCTURE calculations for the dielectronic recombination of Fe^{13+}. We present a detailed comparison of the results with the high-energy resolution measurements reported recently from the Heidelberg Test Storage Ring by Schmidt et al. Many Rydberg series contribute significantly from this initial 3s^2 3p M-shell ion, resulting in a complex recombination `spectrum'. While there is much close agreement between theory and experiment, differences of typically 50% in the summed resonance strengths over 0.1-10 eV result in the experimentally based total Maxwellian recombination rate coefficient being a factor of 1.52-1.38 larger than theory over 10^4-10^5 K, which is a typical temperature range of peak abundance for Fe^{13+} in a photoionized plasma. Nevertheless, this theoretical recombination rate coefficient is an order of magnitude larger than that used by modellers to-date. This may help explain the discrepancy between the iron M-shell ionization balance predicted by photoionization modelling codes such as ION and CLOUDY and that deduced from the iron M-shell unresolved-transition-array absorption feature observed in the X-ray spectrum of many active galactic nuclei. Similar data are required for Fe^{8+} through Fe^{12+} to remove the question mark hanging over the atomic data though.
△ Less
Submitted 18 July, 2006;
originally announced July 2006.
-
Collisional Ionization Equilibrium for Optically Thin Plasmas. I. Updated Recombination Rate Coefficients for Bare though Sodium-like Ions
Authors:
P. Bryans,
N. R. Badnell,
T. W. Gorczyca,
J. M. Laming,
W. Mitthumsiri,
D. W. Savin
Abstract:
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 are often highly suspect. This translates directly into the reliability of the collisional ionization equilibrium (CIE) calculati…
▽ More
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 are often highly suspect. This translates directly into the reliability of the collisional ionization equilibrium (CIE) calculations. We make use of state-of-the-art calculations of dielectronic recombination (DR) rate coefficients for the hydrogenic through Na-like ions of all elements from He up to and including Zn. We also make use of state-of-the-art radiative recombination (RR) rate coefficient calculations for the bare through Na-like ions of all elements from H through to Zn. Here we present improved CIE calculations for temperatures from $10^4$ to $10^9$ K using our data and the recommended electron impact ionization data of \citet{Mazz98a} for elements up to and including Ni and Mazzotta (private communication) for Cu and Zn. DR and RR data for ionization stages that have not been updated are also taken from these two additional sources. We compare our calculated fractional ionic abundances using these data with those presented by Mazzotta et al. for all elements from H to Ni. The differences in peak fractional abundance are up to 60%. We also compare with the fractional ionic abundances for Mg, Si, S, Ar, Ca, Fe, and Ni derived from the modern DR calculations of \citet{Gu03a,Gu04a} for the H-like through Na-like ions, and the RR calculations of \citet{Gu03b} for the bare through F-like ions. These results are in better agreement with our work, with differences in peak fractional abundance of less than 10%.
△ Less
Submitted 1 August, 2006; v1 submitted 17 April, 2006;
originally announced April 2006.
-
Radiative recombination data for modelling dynamic finite-density plasmas
Authors:
N. R. Badnell
Abstract:
We have calculated partial final-state resolved radiative recombination (RR) rate coefficients from the initial ground and metastable levels of all elements up to and including Zn, plus Kr, Mo, and Xe, for all isoelectronic sequences up to Na-like forming Mg-like. The data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf48, which spans a temperature range of…
▽ More
We have calculated partial final-state resolved radiative recombination (RR) rate coefficients from the initial ground and metastable levels of all elements up to and including Zn, plus Kr, Mo, and Xe, for all isoelectronic sequences up to Na-like forming Mg-like. The data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf48, which spans a temperature range of z^2(10^1 - 10^7)K, where z is the initial ion charge. Fits to total rate coefficients have been determined, for both the ground and metastable levels, and those for the ground are presented here. Comparison is made both with previous RR rate coefficients and with (background) R-matrix photoionization cross sections. This RR database complements a DR database already produced and both are being used to produce updated ionization balances for both (electron) collisionally ionized and photoionized plasmas.
△ Less
Submitted 6 April, 2006;
originally announced April 2006.
-
Ionisation state, excited populations and emission of impurities in dynamic finite density plasmas. I: The generalised collisional-radiative model for light elements
Authors:
H. P. Summers,
W. J. Dickson,
M. G. O'Mullane,
N. R. Badnell,
A. D. Whiteford,
D. H. Brooks,
J. Lang,
S. D. Loch,
D. C. Griffin
Abstract:
The paper presents an integrated view of the population structure and its role in establishing the ionisation state of light elements in dynamic, finite density, laboratory and astrophysical plasmas. There are four main issues, the generalised collisional-radiative picture for metastables in dynamic plasmas with Maxwellian free electrons and its particularising to light elements, the methods of…
▽ More
The paper presents an integrated view of the population structure and its role in establishing the ionisation state of light elements in dynamic, finite density, laboratory and astrophysical plasmas. There are four main issues, the generalised collisional-radiative picture for metastables in dynamic plasmas with Maxwellian free electrons and its particularising to light elements, the methods of bundling and projection for manipulating the population equations, the systematic production/use of state selective fundamental collision data in the metastable resolved picture to all levels for collisonal-radiative modelling and the delivery of appropriate derived coefficients for experiment analysis. The ions of carbon, oxygen and neon are used in illustration. The practical implementation of the methods described here is part of the ADAS Project.
△ Less
Submitted 18 November, 2005;
originally announced November 2005.
-
Up-dated opacities from the Opacity Project
Authors:
N. R. Badnell,
M. A. Bautista,
K. Butler,
F. Delahaye,
C. Mendoza,
P. Palmeri,
C. J. Zeippen,
M. J. Seaton
Abstract:
Using the code AUTOSTRUCTURE, extensive calculations of inner-shell atomic data have been made for the chemical elements He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe and Ni. The results are used to obtain up-dated opacities from the Opacity Project, OP. A number of other improvements on earlier work have also been included.
Rosseland-mean opacities from OP are compared with those fro…
▽ More
Using the code AUTOSTRUCTURE, extensive calculations of inner-shell atomic data have been made for the chemical elements He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe and Ni. The results are used to obtain up-dated opacities from the Opacity Project, OP. A number of other improvements on earlier work have also been included.
Rosseland-mean opacities from OP are compared with those from OPAL. Differences of 5 to 10% occur. OP gives the `Z-bump', at log(T)=5.2, to be shifted to slightly higher temperatures. The opacities from OP, as functions of temperature and density, are smoother than those from OPAL. Extensive tests show that the numerical accuracy of the OP opacities following integration over frequency mesh and interpolation on temperature-density mesh is better than 1%.
Prior to a number of recent investigations which have indicated a need for a downward revision in the solar abundances of oxygen and other elements, there was good agreement between properties of the sun deduced from helioseismology and from stellar evolution models calculated using OPAL opacities. The revisions in abundances destroy that agreement. In a recent paper Bahcall et al argue that the agreement would be restored if opacities for the regions of the sun between 0.7 and 0.4 solar radii were larger than those given by OPAL by about 10%. In the region concerned, the present results from OP do not differ from those of OPAL by more than 2.5%.
△ Less
Submitted 18 December, 2004; v1 submitted 29 October, 2004;
originally announced October 2004.
-
A comparison of Rosseland-mean opacities from OP and OPAL
Authors:
M. J. Seaton,
N. R. Badnell
Abstract:
Monochromatic opacities from the Opacity Project (OP) (Seaton et al.) have been augmented by hitherto missing inner-shell contributions (Badnell & Seaton). OP Rosseland-mean opacities are compared with results from OPAL (Iglesias & Rogers) for the elements H, He, C, O, S and Fe. Overall, there is good agreement between OP and OPAL Rosseland-mean opacities for the 6-elements, but there are some d…
▽ More
Monochromatic opacities from the Opacity Project (OP) (Seaton et al.) have been augmented by hitherto missing inner-shell contributions (Badnell & Seaton). OP Rosseland-mean opacities are compared with results from OPAL (Iglesias & Rogers) for the elements H, He, C, O, S and Fe. Overall, there is good agreement between OP and OPAL Rosseland-mean opacities for the 6-elements, but there are some differences. In particular, recent work (Bahcall et al.) has shown that helioseismology measurements give a very accurate value for the depth of the solar convection zone and that solar models give agreement with that value only if opacities there are about 7% larger than OPAL values. For the 6-element mix, here, we obtain Rosseland-mean opacities to be larger than OPAL by 5%.
△ Less
Submitted 22 April, 2004;
originally announced April 2004.
-
On the importance of inner-shell transitions for opacity calculations
Authors:
N. R. Badnell,
M. J. Seaton
Abstract:
For high temperatures and densities, stellar opacities obtained from the Opacity Project (OP) were smaller than those obtained from the OPAL project. Iglesias and Rogers [Ap.J. v443, 469 (1995)] suggested that the discrepancy was due to the omission by OP of important atomic inner-shell processes, and considered in detail results for a mixture of 6 elements: H, He, C, O, S and Fe. Extensive new…
▽ More
For high temperatures and densities, stellar opacities obtained from the Opacity Project (OP) were smaller than those obtained from the OPAL project. Iglesias and Rogers [Ap.J. v443, 469 (1995)] suggested that the discrepancy was due to the omission by OP of important atomic inner-shell processes, and considered in detail results for a mixture of 6 elements: H, He, C, O, S and Fe. Extensive new inner-shell data have now been computed using the code AUTOSTRUCTURE. It is shown that the inclusion of these data in the OP work gives opacities for the 6-element mix which are in much closer agreement with those from OPAL. We also discuss a number of problems relating to the calculation of opacities and of equations-of-state in dense plasmas.
△ Less
Submitted 22 August, 2003;
originally announced August 2003.
-
Dielectronic recombination data for dynamic finite-density plasmas I. Goals and methodology
Authors:
N. R. Badnell,
M. G. O'Mullane,
H. P. Summers,
Z. Altun,
M. A. Bautista,
J. Colgan,
T. W. Gorczyca,
D. M. Mitnik,
M. S. Pindzola,
O. Zatsarinny
Abstract:
A programme is outlined for the assembly of a comprehensive dielectronic recombination database within the generalized collisional--radiative (GCR) framework. It is valid for modelling ions of elements in dynamic finite-density plasmas such as occur in transient astrophysical plasmas such as solar flares and in the divertors and high transport regions of magnetic fusion devices. The resolution a…
▽ More
A programme is outlined for the assembly of a comprehensive dielectronic recombination database within the generalized collisional--radiative (GCR) framework. It is valid for modelling ions of elements in dynamic finite-density plasmas such as occur in transient astrophysical plasmas such as solar flares and in the divertors and high transport regions of magnetic fusion devices. The resolution and precision of the data are tuned to spectral analysis and so are sufficient for prediction of the dielectronic recombination contributions to individual spectral line emissivities. The fundamental data are structured according to the format prescriptions of the Atomic Data and Analysis Structure (ADAS) and the production of relevant GCR derived data for application is described and implemented following ADAS. The requirements on the dielectronic recombination database are reviewed and the new data are placed in context and evaluated with respect to older and more approximate treatments. Illustrative results validate the new high-resolution zero-density dielectronic recombination data in comparison with measurements made in heavy-ion storage rings utilizing an electron cooler. We also exemplify the role of the dielectronic data on GCR coefficient behaviour for some representative light and medium weight elements.
△ Less
Submitted 15 April, 2003;
originally announced April 2003.
-
Assessment of the Fluorescence and Auger Data Base used in Plasma Modeling
Authors:
T. W. Gorczyca,
C. N. Kodituwakku,
K. T. Korista,
O. Zatsarinny,
N. R. Badnell,
E. Behar,
M. H. Chen,
D. W. Savin
Abstract:
We have investigated the accuracy of the 1s-vacancy fluorescence data base of Kaastra & Mewe (1993, A&AS, 97, 443) resulting from the initial atomic physics calculations and the subsequent scaling along isoelectronic sequences. In particular, we have focused on the relatively simple Be-like and F-like 1s-vacancy sequences. We find that the earlier atomic physics calculations for the oscillator s…
▽ More
We have investigated the accuracy of the 1s-vacancy fluorescence data base of Kaastra & Mewe (1993, A&AS, 97, 443) resulting from the initial atomic physics calculations and the subsequent scaling along isoelectronic sequences. In particular, we have focused on the relatively simple Be-like and F-like 1s-vacancy sequences. We find that the earlier atomic physics calculations for the oscillator strengths and autoionization rates of singly-charged B II and Ne II are in sufficient agreement with our present calculations. However, the substantial charge dependence of these quantities along each isoelectronic sequence, the incorrect configuration averaging used for B II, and the neglect of spin-orbit effects (which become important at high-Z) all cast doubt on the reliability of the Kaastra & Mewe data for application to plasma modeling.
△ Less
Submitted 10 April, 2003;
originally announced April 2003.