A New And Improved Version Of HULLAC
|
|
- Gabriella Barton
- 6 years ago
- Views:
Transcription
1 A New And Improved Version Of HULLAC M. Klapisch a, M. Busquet a,b and A. Bar-Shalom a,c a ARTEP #,inc Excelsior Springs Court, Ellicott City, MD 21042,USA b LERMA, Observatoire de Paris, UPMC, CNRS, Place Jules Janssen, Meudon, France c NRCN, 9001 Be er Sheva, Israel. Abstract: We present a new version of the collisional radiative model (CRM) generator code HULLAC. The main features are: (i) input considerably simplified and flexible, (ii) capacity of post-averaging configurations and superconfigurations in mixed mode, (iii) a new fitting formula for cross sections, completely correcting the problems of the classical Sampson and Golden formula, (iv) a new algorithm for solving the rate equations of the CRM, more robust and giving more insight in the quality of the model than the biconjugate gradient method, (v) thanks to thorough comparisons with the LANL code, some errors were corrected, and very good agreement has been obtained on all types of transitions, (vi) finally, most of the code has been re-written according to up-to-date standards. This code is now ready for distribution. Keywords: Non-LTE, collisional radiative model, collisional cross sections, rate equation solver. PACS: 33.20Rm, 52.20Fs, 52.25Dg, 52.25Os. INTRODUCTION The HULLAC code [1] is an integrated collisional radiative model generator and solver. Its efficiency is due to original methods and algorithms. All relevant processes are computed with the same relativistic wavefunctions [2]. The angular momentum recoupling coefficients are computed using graphical methods [3]. The collision processes use the factorization-linearization method [4],[5]. It has been used by a number of groups, and was advertised at the 2001 APIP [6]. However, because the different modules were developed over the years, most of the coding was obsolescent. This is why it was decided to rewrite 1 most of the code in a modern style, more conform to standards of software engineering, in order for the code source to be understood and compiled easily by potential users. By the same token, many improvements have been made, and errors have been corrected. The rate equation solver, which was missing in previous versions, has been developed, based on a new algorithm. It is the purpose of this paper to describe these improvements and additions. The reader is referred to the above references for the theory. to whom correspondence should be addressed: marcel.klapisch@nrl.navy.mil # Contractor to Naval Research Laboratory, Washington, DC, mostly done by M.B. 206
2 NEW INPUT POSSIBILITIES The module anglar acts as the master module, in addition to computing all the angular momentum coefficients of the various matrix elements of the energy and of all the transitions, spontaneous and collisional. It reads an input file in which, among other data, the list of configurations to be computed should be listed. This was inefficient and painstaking for complex spectra. Defining Configurations It is now possible to use the concept of superconfigurations [7] to generate all the desired configurations in a very compact manner. FIGURE 1. Flow chart of the HULLAC package. AME stands for angular matrix elements. For instance, the following lines: superconf nel=15 abc 1s2 2L [5:8] 3L [0:18] (4L,5L) [0:1] endsuper generate all and only those configurations with a total of 15 electrons distributed between the supershells (2s,2p) with all combinations from 5 to 8 electrons, (3s,3p,3d) with zero to 18 electrons, and (4s,4p,4d,4f,5s,5p,5d,5f,5g) supershell with 0 to 1 electron. As a new addition, the symbols Le and Lo would select only even or odd values of l. The sizes of the arrays have been greatly increased to do very complex computations, but if there is an overflow, it will now be detected and the relevant parameter will be clearly designated in the output file. 207
3 Choosing Potentials There was in nearly all versions of HULLAC the possibility of computing different groups of configurations in different potentials. This can greatly increase the accuracy of the energies and wavefunctions. This is now achieved in a very simple way with the keyword: potential new. In the same way, the potential to be chosen as the final one with respect to which the other configurations are shifted to take into account their accurate computations is simply designated as potential new final. Direct- and Post- Averaging In the previous versions, all the input configurations created a list of fine structure levels in jj coupling. For complex spectra, this generates too many detailed levels. Therefore, several modes were introduced in the SETTING chapter of the input file: mode = l(evels) this generates as before fine structure levels for all configurations listed. mode = c(onfigurations) this generates directly, by simple sum rules, relativistic configuration averages as basic entities for the whole chain of program. mode = l+c makes two sets of output files for all the rates, and two different solutions can be obtained at the level of the CRM solver colrad. In addition, one can now perform post averaging, that is, anglar generates one of the two previous sets, but can then perform additional, brute force averaging. Since anglar is the less time consuming module of the chain, it is not a significant penalty. Two additional keywords are possible: nrc - creates non relativistic configurations averages (NRC) of relativistic configurations by statistical weight, and sc - creates superconfigurations as averages of configurations. These two averages can be done starting from mode c or l. In the latter case, the effects of configuration interaction will be included in the averages. These averages can be performed on all the levels or configurations in the list, or on selected levels/configurations, in effect creating a mixed model in which some configurations can be described as fine structure levels, others a configuration averages, and some others as superconfigurations. The reader is referred to the mode of use for more details A NEW FITTING FORMULA Collisional excitation cross sections are the most time-consuming part of setting the model. The temperature dependent rates are obtained by integrating them over the collision energy. Therefore, it is has become standard to represent the cross sections by an easily integrable analytical formula, the parameters of which are transmitted from the cross section module to the CRM equation solver module, so the same cross sections can be used to model several temperatures. The classical, nearly ubiquitous, formula is usually called the Golden and Sampson (GS) formula [8]. It reads: c!(u = " / #E) = A + D ln(u) + 1 (a + u) + c 2 (a + u), 2 208
4 where Ω is the collision strength and u is the collision energy divided by the transition energy. A, D, a, c 1, c 2 are the free parameters that are fitted to the computed cross sections for each transition by a least square method. The rates are then obtained as a combination of exponential integrals [9]. In principle, the 4 coefficients, A, D,c 1 and c 2 of the GS fit are positive and one or two are dominant. a has been introduced to further improve the fit and reproduce a possible slope change near threshold. It has to be larger than -1 to prevent singularity. However, in practice the coefficients are obtained with a least square method as suggested in GS, possibly with the Bethe limit constraint, and from a limited range of energies (because a finite number of sampled cross-sections is computed). But the least square method cannot guarantee accuracy of the fit over the whole sampled range, neither does it allow for a correct asymptotic behavior. Different shapes yield fitting inadequacies of different kind as can be seen in two cases displayed in Fig.2. In Fig. 2b, the reconstructed cross section (dashes) differs from the computed values (diamond) by less than 20% of the maximum cross-section but is a very poor approximation. The origin of disagreement is the presence of a maximum of the crosssection at a distance from the threshold, that cannot be reproduced by the GS formula (except using a less than -1 leading to a singularity) (a) 10-3 (b) Collision Strength Collision Strength E inc /!E E inc /!E FIGURE 2: Examples of problems with the GS formula. Reconstructed GS fit (dashed lines) and BBTC fit (solid line) compared to initial cross sections (diamonds) The consequence of using GS fit in the CRM solver is that a number of rates (sometimes up to 20%) turn out negative, and have to be ignored. This leads to errors in the solution or singularities in the matrices. The BBTC fitting procedure We developed a new fitting algorithm [10] based on the idea of Burgess, Tully and Chidichimo [11]. The idea is to make a mapping of the incident energy (with range from E to! ) into the compact support [0,1]. In reference [11], the actual form of the mapping and the position of energy samples are dependent on the transition, on the selection rules which are in force and on the actual values of the cross-sections. In order to make it tractable in our CRM model, we improved their procedure with a more general 209
5 procedure and a fixed, optimized set of sampling point. We map the function!(u) to a function z(x) with the following transformation: x = 1! (1 + C) / (u + C), z = ("! Alnu).(u + C) n. here n is the integer power of the residue (i.e. after the logarithmic term has been taken out) obtained from the asymptotic behavior deduced from the last 2 or 3 computed crosssections, and C is a free parameter to optimize the fit. The collision strength then reads:! = Alnu + z / (u + C) n. Note that our mapping differs from the one proposed by Burgess et al, as we have one transformation only, regardless of the selection rule relevant to the specific transition, and that we allow any integer asymptotic power expansion n. We then sample the cross sections on a fixed set of x values. We found that the set of 9 values { } l=0:8 = {0,!1,!3,!6,!10,!14,!17,!19,!20} / 20 is optimum for most cases. The cross x l sections are now fully described by 13 coefficients, and can be interpolated by splines. An example of the mapping can be seen on figure 3. FIGURE 3: Example of mapping of a cross section with the BBTC mapping scheme. The rates are easily recovered in the CRM solver by a combination of exponential integrals, and spline coefficients depending only on the set {x}. With this algorithm, we did not find any problematic values of rates in the Xe case with configurations described below. COMPARING WITH LANL CODE In the recent NLTE workshops [12], it became clear that the discrepancy between the various codes for CRM is noticeable. In an effort to understand the source of these differences, we compared HULLAC cross sections for collisional processes with the results of the LANL code 2. These two codes are based on the same physical assumptions, 2 We are thankful to T. Durmaz and R. Mancini for their help in this comparison. 210
6 but use quite different algorithms, so the results should be similar if not identical. We found at first some important discrepancies that forced us to check thoroughly some routines, and eventually we found and corrected previously undetected errors. The DWBA code of LANL was checked at LANL by comparing to another one, based on many body perturbation theory. HULLAC now agrees with both codes, as shown on figures 4 and 5. We are now confident that the cross sections for collisional processes are correct. FIGURE 4: Comparison of cross sections for collisional excitation of He-like Fe XXV 1s 2-1s2p. FIGURE 5: Comparison of collisional ionization cross sections on H-like Fe XXVI 1s nl. Figure 6 shows comparison between the codes on photo-ionization of H-like Fe XXVI. Again the agreement is very good. 211
7 FIGURE 6: Comparison between HULLAC and LANL codes on photo-ionization of H-like FeXXVI. The figures above show that the agreement between the codes is as good as can be for codes that were written differently by different people with different algorithms. THE NEW CRM EQUATIONS SOLVER Generalities HULLAC was until now distributed without a solver for the rate equations. Our aim is to be able to use HULLAC for building opacities databases for cases where the non-lte super transition array approach (SCROLL) [13] is not detailed enough. For this purpose, one needs a robust solver for very large linear systems. Indeed, it has been recognized that many ion stages are necessary, as well as many thousands of states for each ion stage. The matrix involved is sparse because of selection rules, and one popular method is the bi-conjugate gradient [14]. The problem with the latter is that when it does not converge, one cannot extract any useful information. This is a real concern, because the lifetimes of metastable and ground states are orders of magnitude longer than those of resonant states. Thus, the condition of the rate matrix can be in excess of The method adopted here is to solve each ion stage separately, and iterate. Each ion matrix s condition is bound to be smaller than the condition of the whole matrix, and because of the connection to neighboring ions, these ion matrices are not singular. Similar methods [15] involve the separation of ground and metastable states from excited states. This does not converge easily to the LTE limit and ignores transitions from excited states of an ion to excited states of the neighboring ions. Therefore, we factorize the population of individual states as a product of reduced population and total ion population M Q N Q =! N p ; n Q p = N p N Q p=1 212
8 where N p, N Q, n p are, respectively, individual population of state p, total population of charge Q, and reduced population of level p of charge state Q. The advantage of this factorization is that the global effective rates W Q!Q ' i.e. temperature and density dependent contribution from an ion to another, depend only on the reduced population: $ 1 ' W Q!Q ' = # & # N p W p! p ' % N ) = # # n Q p W p! p ' p '"Q ' Q p"q ( p '"Q ' p"q Therefore, if the reduced populations are found at a given iteration, the global rates can be computed. Then a matrix, the order of which is only the number of active charge states, i.e. about a dozen, can be easily set up and solved, giving the charge state distribution, i.e. the N Q of the next iteration. Using a Newton-style development, we find the increments of the reduced populations:!n Q = " 1!N Q ' M -1 Q R Q %n Q ' +!N Q" M -1 Q S Q %n Q" +!N Q %n Q "1 + M Q %G Q N Q { } where Q =Q+1, Q =Q-1, R Q and S Q are, respectively, the rectangular matrices of recombination and ionization into ion Q, and M Q is the matrix of transitions inside ion Q, including losses to other ions on the diagonal. %G Q is the vector of the residues for the states at the previous iteration. More details can be found in reference [16]. Results Using this new version of HULLAC, with the new fit and the new solver, we were able to obtain the solution of the CRM for Xe at T e = 50 ev, for a range of electron densities from to cm -3. The model included 39,683 states (relativistic configurations). Table 1: Details of number of levels and rates for the full Xe model Ion States Rad. Trans. Coll. Excit. Ioniz. & R. R. Autoioniz. Xe IX Xe X Xe XI Xe XII Xe XIII Xe XIV Xe XV Xe XVI Xe XVII Total x At each iteration, one gets automatically the ion charge distribution, therefore this algorithm allows one to follow the convergence of the process, and it is more meaningful 213
9 than a black box algorithm in which the iteration do not have physical meaning. It is then possible to stop at the required accuracy, and create a compact database of charge distribution. Figure 7 shows the evolution of the charge distribution of a smaller model (with 2386 non-relativistic configuration averages). In this case, the convergence criterion was very stringent, (10-4 relative variation of individual populations), to show the stability of the algorithm. Figure 8 displays the spectra obtained from the full model (39,683 states) at two different densities. There is no statistical approximation besides the configuration fraction(z) ion8 ion9 ion10 ion11 ion12 ion13 ion14 ion15 ion16 ion fraction(ion) NQest NQ 100 NQ 200 NQ 500 NQ 1.e_ average. iteration # ion FIGURE 7: Evolution of charge distribution for Xe at T e =50 ev and N e =10 20 cm N e = cm -3 Intensity N e = cm h! (ev) FIGURE 8: Spectra of Xe (39,683 states) at T e =50 ev and N e =10 16 and cm -3. SUMMARY AND CONCLUSION This new version of HULLAC is much more reliable, complete, and easy to use than the previous ones. It was thoroughly checked by comparison to the LANL codes with excellent agreement. The new fitting procedure guarantees a trustworthy evaluation of the rates in the CRM solver. The latter has been shown to be robust and to give insight into 214
10 the atomic model. All the code is in Fortran mostly 77 and 95 compatible, and it was successfully compiled on many platforms. It is now ready for distribution. ACKNOWLEDGEMENTS We thank Dr. S. Obenschain, Head of the Laser Plasma Branch (LPB) of the Naval Research Laboratory for his encouragements and support. This work was done under a contract with LPB through a grant from the USDOE. We are grateful to Prof. R Mancini and T. Durmaz for help with the LANL codes. REFERENCES 1. A. Bar-Shalom, M Klapisch, and J Oreg, J. Quant. Spectros. Rad. Transf. 71, (2001). 2. M Klapisch, Comput. Phys. Comm. 2, (1971); E. Koenig, Physica(Utrecht) 62, 393 (1972). 3. A. Bar-Shalom and M. Klapisch, Comput. Phys. Comm. 50, (1988). 4. A. Bar-Shalom, M Klapisch, and J. Oreg, Phys. Rev. A. 38, (1988). 5. J. Oreg, W. H. Goldstein, M Klapisch et al., Phys. Rev. A 44, 1750 (1991). 6. A. Bar-Shalom, M Klapisch, and J Oreg, presented at the 13 th Int'l Conference on Atomic Processes in Plasmas, Gatlinburg, TN, (2002), D R Schultz, F W Meyer, and F Ownby Eds., AIP Conference Proceedings, 635, A. Bar-Shalom, J. Oreg, W. H. Goldstein et al., Phys. Rev. A. 40, (1989). 8. L B Golden and D H Sampson, Astroph. J. Supp. 38, 19 (1978); L B Golden, R E H Clark, S J Goett et al., Astroph. J. Supp. 45, (1981); S J Goett, R E H Clark, and D H Sampson, Atomic Data Nuclear Data Tables 25, 185 (1980). 9. M Abramowitz and I A Stegun, Handbook of Mathematical Functions. (Dover Publications, New York, 1965). 10. M Busquet, High Ener. Dens. Phys. doi: /j.hedp (2007). 11. A Burgess and J A Tully, Astronom. & Astrophys. 254, (1992); A Burgess, M. C. Chidichimo, and J A Tully, Astronom. & Astrophys. Supp. 131, (1998). 12. Y Ralchenko, R W Lee, and C Bowen, presented at the 14th APS Topical conference on Atomic Processes in Plasmas, Santa Fe, NM, (2004), J. S. Cohen, S. Mazevet, and D. P. Kilcrease Eds., AIP conference proceedings, 730, A Bar-Shalom, J Oreg, and M Klapisch, J. Quant. Spectros. Rad. Transf. 65, 43 (2000). 14. William H Press, Brian P Flannery, Saul A Teukolsky et al., Numerical Recipes in Fortran 77 :The art of scientific computing., 2nd ed. (Cambridge University Press, Cambridge, UK, 1996). 15. S. D. Loch, C. J. Fontes, J. Colgan et al., Phys. Rev. E 69, (2004); M Busquet, J P Raucourt, and J C Gauthier, J. Quant. Spectros. Rad. Transf. 54, (1995). 16. M. Klapisch and M. Busquet, High Ener. Dens. Phys. 3, (2007). 215
A new algorithm for the solution and analysis of the Collisional Radiative Model equations
High Energy Density Physics 3 (2007) 143e148 www.elsevier.com/locate/hedp A new algorithm for the solution and analysis of the Collisional Radiative Model equations Marcel Klapisch a, *, Michel Busquet
More informationConfiguration interaction effect on open M shell Fe and Ni LTE spectral opacities, Rosseland and Planck means
Journal of Physics: Conference Series PAPER OPEN ACCESS Configuration interaction effect on open M shell Fe and Ni LTE spectral opacities, Rosseland and Planck means To cite this article: D Gilles et al
More informationUncertainties on atomic data
Uncertainties on atomic data Connor Ballance (Auburn University) Stuart Loch (Auburn University) Mike Witthoeft & Tim Kallman (NASA-Goddard) Adam Foster & Randall Smith (CfA,Harvard) Joint ITAMP-IAEA Technical
More informationThe LANL atomic kinetics modeling effort and its application to W plasmas
The LANL atomic kinetics modeling effort and its application to W plasmas James Colgan, Joseph Abdallah, Jr., Christopher Fontes, Honglin Zhang Los Alamos National Laboratory IAEA CRP December 2010 jcolgan@lanl.gov
More informationThe validity of classical trajectory and perturbative quantal methods for electron-impact ionization from excited states in H-like ions
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS B: ATOMIC, MOLECULAR AND OPTICAL PHYSICS J. Phys. B: At. Mol. Opt. Phys. 38 (25) L199 L26 doi:1.188/953-475/38/12/l1 LETTER TO THE EDITOR The validity
More informationINTERNATIONAL BULLETIN ON ATOMIC AND MOLECULAR DATA FOR FUSION
INTERNATIONAL BULLETIN ON ATOMIC AND MOLECULAR DATA FOR FUSION Number 66 August 2007 Contributors: M. E. Bannister, J. Bretagne, J. Fuhr, H. B. Gilbody, C. C. Havener, T. Kato, A. Kramida, P. S. Krstic,
More informationAtomic Data for Astrophysics. II. New Analytic Fits for Photoionization Cross Sections of Atoms and Ions
Atomic Data for Astrophysics. II. New Analytic Fits for Photoionization Cross Sections of Atoms and Ions D. A. Verner 1, G. J. Ferland, and K. T. Korista Department of Physics and Astronomy, University
More informationA PROGRAM TO COMPUTE EXACT HYDROGENIC OSCILLATOR STRENGTHS, EINSTEIN COEFFICIENTS, AND LIFETIMES
A PROGRAM TO COMPUTE EXACT HYDROGENIC OSCILLATOR STRENGTHS, EINSTEIN COEFFICIENTS, AND LIFETIMES D. Hoang Binh LERMA, BAT 18, Observatoire de Paris Meudon, F 92195 Meudon Principal Cedex, France Abstract
More informationEffective collision strengths for fine-structure transitions for the electron impact excitation of N II. C. E. Hudson and K. L.
A&A 43, 725 729 (25) DOI: 1.151/4-6361:241969 c ESO 25 Astronomy & Astrophysics Effective collision strengths for fine-structure transitions for the electron impact excitation of N II C. E. Hudson and
More informationAtomic data from the IRON project ABSTRACT
A&A 446, 361 366 (26) DOI: 1.151/4-6361:253631 c ESO 26 Astronomy & Astrophysics Atomic data from the IRON project LX. Electron-impact excitation of n = 3, 4 levels of Fe 17+ M. C. Witthoeft 1,N.R.Badnell
More informationR. Clark, D. Humbert, K. Sheikh Nuclear Data Section
Calculation of Atomic Data for Plasma Modeling: Introduction and Atomic Structure Part 1 R. Clark, D. Humbert, K. Sheikh Nuclear Data Section Overview Plasmas in fusion research Data needs for plasma modeling
More informationEUV spectra from the NIST EBIT
EUV spectra from the NIST EBIT D. Kilbane and G. O Sullivan Atomic and Molecular Plasma Physics group, UCD, Ireland J. D. Gillaspy, Yu. Ralchenko and J. Reader National Institute of Standards and Technology,
More informationSolutions of the Schrödinger equation for an attractive 1/r 6 potential
PHYSICAL REVIEW A VOLUME 58, NUMBER 3 SEPTEMBER 998 Solutions of the Schrödinger equation for an attractive /r 6 potential Bo Gao Department of Physics and Astronomy, University of Toledo, Toledo, Ohio
More informationNumerical Modeling of Radiative Kinetic Plasmas
2014 US-Japan JIFT Workshop on Progress in kinetic plasma simulations Oct.31-Nov.1, 2014, Salon F, New Orleans Marriott, New Orleans, LA, U.S.A Numerical Modeling of Radiative Kinetic Plasmas T. Johzaki
More informationTheoretical approaches to electronimpact
Theoretical approaches to electronimpact ionization James Colgan, C. J. Fontes, D. P. Kilcrease, J. Abdallah, Jr., M. E. Sherrill, P. Hakel, and M. S. Pindzola 2 Los Alamos National Laboratory, Los Alamos,
More informationABSTRACT Electron collision excitation strengths of inelastic transitions among the 52 Ðne-structure levels of the
THE ASTROPHYSICAL JOURNAL, 530:1091È1104, 2000 February 20 ( 2000. The American Astronomical Society. All rights reserved. Printed in U.S.A. ELECTRON COLLISION EXCITATION OF FINE-STRUCTURE LEVELS IN S
More informationTransition Probabilities for the Dipole Allowed Fine Structure Transitions in Sll
Physica Scripta., Vol.55, 200-238, 1997 ISSN: 1402-4896/55/2/012 (print) doi: 10.1088/0031-8949/55/2/012 IOP Publishing http://journals.iop.org/ http://www.iop.org/ej/journal/physscr Transition Probabilities
More informationTime-dependent kinetics model for a helium discharge plasma
J. Phys. B: At. Mol. Opt. Phys. 32 (1999) 1001 1008. Printed in the UK PII: S0953-4075(99)97893-8 Time-dependent kinetics model for a helium discharge plasma J Abdallah Jr, N Palmer, W Gekelman, J Maggs
More informationAtomic data from the IRON Project
A&A 372, 083 087 (200) DOI: 05/0004-636:2000562 c ESO 200 Astronomy & Astrophysics Atomic data from the IRON Project L. Electron impact excitation of Fe XIX K. Butler and C. J. Zeippen 2 Institut für Astronomie
More informationCalculating Radiative Recombination Continuum From a Hot Plasma
Calculating Radiative Recombination Continuum From a Hot Plasma Randall Smith August 2, 2008 When an electron collides with an atom or ion, it may excite the atom/ion (I + e I + e) ionize the atom/ion
More informationThe CHIANTI Atomic Database
The CHIANTI Atomic Database An Overview of Data, Software and Applications Dr Peter Young George Mason University, USA NASA Goddard Space Flight Center, USA Overview 1. Quick guide 2. History of project
More informationEmpirical formula for cross section of direct electron-impact ionization of ions
J. Phys. B: At. Mol. Opt. Phys. 33 (2000) 5025 5032. Printed in the UK PII: S0953-4075(00)16500-9 Empirical formula for cross section of direct electron-impact ionization of ions V A Bernshtam, Yu V Ralchenko
More informationNIST Research on Spectroscopy and Collisional-Radiative Modeling of Highly-Charged Ions of Tungsten
NIST Research on Spectroscopy and Collisional-Radiative Modeling of Highly-Charged Ions of Tungsten Yuri Ralchenko National Institute of Standards and Technology Gaithersburg, USA Vienna, Austria, Dec
More informationDetailed diagnostics for a hot bromine plasma by the open M-shell opacity
PHYSICAL REVIEW E 72, 016404 2005 Detailed diagnostics for a hot bromine plasma by the open M-shell opacity Fengtao Jin and Jianmin Yuan* Department of Physics, National University of Defense Technology,
More informationCHIANTI An atomic database for astrophysical plasmas
CHIANTI An atomic database for astrophysical plasmas Enrico Landi University of Michigan On behalf of the CHIANTI team: Enrico Landi Ken Dere Peter Young Giulio Del Zanna Helen Mason University of Michigan,
More informationAtomic structure and dynamics
Atomic structure and dynamics -- need and requirements for accurate atomic calculations Analysis and interpretation of optical and x-ray spectra (astro physics) Isotope shifts and hyperfine structures
More informationA New Computational Method for non-lte, the Linear Response Matrix
UCRL-JC-3407 Rev PREPRINT A New Computational Method for non-lte, the Linear Response Matrix J. A. Harte, R. M. More, G. B. Zimmerman, S. B. Libby, F. R. Graziani, K. B. Fournier This paper was prepared
More informationJoint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012
2327-4 Joint ICTP- Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Atomic Processes Modeling in Plasmas Modeling Spectroscopic Observables from Plasmas Hyun-Kyung
More informationExtreme-ultraviolet emissivity from Xe 8+ to Xe 12+ by using a detailed line-by-line method
Eur. Phys. J. D (21) DOI: 1.114/epjd/e21-192-6 Regular Article THE EUROPEAN PHYSICAL JOURNAL D Extreme-ultraviolet emissivity from Xe 8+ to Xe 12+ by using a detailed line-by-line method J. Zeng a, C.
More informationAnalysis of the x-ray spectrum emitted by laser-produced plasma of dysprosium
Marcus et al. Vol. 24, No. 5/ May 2007/J. Opt. Soc. Am. B 1187 Analysis of the x-ray spectrum emitted by laser-produced plasma of dysprosium Gilad Marcus, Einat Louzon, Zohar Henis, and Shlomo Maman Soreq
More informationOne and Two Photon Ionization along the Fe Isonuclear Sequence
I R A M P 8(), December 017, pp. 81-91 One and Two Photon Ionization along the Fe Isonuclear Sequence International Science Press ISSN: 9-3159 One and Two Photon Ionization along the Fe Isonuclear Sequence
More information1. Introduction ABSTRACT
Determination of the average ionization and thermodynamic regimes of xenon plasmas with an application to the characterization of blast waves launched in xenon clusters R. Rodriguez a ' b '*, J.M. Gil
More information7. Non-LTE basic concepts
7. Non-LTE basic concepts LTE vs NLTE occupation numbers rate equation transition probabilities: collisional and radiative examples: hot stars, A supergiants 10/13/2003 Spring 2016 LTE LTE vs NLTE each
More informationAtomic Structure and Processes
Chapter 5 Atomic Structure and Processes 5.1 Elementary atomic structure Bohr Orbits correspond to principal quantum number n. Hydrogen atom energy levels where the Rydberg energy is R y = m e ( e E n
More informationATOMIC AND PLASMA-MATERIAL INTERACTION DATA FOR FUSION VOLUME 11
ATOMIC AND PLASMA-MATERIAL INTERACTION DATA FOR FUSION VOLUME 11 INTERNATIONAL ATOMIC ENERGY AGENCY, VIENNA, 2003 The volumes of ATOMIC AND PLASMA MATERIAL INTERACTION DATA FOR FUSION are published by
More informationAtomic data from the IRON Project
ASTRONOMY & ASTROPHYSICS JANUARY II 000, PAGE 85 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 141, 85 96 (000) Atomic data from the IRON Project XL. Electron impact excitation of the Fe XIV EUV transitions
More informationComputations of Scattering Cross Sections for He, Ne, Ar, Kr, Xe and Rn
Computations of Scattering Cross Sections for He, Ne, Ar, Kr, Xe and Rn Full Length Research Article * ABDU, S. G. 1 1 Department of Physics, Kaduna State University, Kaduna-Nigeria *(sgabdul@kasu.edu.ng)
More informationDielectronic recombination data for dynamic finite-density plasmas. IV. The carbon isoelectronic sequence
A&A 417, 1173 1181 (2004) DOI: 10.1051/0004-6361:20034174 c ESO 2004 Astronomy & Astrophysics Dielectronic recombination data for dynamic finite-density plasmas IV. The carbon isoelectronic sequence O.
More informationarxiv: v1 [physics.atom-ph] 2 Dec 2015
J. Phys. B: At. Mol. Opt. Phys. arxiv:1512.657v1 [physics.atom-ph] 2 Dec 215 Theoretical investigation of spectroscopic properties of W 26+ in EBIT plasma V. Jonauskas, A. Kynienė, P. Rynkun, S. Kučas,
More informationWhy write proofs? Why not just test and repeat enough examples to confirm a theory?
P R E F A C E T O T H E S T U D E N T Welcome to the study of mathematical reasoning. The authors know that many students approach this material with some apprehension and uncertainty. Some students feel
More informationOptimization Methods
Optimization Methods Decision making Examples: determining which ingredients and in what quantities to add to a mixture being made so that it will meet specifications on its composition allocating available
More informationThe Extreme Ultraviolet Emissions of W 23+ (4f 5 )
The Extreme Ultraviolet Emissions of W 23+ (4f 5 ) T Pütterich, V. Jonauskas, R Neu, R Dux and ASDEX Upgrade Team Max Planck Institut für Plasmaphysik, EURATOM Association, D 85748 Garching, Germany Institute
More informationIso-nuclear tungsten dielectronic recombination rates for use in magnetically-confined fusion plasmas
Iso-nuclear tungsten dielectronic recombination rates for use in magnetically-confined fusion plasmas D.-H. Kwon a, W. Lee a, S. Preval b, C. P. Ballance c,, E. Behar d, J. Colgan e, C. J. Fontes e, T.
More informationPHYS852 Quantum Mechanics II, Spring 2010 HOMEWORK ASSIGNMENT 8: Solutions. Topics covered: hydrogen fine structure
PHYS85 Quantum Mechanics II, Spring HOMEWORK ASSIGNMENT 8: Solutions Topics covered: hydrogen fine structure. [ pts] Let the Hamiltonian H depend on the parameter λ, so that H = H(λ). The eigenstates and
More informationDiode Lasers and Photonic Integrated Circuits
Diode Lasers and Photonic Integrated Circuits L. A. COLDREN S. W. CORZINE University of California Santa Barbara, California A WILEY-INTERSCIENCE PUBLICATION JOHN WILEY & SONS, INC. NEW YORK / CHICHESTER
More informationChapter 6. Atomic Physics and Process in Partially Ionized Plasma
Chapter 6 Atomic Physics and Process in Partially Ionized Plasma 6.1 Fundamentals of Atomic Physics Hydrogen Atom E 2 1 1 2 2 2m a n Photon(Radiation) Emission from H Atom Opacity Photon Energy Energy
More informationFollow links Class Use and other Permissions. For more information, send to:
COPYRIGHT NOTICE: Stephen L. Campbell & Richard Haberman: Introduction to Differential Equations with Dynamical Systems is published by Princeton University Press and copyrighted, 2008, by Princeton University
More informationPlasma EUV source has been studied to achieve 180W of power at λ=13.5nm, which is required for the next generation microlithography
Acknowledgement K. Nishihara, H. Nishimura, S. Fujioka Institute for Laser Engineering, Osaka University A. Sunahara, H. Furukawa Institute for Laser Technology T. Nishikawa, Okayama University F. Koike,
More informationAccelerator Physics. Tip World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI HONG KONG TAIPEI BANGALORE. Second Edition. S. Y.
Accelerator Physics Second Edition S. Y. Lee Department of Physics, Indiana University Tip World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI HONG KONG TAIPEI BANGALORE Contents Preface Preface
More informationCONTRIBUTION FUNCTION OF MOLYBDENUM ATOM AND IONS IN ASTROPHYSICAL AND LABORATORY PLASMA AS A FUNCTION OF ELECTRON TEMPERATURE
COTRIBUTIO FUCTIO OF MOLYBDEUM ATOM AD IOS I ASTROPHYSICAL AD LABORATORY PLASMA AS A FUCTIO OF ELECTRO TEMPERATURE A.. Jadhav Department of Electronics, Yeshwant Mahavidyalaya, anded. Affiliated to Swami
More informationExcitation of the Fe XIV spectrum in the Sun, stars and Seyfert galaxies: reconciling theory, observations and experiment.
Excitation of the Fe XIV spectrum in the Sun, stars and Seyfert galaxies: reconciling theory, observations and experiment Peter Storey, UCL 1 Introduction Observational background Modelling Seyfert spectra
More informationJ S Parker (QUB), Martin Plummer (STFC), H W van der Hart (QUB) Version 1.0, September 29, 2015
Report on ecse project Performance enhancement in R-matrix with time-dependence (RMT) codes in preparation for application to circular polarised light fields J S Parker (QUB), Martin Plummer (STFC), H
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationStudy of Electron Energy and Angular Distributions and Calculations of X-ray, EUV Line Flux and Rise Times
J. Astrophys. Astr. (1987) 8, 263 270 Study of Electron Energy and Angular Distributions and Calculations of X-ray, EUV Line Flux and Rise Times Ranjna Bakaya, Sunil Peshin, R. R. Rausaria & P. N. Khosa
More informationATOMIC AND LASER SPECTROSCOPY
ALAN CORNEY ATOMIC AND LASER SPECTROSCOPY CLARENDON PRESS OXFORD 1977 Contents 1. INTRODUCTION 1.1. Planck's radiation law. 1 1.2. The photoelectric effect 4 1.3. Early atomic spectroscopy 5 1.4. The postulates
More informationCHIANTI: an atomic database for emission lines
Astron. Astrophys. 329, 291 314 (1998) ASTRONOMY AND ASTROPHYSICS CHIANTI: an atomic database for emission lines II. Comparison with the SERTS-89 active region spectrum P.R. Young 1, E. Landi 2, and R.J.
More informationOscillator strengths and E1 radiative rates for Ca-like titanium, Ti III
Int. J. New. Hor. Phys. 2, No. 1, 25-31 (2015) 25 International Journal of New Horizons in Physics http://dx.doi.org/10.12785/ijnhp/020105 Oscillator strengths and E1 radiative rates for Ca-like titanium,
More informationCONTENTS. vii. CHAPTER 2 Operators 15
CHAPTER 1 Why Quantum Mechanics? 1 1.1 Newtonian Mechanics and Classical Electromagnetism 1 (a) Newtonian Mechanics 1 (b) Electromagnetism 2 1.2 Black Body Radiation 3 1.3 The Heat Capacity of Solids and
More informationNeutron Generation from 10MeV Electron Beam to Produce Mo99
From the SelectedWorks of Innovative Research Publications IRP India Winter January 1, 2015 Neutron Generation from 10MeV Electron Beam to Produce Mo99 Innovative Research Publications, IRP India, Innovative
More informationOPSIAL Manual. v Xiaofeng Tan. All Rights Reserved
OPSIAL Manual v1.0 2016 Xiaofeng Tan. All Rights Reserved 1. Introduction... 3 1.1 Spectral Calculator & Fitter (SCF)... 3 1.2 Automated Analyzer (AA)... 3 2. Working Principles and Workflows of OPSIAL...
More informationNuclear Effects in Electron Capture into Highly Charged Heavy Ions
Nuclear Effects in Electron Capture into Highly Charged Heavy Ions W. Scheid 1,A.Pálffy 2,Z.Harman 2, C. Kozhuharov 3, and C. Brandau 3 1 Institut für Theoretische Physik der Justus-Liebig-Universität
More informationLINE INTENSITY RATIOS IN THE EIS RANGE SENSITIVE TO ELECTRON DENSITIES IN 10 7 K PLASMAS
The Astrophysical Journal, 679:843Y847, 2008 May 20 # 2008. The American Astronomical Society. All rights reserved. Printed in U.S.A. LINE INTENSITY RATIOS IN THE EIS RANGE SENSITIVE TO ELECTRON DENSITIES
More informationSolving the Time-Dependent Schrödinger Equation: Instructor Notes a. Abstract
Solving the Time-Dependent Schrödinger Equation: Instructor Notes a David G. Robertson Department of Physics, Otterbein University, Westerville, OH 43081 (Dated: October 10, 2011) Abstract Some notes on
More informationVibrationally resolved ion-molecule collisions
Vibrationally resolved ion-molecule collisions CRP: Atomic and Molecular Data for State-Resolved Modelling of Hydrogen and Helium and Their Isotopes in Fusion plasma Predrag Krstic Physics Division, Oak
More informationProper initial conditions for long-term integrations of the solar system
Proper initial conditions for long-term integrations of the solar system M. Arminjon Laboratoire Sols, Solides, Structures [CNRS / Université Joseph Fourier / Institut National Polytechnique de Grenoble]
More informationA. N. Jadhav 1. Yeshwant Mahavidyalaya, Nanded. Maharashtra, India.
Departure from Maxwellian Velocity Distribution of High Temperature Plasma in Solar Flares where Almost All Electrons have been Stripped out from Iron Atoms A. N. Jadhav 1 1 Department of Electronics,
More informationGA A25842 STUDY OF NON-LTE SPECTRA DEPENDENCE ON TARGET MASS IN SHORT PULSE LASER EXPERIMENTS
GA A25842 STUDY OF NON-LTE SPECTRA DEPENDENCE ON TARGET MASS IN SHORT PULSE LASER EXPERIMENTS by C.A. BACK, P. AUDBERT, S.D. BATON, S.BASTIANI-CECCOTTI, P. GUILLOU, L. LECHERBOURG, B. BARBREL, E. GAUCI,
More informationElectron impact ionization of diatomic molecules
Eur. Phys. J. D 8, 5 5 (8) DOI:./epjd/e8-- Electron impact ionization of diatomic molecules I. Tóth, R.I. Campeanu, V. Chiş and L. Nagy Eur. Phys. J. D 8, 5 5 (8) DOI:./epjd/e8-- THE EUROPEAN PHYSICAL
More informationIntroduction to Modern Quantum Optics
Introduction to Modern Quantum Optics Jin-Sheng Peng Gao-Xiang Li Huazhong Normal University, China Vfe World Scientific» Singapore* * NewJerseyL Jersey* London* Hong Kong IX CONTENTS Preface PART I. Theory
More informationAnalytic Benchmarking of the 2DX eigenvalue code
Analytic Benchmarking of the 2DX eigenvalue code D. A. Baver, J. R. Myra Lodestar Research Corporation M. Umansky Lawrence Livermore National Laboratory Analytic benchmarking of the 2DX eigenvalue code
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationTheoretical Analysis of Neutron Double-Differential Cross Section of n + 19 F at 14.2 MeV
Commun. Theor. Phys. (Beijing, China) 47 (2007) pp. 102 106 c International Academic Publishers Vol. 47, No. 1, January 15, 2007 Theoretical Analysis of Neutron Double-Differential Cross Section of n +
More informationPlasma shielding during ITER disruptions
Plasma shielding during ITER disruptions Sergey Pestchanyi and Richard Pitts 1 Integrated tokamak code TOKES is a workshop with various tools objects Radiation bremsstrahlung recombination s line s cyclotron
More informationInteractive system for the interpretation of atomic spectra
J. Straus, K. Kolacek, V. Bohacek, O. Frolov, V. Prukner, M. Ripa, V. Sember, J. Schmidt, P. Vrba Institute of Plasma Physics AS CR, Za Slovankou 3, P.O. Box 17, 182 21 Prague 8 D. Klir Faculty of Electrical
More informationSTEAMEST: A Software Tool for Estimation of Physical Properties of Water and Steam
226 JOURNAL OF SOFTWARE, VOL. 4, NO. 3, MAY 2009 STEAMEST: A Software Tool for Estimation of Physical Properties of Water and Steam Muhammad Faheem Department of Chemical Engineering, University of Engineering
More informationxii ASTROPHYSICAL BACKGROUND
Preface The purpose of this book is to provide an account of the major advances achieved in making quantitative spectroscopic analyses of the observable outer layers of stars. These analyses can now be
More informationDual Laser Plasma Photoabsorption Studies Of Gadolinium In The Extreme Ultraviolet Region
Dual Laser Plasma Photoabsorption Studies Of Gadolinium In The Extreme Ultraviolet Region Paddy Hayden National Centre for Plasma Science & Technology (NCPST) / School of Physical Sciences, Dublin City
More informationUse of intermediate coupling relationships to test measured branching fraction data
J. Phys. B: At. Mol. Opt. Phys. 31 (1998) L769 L774. Printed in the UK PII: S0953-4075(98)95845-X LETTER TO THE EDITOR Use of intermediate coupling relationships to test measured branching fraction data
More informationTHE OBSERVATION AND ANALYSIS OF STELLAR PHOTOSPHERES
THE OBSERVATION AND ANALYSIS OF STELLAR PHOTOSPHERES DAVID F. GRAY University of Western Ontario, London, Ontario, Canada CAMBRIDGE UNIVERSITY PRESS Contents Preface to the first edition Preface to the
More informationHe + and Xe + velocities near the presheath-sheath boundary in a He/Xe discharge
University of Iceland Faculty of Physical Sciences He + and Xe + velocities near the presheath-sheath boundary in a He/Xe discharge by Davíð Ingvi Snorrason 16 ECTS thesis submitted in partial fulfillment
More informationAtomic Physics 3 rd year B1
Atomic Physics 3 rd year B1 P. Ewart Lecture notes Lecture slides Problem sets All available on Physics web site: http:www.physics.ox.ac.uk/users/ewart/index.htm Atomic Physics: Astrophysics Plasma Physics
More informationNUMERICALLY COMPUTING THE LYAPUNOV EXPONENTS OF MATRIX-VALUED COCYCLES
NUMERICALLY COMPUTING THE LYAPUNOV EXPONENTS OF MATRIX-VALUED COCYCLES RODRIGO TREVIÑO This short note is based on a talk I gave at the student dynamical systems seminar about using your computer to figure
More informationarxiv:quant-ph/ v1 20 Apr 1995
Combinatorial Computation of Clebsch-Gordan Coefficients Klaus Schertler and Markus H. Thoma Institut für Theoretische Physik, Universität Giessen, 3539 Giessen, Germany (February, 008 The addition of
More informationA new class of sum rules for products of Bessel functions. G. Bevilacqua, 1, a) V. Biancalana, 1 Y. Dancheva, 1 L. Moi, 1 2, b) and T.
A new class of sum rules for products of Bessel functions G. Bevilacqua,, a V. Biancalana, Y. Dancheva, L. Moi,, b and T. Mansour CNISM and Dipartimento di Fisica, Università di Siena, Via Roma 56, 5300
More informationADAS404: Isonuclear master data - extract from isoelectronic master data
ADAS404: Isonuclear master data - extract from isoelectronic master data The program implements a conversion of collisional-dielectronic data for isoelectronic sequences of ADAS data format adf0 (iso-electronic
More informationA Collisional Radiative Model for Caesium and Its Application to an RF Source for Negative Hydrogen Ions
A Collisional Radiative Model for Caesium and Its Application to an RF Source for Negative Hydrogen Ions D. Wünderlich 1,a), C. Wimmer 1) and R. Friedl 2) 1 Max-Planck-Institut für Plasmaphysik, Boltzmannstr.
More informationThe Performance Evolution of the Parallel Ocean Program on the Cray X1
The Performance Evolution of the Parallel Ocean Program on the Cray X1 Patrick H. Worley Oak Ridge National Laboratory John Levesque Cray Inc. 46th Cray User Group Conference May 18, 2003 Knoxville Marriott
More informationxvi xxiii xxvi Construction of the Real Line 2 Is Every Real Number Rational? 3 Problems Algebra of the Real Numbers 7
About the Author v Preface to the Instructor xvi WileyPLUS xxii Acknowledgments xxiii Preface to the Student xxvi 1 The Real Numbers 1 1.1 The Real Line 2 Construction of the Real Line 2 Is Every Real
More informationP. W. Atkins and R. S. Friedman. Molecular Quantum Mechanics THIRD EDITION
P. W. Atkins and R. S. Friedman Molecular Quantum Mechanics THIRD EDITION Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1997 Introduction and orientation 1 Black-body radiation 1 Heat capacities 2 The
More informationAtomic and Molecular Databases and Data Evaluation Activities at NIFS
Joint IAEA-NFRI Technical Meeting on Data Evaluation for Atomic, Molecular and Plasma Material Interaction Processes in Fusion Daejon, 4-7 Sep., 2012 Atomic and Molecular Databases and Data Evaluation
More informationNotes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Multidimensional Unconstrained Optimization Suppose we have a function f() of more than one
More informationHS AP Physics 1 Science
Scope And Sequence Timeframe Unit Instructional Topics 5 Day(s) 20 Day(s) 5 Day(s) Kinematics Course AP Physics 1 is an introductory first-year, algebra-based, college level course for the student interested
More informationarxiv: v1 [astro-ph.sr] 6 Jan 2011
Theoretical and experimental activities on opacities for a good interpretation of seismic stellar probes arxiv:1101.1182v1 [astro-ph.sr] 6 Jan 2011 S Turck-Chièze 1, G Loisel 1,2, D Gilles 1, F Thais 2,
More informationIMPLEMENTATION OF THE MONTE CARLO-LIBRARY LEAST- SQUARES APPROACH TO ENERGY DISPERSIVE X-RAY FLUORESCENCE ANALYSIS
227 IMPLEMENTATION OF THE MONTE CARLO-LIBRARY LEAST- SQUARES APPROACH TO ENERGY DISPERSIVE X-RAY FLUORESCENCE ANALYSIS Fusheng Li, Weijun Guo, and Robin P. Gardner Center for Engineering Applications of
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 5.76 Modern Topics in Physical Chemistry Spring, Problem Set #2
Reading Assignment: Bernath Chapter 5 MASSACHUSETTS INSTITUTE O TECHNOLOGY 5.76 Modern Topics in Physical Chemistry Spring 994 Problem Set # The following handouts also contain useful information: C &
More informationPhotoionization of Fe
J. Phys. B: At Mol. Opt. Phys. 26 (1993) LlLL6. Printed in the UK LETIXR TO THE EDITOR Photoionization of Fe Maryvonne Le Doumeuf?, Sultana N NaharS and Ani1 K F radhans t Laboratoire de ThCorie des Processus
More informationNew theoretical insights on the physics of compound nuclei from laser-nucleus reactions
New theoretical insights on the physics of compound nuclei from laser-nucleus reactions Adriana Pálffy Max Planck Institute for Nuclear Physics, Heidelberg, Germany Laser-Driven Radiation Sources for Nuclear
More informationComputer Simulation Using Particles
Computer Simulation Using Particles R W Hockney Emeritus Professor, University ofreading and J W Eastwood Culham Laboratory, Abingdon Institute of Physics Publishing Bristol and Philadelphia CONTENTS Foreword
More informationUNIVERSITY OF CAMBRIDGE
UNIVERSITY OF CAMBRIDGE DOWNING COLLEGE MATHEMATICS FOR ECONOMISTS WORKBOOK This workbook is intended for students coming to Downing College Cambridge to study Economics 2018/ 19 1 Introduction Mathematics
More information