BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

N. R. Badnell

N. R. Badnell appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.atom-phastro-ph.SRastro-ph.GAastro-ph.HEastro-ph.IMphysics.plasm-phastro-phastro-ph.COphysics.comp-phphysics.space-ph

Catalog footprint

What is connected

29works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2021arXiv

The influence of photo-induced processes and charge transfer on carbon and oxygen in the lower solar atmosphere

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of density on the oxygen ionisation equilibrium in collisional plasmas

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fe III emission in quasars: evidence for a dense turbulent medium

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Helium line emissivities in the solar corona

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PSTGF : time-independent R-Matrix atomic electron-impact code

STGF is a community code employed for outer-region R-matrix calculations, describing electron-impact collisional processes. It is widely recognised that the original version of STGF was written by M. J. Seaton in 1983, but through constant refinement over the next decades by worldwide contributors has evolved into its current form that more reflects modern coding practice and current computer architectures. Despite its current wide acceptance, it was never formally published. Therefore, we present an updated high-performance parallel version of PSTGF, that balances the requirements of small university clusters, yet can exploit the computational power of cutting edge supercomputers. There are many improvements over the original STGF, but most noticeably, the full introduction of MQDT options that provide subsequent integration with ICFT (Intermediate Coupling Frame Transformation) codes, and for either Breit-Pauli / DARC (Dirac Atomic R-matrix Codes), better load balancing, high levels of vectorisation and simplified output. Semantically, the program is full fortran 90 in conjunction with MPI (Message Passing Interface) though has CUDA fortran options for the most numerically intensive code sections.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

R-matrix electron-impact excitation data for the C-like iso-electronic sequence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

H, He-like recombination spectra I: $l$-changing collisions for hydrogen

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

H, He-like recombination spectra II: $l$-changing collisions for He Rydberg states

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hot DA white dwarf model atmosphere calculations: Including improved Ni PI cross sections

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Partial and Total Dielectronic Recombination rate coefficients for W73+ to W56+

Dielectronic recombination (DR) is a key atomic process which affects the spectroscopic diagnostic modelling of tungsten, most of whose ionization stages will be found somewhere in the ITER fusion reactor: in the edge, divertor, or core plasma. Accurate DR data is sparse while complete DR coverage is unsophisticated (e.g. average-atom or Burgess General Formula) as illustrated by the large uncertainties which currently exist in the tungsten ionization balance. To this end, we present a series of partial final state-resolved and total DR rate coefficients for W73+ to W56+ Tungsten ions. This is part of a wider effort within The Tungsten Project to calculate accurate dielectronic recombination rate coefficients for the tungsten isonuclear sequence for use in collisional-radiative modelling of finite-density tokamak plasmas. The recombination rate coefficients have been calculated with autostructure using kappa-averaged relativistic wavefunctions in level resolution (intermediate coupling) and configuration resolution (configuration average). The results are available from OPEN-ADAS according to the adf09 and adf48 standard formats. Comparison with previous calculations of total DR rate coefficients for W63+ and W56+ yield agreement to within 20% and 10%, respectively, at peak temperature. It is also seen that the Juttner correction to the Maxwell distribution has a significant effect on the ionization balance of tungsten at the highest charge states, changing both the peak abundance temperatures and the ionization fractions of several ions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scaling of collision strengths for highly-excited states of ions of the H- and He-like sequences

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bootstrapping dielectronic recombination from second-row elements and the Orion Nebula

Dielectronic recombination (DR) is the dominant recombination process for most heavy elements in photoionized clouds. Accurate DR rates for a species can be predicted when the positions of autoionizing states are known. Unfortunately such data are not available for most third and higher-row elements. This introduces an uncertainty that is especially acute for photoionized clouds, where the low temperatures mean that DR occurs energetically through very low-lying autoionizing states. This paper discusses S$^{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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dirac $R$-matrix and Breit-Pauli distorted wave calculations of the electron-impact excitation of W$^{44+}$

With construction of ITER progressing and existing tokamaks carrying out ITER-relevant experiments, accurate fundamental and derived atomic data for numerous ionization stages of tungsten (W) is required to assess the potential effect of this species upon fusion plasmas. The results of fully relativistic, partially radiation damped, Dirac $R$-matrix electron-impact excitation calculations for the W$^{44+}$ ion are presented. These calculations use a configuration interaction and close-coupling expansion that opens-up the 3d-subshell, which does not appear to have been considered before in a collision calculation. As a result, it is possible to investigate the arrays, [3d$^{10}$4s$^2-$3d$^9$4s$^2$4f] and [3d$^{10}$4s$^2-$3d$^9$4s4p4d], which are predicted to contain transitions of diagnostic importance for the soft x-ray region. Our $R$-matrix collision data are compared with previous $R$-matrix results by Ballance and Griffin as well as our own relativistically corrected, Breit-Pauli distorted wave and plane-wave Born calculations. All relevant data are applied to the collisional-radiative modelling of atomic populations, for further comparison. This reveals the paramount nature of the 3d-subshell transitions from the perspectives of radiated power loss and detailed spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Validity of the ICFT R-matrix method: Be-like Al$^{9+}$ a case study

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A possible solution to the Lyman/Balmer line problem in hot DA white dwarfs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

Modeling the spectral emission of low-charge iron group ions enables the diagnostic determination of the local physical conditions of many cool plasma environments such as those found in H II regions, planetary nebulae, active galactic nuclei etc. Electron-impact excitation drives the population of the emitting levels and, hence, their emissivities. By carrying-out Breit-Pauli and intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of Fe$^{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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recombination of W18+ ions with electrons: Absolute rate coefficients from a storage-ring experiment and from theoretical calculations

We present new experimentally measured and theoretically calculated rate coefficients for the electron-ion recombination of W$^{18+}$([Kr] $4d^{10}$ $4f^{10}$) forming W$^{17+}$. At low electron-ion collision energies, the merged-beam rate coefficient is dominated by strong, mutually overlapping, recombination resonances. In the temperature range where the fractional abundance of W$^{18+}$ is expected to peak in a fusion plasma, the experimentally derived Maxwellian recombination rate coefficient is 5 to 10 times larger than that which is currently recommended for plasma modeling. The complexity of the atomic structure of the open-$4f$-system under study makes the theoretical calculations extremely demanding. Nevertheless, the results of new Breit-Wigner partitioned dielectronic recombination calculations agree reasonably well with the experimental findings. This also gives confidence in the ability of the theory to generate sufficiently accurate atomic data for the plasma modeling of other complex ions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Collision Strengths for [O III] Optical and Infrared Lines

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Opacities and Spectra of the r-process Ejecta from Neutron Star Mergers

Material ejected during (or immediately following) the merger of two neutron stars may assemble into heavy elements 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Suppression of Dielectronic Recombination Due to Finite Density Effects

We have developed a general model for determining density-dependent effective dielectronic recombination (DR) rate coefficients in order to explore finite-density effects on the ionization balance of plasmas. Our model consists of multiplying by a suppression factor those highly-accurate total zero-density DR rate coefficients which have been produced from state-of-the-art theoretical calculations and which have been benchmarked by experiment. The suppression factor is based-upon earlier detailed collision-radiative calculations which were made for a wide range of ions at various densities and temperatures, but used a simplified treatment of DR. A general suppression formula is then developed as a function of isoelectronic sequence, charge, density, and temperature. These density-dependent effective DR rate coefficients are then used in the plasma simulation code Cloudy to compute ionization balance curves for both collisionally ionized and photoionized plasmas at very low (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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dielectronic recombination of W^20+ (4d^10 4f^8): addressing the half-open f-shell

A recent measurement of the dielectronic recombination (DR) of W^20+ [Schippers et al Phys. Rev. A83, 012711 (2011)] found an exceptionally large contribution from near threshold resonances (<1eV). This still affected the Maxwellian rate coefficient at much higher temperatures. The experimental result was found to be a factor 4 or more than that currently in use in the 100-300eV range which is of relevance for modeling magnetic fusion plasmas. We have carried-out DR calculations with AUTOSTRUCTURE which include all significant single electron promotions. Our intermediate coupling (IC) results are more than a factor of 4 larger than our LS-coupling ones at 1eV but still lie a factor 3 below experiment here. If we assume complete (chaotic) mixing of near-threshold autoionizing states then our results come into agreement (to within 20%) with experiment below about 2eV. Our total IC Maxwellian rate coefficients are 50-30% smaller than those based-on experiment over 100-300eV.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electron-ion Recombination of Fe XII forming Fe XI: Laboratory Measurements and Theoretical Calculations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

R-matrix electron-impact excitation data for astrophysically abundant sulphur ions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dielectronic Recombination of Argon-Like Ions

We present a theoretical investigation of dielectronic recombination (DR) of Ar-like ions that sheds new light on the behavior of the rate coefficient at low-temperatures where these ions form in photoionized plasmas. We provide results for the total and partial Maxwellian-averaged DR rate coefficients from the initial ground level of K II -- Zn XIII ions. It is expected that these new results will advance the accuracy of the ionization balance for Ar-like M-shell ions and pave the way towards a detailed modeling of astrophysically relevant X-ray absorption features. We utilize the AUTOSTRUCTURE computer code to obtain the accurate core-excitation thresholds in target ions and carry out multiconfiguration Breit-Pauli (MCBP) calculations of the DR cross section in the independent-processes, isolated-resonance, distorted-wave (IPIRDW) approximation. Our results mediate the complete absence of direct DR calculations for certain Ar-like ions and question the reliability of the existing empirical rate formulas, often inferred from renormalized data within this isoelectronic sequence.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Revised Results for Non-thermal Recombination Flare Hard X-Ray Emission

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Electron-ion Recombination of Fe X forming Fe IX and of Fe XI forming Fe X: Laboratory Measurements and Theoretical Calculations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Radiative transition rates and collision strengths for Si II

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Dielectronic recombination of Fe 3p^q ions: a key ingredient for describing X-ray absorption in active galactic nuclei

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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Ionisation state, excited populations and emission of impurities in dynamic finite density plasmas. I: The generalised collisional-radiative model for light elements

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.

0 citations0 reviewsNo code linkedOpen access