Документ взят из кэша поисковой машины. Адрес оригинального документа : http://star.arm.ac.uk/preprints/2013/637.pdf Дата изменения: Tue Jul 30 13:47:49 2013 Дата индексирования: Fri Feb 28 00:54:32 2014 Кодировка: Поисковые слова: mercury program

July 30, 2013

Diagnosing transient ionization in dynamic events
J.G. Doyle1 , A. Giunta2,3 , M.S. Madjarska1 , H. Summers2,3 , M. O'Mullane2,3 , and A. Singh1
1 2 3 4

,4

Armagh Observatory, College Hill, Armagh BT61 9DG, N. Ireland Department of Physics, University of Strathclyde, 107 Rottenrow, Glasgow, G4 0NG, Scotland Space Science and Technology Department, STFC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire, OX11 0QX, UK. Dept.of Physics and Electronics, Deen Dayal Upadhyaya College, University of Delhi, India

Received date, accepted date
ABSTRACT

Aims. The present study aims to provide a diagnostic line ratio that will enable the observer to determine whether a plasma is in a state of transient ionization. Methods. We use the Atomic Data and Analysis Structure (ADAS) to calculate line contribution functions for two lines, Si IV 1394 е and O IV 1401 е, formed in the solar transition region. The generalized collisional-radiative theory is used. It includes all radiative and electron collisional processes, except for photon-induced processes. State-resolved direct ionization and recombination to and from the next ionization stage are also taken into account. Results. For dynamic bursts with a decay time of a few seconds, the Si IV 1394 е line can be enhanced by a factor of 2­4 in the first fraction of a second with the peak in the line contribution function occurring initially at a higher electron temperature due to transient ionization compared to ionization equilibrium conditions. On the other hand, the O IV 1401 е does not show such any enhancement. Thus the ratio of these two lines, which can be observed with the Interface Region Imaging Spectrograph, can be used as a diagnostic of transient ionization. Conclusions. We show that simultaneous high-cadence observations of two lines formed in the solar transition region may be used as a direct diagnostic of whether the observed plasma is in transient ionization. The ratio of these two lines can change by a factor of four in a few seconds owing to transient ionization alone.
Key words. Sun: corona - Sun: transition region - Line: profiles - Atomic processes ­ Line: formation

1. Introduction
All spectral lines have encoded information that allows the observer to diagnose important physical parameters of the underlying plasma. For example, forbidden or inter-system lines allow an evaluation of the plasma's electron density, while resonance lines allow the observer to diagnose the electron temperature. With sufficient spectral resolution, all spectral lines can be used to give the observer information on plasma flows and/or turbulent motions. To diagnose a plasma's state of ionization however, requires both high-cadence data and suitable spectral lines. Various authors have presented simulated spectra based on transient ionization (Mewe & Schrijver, 1980; Reale & Orlando, 2008) for comparison with observational data. In a recent paper, Doyle et al. (2012) discussed the diagnostic potential of high-cadence ultraviolet spectral data when transient ionization is considered. The above paper used high-cadence spectral line data from the Solar Maximum Mission (SMM) observed in O V 1371 е which allowed the authors to measure electron densities and temperatures during the early stages of a feature's evolution, something that is not currently possible. The high cadence UV spectrometer on SMM allowed observations of selected spectral lines with a subsecond cadence. The forthcoming observations from the Interface Region Imaging Spectrograph1 will once again enable high quality observations of lines formed in the solar transition region. What
Send offprint requests to: madj@arm.ac.uk 1 http://iris.lmsal.com/

we would like to have is a simple line ratio consisting of one line that shows a response to transient ionization and another line that is not responsive to transient ionization. Then, provided we can monitor these lines with sufficient cadence, any fast ( 0.5 s) increase in the ratio must be due to non-equilibrium conditions. Here, we look at the response of two such lines, Si IV 1394 е and O IV 1401 е. The Si IV 1394 е line is in the ground spin system (3s 2 S1/2 ­ 3p 2 P3/2 ). The character of the excitation cross-sections (dipole and non-spin change) and access to higher n-shell/cascade will result in a line enhancement. However, the O IV 1401 line, by contrast, has an upper state in the quartet spin system (2s2 2p 2 P3/2 ­ 2s2p2 4 P5/2 ) which means that access to this spin system is driven by cross-sections that decrease with temperature, hence little or no line enhancement, therefore, transient ionization under-fills the population structure of the alternate spin systems from the ground. By contrast, the ground spin system population structure gets enhanced (due to the exponential factor in the rate coefficient) with increasing electron temperature.

2. Transient Ionization
The atomic structure of atoms and ions is in principle an infinite assembly of levels with an infinite number of reactions between them, however, simplifying assumptions about the nature of the plasma, its dynamic character and the relative importance of the various reactions have to be made. As shown in the recent paper by Doyle et al. (2012), the generalized collisional-radiative theory (GCR) is the most appropriate approach. The above
1


Si IV 1394 contribution function 5e-08

contribution function

4e-08

0.2 sec 0.5 sec

3e-08

1 sec 2 sec

2e-08 10 sec 1e-08

EQ

0

4

4.5

5 log Te

5.5

6

6e-10 O IV 1401 contribution function 5e-10

i-i >> e-e . More details on the above may be found in Summers et al. (2006). From a dynamic point of view, the intrinsic and extrinsic groups are to be compared with each other and with timescales, pla sma , representing plasma ion diffusion across temperature or density scale lengths, the relaxation times of transient phenomena, and observation times. For astrophysics plasmas, pla sma g m >> o >> e-e , where g represents the relaxation time of the ground-state populations of ions. These time scales imply that the dominant population of impurities in the plasma are those of the ground and metastable states of the various ions. In a transient event, the dominant populations evolve on time scales similar to the plasma diffusion time scales and so should be modelled dynamically. From a theoretical point-of-view, non-equilibrium occurs when pla sma ion,rec = 1/[Ne (S (zz+1) + (z+1z) )] (McWhirter, 1965), where Ne is the electron density, and S (zz+1) and (z+1z) are the effective ionization and recombination coefficients (see Section 2.2). This means that in a transient event, the plasma time scale ( pla sma ) is shorter than the atomic ionization/recombination time scale (ion,rec ).
2.2. GCR

contribution function

4e-10 0.2 sec 3e-10 0.5 sec 1 sec 2 sec 10 sec

2e-10

1e-10

EQ

0

4

4.5

5 log Te

5.5

6

Here, we calculate the contribution functions for Si IV 1394 е and O IV 1401 е using the GCR theory. Each ion in an optically thin plasma is described by a complete set of levels with collisional and radiative couplings between them. All radiative and electron collisional processes are included, except for photon-induced processes. In addition, state-resolved direct ionization and recombination to and from the next ionization stage are also taken into account. The emissivity of a spectral line is given by
ji

= Ael

Fig. 1. Line contribution function for Si IV 1394 е and O IV 1401 е in cm-3 s-1 at ionization equilibrium and for transient ionization. The electron density was fixed at 1010 cm-3 , with an initial starting temperature of 30 000K.

NH 2 NG Ne e

(z) ji

(T e , Ne , t)

(1)

where Ael = N (Z ) /NH is the abundance of the element Z relative to hydrogen, NH /Ne , tabulated by McWhirter et al. (1975), and G(jz) i (T e , Ne , t) is the time-dependent contribution function defined as G
(z) ji

theory, plus ionization, recombination, excitation and radiative decay data are all as implemented in ADAS, the Atomic Data and Analysis Structure2 . A good description of this is given in the Appendix of Lanza et al. (2001). Below we first discuss the various time constants and their relative importance.
2.1. Time constants

(T e , Ne , t) = PEC(exc, ji
(exc,z) ji

z)

N (z+1) (t) N (z) (t) + PEC(reci,z) j (Z ) N N (Z )
(exc,z) (rec,z) ji

(2)

where PEC

=

A

ji

F

is the excitation photon

The lifetimes of the various states of ions and electrons in a plasma too radiative or collisional processes divide into two groups. The first is the intrinsic group, comprising purely of atomic parameters, and includes the metastable radiative decay time, m , the ordinary excited-state radiative decay time o , and auto-ionizing-state decay time, a . These are generally ordered as a << o << m . The second is the extrinsic group, which depends on plasma conditions, especially density. It includes free particle thermalization including electron-electron e-e , ion-ion i-i , and ion-electron i-e , charge-state change (ionization ion and recombination rec ) and redistribution of population amongst excited ion states (red ). Their order is ion,rec >> i-e >>
2

emissivity coefficient, and PEC = A ji F (rec,z) the recombination photon emissivity coefficient. In the above, A ji is the radiative transition probability, N (z) and N (z+1) are the population densities of the ground states of the ions of charge z and z + 1 and F (exc,z) , and F (rec,z+1) are the effective contributions to the population of the upper excited state i. In a time-dependent plasma model, the line emissivity is no longer a unique function of the local temperature and density conditions, but it depends on the past history of the temperature, density, and state of ionization of the plasma. Therefore, the assumption of ionization equilibrium in calculating the ionization balance is not appropriate and time-dependent fractional abundances must be determined. The time dependence of ionization stage populations, N (z) , leads to the following equation: d N (z) = Ne [S dt
(z-1)

http://www.adas.ac.uk/manual.php

N

(z-1)

+ (S

(z)

+ (z) )N

(z)

+ (z+1) N

(z+1)

] (3)

2


where S and are the collisional-dielectronic ionization and recombination coefficients. They give the contribution to the growth rates for the ground state population due to the effective ionization, which includes direct and excitation/auto-ionization contributions, and the effective recombination, which includes radiative, dielectronic, and three-body contributions. The values of these coefficients, currently within the ADAS database, have been obtained following the GCR approach as described in Summers et al. (2006). The solution of Equation (3) is such that the number density of the element of a nuclear Z charge Z , N (Z ) , is equal to z=0 N (z) . The time-dependent (z) (Z ) are calculated as follows: fractional abundances N (t)/N the code derives the solution for a range of fixed plasma electron temperature and density pairs, starting from an initial population distribution N (z) (t = 0)/N (Z ) using an eigenvalue approach. At temperatures close to ionization equilibrium, the contribution of recombination to the emissivity is small; however, this contribution becomes large in non-equilibrium conditions. Further details on the above may be found in Doyle et al. (2012) and references therein. Here we use an initial starting temperature of 30 000K and a fixed electron density of 1010 cm-3 . We then calculate the ionization fraction and the resulting line contribution function for different relaxation times, ranging from 0.2 sec after the initial start (i.e., in a highly transient state) to 10 sec (close to ionization balance) and finally when the plasma is in ionization equilibrium.

References
Curdt, W., Brekke, P., Feldman, U., Wilhelm, K., Dwivedi, B.N., Schuhle, U. & Lemaire, P., 2001, A&A 375, 591 Doyle, J.G., Giunta, A., Singh, A., Madjarska, M.S., Summers, H., Kellet, B.J. & O'Mullane, M., 2012, Sol Phys. 280, 111 Lanza, A.F., Spadaro, D., Lanzafame, A.C., Antiochos, S.K., MacNeice, P.J., Spicer, D.S. & O'Mullane, M.G., 2001, ApJ 547, 1116 McWhirter, R.W.P., Thonemann, P.C. & Wilson, R., 1975, A&A 40, 63 McWhirter, R.W.P., 1965, in Plasma diagnostic techniques, ed. R. H. Huddlestone & S. L. Leonard (Academic Press New York London) Mewe, R. & Schrijver, J., 1980, A&A 87, 261. Phillips, K. J. H., Feldman, U., & Landi, E. 2008, Ultraviolet and X-ray spectroscopy of the solar atmosphere (Cambridge University Press) Reale, F. & Orlando, S., 2008, Sol Phys. 684. 715. Spitzer, L., 1962, Physics of Fully Ionized Gases, Interscience, New York, p120 Summers, H.P., Dickson, W.J., O'Mullane, M.G., Badnell, N.R., Whiteford, A.D., Brooks, D.H., Lang, J., Loch, S.D. & Griffin, D.C., 2006, Plasma Phys. Controlled Fusion, 48, 263 (arXiv:astro-ph/0511561)

3. Results & discussion
In Fig. 1, we give the line contribution function for different relaxation times, that are shorter than the time for the plasma to reach ionization equilibrium. This figure shows that the peak of the line contribution function for both lines occurs at a higher temperature for short relaxation times. However, for O IV 1401 е, transient ionization does not produce a flux increase over the one for ionization equilibrium, while the Si IV 1394 е line increases over a factor of three. The reason these two lines behave in a different manner under transient ionization is due to their formation process as outlined in Section 1. For O IV 1401 е at an electron density of 1010 cm-3 , it takes 10 sec for the line to reach ionization equilibrium, while Si IV 1394 е takes 100 sec to reach equilibrium. Figure 1 shows that the ratio of Si IV 1394 to O IV 1401 ranges from 30 in ionization equilibrium to 120 for highly transient conditions (0.2 s). Multiplying by the relative elemental Si/O coronal abundance (Phillips et al. 2008) gives a ratio ranging from 3.3 to 12.9. The Curdt et al. (2001) SUMER spectral atlas gives the observed ratio for these lines as 4.6 in a sunspot and 6.9 in the quiet Sun. The above observational ratios are in general agreement (to within a factor of two) with the present Si IV 1394/O IV 1401 ratio when assuming ionization equilibrium. However, with a different DEM distribution, different abundances will change the ratio but the important point to note by the observer is not the absolute ratio but whether the observed ratio would change quickly. As a result, if observing these two lines simultaneously with, say, a 0.5 s cadence, then any sudden increase by a factor of 2­3 in the line ratio must be due to transient ionization. Such observations will therefore enable an evaluation of the plasma's electron density as outlined in Doyle et al. (2012).
Acknowledgements. Research at the Armagh Observatory is grant-aided by the N. Ireland Dept. of Culture, Arts and Leisure. We thank UK STFC for support via ST/J001082/1 and the Leverhulme Trust.

3