the broad mechanism for radiation emission from a hot tenuous plasma is simple .
thermal kinetic energy of free electrons in the plasma is transferred by collisions to the internal energy of impurity ions , @xmath0 where @xmath1 denotes an excited state and @xmath2 the ground state of the impurity ion .
this energy is then radiated as spectrum line photons which escape from the plasma volume @xmath3 where @xmath4 is the emitted photon energy and @xmath5 its frequency .
similarly ions in general increase or decrease their charge state by collisions with electrons @xmath6 where @xmath7 denotes the next ionisation stage of impurity ion @xmath2 .
the situation is often referred to as the _ coronal picture_. the coronal picture has been the basis for the description of impurities in fusion plasmas for many years .
however , the progress towards ignition of fusion plasmas and to higher density plasmas requires a description beyond the coronal approximation . models of finite density plasmas which include some parts of the competition between radiative and collisional processes are loosely called _ collisional - radiative_. however , collisional - radiative theory in its origins ( bates , 1962 ) was designed for the description of dynamic plasmas and this aspect is essential for the present situations of divertors , heavy species , transport barriers and transient events . the present work is centred on _ generalised collisional radiative _ ( @xmath8 ) theory ( mcwhirter and summers , 1984 ) which is developed in the following sections .
it is shown that consideration of relaxation time - scales , metastable states , secondary collisions etc . -
aspects rigorously specified in collisional - radiative theory - allow an atomic description suitable for modelling the newer areas above .
the detailed quantitative description is complicated because of the need to evaluate individually the many controlling collisional and radiative processes , a task which is compounded by the variety of atoms and ions which participate .
the focus is restricted to plasmas which are optically thin and not influenced by external radiation fields , and for which ground and metastable populations of ions dominate other excited ion populations .
the paper provides an overview of key methods used to expedite this for light elements and draws illustrative results from the ions of carbon , oxygen and neon .
the paper is intended as the first of a series of papers on the application of collisional - radiative modelling in more advanced plasma scenarios and to specific important species .
the practical implementation of the methods described here is part of the adas ( atomic data and analysis structure ) project ( summers , 1993 , 2004 ) .
illustrations are drawn from adas codes and the adas fundamental and derived databases .
the lifetimes of the various states of atoms , ions and electrons in a plasma to radiative or collisional processes vary enormously .
of particular concern for spectroscopic studies of dynamic finite density plasmas are those of translational states of free electrons , atoms and ions and internal excited states ( including states of ionisation ) of atoms and ions .
these lifetimes determine the relaxation times of the various populations , the rank order of which , together with their values relative to observation times and plasma development times determines the modelling approach .
the key lifetimes divide into two groups .
the first is the _ intrinsic _ group , comprising purely atomic parameters , and includes metastable radiative decay , @xmath9 , ordinary excited state radiative decay @xmath10 and auto - ionising state decay ( radiative and auger ) , @xmath11 .
the intrinsic group for a particular ion is generally ordered as @xmath12 with typical values @xmath13 where @xmath14 is the ion charge .
the second is the _ extrinsic _ group , which depends on plasma conditions - especially particle density .
it includes free particle thermalisation ( including electron - electron @xmath15 , ion - ion @xmath16 and ion - electron @xmath17 ) , charge - state change ( ionisation @xmath18 and recombination @xmath19 ) and redistribution of population amongst excited ion states ( @xmath20 ) .
the extrinsic group is ordered in general as @xmath21 with approximate expressions for the time constants given by @xmath22(1/(z+1)^2)(kt_e / i_h)^{1/2}(\rm{cm}^{-3}/n_e)~~s    \nonumber\\   \tau_{ion } & \sim & [ 10 ^ 5 - 10 ^ 7](z+1)^4(i_h / kt_e)^{1/2}e^{\chi / kt_e}(\rm{cm}^{-3}/n_e)~~s     \nonumber\\   \tau_{i - i } & \sim & [ 7.0\times10 ^ 7](m_i / m_p)^{1/2 }   \nonumber\\ & & ( kt_e / i_h)^{3/2}(1/z^4)(\rm{cm}^{-3}/n_i)~~s      \\   \tau_{i - e } & \sim & [ 1.4\times10 ^ 9](m_i / m_p)^{1/2 }       \nonumber\\ & & ( ( kt_e / i_h)+5.4\times10^{-4}(kt_i / i_h)(m_p / m_i))^{3/2 }   \nonumber\\ & & ( 1/z^2)(\rm{cm}^{-3}/n_i)~~s                     \nonumber\\   \tau_{e - e } & \sim & [ 1.6\times10 ^ 6](kt_e / i_h)^{3/2}(\rm{cm}^{-3}/n_i)~~s .
\nonumber \end{aligned}\ ] ] the ion mass is @xmath23 , the proton mass @xmath24 , the ionisation potential @xmath25 , the ion density @xmath26 , the electron density @xmath27 , the ion temperature @xmath28 , the electron temperature @xmath29 and the ionisation energy of hydrogen is @xmath30 . @xmath20 may span across the inequalities of equation [ eqn : eqn6 ] and is discussed in a later paragraph .    from a dynamic point of view , the intrinsic and extrinsic groups are to be compared with each other and with timescales , @xmath31 , representing plasma ion diffusion across temperature or density scale lengths , relaxation times of transient phenomena and observation times . for most plasmas in magnetic confinement fusion and astrophysics @xmath32 where @xmath33 represents the relaxation time of ground state populations of ions ( a composite of @xmath19 and @xmath18 ) and
it is such plasmas which are addressed in this paper .
these time - scales imply that the dominant populations of impurities in the plasma are those of the ground and metastable states of the various ions .
the dominant populations evolve on time - scales of the order of plasma diffusion time - scales and so should be modelled dynamically , that is in the time - dependent , spatially varying , particle number continuity equations , along with the momentum and energy equations of plasma transport theory .
illustrative results are shown in figure [ fig : fig1]a .
the excited populations of impurities and the free electrons on the other hand may be assumed relaxed with respect to the instantaneous dominant populations , that is they are in a _ quasi - equilibrium_. the quasi - equilibrium is determined by local conditions of electron temperature and electron density .
so , the atomic modelling may be partially de - coupled from the impurity transport problem into local calculations which provide quasi - equilibrium excited ion populations and emissivities and then effective source coefficients ( collisional - radiative coefficients ) for dominant populations which must be entered into the plasma transport equations .
the solution of the transport equations establishes the spatial and temporal behaviour of the dominant populations which may then be re - associated with the local emissivity calculations for matching to and analysis of observations .    for excited populations
, @xmath20 plays a special and complicated role due to the very large variation in collisional excitation / de - excitation reaction rates with the quantum numbers of the participating states . in the low density coronal picture @xmath34 and redistribution
plays no part .
critical densities occur for @xmath35 and for @xmath36 and allow division of the ( in principle ) infinite number of excited populations into categories including _ low levels _ , _ high singly excited levels _ and _ doubly excited levels _ for which important simplifications are possible .
these are examined in section [ sec : sec2 ] .
light element ions in fusion plasmas are generally in the singly excited state redistibutive case , approaching the doubly excited redistribution case at the higher densities .
highly ionised ions of heavy species in fusion plasmas approach the coronal picture .
illustrative results on critical densities are shown in figure [ fig : fig1]b .    finally , because of the generally short @xmath15 compared with other timescales ( including those of free - free and free - bound emission ) , it is usually the case that the free electrons have close to a maxwellian distribution .
this assumption is made throughout the present paper , but is relaxed in the next paper of the series ( bryans , 2005 ) .
the basic model was established by bates ( 1962 ) .
the ion in a plasma is viewed as composed of a complete set of levels indexed by @xmath37 and @xmath38 and a set of radiative and collisional couplings between them denoted by @xmath39 ( an element of the _ collisional - radiative matrix _ representing transition from @xmath38 to @xmath37 ) to which are added direct ionisations from each level of the ion to the next ionisation stage ( coefficient @xmath40 ) and direct recombinations to each level of the ion from the next ionisation stage ( coefficient @xmath41 ) .
thus , for each level , there is a total loss rate coefficient for its population number density , @xmath26 , given by @xmath42 following the discussion in the introduction , it is noted that populated metastable states can exist and there is no real distinction between them and ground states .
we use the term _ metastables _ to denote both ground and metastables states .
metastables are the dominant populations and so only recombination events which start with a metastable as a collision partner matter .
we condider the population structure of the @xmath14-times ionised ion , called the recombined or child ion .
the @xmath43-times ionised ion is called the recombining or parent ion and the @xmath44-times ionised ion is called the grandchild .
the metastables of the recombined ion are indexed by @xmath45 and @xmath46 , those of the recombining ion by @xmath47 and @xmath48 and those of the grandchild by @xmath49 and @xmath50 .
therefore the ion of charge state @xmath14 has metastable populations @xmath51 , the recombining ion of charge @xmath43 has metastable populations @xmath52 and the grandchild ion of charge @xmath44 has metastable populations @xmath53 . we designate the remaining excited states of the @xmath14-times ionised ion , with the metastables separated , as _ ordinary _ levels for which we reserve the indices @xmath37 and @xmath38 and populations @xmath26 and @xmath54 .
there are then , for example , direct recombination coefficients @xmath55 from each parent metastable into each child ordinary level and direct ionisation coefficients from each child ordinary level to each parent metastable @xmath56 such that @xmath57 .
also there are direct ionisation coefficients @xmath58to the metastables of the child from the metastables of the grandchild .
then the continuity equations for population number densities are @xmath59       =      \left[\begin{array}{llll }          \scriptnew{c}_{\mu \mu ' } & n_e~\scriptnew{r}_{\mu,\sigma } & 0 & 0                \\
n_e\scriptnew{s}_{\rho \mu ' } & c_{\rho \sigma } & c_{\rho j } & n_e~r_{\rho \nu ' } \\          0   & c_{i \sigma } & c_{i j } & n_e~r_{i \nu ' }         \\          0 & n_es_{\nu \sigma } & n_es_{\nu j } & \scriptnew{c}_{\nu \nu ' }             \end{array }       \right ]      \left[\begin{array}{l }          n^-_{\mu ' }           \\
n_{\sigma }           \\
n_{j }            \\          n^+_{\nu ' }               \end{array }       \right ]         \nonumber \\\end{aligned}\ ] ] where the equations for the @xmath44-times and @xmath43-times ionised ions have been simplified by incorporating their ordinary population contributions in their metastable contributions ( shown as script capital symbols ) as the immediate focus is on the @xmath14-times ionised ion .
this incorporation procedure is shown explicitly in the following equations for the @xmath14-times ionised ion through to equations [ eqn : eqn16 ] and may be done for each ionisation stage separately . note additionally the assumption ( made by omission of the ( 3,1 ) partition element , where 3 denotes the row and 1 the column ) that state - selective ionisation from the stage @xmath44 takes place only into the metastable manifold of the stage @xmath14 .      from the quasi - static assumption , we set @xmath60 and then the matrix equation for the ordinary levels of the @xmath14-times ionised ion gives @xmath61 where we have used summation convention on repeated indices .
substitution in equations [ eqn : eqn10 ] for the metastables of the @xmath14-times ionised ion gives @xmath62n^-_{\mu ' }     \nonumber \\   & & + \left [ c_{\rho \sigma}-c_{\rho j}c^{-1}_{ji}c_{i \sigma}\right ] n_{\sigma } \nonumber \\   & & + n_e\left [ r_{\rho \nu'}- c_{\rho j}c^{-1}_{ji}r_{i \nu'}\right]n^+_{\nu'}.    \end{aligned}\ ] ] the left - hand - side is interpreted as a total derivative with time - dependent and convective parts and the right - hand - side comprises the source terms .
the terms in square brackets in equations [ eqn : eqn12 ] give the effective growth rates of each metastable population of the @xmath14-times ionised ion driven by excitation ( or de - excitation ) from other metastables of the @xmath14-times ionised ion , by ionisation to the @xmath43-times ionised ion and excitation to other metastables of the @xmath14-times ionised ion ( a negatively signed growth ) and by recombination from the metastables of the @xmath43-times ionised ion .
these are called the @xmath8 coefficients . following burgess and summers ( 1969 ) , who used the name ` collisional - dielectronic ' for ` collisional - radiative ' when dielectronic recombination is active , we use the nomenclature @xmath63 for the @xmath8 recombination coefficients which become @xmath64 the @xmath8 metastable cross - coupling coefficients ( for @xmath65 ) are @xmath66 /n_e.\ ] ] note that the on - diagonal element @xmath67 /n_e$ ] with @xmath68 is a total loss rate coefficient from the metastable @xmath45 .
substitution from equation [ eqn : eqn11 ] in the equations [ eqn : eqn10 ] for the metastables of the @xmath69-times ionised ion gives @xmath70 n_{\sigma } \nonumber \\   & & + \left [ c_{\nu \nu'}- n_e^2s_{\nu j}c^{-1}_{ji}r_{i \nu'}\right]n^+_{\nu'}.     \end{aligned}\ ] ] the @xmath8 ionisation coefficients resolved by initial and final metastable state are @xmath71\ ] ] and note that there is contribution to cross - coupling between parents via recombination to excited states of the @xmath14-times ionised ion followed by re - ionisation to a different metastable , @xmath72.\ ] ] consider the sub - matrix comprising the ( 2,2 ) , ( 2,3 ) , ( 3,2 ) and ( 3,3 ) partitions of equations [ eqn : eqn10 ] . introduce the inverse of this sub - matrix as @xmath73       =      \left[\begin{array}{ll }              c_{\rho \sigma } & c_{\rho j }                   \\
c_{i \sigma }     & c_{i j }                      \\
\end{array }       \right]^{-1}\end{aligned}\ ] ] and note that the inverse of the ( 1,1 ) partition @xmath74^{-1 } \equiv \scriptnew{c}_{\rho \sigma}= \left [ c_{\rho \sigma}-c_{\rho j}c^{-1}_{ji}c_{i \sigma}\right].\ ] ] this compact representation illustrates , that the imposition of the quasi - static assumption leading to elimination of the ordinary level populations in favour of the metastable populations , may be viewed as a _ condensation _ in which the influence of the ordinary levels is _ projected _ onto the metastable levels .
the metastables can be condensed in a similar manner onto the ground restoring the original ( ground states only ) collisional - radiative picture .
the additive character of the direct metastable couplings @xmath75 means that these elements may be adjusted retrospectively after the main condensations .
there are two kinds of derived coefficients associated with individual spectrum line emission in common use in fusion plasma diagnosis
. these are _ photon emissivity coefficients _ ( @xmath76 ) and _ ionisation per photon ratio _ ( @xmath77 ) .
the reciprocals of the latter are also known as _
photon efficiencies_. from equations [ eqn : eqn11 ] , the emissivity in the spectrum line @xmath78 may be written as @xmath79 this allows specification of the _ excitation _ photon emissivity coefficient @xmath80 and the _ recombination _ photon emissivity coefficient @xmath81 the ionisation per photon ratios are most meaningful for the excitation part of the emissivity and are @xmath82 each of these coefficients is associated with a particular metastable @xmath46 , @xmath83 or @xmath84 of the @xmath85 , @xmath86 or @xmath87 ions respectively .
the radiated power in a similar manner separates into parts driven by excitation and by recombination as @xmath88 called the _ low - level line power coefficient _ and @xmath89 called the _ recombination - bremsstrahlung power coefficient _ where it is convenient to include bremsstrahlung with @xmath90 . note that in the generalised picture , additional power for the @xmath14-times ionised ions occurs in forbidden transitions between metastables as @xmath91 for the @xmath14-times ionised ion . in the fusion context , this is usually small .
radiated power is the most relevant quantity for experimental detection . for modelling , it is the _ electron energy loss function _ which enters the fluid energy equation .
the contribution to the total electron energy loss rate for the @xmath14-times ionised ion associated with ionisation and recombination from the @xmath43-times ionised ion is @xmath92 where the @xmath93 and @xmath94 are absolute energies of the metastables @xmath45 of the @xmath14-times ionised ion and of the metastables @xmath47 of the @xmath43-times ionised ion respectively .
thus the electron energy loss is a derived quantity from the radiative power coefficients and the other generalised collisonal - radiative coefficients .
note that cancellations in the summations cause the reduction to relative energies and that in ionisation equilibrium , the electron energy loss equals the radiative power loss .    in the following sections
, it is shown how the quasi - static assumption and population categorisation allow us to solve the infinite level population structure of each ion in a manageable and efficient way .
such solution is necessary for low and medium density astrophysical and magnetic confinement fusion plasmas .
this is unlike the simpler situation of very dense plasmas where heavy level truncation is used , because of continuum merging .
the handling of metastables in a generalised collisonal - radiative framework requires a detailed specific classification of level structure compatible with both recombining and recombined ions . for light element ions , russell - saunders
( l - s ) coupling is appropriate and it is sufficient to consider only terms , since although fine structure energy separations may be required for high resolution spectroscopy , relative populations of levels of a term are close to statistical .
so the parent ion metastables are of the form @xmath95 with @xmath96 the configuration and the recombined ion metastables are of the form @xmath97 with the configuration @xmath98 and the excited ( including highly excited ) terms are @xmath99 .
@xmath100 and @xmath101 denote individual electron principal quantum number and orbital angular momentum , @xmath102 and @xmath103 denote total orbital and total spin angular momenta of the multi - electron ion respectively in the specification of an ion state in russell - saunders coupling .
the configuration specifies the orbital occupancies of the ion state . for ions of heavy elements , relative populations of fine structure levels can differ markedly from statistical and it is necessary to work in intermediate coupling with parent metastables of the form @xmath104 and recombined metastables of the form @xmath105 with the configuration @xmath98 and the excited ( including highly excited ) levels @xmath106 .
there is a problem . to cope with the very many principal quantum shells participating in the calculations of collisional - dielectronic coefficients at finite density
necessitates a grosser viewpoint ( in which populations are _ bundled _ ) , whereas for modelling detailed spectral line emission , the finer viewpoint ( in which populations are fully _ resolved _ ) is required . in practice
, each ion tends to have a limited set of low levels principally responsible for the dominant spectrum line power emission for which a bundled approach is too imprecise , that is , averaged energies , oscillator strengths and collision strengths do not provide a good representation . note also that key parent transitions for dielectronic recombination span a few ( generally the same ) low levels for which precise atomic data are necessary . in the recombined ion ,
parentage gives approximate quantum numbers , that is , levels of the same @xmath100 ( and @xmath101 ) divide into those based on different parents .
lifetimes of levels of the same @xmath100 but different parents can vary strongly ( for example through secondary autoionisation ) .
also the recombination population of such levels is generally from the parent with which they are classified .       we therefore recognise three sets of non - exclusive levels of the recombined ion :    * metastable levels - indexed by @xmath45 , @xmath46 .
* low levels - indexed by @xmath37 , @xmath38 in a resolved coupling scheme , being the complete set of levels of a principal quantum shell range @xmath107 , including relevant metastables and spanning transitions contributing substantially to radiative power or of interest for specific observations . *
bundled levels - segregated according to the parent metastable upon which they are built and possibly also by spin system - which can include _ bundle - nl _ and _ bundle - n_.    viewed as a recombining ion , the set ( i ) must include relevant parents and set ( ii ) must span transitions which are dielectronic parent transitions .
time dependence matters only for the populations of ( i ) , high precision matters only for groups ( i ) and ( ii ) and special very many level handling techniques matter only for group ( iii ) .    to satisfy the various requirements and to allow linking of population sets at different resolutions , a series of manipulations on the collisional - radiative matrices
are performed ( summers and hooper , 1983 ) . to illustrate this , suppose there is a single parent metastable state .
consider the collisional - radiative matrix for the recombined ion and the right hand side in the bundle - n picture , and a partition of the populations as @xmath108 $ ] with @xmath107 and @xmath109 .
elimination of the @xmath110 yields a set of equations for the @xmath111 .
we call this a ` condensation ' of the whole set of populations onto the @xmath100 populations .
the coefficients are the effective ionisation coefficients from the @xmath100 , the effective cross - coupling coefficients between the @xmath100 and the effective recombination coefficients into the @xmath100 , which now include direct parts and indirect parts through the levels @xmath112 .
exclusion of the direct terms prior to the manipulations yields only the the indirect parts . call these @xmath113 and @xmath114 .
we make the assumption that @xmath113 and @xmath114 may be expanded over the resolved low level set ( see section [ sec : sec2.2 ] ) to give the expanded indirect matrix @xmath115 and @xmath116 where @xmath37 and @xmath38 span the resolved low level set ( ij ) .
these indirect couplings are then combined with higher precision direct couplings @xmath117 and @xmath118 so that    @xmath119    and @xmath120 the procedure is shown schematically in figure [ fig : fig2 ] .
the process may be continued , condensing the low level set onto the metastable set .
the generalised collisional - dielectronic coefficients are the result .
the time dependent and/or spatial non - equilibrium transport equations which describe the evolution of the ground and metastable populations of ions in a plasma use these generalised coefficients . following solution
, the condensations can be reversed to recover the complete set of excited populations and hence any required spectral emission .
the progressive condensation described above can be viewed as simply one of a number of possible paths which might be preferred because of special physical conditions or observations .    for light element ions ,
four types of bundling and condensation are distinguished in this work :    * ground parent , spin summed bundle - n @xmath121 lowest n - shell . * parent and spin separated bundle - n @xmath121 lowest spin system n - shell ( the _ bundle - ns population model _ ) .
* low ls resolved @xmath121 metastable states ( the _ low - level population model _ ) .
* parent and spin separated bundle - n @xmath121 low ls resolved @xmath121 metastable states .
type ( a ) corresponds to the approach used in summers ( 1974 ) .
type ( c ) corresponds to the usual population calculation for low levels in which ( consistent ) recombination and ionisation involving excited states are ignored .
it establishes the dependence of each population on excitation for the various metastables only .
type ( d ) , effectively the merging of ( b ) and ( c ) in the manner described earlier , is the principal procedure to be exploited in this work for first quality studies .
details are in the following sub - sections .
.[table : tab1]bundle-@xmath122 calculation pathways .
the parent / spin system weight factor is defined in equation [ eqn : eqn32 ] and [ eqn : eqn33 ] .
@xmath123 indicates the number of metastables of the recombined ion associated with the parent / spin system .
the _ parent index _
, sequentially numbering the different parents shown in brackets , is used as the reference for tabulation of coefficients . [ cols="^,^,^,^,^,^ " , ]     figure [ fig : fig5 ] illustrates the behaviour of low level populations .
the graph is of the parameter @xmath124 from equation [ eqn : eqn20 ] .
energy levels for high bundle-@xmath122 levels and their a - values , maxwell averaged collision strengths , radiative recombination coefficients , dielectronic recombination coefficients and ionisation coefficients are generated from a range of parametric formulae and approximations described in earlier works ( burgess and summers , 1976 ; summers , 1977 ; summers and hooper , 1983 ; burgess and summers , 1987 ; summers , 2004 ) .
high quality specific data when available are acccessed from archives and substituted for the default values .
this is a systematic procedure for dielectronic coefficients ( see sections [ sec : sec2.3.1 ] and [ sec : sec5 ] below ) .    for low levels , a complete basis of intermediate coupling energy level , a - value and scaled born approximation collision data , spanning the principal quantum shell range @xmath125
is generated automatically using the cowan ( 1981 ) or autostructure ( badnell , 1986 ) procedures .
this is called our _ baseline _ calculation .
these data are merged with more restricted ( in level coverage ) but similarly organised highest quality data from archives where available ( e.g. r - matrix data such as ramsbottom ( 1995 ) ) .
the data collection is compressed by appropriate summing and averaging to form a complete ls term basis and augmented with comprehensive high quality ls resolved dielectronic recombination , radiative recombination and collisional ionisation coefficients mapped from archives ( see sections [ sec : sec2.3.1 ] , [ sec : sec2.3.2 ] and [ sec : sec5 ] below ) .
the radiative data has its origin in the work of burgess and summers ( 1987 ) .
the final data collection is called a _ specific ion file _ ( adas data format _ adf04 _ ) .
the detailed content is examined in section [ sec : sec3 ] .
details of state selective dielectronic recombination and ionisation coefficients are given in the following sub - sections .
state selective dielectronic recombination coefficients are required to all resolved low levels and to all bundle-@xmath122 shells for the various initial and intermediate state metastable parents @xmath47 for @xmath8 modelling .
these are very extensive data and have been prepared for the present @xmath8 work through an associated international ` dr project ' summarised in badnell ( 2003 ) [ hereafter called _ dr - paper i _ ] and detailed in subsequent papers of the series .
based on the independent particle , isolated resonance , distorted wave ( ipirdw ) approximation , the partial dielectronic recombination rate coefficient @xmath126 from an initial metastable state @xmath47 into a resolved final state @xmath37 of an ion @xmath127 is given by @xmath128 where @xmath129 is the statistical weight of the @xmath130-electron doubly - excited resonance state @xmath38 , @xmath131 is the statistical weight of the @xmath132-electron target state and the autoionization ( @xmath133 ) and radiative ( @xmath134 ) rates are in inverse seconds .
the suffix @xmath135 is used here to denote a parent ion state .
@xmath136 is the energy of the continuum electron ( with angular momentum @xmath101 ) , which is fixed by the position of the resonances .
note that the parent states @xmath135 are excited , that is they exclude the metastable parents @xmath83 .
the code autostructure ( badnell , 1986 ; badnell and pindzola , 1989 ; badnell , 1997 ) was used to calculate multi - configuration ls and intermediate coupling energy levels and rates within the ipirdw approximation .
the code makes use both of non - relativistic and semi - relativistic wavefunctions ( pindzola and badnell , 1990 ) and is efficient and accurate for both the resolved low level and high-@xmath100 shell problems .
lookup tables ( see section [ sec : sec5 ] ; _ adf09 _ ) are prepared comprising state selective recombination coefficients at a standard set of @xmath14-scaled temperatures , for each metastable parent , to all ls resolved terms of the recombined ion with valence electron up to @xmath100-shell @xmath137 ( normally @xmath138 ) and to bundle-@xmath122 levels of a representative set of @xmath100-shells ( usually spanning @xmath139 to 999 ) . these bundle-@xmath122 coefficients are simple sums over orbital states @xmath140 and so apply at zero density .
this provides an extensive , but still economical , tabulation .       in _ dr paper - i _
, we introduced an associated code , the ` burgess
bethe general program ' bbgp .
bbgp is used here as a support function in a model for the @xmath101-redistribution of doubly - excited states which provide a correction to the accurate , but unredistributed , dielectronic data .
the redistributed data ( regenerated in the same _ adf09 _ format ) are normalised to ipirdw at zero density .
the procedure is similar to that for singly excited systems . for @xmath141-averaged levels ,
the number densities expressed in terms of their deviations , @xmath142 from saha - boltzmann , and referred to the initial parent @xmath83 , are given by @xmath143^{3/2 }   \frac{\omega_{p , nl}}{\omega_{\nu'}}e^{-e / kt_\rme}b_{p , nl}\ , .
\label{eq22}\end{aligned}\ ] ] thus , in the bbgp zero - density limit , with only resonant capture from the @xmath47 parent balanced by auger breakup and radiative stabilization back to the same parent , we have @xmath144 in the extended bbgp model , we include resonant - capture from initial metastables other than the ground , dipole - allowed collisional redistribution between adjacent doubly - excited @xmath101-substates of the same @xmath100 by secondary ion- and electron - impact , and losses by ` alternate ' auger break - up and parent radiative transition pathways . the population equations for the @xmath101-substates of
a doubly - excited @xmath100-shell become @xmath145 @xmath146 denotes resonance capture coefficients .
@xmath147 denotes the number of parent metastables which are the starting point for resonance capture whereas @xmath148 denotes the number of true excited parent states .
ion impact redistributive collisions are effective and are represented in the equations as @xmath149 ion contributions .
the density corrected bundle-@xmath122 recombination coefficients are then @xmath150 details are given in _
dr paper - i_. figure [ fig : fig6 ] shows the effects of redistribution of the doubly excited @xmath142 factors in the recombination of @xmath151 .
autoionisation rate calculations in ls coupling exclude breakup of non - ground metastable parent based ` singly excited ' systems by spin change .
such breakup only occurs in intermediate coupling but can not be ignored for a correct population assessment even in light element ions .
these rates are computed in separate intermediate coupling dielectronic calculations .
these rates are main contributors to the parent metastable cross - coupling identified earlier ( cf . the @xmath152 $ ] intermediate states in figure [ fig : fig4]b leading back to the @xmath153 parent ) .
direct ionisation coefficients for excited n - shells in bundle-@xmath122 population modelling are evaluated in the exchange - classical - impact parameter ( ecip ) approximation ( burgess and summers , 1976 ) .
the method which merges a symmetrised classical binary encounter model with an impact parameter model for distant encounters is relatively simple and of moderate precision , but has a demonstrated consistency and reliability for excited n - shells in comparison with more eleborate methods .
there is a higher precision requirement for state selective ionisation from metastable and low levels .
such ionisation includes direct and excitation / autoionisation parts , but the latter only through true excited parents .
stepwise ionisation is handled in the collisional - radiative population models .
there are extensive ionisation cross - section measurements , but these are in general unresolved and so are of value principally for renormalising of resolved theoretical methods .
the most powerful theoretical methods used for excitation , namely r - matrix with pseudostates ( rmps ) ( bartschat and bray , 1996 ; ballance , 2003 ) and convergent - close - coupling ( ccc ) ( bray and stelbovics , 1993 ) , have some capability for ionisation , but at this stage are limited to very few electron systems as is the time - dependent - close - coupling ( tdcc ) ( pindzola and robicheaux , 1996 ; colgan , 2003 ) ionisation method . for @xmath8 calculations , we have relied on the procedures of summers and hooper ( 1983 ) for initial and final state resolution of total ionisation rate coefficients and on the distorted wave approximation .
the distorted wave method is our main method for extended @xmath8 studies for many elements and most stages of ionisation .
rmps and tdcc studies are directed mostly at the neutral and near neutral systems where the distorted wave method is least reliable ( loch , 2005 ) .
we use the configuration average distorted ( _ cadw _ ) wave approach of pindzola , griffin and bottcher ( 1986 ) .
it has reasonable economy of computation while allowing access to complex , multi - electron ions , highly excited states , excitation / autoionisation and radiation damping .
it expresses the configuration averaged excitation cross - section for @xmath154 as @xmath155 and the configuration averaged ionisation cross - section for @xmath156 as @xmath157 where @xmath158 are squared two - body coulomb matrix elements , @xmath159 denotes the average cross - section and @xmath160 , @xmath161 and @xmath162 denote average initial , final and ejected electron momenta respectively in the configuration average picture .
the configuration average direct ionisation cross - sections and rate coefficients may be unbundled back to resolved form using angular factors obtained by sampson and zhang ( 1988 ) .
note that the excitations described here are to auto - ionising levels and so resolved auger yields may be used which are the same as those in the dielectronic calculations of section [ sec : sec2.3.1 ] above .
in fact the ionisation coefficient calculation results are structured and archived in adas ( format _ adf23 _ ) in a manner very similar to state selective dielectronic recombination .
extensive studies have been carried out on light and heavy elements ( colgan , 2003 ; loch , 2003 ) .
for bundle-@xmath122 modelling , the expected fundamental data precision is @xmath163 30% for excitation and ionising collisional rate coefficients , @xmath163 10% for a - values and state selective recombination coefficients and
@xmath163 1% for energies .
cccccccc ion & e & a & & @xmath164 & @xmath165 & @xmath103 + & . & & @xmath166 & @xmath167 & & & + @xmath168 & a & b & b & b , c & b & b & b + @xmath169 & a & b & b & b , c & b & b & b + @xmath170 & a & b & a , b & b , c & b & b & b + @xmath171 & a & a , b & b & b , c & b & b & b + @xmath172 & a & a & b & b , c & b & b & b + @xmath173 & a & a & b & b & b & b & b + @xmath174 & a & b , c & b & b & b & b & b + @xmath175 & a & b , c & b & b & b & b & b + @xmath176 & a & b , c & a , b & b & b & b & b + @xmath177 & a & b & a & a & b & b & b + @xmath178 & a & b & a & a & b & b & b + @xmath179 & a & b & b & b & b & b & b + @xmath180 & a & a & b & b , c & b & b & b + @xmath181 & a & a & b & b , c & b & b & b + @xmath182 & a & b , c & b & c & b & b & b + @xmath183 & a & c & b & b & b & b & b + @xmath184 & a & b & a & b & b & b & b + @xmath185 & a & b & a & b & b & b & b + @xmath186 & a & b & a & b & b & b & b + @xmath187 & a & b & a & b & b & b & b + @xmath188 & a & b & a & b & b & b & b + @xmath189 & a & b & b & b & b & b & b + @xmath190 & a & a & b & c & b & b & b + @xmath191 & a & a & b & c & b & b & b +    for low level modelling many sources are used .
a rating is given for the classes of fundamental data for the ions of carbon , oxygen and neon in table [ table : tab4 ] which is based on the following considerations .
for energy levels , categories are _ a _ spectroscopic , _ b _ @xmath163 0.5% and _ c _ @xmath163 1.0% .
c _ is anticipated from ab initio multi - configuration structure calculations , _ b _ from such calculations with extended optimising and _ a _ reflects direct inclusion of experimental energies from reference sources . for a - values , categories are _ a _ @xmath163 5% , _ b _ @xmath163 10% and _ c _ @xmath163 25% .
c _ is anticipated from our baseline calculations , _ b _ from optimised multi - configuration structure calculations with extended optimising and _ a _ from specific studies in the literature . for electron impact maxwell averaged collision strengths , @xmath192 , the categories are _ a _
@xmath163 10% , _ b _ @xmath163 20% and _ c _ @xmath163 35% .
c _ is from our baseline calculations , _ b _ from distorted wave calculations and _ a _ from specific r - matrix calculations , equivalent methods or experiment . for radiative recombination , the categories are _ a _ @xmath163 5% , _ b _ @xmath163 10% and _ c _ @xmath163 50% .
c _ is from scaled methods using hydrogenic matrix elements , _ b _ from distorted wave one - electron wave functions in an optimised potential using spectroscopic energies and _ a _ from specific r - matrix calculations and experiment .
b _ is the baseline in adas . for dielectronic recombination ,
the categories are _ a _
@xmath163 20% , _ b _ @xmath163 30% and _ c _ @xmath163 45% .
c _ is from bbgp approximations , _ b _ from ls - coupled calculations using autostructure from the dr project and _ a _ from ic - coupled calculations using autostructure with parent and lowest resonance energy level adjustments from the dr project .
b _ is the baseline in adas .
it is to be noted that the theoretical relative precision which is consistent with the variation between the three categories is 15% better but dielectronic recombination comparisons with experiment still show unexplained discrepancies at the 20% level so the present categories are safe . for ionisation , the categories are _ a _ @xmath163 10% , _ b _ @xmath163 25% and _ c _ @xmath163 40% .
c _ is from ecip approximations , _ b _ from configuration average with angular factor term resolution and _ a _ from rmps , tdcc calculations and experiment .
b _ is the baseline in adas .    as indicated in section [ sec : sec2.3 ] , these various data are assembled in an _
adf04 _ file which is sufficient to support the primary low - level population calculation .
the tabulations are at a set of temperatures arising from a fixed set of z - scaled temperatures ( see section [ sec : sec2.1.1 ] ) which spans the full range to asymptotic regions of reaction data .
collision data are converted to these ranges using _ c - plots _ ( burgess and tully ( 1992 ) and this procedure also flags data errors or queries .
the precision of the specific ion file determines the achievable precision of all derived populations , emissivities and collisional - radiative coefficients ( whiteford , 2005 ; omullane , 2005 ) .
the assembling of data in the _ adf04 _ is systematic and orderly . a comment section at the end of the file details the assembly steps , implementer , codes used and dates .
this includes baseline and supplementation files , merging , ls compression , dielectronic recombination data inclusion etc .
also there is extended detail of orginal data sources and a history of updates . a given _ adf04 _ file represents a ` snapshot ' of the state of available knowledge at the time .
it is subject to periodic review and adas codes ( see section [ sec : sec5 ] ) are designed to enable easy reprocessing of all derived data following fundamental data update .
the grading for each ion given in table [ table : tab4 ] is justified and supported by the comments from its _ adf04 _ file .
a summary from @xmath193 is given in the following paragraphs in illustration .
for @xmath188 , the low levels span 44 terms , including up to the @xmath194 shell built on the parent @xmath195 and @xmath196 shell built on the @xmath197 .
only the @xmath195 parent is treated as a metastable from the @xmath8 point of view .
the intermediate coupled baseline data set is _copmm@xmath19810-ic@xmath198ne6.dat_ with preferred supplementary energy a - value and @xmath192 data merged from _
copjl#be - ic@xmath198ne6.dat_. ionisation potentials and energy levels are from the nist standard reference database apart from @xmath199 ( kelly , 1987 ) and @xmath200 levels by ( ramsbottom , 1995 ) .
the categorisation is _
a_. a - values were drawn from the opacity project ( opacity project team , 1996 ; tully , 1991 ) as justified ( for n iv and o v ) by wiese , fuhr and deters ( 1996 ) and supported / adjusted by @xmath163 3% by fleming ( 1996a , b ) , jnsson ( 1998 ) , froese fischer , godefeid and olsen ( 1997 ) , froese fischer , gaigalas and godefried(1997 ) , nussbaumer and storey ( 1979 ) , sampson , goett and clark ( 1984 ) and ramsbottom ( 1995 ) .
b _ is safe .
@xmath192s are taken from ramsbottom ( 1995 ) .
these are for a 26 ls eigenstate multichannel r - matrix calculation .
the categories assigned are _ a _ and _ b_. radiative and dielectronic recombination and ionisation all follow the baselines in adas , that is categories _ b _ although the summed and averaged ionisation rate coefficient ( over metastables ) is normalised to experiment .
figure [ fig : fig7 ] shows @xmath201 coefficients for the c ii 858    spectrum line .
the coefficients depend on both electron temperature and electron density in general .
a common practice in spectral analysis is to observe principally the strongest resonance line of an ion .
such emission is driven largely from the ground state ( figure [ fig : fig7]a ) and because of the large a - value , the density sensitivity occurs at relatively high density .
thus such resonance line emission at moderate to low densities mostly reflects temperature and the distribution of the ionisation stages .
comparison in near equilibrium ionisation balance plasmas of line ratios from the same ion , by contrast , is mostly directed at electron density and relies on the presence of metastables and spin changing collisional processes to confer the sensitivity . as shown earlier , the balance of the dominant ground and metastable populations is disturbed in dynamic plasmas and so density sensitivity may be modified by the dynamic state .
the distinction of the metastable driven @xmath201 in @xmath8 modelling , as illustrated in figure [ fig : fig7]b , allows more complete diagnostic study and the possibility of separation of the two effects . in strongly recombining plasmas
( most commonly photo - ionised astrophysical plasmas ) the direct contribution of recombination to the emission may dominate the excitation part .
the @xmath202 , as illustrated in figure [ fig : fig7]c , is then required . in the fusion context ,
multi - chordal spectral observations are important for the study of impurity transport especially near sources .
visible and quartz uv observations are convenient and this places a requirement for @xmath201s from higher quantum shells .
it is this primarily which defines the span of our _ low - levels _ for population modelling .
emission from higher @xmath100-shells is significantly affected by cascading from yet higher levels .
the full machinery of projection as described in section [ sec : sec2.1 ] is necessary for our global ambition of 20% precision for emissivity coefficients .
figure [ fig : fig8 ] illustrates the @xmath8 recombination coefficients . at low electron density , radiative and dielectronic recombination dominate . for capture from metastables ,
alternate auger branching can largely suppress the dielectronic part of the surfaces .
thus figure [ fig : fig8]a shows the characteristic exponential rise at the temperature for excitation of the main parent ( core ) transition of dielectronic recombination and then the subsequent fall - off ( as @xmath203 ) , in contrast with [ fig : fig8]c . at moderate densities ,
suppression of the high @xmath100-shell populations , principally populated by dielectronic recombination , through re - ionisation occurs and the coefficient falls in the dielectronic recombination region . at very high density ,
three - body recombination becomes effective , preferentially beginning at the lower electron temperatures .
it is evident that recombination is less effective from the metastable at relevant ionisation balance temperatures .
models which ignore the role of capture from the metastable parent ( which may be the dominant population ) can lead to substantial errors in recombination coefficients , while simply excluding all capture from the metastable can not deliver the precision sought for current modelling . for light elements in astrophysical plasmas , the zero - density coronal assumption for the recombination coefficient is still frequently made .
this can not be justified even at solar coronal densities .
figure [ fig : fig9 ] illustrates the @xmath8 ionisation coefficients . at low electron density , the coefficient is dominated by direct ionisation , including excitation / auto - ionisation , from the driving metastable .
relatively high electron densities are required before the stepwise contribution begins .
it is primarily the excitation to , and then further excitation and ionisation from , the first excited levels which controls this .
the ground and metastable resolved coefficients both show the same broad behaviour as the usual ( stage to stage ) collisonal - radiative coefficient , tending to a finite limit at very high density .
it is to be noted that the @xmath204 coefficients are required to be able to construct a meaningful stage to stage collisional - radiative ionisation coefficient from the generalised progenitors .
it remains the case that most plasma modelling ( certainly in the fusion area ) is not adjusted to the use of the generalised coefficients as source terms .
reconstruction of stage to stage source terms ( at the price of a reduction in modelling accuracy ) from the generalised coefficients is still is a requirement and is addressed more fully in section [ sec : sec5 ] .
the generalised coefficients may be used to establish the equilibrium ionisation balance for an element in which the dominant ground and metastable populations are distinguished , that is the fractional abundances @xmath205}}{n^{[tot ] } } \right ) : \sigma=1,\cdots , m_z ; z=0 , \cdots , z_0\ ] ] where @xmath206 is the number of metastables for ionisation stage @xmath14 and @xmath207}=\sum_{z=0}^{z_0 } n^{[z]}=\sum_{z=0}^{z_0}\sum_{\sigma=0}^{m_z } n_{\sigma}^{[z]}\ ] ] in equilibrium .
writing @xmath208}$ ] for the vector of populations @xmath209}$ ] , the equilibrium population fractions are obtained from solution of the matrix equations . @xmath210 } & n_e~\underline{\scriptnew{r}}^{[1 \rightarrow 0 ] } & 0 & .    \\
n_e\underline{\scriptnew{s}}^{[0 \rightarrow 1 ] } & \underline{\scriptnew{c}}^{[1,1 ] } &          n_e~\underline{\scriptnew{r}}^{[2 \rightarrow 1 ] } & .    \\
0 & n_e\underline{\scriptnew{s}}^{[1 \rightarrow 2 ] } & \underline{\scriptnew{c}}^{[2,2 ] } & .     \\
0 & 0                                     & n_e\underline{\scriptnew{s}}^{[2 \rightarrow 3 ] } & .
\\           .
& .
\end{array }       \right ]      \left[\begin{array}{l }          \underline{n}^{[0 ] }          \\          \underline{n}^{[1 ] }          \\          \underline{n}^{[2 ] }          \\          \underline{n}^{[3 ] }          \\          .
\end{array }       \right]_{equil}=0.\end{aligned}\ ] ] these in turn may be combined with the @xmath211 and @xmath90 to obtain the equilibrium radiated power loss function for the element as @xmath212 } & = & \sum_{z=0}^{z_0 } p^{[z]}\left ( \frac{n^{[z]}}{n^{[tot ] } } \right ) _ { equil } \nonumber \\ &   = & \sum_{z=0}^{z_0}\sum_{\sigma=0}^{m_z } ( \scriptnew{plt}^{[z]}_{\sigma}+\scriptnew{prb}^{[z]}_{\sigma})\left ( \frac{n^{[z]}_{\sigma}}{n^{[tot]}}\right ) _ { equil}.\end{aligned}\ ] ] the equilibrium fractional abundances and equilibrium radiated power function are illustrated in figure [ fig : fig10 ] .
it is useful at this point to draw attention to emission functions which combine emission coefficients with equilibrium fractional abundances .
they are commonly used in differential emission measure analysis of the solar atmosphere ( lanzafame , 2002 ) where they are called @xmath213 functions .
in solar astrophysics , it is assumed that the @xmath213 are functions of the single parameter @xmath29 ( either from a zero - density coronal approximation or by specification at fixed density or pressure ) and usually the abundance of hydrogen relative to electrons @xmath214 is incorporated in the definition . for finite density plasmas , in the generalised collisonal - radiative picture , we define @xmath215 functions , parameterised by @xmath29 and @xmath27 as @xmath216}_{j \rightarrow k } &   = &   \sum_{\sigma}\scriptnew{pec}^{[z](exc)}_{\sigma , j \rightarrow k}\left ( \frac{n_{\sigma}^{[z]}}{n^{[tot ] } } \right ) _ { equil } \nonumber \\ & & + \sum_{\nu'}\scriptnew{pec}^{[z](rec)}_{\nu',j \rightarrow k}\left ( \frac{n_{\nu'}^{[z+1]}}{n^{[tot ] } } \right ) _ { equil}. \end{aligned}\ ] ] these are strongly peaked functions in @xmath29 .
the precision of the present modelling and data , including the full density dependence , is in principle sufficient to allow bivariate differential emission measure analysis ( judge , 1997 ) although this remains to be carried out .
the general behaviour of a @xmath215 function is illustrated in figure [ fig : fig10]c .
the organisation of the main calculations and data flow is shown in figure [ fig : fig11 ] .
the calculations executed for the paper have been implemented in general purpose codes and attention has been given to the precise specification of all data sets and the machinery for accessing and manipulating them .
this includes initial data , intermediate and driver data as well as the final @xmath8 products and follows adas project practice .
population structure modelling requires an initial input dataset of energy levels , transition probabilities and collisional rates of adas data format _ adf04 _ which is complete for an appropriate designated set of low levels . for the light elements ,
best available data were assembled and verified as described in section [ sec : sec3 ] .
these data are substitutes for more moderate quality , but complete , _ baseline _ data prepared automatically ( see _ gcr - paper iii _ ) .
the _ adf04 _ files are required for ls - terms . in practice , we find it most suitable to prepare the data for lsj levels and then bundle back to terms . for @xmath8 modelling ,
state selective recombination and ionisation coefficients must be added to form the fully specified _
the ls - coupled dielectronic coefficients are mapped in from very large comprehensive tabulations of data format _
_ prepared as part of the _ dr project _
( badnell , 2003 ) .
the resulting _ adf04 _
files are the adas preferred data sets and are available for elements helium to neon .    the two population codes , called adas204 and adas208 , work together .
adas204 is the bundle-@xmath122 model .
data on the metastable parent structure , quantum defects , auto - ionisation thresholds and autoionsation rates are required .
these data may be extracted from the _ adf04 _ files and it has been helpful to make the driver dataset preparation automatic .
high quality shell selective dielectronic data are essential and this is part of the provision in the _ adf09 _
files described in the previous paragraph .
adas204 provides complete population solutions , but extracts from these solutions the condensed influence on the low n - shells , as _ projection matrices _
, for connection with the calculations of adas208 .
it is to be noted that the main on - going development is refinement of atomic collision rates between the key low levels .
the projection matrices are not subject to frequent change and so are suitable for long - term archiving ( _ adf17 _ ) .
adas208 is the low - level resolved population model which delivers the final data for application .
it draws its key data from the fully configured _
_ file and supplements these with projection data .
the evaluation of the population structure takes place at an extended set of @xmath14-scaled electron temperatures and densities ( see section [ sec : sec2 ] ) and this means that the resulting @xmath8 coefficient data are suited to interpolation along iso - electronic sequences .
thus the initial tabulation of @xmath8 coefficients is in iso - electronic datasets .
it is convenient to implement the gathering and mapping from iso - electronic to iso - nuclear in a separate step , which also supports the merging back to the unresolved stage - to - stage picture if required .
the production of @xmath217 and @xmath218 coefficients is directly to iso - nuclear oriented collections .
the number of @xmath217s from the population calculations can , in principle , be very large .
we restrict these by a threshold magnitude and to particular important spectral regions .
it is straightforward to rerun adas208 to generate @xmath217 or @xmath218 coefficients alone in spectral intervals of one s choice .
the separation of the adas204 and adas208 tasks and the ease of modifying data within an _
adf04 _ file means that ` what - if ' studies on the sensitivity of the derived data to fundamental data uncertainly can readily carried out .
special adas codes enable detailed study of cumulative error and dominating sources of uncertainty such that error surfaces , for example for a @xmath217 as a function of electron temperature and density , may be generated .
such error ( uncertainty ) analysis for theoretical derived coefficients and its utilisation in the confrontation with diagnostic experiments is the subject of a separate work ( see omullane , 2005 ) .
the requirements for precise modelling of spectral emission and the relating of ionisation stages in thermal plasmas have been considered .
collisional - radiative methodologies have been developed and extended to enable the full role of metastables to be realised , so that this generalised ( @xmath8 ) picture applies to most dynamically evolving plasmas occurring in magnetic confinement fusion and astrophysics .
the procedures are valid up to high densities .
the studies presented in the paper explore the density effects in detail within the @xmath8 picture and show that the density dependencies of excited ion populations and of effective rate coefficients can not be ignored .
specific results are presented for light elements up to neon and the computations are carried out in an atomic basis of terms ( ls - coupled ) .
such modelling will remain sufficient up to about the element argon , beyond which a level basis ( intermediate coupling ) becomes necessary .
heavier elements will be examined in further papers of this series .
considerable attention has been given to the generation and assembly of high quality fundamental data in support of the @xmath8 modelling .
also datasets of fundamental and derived data have been specified precisely and codes have been organised following the principles of the adas project .
the product of the study is the preferred adas data for the light element ions at this time .
the atomic data and analysis structure , adas , was originally developed at jet joint undertaking .
bates d r , kingston a e and mcwhirter r w p 1962 _ proc .
soc . a _ * 267 * 297 .