Calculation of Tin Atomic Data and Plasma Properties

Similar documents
k p theory for bulk semiconductors

LECTURE 21 Mohr s Method for Calculation of General Displacements. 1 The Reciprocal Theorem

22.51 Quantum Theory of Radiation Interactions

Research on Complex Networks Control Based on Fuzzy Integral Sliding Theory

MARKOV CHAIN AND HIDDEN MARKOV MODEL

Chapter 6. Rotations and Tensors

( ) r! t. Equation (1.1) is the result of the following two definitions. First, the bracket is by definition a scalar product.

Note 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2

G : Statistical Mechanics

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Neural network-based athletics performance prediction optimization model applied research

Numerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes

Quantum Runge-Lenz Vector and the Hydrogen Atom, the hidden SO(4) symmetry

Image Classification Using EM And JE algorithms

Associative Memories

Lowest-Order e + e l + l Processes in Quantum Electrodynamics. Sanha Cheong

Atomic Scattering Factor for a Spherical Wave and the Near Field Effects in X-ray Fluorescence Holography

Structure and Drive Paul A. Jensen Copyright July 20, 2003

3. Stress-strain relationships of a composite layer

Predicting Model of Traffic Volume Based on Grey-Markov

Module III, Lecture 02: Correlation Functions and Spectral Densities

A finite difference method for heat equation in the unbounded domain

Semiclassical theory of molecular nonlinear optical polarization

Title: Radiative transitions and spectral broadening

Development of whole CORe Thermal Hydraulic analysis code CORTH Pan JunJie, Tang QiFen, Chai XiaoMing, Lu Wei, Liu Dong

Nested case-control and case-cohort studies

The Application of BP Neural Network principal component analysis in the Forecasting the Road Traffic Accident

Cyclic Codes BCH Codes

SUPPLEMENTARY INFORMATION

Multispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory

Calculation method of electrical conductivity, thermal conductivity and viscosity of a partially ionized gas. Ilona Lázniková

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

arxiv: v1 [physics.comp-ph] 17 Dec 2018

The line method combined with spectral chebyshev for space-time fractional diffusion equation

QUARTERLY OF APPLIED MATHEMATICS

Inductance Calculation for Conductors of Arbitrary Shape

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.

arxiv: v3 [cond-mat.str-el] 15 Oct 2009

ON AUTOMATIC CONTINUITY OF DERIVATIONS FOR BANACH ALGEBRAS WITH INVOLUTION

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

COXREG. Estimation (1)

Strain Energy in Linear Elastic Solids

COMBINING SPATIAL COMPONENTS IN SEISMIC DESIGN

D hh ν. Four-body charm semileptonic decay. Jim Wiss University of Illinois

Supplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Example: Suppose we want to build a classifier that recognizes WebPages of graduate students.

The Feynman path integral

Applied Nuclear Physics (Fall 2004) Lecture 23 (12/3/04) Nuclear Reactions: Energetics and Compound Nucleus

THEOREMS OF QUANTUM MECHANICS

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

5.76 Lecture #5 2/07/94 Page 1 of 10 pages. Lecture #5: Atoms: 1e and Alkali. centrifugal term ( +1)

MULTIVARIABLE FUZZY CONTROL WITH ITS APPLICATIONS IN MULTI EVAPORATOR REFRIGERATION SYSTEMS

EXPERIMENT AND THEORISATION: AN APPLICATION OF THE HYDROSTATIC EQUATION AND ARCHIMEDES THEOREM

NUMERICAL DIFFERENTIATION

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Numerical integration in more dimensions part 2. Remo Minero

MODEL TUNING WITH THE USE OF HEURISTIC-FREE GMDH (GROUP METHOD OF DATA HANDLING) NETWORKS

STAT 511 FINAL EXAM NAME Spring 2001

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Physics 2D Lecture Slides Lecture 30: Mar 12th

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Absorption and Recurrence Spectra of Li Rydberg Atom in Perpendicular Electric and Magnetic Fields

On the Power Function of the Likelihood Ratio Test for MANOVA

An Effective Space Charge Solver. for DYNAMION Code

L-Edge Chromatic Number Of A Graph

Lecture 12: Discrete Laplacian

Canonical transformations

[WAVES] 1. Waves and wave forces. Definition of waves

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Probabilistic method to determine electron correlation energy

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Dynamic Analysis Of An Off-Road Vehicle Frame

CHAPTER 2. PHOTOABSORPTION AND PHOTOELECTRIC PROCESS IN Si QUANTUM DOTS

Delay tomography for large scale networks

The Order Relation and Trace Inequalities for. Hermitian Operators

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Application of Particle Swarm Optimization to Economic Dispatch Problem: Advantages and Disadvantages

Lecture 14: Forces and Stresses

A General Column Generation Algorithm Applied to System Reliability Optimization Problems

) is the unite step-function, which signifies that the second term of the right-hand side of the

Xin Li Department of Information Systems, College of Business, City University of Hong Kong, Hong Kong, CHINA

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Quantitative Evaluation Method of Each Generation Margin for Power System Planning

Linear Approximation with Regularization and Moving Least Squares

APPENDIX A Some Linear Algebra

CHAPTER 7 ENERGY BALANCES SYSTEM SYSTEM. * What is energy? * Forms of Energy. - Kinetic energy (KE) - Potential energy (PE) PE = mgz

MODELING OF MULTICOMPONET FUEL EFFECTS IN INTERNAL COMBUSTION ENGINES

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Andre Schneider P622

REAL-TIME IMPACT FORCE IDENTIFICATION OF CFRP LAMINATED PLATES USING SOUND WAVES

Kernel Methods and SVMs Extension

DISCRETE VARIABLE REPRESENTATION OF THE ANGULAR VARIABLES IN QUANTUM THREE-BODY SCATTERING DAVID CABALLERO

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Unified spin-wave theory for quantum spin systems with single-ion anisotropies

Transcription:

AN-ET-04/4 Cacuaton of Tn Atomc Data and Pasma Propertes prepared by Energy Technoogy Dvson Argonne Natona aboratory Argonne Natona aboratory s managed by The Unversty of Chcago for the U. S. Department of Energy

About Argonne Natona aboratory Argonne s managed by The Unversty of Chcago for the U.S. Department of Energy under contract W-3-09-Eng-38. The aboratory s man facty s outsde Chcago, at 9700 South Cass Avenue, Argonne, Inos 60439. For nformaton about Argonne and ts poneerng scence and technoogy programs, see www.an.gov. Avaabty of Ths Report Ths report s avaabe, at no cost, at http://www.ost.gov/brdge. It s aso avaabe on paper to U.S. Department of Energy and ts contractors, for a processng fee, from: U.S. Department of Energy Offce of Scentfc and Technca Informaton P.O. Box 6 Oak Rdge, TN 3783-006 phone (865) 576-840 fax (865) 576-578 reports@adons.ost.gov Dscamer Ths report was prepared as an account of work sponsored by an agency of the Unted States Government. Nether the Unted States Government nor any agency thereof, nor The Unversty of Chcago, nor any of ther empoyees or offcers, makes any warranty, express or mped, or assumes any ega abty or responsbty for the accuracy, competeness, or usefuness of any nformaton, apparatus, product, or process dscosed, or represents that ts use woud not nfrnge prvatey owned rghts. Reference heren to any specfc commerca product, process, or servce by trade name, trademark, manufacturer, or otherwse, does not necessary consttute or mpy ts endorsement, recommendaton, or favorng by the Unted States Government or any agency thereof. The vews and opnons of document authors expressed heren do not necessary state or refect those of the Unted States Government or any agency thereof, Argonne Natona aboratory, or The Unversty of Chcago.

AN-ET-04/4 Cacuaton of Tn Atomc Data and Pasma Propertes by V. Morozov, V. Tokach, and A. Hassanen Energy Technoogy Dvson, Argonne Natona aboratory September 004 Argonne Natona aboratory s managed by The Unversty of Chcago for the U. S. Department of Energy

CONTENTS ABSTRACT... INTRODUCTION.... HFS ENERGY EVES AND OTHER DATA.... ION STATES AND POPUATIONS OF ATOMIC EVES...4 3. THERMODYNAMIC PASMA PROPERTIES AND EQUATION OF STATE...7 4. DETAIED SET OPACITIES IN WIDE ENERGY RANGE...8 5. PANCK GROUP AVERAGED OPACTITES IN WIDE ENERGY RANGE...9 6. HF ENERGY EVES AND OTHER DATA FOR EUV EMITTING TIN ION SPECIES... 7. DETAIED TIN EUV TRANSITIONS...7 8. DETAIED TIN EUV OPACITIES... CONCUSION...3 ACKNOWEDGMENTS...4 REFERENCES...4 APPENDIX A...6 APPENDIX B...9 APPENDIX C...3

FIGURES Fg. : Concentratons of tn ons for varous denstes...5 Fg. : Cacuatons from equaton of state for tn ons at varous denstes...7 Fg. 3: Optca coeffcents for tn pasma at T = 6 ev and N = 0 6 cm -3...0 Fg. 4: Panck mean optca coeffcents for tn pasma...0 Fg. 5: Expermenta data vs. cacuaton for Sn VIII 4 p 7 6 d 4d 5 transton...9 Fg. 6: Cacuaton of tn EUV transtons by the HEIGHTS-ATOM code...0 Fg. 7: Cacuaton of tn EUV transtons by the Cowan code... Fg. 8: Tn pasma optca coeffcents n the EUV range...3 TABES Tabe : Concentraton (n percent) of xenon ons at 0 7 cm -3...6 Tabe : Concentraton (n percent) of tn ons at 0 7 cm -3...6 Tabe 3: Energes of Xe XI nner shes by SD HF and average term HF...3 Tabe 4: Energes of Sn VIII nner shes by SD HF and average term HF...3 Tabe 5: O I p 3 3s excted state energy eves...4 Tabe 6: Ar II 3p 4 3d excted state energy eves...4 Tabe 7: Sn VIII 4d 7 ground state energy eves...5 Tabe 8: Xe XI 4d 8 ground state energy eves (n cm - )...6 Tabe 9: Xe XI 4d 7 5p excted state energy eves (n cm - )...6 Tabe 0: Summary of the tn EUV transtons...

ABSTRACT Ths report revews the major methods and technques we use n generatng basc atomc and pasma propertes reevant to extreme utravoet (EUV) thography appcatons. The bass of the work s the cacuaton of the atomc energy eves, transtons probabtes, and other atomc data by varous methods, whch dffer n accuracy, competeness, and compcaton. ater on, we cacuate the popuatons of atomc eves and on states n pasmas by means of the coson-radaton equbrum (CRE) mode. The resuts of the CRE mode are used as nput to the thermodynamc functons, such as pressure and temperature from the nterna energy and densty (equaton of state), eectrc resstance, therma conducton, and other pasma propertes. In addton, optca coeffcents, such as emsson and absorpton coeffcents, are generated to resove a radaton transport equaton (RTE). The capabtes of our approach are demonstrated by generatng the requred atomc and pasma propertes for tn ons and pasma wthn the EUV regon near 3.5 nm. INTRODUCTION To meet the requrements of the Inte thography Roadmap goas for hgh voume manufacturng n the future [] and Internatona SEMATECH s EUV Source Program goas [], the EUV source s requred to have a power of 80-0 W at 3.5 nm (% bandwdth). Varous aser produced pasma (PP) and gas dscharge produced pasma (DPP) devces are under deveopment and nvestgaton by dfferent research groups. Both types of EUV sources have advantages and dsadvantages. At present, none of the current EUV sources can dever enough power eves demanded by commerca chp manufacturers. The effcency of generatng EUV radaton s the key factor n successfu deveopment of the source. The eadng EUV groups are usng a varety of scentfc and engneerng approaches to maxmze EUV brghtness from ther devces. A common technque s the ncrease of the optmzed overa converson effcency of a devce, because t mnmzes the requred nput power for a requred EUV output. Because many physca processes are nvoved, and many technca probems need to be soved, when optmzng a partcuar EUV devce, ony computer modeng can generate a compete pcture at a reasonabe pace. We are deveopng an ntegrated mode to smuate the envronment of the EUV source and optmze the output of the source. The mode descrbes the hydrodynamc and optca processes that occur n EUV devces. It takes nto account pasma evouton and magnetohydrodynamc (MHD) processes as we as photon radaton transport. It uses the tota varaton dmnshng scheme n the ax-fredrch formuaton for the descrpton of magnetc compresson and dffuson n a cyndrca geometry. Aso under deveopment are modes for opacty cacuatons: a cosona radaton equbrum mode, a sef-consstent fed mode wth Auger processes, and a nonstatonary knetc mode. Radaton transport for both contnuum and nes wth detaed spectra profes s taken nto account. The deveoped modes are beng ntegrated nto the HEIGHTS-EUV computer smuaton package [3, 4]. Beng sef-consstent,

the HEIGHTS-EUV package can generate a requred nformaton and s competey ndependent of the data from other externa packages. Furthermore, expermentay (or numercay) obtaned reabe data can be ncorporated nto HEIGHTS-EUV to ncrease the overa accuracy and effcency of the smuaton resuts. The focus of ths report s the major methods and technques we use n our HEIGHTS- ATOM code to generate basc atomc and pasma propertes. Based upon accuratey generated atomc energy eves, transton probabtes, and other atomc data, we can cacuate the popuatons of atomc eves and on states n pasma by means of the CRE mode. In turn, the resuts of the CRE mode are used as nput to thermodynamc functons, such as those for pressure and temperature from the nterna energy and densty (equaton of state), eectrc resstance, therma conducton, and other pasma propertes. Optca coeffcents are aso generated to resove the RTE. Combnng the generated propertes wth an approprate descrpton of an EUV devce n terms of the MHD boundary condtons, we are abe to smuate the dynamcs of the devce and ts tota radaton output. The accuracy and competeness of atomc data are key factors n the successfu numerca smuaton of the EUV devce. Ths report contans the resuts of our computer smuaton of the tn atomc data and pasma propertes n tabuated form. The cacuaton of atomc propertes by the Hartree-Fock- Sater (HFS) approxmaton s descrbed n the frst secton. Secton s devoted to the descrpton of the CRE mode, whch s used to obtan the popuatons of the atomc eves and the on composton of the pasma. The cacuaton of thermodynamc pasma propertes and the equaton of state are descrbed n the thrd secton. The resouton of the radaton transport equaton n a wde spectra range depends upon the competeness of the optca coeffcents, whch are dscussed n Secton 4. The Panck mean opactes are presented n Secton 5. The next three sectons are prmary dedcated to cacuatng hghy accurate EUV nformaton for tn. Secton 6 dscusses the mtatons of usng the smpfed methods for generatng EUV atomc data and descrbes the advanced S-dependent Hartree-Fock (SD HF) atomc method. Detaed EUV transtons are presented n Secton 7. Fnay, Secton 8 presents cacuatons of tn EUV opactes wth very hgh resouton.. HFS ENERGY EVES AND OTHER DATA In modeng the dynamcs and the output of an EUV source, one needs to dstngush two aspects. Frst, n sovng the hydrodynamc part of the probem, the pasma nterna energy must be corrected and re-dstrbuted accordngy to the radaton transport. A key eement n ths process s the competeness of the radant energy redstrbuton n the whoe pasma doman wthn the very broad spectra range of partcpatng photons. The second aspect s the detaed cacuaton of the effectve radaton of the EUV source wthn the operatng energy range of 3.5 ± % nm. In ths case, the cacuaton of the radaton transport must be orented to the accurate accountng of ony those photons, whose energes are wthn the narrow EUV range. The smuaton of the dynamcs of the pasma evouton typcay nvoves a wde range of temperature and densty vaues and a very compcated onc structure. From the atomc physcs vewpont, detaed resouton of each possbe eve for each possbe on (and

consequenty, each eectronc transton n the on) n a wde range of temperatures, denstes, and energes s enormousy aborous, especay when one s accountng for the fne structure of each spt eve. The dffcutes come from the very arge number of atomc terms and spt 5 eves. For exampe, to spt a d she, one needs to account for 6 terms and 37 eves. The 7 f she w have 9 terms and 37 eves. In approxmaton of confguratons, the transton n n 4d 4d 4 f w ony have one strong ne, whe spttng the shes nto eves repaces ths ne by severa hundred weak nes. Knowng that the tota transton strength s unchanged, the radaton n strong nes may be coapsed for a dense pasma, whe the radaton n weak nes stays optcay transparent. The ne spttng may dramatcay nfuence the tota hydrodynamc behavor of a pasma through the radaton transport mechansm [5]. The resuts of theoretca approxmatons of atomc data strongy depend on the chosen theoretca modes [6-8], and ths s partcuary true for ntermedate- and hgh-z materas. To descrbe a sphercay symmetrca quantum system, the sef-consstent fed methods, such as Drac-Fock (DF) or Hartree-Fock (HF) methods, are beeved to be the most effectve. However, the very arge number of ons and ther atomc eves and transtons nvoved nto the pasma computaton, mts ther appcabty. In consstenty sovng the MHD equatons and the RTE, the pasma propertes, opactes, equaton of state, and on concentratons can be obtaned from the structures of atomc eves and transton probabtes, cacuated by smpfed methods. One such method s the HFS approxmaton [9], whch aows the determnaton of the energes and other atomc characterstcs for each n-confguraton of each on that mght exst n the pasma. The smpfcaton of the method normay resuts n the shft of some spectra nes from ther true paces by severa percent, whch s not probematc n determnng the ntegra radaton fux. Nevertheess, we cannot negect the fne structure of open shes, and have thus combned the HFS methods wth the spttng of the energy eves and spectra nes by means of the Racah theory of anguar momentum wthn the framework of perturbaton theory [0,, 5]. The accuracy of the HFS mode s typcay wthn severa percents for the spt eves. However, ths accuracy s not good enough for the second part of our project, whch nvoves determnng the EUV output of the source wthn the very narrow % bandwdth. [To obtan a spectroscopc accuracy of the EUV optca coeffcents, the atomc data are cacuated by the Sdependent (SD) HF method, as dscussed n Secton 6.] In the HFS method the potenta of drect eectron nteracton s cacuated from the rada wavefunctons of partcpatng eectrons, and exchange nteractons are averaged n the form of the exchange potenta V ex. The rada wavefuncton of an atom can be represented as the product of rada wavefunctons P n (r) of the eectrons. It s assumed that a q n equvaent eectrons have the same rada wavefuncton: d dr Z( r) + + V r r ( + ) ( r) + ε n P r ( r) = 0, Z0 ( ) ρ ( x) dx Z r 3 3 0 = dr, Vex ( r) = ρ, ρ ( r) = r q r π n r ex n 3 n P n (). r (.) 3

Tradtonay, n atomc physcs, the energes are expressed n Ry, and dstances are n a 0 or Bohr unts. In the above equatons, ε n s bndng energy of the eectron, Z (r) s the effectve charge of the on fed, Z 0 s the nucear charge, and ρ (r) s the eectron densty. From souton of the HFS equatons, one may fnd one rada wavefuncton for each she of the atom or on. We normay cacuate the functons for each on, startng from the neutra and endng wth the totay onzed core. Obtaned wavefunctons are used to ater cacuate reatvstc correctons, transton waveengths, and dpoe transtons from the ground state to the hghy excted state (prncpa quantum number may reach up to 0), aowed by the seecton rues. Dscrete transtons from nner shes wth energes mted to the vaue dependent on onzaton potenta were aso taken nto account. Photoonzaton crosssectons were cacuated for a nner and excted states. Detaed spn-orbt spttng of non-fed shes was mpemented k k outsde the HFS, but wthn the CRE mode, by usng Sater ntegras F, G, and constants of k spn-orbt spttng ξ as descrbed n Ref. 5.. ION STATES AND POPUATIONS OF ATOMIC EVES To descrbe the popuatons of atomc eves, we utze the CRE mode, whch s equay appcabe for ow, medum, and hgh temperature ranges. The CRE mode accounts for the exctaton and onzaton processes that can take pace n a pasma. The fact that the CRE mode consders the transtons between a atomc eves s of partcuar mportance. Nonoca effects are accounted for n the form of an escape factor [, 3], whch negects by photoexctaton n contnuum and reduces the strength of spontaneous transtons. Such an approxmaton fary descrbes the pasma behavor n the condton of absence of an externa source of hard radaton. The pasma onzaton state and popuaton eves n for a prescrbed set of temperatures and denstes are cacuated accordng to the system of knetc equatons n statonary form: d n = n Kj + n jk dt j j j = 0. (.) The popuaton of atomc eve s determned by the set of transtons from ths eve to other eves j wth transton rates K j, as we as transtons from other eves j to ths eve wth transton rates K j. One equaton s wrtten for each atomc eve. If eve defnes the ground state, then the popuaton of ths atomc eve gves the concentraton of the on n the pasma. Impact-eectron exctaton and de-exctaton, mpact eectron onzaton, three-body recombnaton, spontaneous transton, and photo- and d-eectronc recombnaton are ncuded n cacuaton of tota rates of eectron transton. From the known popuaton of eves, on and eectron concentratons N and N e are defned for a gven temperature. A numerca smuaton of concentraton of tn ons cacuated by the CRE mode wth spttng HFS atomc eves s presented n Fg. for varous pasma denstes. Due to onzaton of the d-she, the ons change very qucky when temperature ncreases. Ths condton resuts n a ow concentraton of the major on, hardy arger than 50% at best. At the same tme, the cear advantage of tn as an EUV workng eement s that severa ons are EUV productve, 4

whch notceaby wden the range of the requred temperature n the devce. Ths stuaton s especay true at ow densty, when EUV emttng ons are appearng at sghty hgher temperature and not changng as fast. Fg. : Concentratons of tn ons for varous denstes Our cacuatons can be vadated and benchmarked by comparson to the resuts of other authors. Accordngy to the resuts cacuated by the CRE mode and reported n [4], the 3 concentratons of xenon ons at n = 0 7 cm e and T e = 3 ev were 0% Xe+9, 5% Xe+0, 0% Xe+, and 3% Xe+, as shown by gray shadng n the second coumn of Tabe. As seen from the other data n the tabe, our cacuatons generay agree wth these fgures, except that we obtaned a twce arger vaue at T e = 3 ev for the concentraton of Xe+0 ons. We tend to rey on neary 45% concentraton of Xe+0, whch woud fx the tota concentraton of a ons to 00% and be more reastc n terms of the reatve concentraton of the major on to the others. The authors of Ref. 4 aso reported on concentratons at smar temperature and densty 5

for a tn pasma to be 5% Sn+8, 40% Sn+9, 30% Sn+0, and 5% Sn+. As shown n coumn of Tabe, these resuts are smar to our computatons for T e = 9 ev. Resuts of the tabes generay confrm the very compex composte nature of xenon and tn pasmas, where at the same condtons one may expect to fnd up to seven on speces. Tabe : Concentraton (n percent) of xenon ons at 0 7 cm -3. Top two rows gve temperature and pasma densty. Shaded data from Ref. 4. T e, ev 3 30.8 3.0 3. 3.4 3.6 3.8 3.0 3. 3.4 ρ e, g/cm 3.9 0-6.8 0-6.7 0-6.6 0-6.6 0-6.5 0-6.4 0-6.3 0-6.3 0-6 Xe+7 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 Xe+8.49.5.04.84.66.50.35..0 Xe+9 0 4.5.90.69 0.5 9.37 8.9 7.3 6. 5.5 Xe+0 5 49.08 49.03 48.88 48.63 48.30 47.9 47.4 46.86 46. Xe+ 0.07 3.34 4.6 5.9 7.9 8.49 9.76 3.0 3.6 Xe+ 3.4.40.70 3.0 3.36 3.70 4.0 4.53 5.00 Xe+3 0.04 0.05 0.06 0.07 0.09 0.0 0. 0.5 0.7 Tabe : Concentraton (n percent) of tn ons at 0 7 cm -3. Top two rows gve temperature and pasma densty. Shaded data from Ref. 4. T e, ev 3 5 8 9 30 3 3 35 ρ e, g/cm 3.37 0-6.9 0-6.4 0-6.09 0-6.05 0-6.0 0-6.9 0-6 Sn+7.4.85.64 0.9 0.5 0.8 0.04 Sn+8 5 40.58 0.60 4.98 0.50 7.5 4.74.3 Sn+9 40 37.07 45.50 4.80 38.5 3.73 6.96.58 Sn+0 30 8.7 6.94 33.69 39.35 43.3 45.30 40.8 Sn+ 5 0.40 3.85 6.59 0.34 5.0 0.4 36.66 Sn+ 0.00 0. 0.9 0.63.5.6 8.80 6

3. THERMODYNAMIC PASMA PROPERTIES AND EQUATION OF STATE After eectron N e and on N concentratons are found by the CRE mode, they can be substtuted nto a set of equatons of state. Generay, the two-temperature approxmaton for pressure p contans the correspondent terms for knetc energy of ons and eectrons. The equaton for nterna energy e nt contans the terms for onzaton and exctaton of eectrons: P = kt e nt 3kTe N = ρ N + kt N, e e e 3kT N + ρ + ρ N j j I + j N ε N Here we use tradtona notaton for Botzmann s constant k, onzaton potenta I of the atomc eve, pasma densty ρ, exctaton energy of eve j to eve ε, eectron and on temperatures T e, T. Resuts of these cacuatons are presented n Fg... j (3.) Fg. : Cacuatons from equaton of state for tn ons at varous denstes 7

Recproca to the resstvty η, the eectrca conductvty σ s found as the sum of conductvtes defned by the eectron scatterng σ and σ on charged and neutra partces. ( kt ) 3 c η 4 e β 3 N ee = = + ; σ c = ; σ. n σ σ σ π π e Z = (3.) Λ π m kt N s c n Here we use standard notaton for Couomb ogarthm Λ, concentraton of neutra atom N, eectron charge e, and mean on charge Z. Parameter β s an eectron-eectron scatterng 0 correcton [5]. Emprca vaues of transport crosssectons s 0 are taken from [6]. n e e 0 0 4. DETAIED SET OPACITIES IN WIDE ENERGY RANGE The eectronc transtons and ther accompanyng absorpton and emsson of photons are subdvded nto three types: bremsstrahung; photoonzaton from ground, excted, and nner eves; and dscrete transtons. The atter s approxmated n the form of dpoe transtons and ncudes transtons between ground and excted states, transtons between excted states, and party the transtons from nner shes. Because of ts mportance, the profes of spectra nes are processed very carefuy by means of a major broadenng mechansms, such as radaton, Stark, Dopper, and resonance broadenngs [7]. The tota absorpton coeffcent κ abs s cacuated as a sum of absorpton coeffcents for free-free κ ff, bound-free κ bf the popuaton eves: κ ( T, ρ, hω) = ff κ ( T, ρ, hω) = bf κ ( T, ρ, hω) = bb σ j j, k, and bound-bound κ bb radaton transtons weghed to the vaue of ( T, hω ) n ( T, ρ ) n ( T, ρ )( exp( hω kt )) σ ( hω) N j π e f m c e jk j ( T, ρ), Φ( T, ρ, ω) N e j ( T, ρ ). Knowng the crosssectons of nverse processes, the tota emsson coeffcent κ em s cacuated by smar formuae. The vaues of eectron densty n e and on densty n are cacuated as descrbed n prevous secton. The oscator strengths f jk, crosssectons for photo absorpton σ, crosssectons for photoonzaton σ j, and ne profe Φ are gven esewhere [5, 6, 8]. As shown earer [5, 8], accountng for eectrostatc and spn-orbt spttng of shes and spectra nes consderaby nfuences the dynamcs of the energy baance n the pasma. A method of cacuatng the optca coeffcents for hgh-z pasma was deveoped and mpemented [5, 9] by means of the CRE mode. The mportant feature of ths method s the jont use of HFS atomc data and Racah technques of anguar moments. Ths feature makes possbe the use of comparatvey unsophstcated methods to consder the compex eectronc structure of each, (4.) 8

partcpatng on, and the compcated spttng of each confguraton nto terms (usng the Scoupng approxmaton), over a wde range of spectra frequences, and n the expected range of temperatures and denstes. By ths technque, the detaed emsson and absorpton spectra are ntay cacuated over a very compete spectra frequency scae (up to 00,000 ponts) for the expected range of MHD vaues. The wdth of the frequency nterva was comparabe to the Dopper wdth of the strongest spectra nes, whch provdes a satsfactory resouton of the ne profes. In Fg. 3 we present cacuated emsson and absorpton coeffcents for the tn pasma at 6 typca temperature and densty (6 ev and 0 cm -3 ). 5. PANCK GROUP AVERAGED OPACTITES IN WIDE ENERGY RANGE Because of the arge sze of the generated opacty tabes, the practca use of such detaed data s not convenent, and the emsson and absorpton coeffcents are thus averaged n spectra groups. A rgorous theory for averagng the opactes wthn a group of frequences does not exst. Such an averagng procedure s consdered correct ony when the absorpton coeffcent s constant wthn the group, or the optca thckness of each ne s very sma, and the absorpton becomes a near-ke functon from the frequency. Ths stuaton s ony possbe for a contnuum, but even n that case, every photoonzaton threshod must become a boundary of a group. It becomes even more compcated for the ned spectrum, when the absorpton coeffcent often drops severa orders of magntude wthn a very mted frequency nterva, say, from the center of a very strong spectra ne to the wngs of the same ne, and the center of the ne s optcay thck. Moreover, the temperature and densty vaues may vary. Ths change eads to spectra nes for the other ons, aong wth changes n the wdth of the exstng nes. From a practca vewpont, an organzed seecton of the strongest nes s a reasonabe way to descrbe the optca coeffcents wthn the most mportant hydrodynamc areas for typca temperature and densty vaues. The other nes are averaged wthn broad groups. Unfortunatey, the prmary goa of the numerca smuaton s the determnaton of the typca hydrodynamc parameters wthn the mportant areas of the pasma doman! For a unform sotherma pasma, the optca thckness of a spectra ne s determned by mutpyng the absorpton coeffcent of the ne by the near dmenson of the pasma τ ε = κ abs ε. In the case of a nonunform nonsotherma pasma, ths defnton s generazed ( ) ( ) as ( ε ) κ ( τ = T, ρ, ε ) d over the nterva, where the on exsts emttng wth frequency ε. abs The borders of the groups are cacuated from the foowng consderatons. Frst, the wdth of a group ε shoud never exceed α T, where T s the pasma temperature, and α s a chosen parameter for averagng. Ths condton provdes the smooth averagng of opactes n the contnuum regon, where the averagng s normay performed wthn the broad groups, and the optca thckness s much ess than one. Second, nvarabty of the optca coeffcent s requred for a chosen vaue β wthn the group. Ths condton provdes a very specfed resouton of the spectra nes wth the optca thckness neary at unty, and the wngs of the nes wth the optca thckness greater than unty. A fna consderaton s that a those 9

frequences whch beong n the detaed spectrum to a ne wth ( ε ) τ are aso ncuded n the doman of the groups. Ths condton provdes a very thorough resouton of the strong nes n the averagng spectrum. By varaton of the parameters α and β, severa group mean opactes are generated wth dfferent eves of competeness and deta. Fg. 3: Optca coeffcents for tn pasma at T = 6 ev and N = 0 6 cm -3 Fg. 4: Panck mean optca coeffcents for tn pasma, averaged over 69 (eft) and 340 (rght) groups 0

Based upon severa recent studes [3, 4], t s supposed that the maxma radaton fux corresponds to the moment of pnch formaton. Typca tn pnch parameters are the foowng: 7 the temperature s cose to 5 ev, the densty s 3 0 cm -3, the spectra range n radaton energes vares from 5 ev to 50 ev, and the average optca pasma thckness s cm. Usng these parameters, we have generated a basc set of opactes, averaged wthn 69 spectra groups. Takng a wder set of temperature and densty pars, we have generated the optma scae of spectra groups. Combnng a scaes, the resutng energy scae has a tota of 340 energy groups. Resuts of the Panck mean absorpton and emsson coeffcents for the mentoned temperature and densty vaues are shown n Fg. 4 above. As the next step, the output of the sef-consstent hydrodynamc/radaton transport cacuatons wth these opactes can be used for further mprovements n the quaty of the coeffcents. 6. HF ENERGY EVES AND OTHER DATA FOR EUV EMITTING TIN ION SPECIES The accuracy of the HFS mode s typcay wthn severa percents for the spt eves, whch s nsuffcent for the narrow % bandwdth. To obtan a hgher accuracy of the EUV optca coeffcents, we use the SD HF method, whch s more accurate but sgnfcanty more advanced and dffcut to mpement. Strcty speakng, the HF method encompasses severa methods for cacuaton of varous atomc structures. Beng apped to the same atomc system, each modfcaton of the HF method can produce dfferent resuts, so the choce of the approprate method s very mportant n terms of the overa accuracy of the cacuatons. In the HF method the wavefunctons, egenvaues, and tota energy of an atom are found from the varatona prncpe for the we-known Schrödnger equaton: Ĥ Ψ = EΨ, (6.) where Ĥ s a Hamtonan of the atom, ncudng the nteracton of each eectron wth the nuceus and the other eectrons; E s the tota energy of an atom; and Ψ s the atomc wavefuncton, whch s expressed through the atomc orbtas ϕ r ) ϕ ( r ) as ϕ ϕ n ( r ) ϕ( r ) ϕ( rn ) ( r ) ϕ ( r ) ϕ ( r ) ( r ) ϕ ( r ) ϕ ( r ) n n ( j n j n Ψ =. (6.) ϕ The number of the atomc orbtas s determned by the number of shes n the atom. For nstance, the tn atom has s s p 6 3s 3p 6 3d 0 4s 4p 6 4d 0 5s 5p ground state, or orbtas. Each orbta s found from the system of ntegra-dfferenta equatons for functons ϕ n :

d dr N n n ( + ) + r x β x Z r + N N x x ( n) y ( r) + α ( n, n ) y ( r) ε ϕ ( r) x ( n, n ) y () r ϕ () r ε ϕ () r = 0. n, n n x f n x n, n n n, n n n n x x n, n The expressons for the rada ntegras y x x n, n ( r), as we as the potentas f yn n ( r) α y x () r, and β y x () r x n, n x n, n n n x,,, are omtted here and can be found esewhere [6]. Here, we just note that they defne the nteracton of the eectrons of the n she wth the other eectrons of the same she, averaged over a anges, and wth the eectrons of the other shes, ncudng both the usua and the exchange nteractons. In genera, the coeffcents f x, α x, and β x depend on the whoe set of quantum numbers, defnng the atomc eve under consderaton, and partcuary on and S. Consequenty, dfferent orbtas (or rada wavefunctons) ϕ, ϕ n n, K correspond to dfferent terms of varous eectron confguratons. In a we-known monograph [8], Cowan ntroduces the concept of the S-dependent Hartree-Fock (SD HF) cacuatons. Hs approach uses a dfferent set of rada wavefunctons for each S term, assumng the pure S-coupng scheme. To be precse, the method shoud be caed SνD HF, because the tota energy for d-shes aso depends on the senorty number ν. The number of terms can be arge, especay for d- and f-shes. Sgnfcant smpfcaton can be ganed by assumng that for a terms the wavefunctons and egenvaues are dentca and cacuated for the center of gravty of the she. The so-caed average term approxmaton works very we for the hghy excted states, but s qute unpredctabe for the outer and nner shes. In contrast to the average term HF, a SD HF wavefunctons, egenfunctons, and Sater ntegras depend on the and S quantum numbers for nner and outer shes. Ths condton requres onger, but more accurate computatons, and very mportant for cacuaton the tota energes of the atomc eves. Tabe 3 presents our cacuatons of the energes of the Xe XI nner shes by the SD HF method. These are compared to the resuts obtaned by the weknown Cowan HF average term code, whch s wdey used n spectroscopc research for dentfcaton of atomc eves. The correspondence s wthn -5%, athough our resuts are consstenty ower than those from the Cowan code for the nner shes, and sghty hgher for the open she. Smar resuts for Sn VIII are presented n Tabe 4, and the energes for the Sn VIII Sn XIII, cacuated by varous methods are shown n Appendx A. The anguar wavefunctons are normay cacuated separatey to the rada wavefunctons by the summaton of the eectron momentums. As a rue, the two mted cases are rarey reazed n practce, that s, when eectrostatc nteracton s consdered predomnant (Scoupng) or spn-orbt nteracton exceeds the eectrostatc nteracton (jj-coupng). In such ntermedate cases, the Hamtonan matrx cannot be wrtten dagona n any coupng schemes. Therefore, the compete matrx s wrtten by transformng the Couomb matrx from Srepresentaton to jj-representaton. The egenvaues of the Hamtonan matrx are found ater by numerca dagonazaton, and the egenvector (purty vector) defnes the composton of the eve, correspondng to ths egenvaue. The eve s normay assgned accordng to the hghest

contrbuton of the bass term from the purty vector. Note that the energy eves, found wthn the ntermedate coupng, never have 00 percent pure S- or jj-coupng, and the dfference aways exsts between the eve assgned wthn the ntermedate coupng or the eve cacuated wthn the pure coupng scheme and accordngy assgned as requred by the scheme. The reatvstc effects can be neggbe for the ow-z eements, but apparenty become evdent for ntermedate-z and hgh-z eements. The most wdespread and reatvey easy way to account for them wthn the HF method s by one-eectron reatvstc correctons wthn the framework of perturbaton theory. The more strct and accurate way to account for reatvstc effects s to use the Drac-Fock approxmaton, but ths sgnfcanty compcates the probem, whe the gan from t woud be pronounced ony for the hgh-z eements. Tabe 3: Energes of Xe XI nner shes by SD HF and average term HF. HEIGHTS-ATOM Cowan Conf S P D 3 F G E av s 46.40 46.73 46.65 46.075 46. 565.4407 s 39.7600 39.7356 39.7349 39.760 39.7394 47.8558 p 6 365.300 365.854 365.847 365.757 365.89 376.567 3s 9.9495 9.993 9.989 9.96 9.935 98.588 3p 6 8.8833 8.869 8.865 8.855 8.866 84.807 3d 0 63.385 63.3609 63.3605 63.3530 63.364 64.559 4s 7.480 7.406 7.406 7.408 7.4080 8.6868 4p 6 3.45 3.33 3.3 3.09 3.5 3.9933 4d 8 7.087 7.508 7.568 7.0 7.39 7.07 Tabe 4: Energes of Sn VIII nner shes by SD HF and average term HF. Conf P 4 P D HEIGHTS-ATOM 3 D F 4 F G Cowan H E av s 09.3064 09.303 09.379 09.37 09.369 09.909 09.309 09.3064 65.7804 s 3.6554 3.6503 3.6769 3.667 3.6659 3.6399 3.650 3.6554 34.0976 p 6 98.609 98.5977 98.644 98.609 98.634 98.5873 98.5994 98.609 307.0733 3s 7.5397 7.535 7.5580 7.5449 7.5484 7.564 7.5367 7.5397 75.3565 3p 6 6.7069 6.703 6.753 6.7 6.757 6.6935 6.7039 6.7069 63.887 3d 0 45.3449 45.3404 45.3633 45.350 45.3537 45.335 45.349 45.3449 46.075 4s 8.6978 8.6948 8.709 8.7009 8.7030 8.6896 8.696 8.6978 9.4893 4p 6 5.0560 5.053 5.0665 5.0589 5.0608 5.0484 5.0544 5.0560 5.699 4d 7 0.634 0.836 0.689 0.343 0.49 0.33 0.796 0.634 0.565 As foows from Tabes 3 and 4, t s very hard to benchmark and verfy the accuracy of the computaton of such compcated eements as xenon or tn wthout expermenta resuts. Before proceedng wth our tn cacuaton, we checked the accuracy of our code on we-known eements wth avaabe expermenta measurements of the eves for oxygen (Tabe 5) and argon (Tabe 6). The NIST tabes are pubshed n Ref. 0. The smpcty of oxygen and argon comes from the pronounced S-coupng approxmaton, whch aows us to unquey dentfy and assgn the eves once the ntermedate coupng s apped. The average term HF method presents smar, but sghty ess accurate resuts, despte the fact that the average term approxmaton s expected to work very we for the excted states [8]. As mentoned above, the ntermedate-z eements, such as tn, do not have any pure 3

coupng scheme, and the ntermedate coupng approxmaton s the ony reasonabe choce n computer cacuaton of atomc eves. However, the assgnment of the cacuated eves s arbtrary, preservng ony the exact quantum number of tota momentum J. In coumn 3 of the Tabe 7, we present our resuts for the Sn VIII energy eves, and n the coumn 7 of the same tabe are the resuts of the Cowan code. We have ordered the eves by ncreasng the tota moment and arranged the Cowan code eves accordngy to our dentfcaton. Snce dfferent ntermedate coupng codes mght have dfferent namng schemes, the reasonabe way of presentng the resuts s arrangng the eves accordngy to the tota moment of the eve. The expermenta energy eves are avaabe thanks to the Ref. 3, and we can benchmark the accuracy of these cacuatons. As one can see, the energy eves computed by us are cose to the Cowan code resuts and dffer from the expermentay defned vaues from 3% to %. Tabe 5: O I p 3 3s excted state energy eves (n cm - ). Excted eve HEIGHTS-ATOM Cowan NIST eve Acc,% eve Acc Tabes p 3 ( P)3s( S) 3 P 0 507.0 5645.60 p 3 ( 4 S)3s( S) 3 S 693.7.00 9454.00.35 306.78 p 3 ( P)3s( S) P 5458.8 60975.40 p 3 ( P)3s( S) 3 P 507.0 5647.0 p 3 ( D)3s( S) 3 D 37.88 3.6 34690.0 6.67 7387. p 3 ( 4 S)3s( S) 5 S 0.00 0.00 0.00 0.00 0.00 p 3 ( P)3s( S) 3 P 504.33 5650.60 p 3 ( D)3s( S) D 3506.93 3948.50 p 3 ( D)3s( S) 3 D 34.9 3.64 3469.0 6.7 7379.33 p 3 ( D)3s( S) 3 D 3 37.88 3.70 34693.0 6.77 7367. Tabe 6: Ar II 3p 4 3d excted state energy eves (n cm - ). Excted eve HEIGHTS-ATOM Cowan NIST eve Acc,% eve Acc Tabes 3p 4 ( D)3d( D) S 0.5 40846.478.09 5850.. 5765.770 3p 4 ( D)3d( D) P 0.5 448.845.93 868.8 56.0 4493.6398 3p 4 ( 3 P)3d( D) P 0.5 35.034 83.63 9307.9 65. 38.69 3p 4 ( 3 P)3d( D) 4 P 0.5 809.789.4 836.6 53.3 4900.690 3p 4 ( 3 P)3d( D) 4 D 0.5 394.0666 3.97 403.8.6 40.340 3p 4 ( S)3d( D) D.5 5884.3505 3.57 63587.7 33.6 47604.4745 3p 4 ( D)3d( D) P.5 449.8.98 9756.3 53. 408.57 3p 4 ( D)3d( D) D.5 40839.64 0.83 4700.9 39.0 4050.899 3p 4 ( 3 P)3d( D) P.5 35483.3365 65.96 9456. 593.0 334.5 3p 4 ( 3 P)3d( D) 4 P.5 8366.83.0 396.9 5.9 575.754 3p 4 ( 3 P)3d( D) D.5 37436.770 06.9 035.9 470.4 847.679 3p 4 ( 3 P)3d( D) 4 D.5 9.355 3.96 99.0.4 303.3660 3p 4 ( 3 P)3d( D) 4 F.5 53.596 0.95 5783.3 4.9 044.0744 3p 4 ( S)3d( D) D.5 5884.407 4.46 63853.6 35. 4764.869 3p 4 ( D)3d( D) D.5 4036.0439.57 5648. 35.9 40008.35 3p 4 ( D)3d( D) F.5 377.7595.76 4606. 49. 3097.990 3p 4 ( 3 P)3d( D) 4 P.5 8690.6764 0. 3679.7 5.3 5548.5860 3p 4 ( 3 P)3d( D) D.5 3640.756 94.04 0099.4 444. 8759.9507 3p 4 ( 3 P)3d( D) 4 D.5 0.3957 8.4 5.4 0.9 53.8450 3p 4 ( 3 P)3d( D) F.5 5.5343 4.7 909.4 63.3 780.340 3p 4 ( 3 P)3d( D) 4 F.5 9.9700 0.59 5500. 43.8 0780.383 3p 4 ( D)3d( D) F 3.5 3759.445 9.8 4653.6 49. 379.749 3p 4 ( D)3d( D) G 3.5 6396.399 0.66 945.5 33.7 876.667 3p 4 ( 3 P)3d( D) 4 D 3.5 0.0000 0.00 0.0 0.0 0.0000 3p 4 ( 3 P)3d( D) F 3.5 387.5803 6.9 7647.0 64. 685.885 3p 4 ( 3 P)3d( D) 4 F 3.5 63.0456.87 508.4 45. 0389.7346 3p 4 ( D)3d( D) G 4.5 6396.399 0.66 9074. 3.9 876.667 3p 4 ( 3 P)3d( D) 4 F 3.5 0.96.79 4509.5 47. 9858.9536 4

From a theoretca standpont, the agreement of the cacuated atomc energy eves wth the expermentay measured vaues wthn 3% to % may be consdered encouragng, but n reaty, the requred EUV bandwdth has much greater restrctons. The we-known effect of overestmaton of the theoretca energy-eve spttng resuts n sghty arger vaues of Couomb eectron-eectron nteractons [8, ]. The eectron correaton cannot be determned scaed-down theoretca vaues of the snge-confguraton Sater ntegras. Any expermenta data can actuay be a great hep n determnng the exact vaues of the scang factors, as demonstrated n Tabe 8. The expermenta vaues n the second coumn are taken from Ref.. As expected, the purey theoretca vaues for the ground energy eves of Xe XI, shown n coumn 3, are sghty overestmated, by approxmatey 0%, but scang the Sater ntegras to a factor of 0.8 yeds a very accurate resut, wthn ess than % for most ground eves. Thanks to the Ref., we can aso benchmark the accuracy of our cacuaton for the Xe XI 4d 7 5p excted eves. Addtonay, we have cacuated the same energy eves by the Cowan code. The cacuated resuts for Xe XI 4d 7 5p excted eves n Tabe 9 show exceent agreement wth the expermenta vaues, as accurate as -% for most eves. As evdent from the tabe, the Cowan code produces sghty hgher energy vaues, whch woud shft the 4d 8 4d 7 5p transton array toward onger waveengths. Tabe 7: Sn VIII 4d 7 ground state energy eves. Exp HEIGHTS-ATOM (th) HEIGHTS-ATOM (ft) Cowan cm - Energy, cm - Acc Energy, cm - Acc Energy, cm - Acc 4d 7 3P 0.5 35458 36645 3.35% 3645.80% 37684 6.8% 4d 74 3P 0.5 3946 4936 4.3% 336.59% 5937 8.3% 4d 7 3P.5 30657 3889 7.8% 337 4.83% 3347 5.5% 4d 74 3P.5 880 907 4.33% 859 0.% 047.99% 4d 7 D.5-69434 7083 75845 4d 7 3D.5 4477 43707.06% 43468.60% 45637 3.3% 4d 74 3F.5 53 3375 0.06% 83 5.5% 083 0.58% 4d 74 3P.5 0373 448 5.8% 955 4.6% 7 8.83% 4d 7 D.5 75377 6548 3.% 6689.35% 795 4.96% 4d 7 3D.5 33670 4009 4.77% 457 3.30% 34785 3.3% 4d 7 3F.5 4545 453 0.48% 44567.95% 48749 7.5% 4d 74 3F.5 034 57 7.89% 0605.49% 00.3% 4d 7 3F 3.5 49476 49870 0.80% 48753.46% 5498 6.% 4d 74 3F 3.5 6986 730.05% 676 3.87% 6704 4.03% 4d 7 3G 3.5 900 8378.5% 9338.6% 9358.3% 4d 74 3F 4.5 0 0 0 0 4d 7 3G 4.5 636 884 3.3% 976.50% 308.97% 4d 7 3H 4.5 3775 36996.00% 3743 0.89% 3807 0.85% 4d 7 3H 5.5 303 9369 3.% 3004 0.89% 30933.05% 5

Tabe 8: Xe XI 4d 8 ground state energy eves (n cm - ). Exp Theor Acc,% Modfed Acc,% 4d 8 3 F 4 0 0 0 4d 8 3 F 3 340 5069 4.68 384 0.34 4d 8 3 F 505 5456.65 59 0.09 4d 8 3 P 6670 940 0.4 705.33 4d 8 3 P 0 30 37699 7.04 340 0.6 4d 8 3 P 3460 39394 3.8 3545.54 4d 8 G 4 40835 4555.55 3997.5 4d 8 D 4900 4604 7.7 43006 0.5 4d 8 S 0 8830 99330.7 85587.89 Tabe 9: Xe XI 4d 7 5p excted state energy eves (n cm - ). J HEIGHTS Exp Acc,% Cowan Exp Acc,% J HEIGHTS Exp Acc,% Cowan Exp Acc,% 75368 704730 748359 7393. 77807 733755. 76489 746 75695 74800.45 70466 7393.6 767530 777 758000 744955.7 73565 74800.5 77638 73749 763470 746445.3 7478 744955.78 78988 793 0.3 7575 77677 74935.89 78734 746445.43 796489 758364 773934 75054.83 733 74935.47 835707 7796 793.74 774880 754860.58 73356 75054.58 5435 74594 36.95 693853 777965 7590.4 7380 754860.6 75737 745470.56 707 7804 7666.53 739969 7590.59 76546 7555.74 70994 785056 766860.3 74543 7666.3 77635 754745.9 79 79066 768773.7 74769 766860.56 77788 758337.5 7444 74594.5 79063 77375.3 75379 768773.8 779903 760950.43 7306 745470.03 79759 775570.7 754544 77375.54 7848 765770.98 735493 7555.7 797858 780503.8 759955 775570.05 7854 767369.7 7383 754745.5 8099 788465.60 767977 780503.63 788088 775030.66 73935 758337.57 8043 79535.89 77407 788465.86 798 778350.74 745760 760950.04 83769 805.54 7865 79535.73 7938 784035.6 746488 765770.58 8675 84474 0.7 78970 805.46 795405 788396 0.88 75477 767369.68 83039 83889 0.95 80478 84474.5 80074 79805.8 763654 775030.49 85948 88639 83889.40 80734 80830 0.0 7684 778350.3 70948 695376 3.55 674946 8504 7754 784035.6 79055 73.3 6930 695376 0.6 894 788076 788396 0.04 733795 75053.9 703565 8450 83060.34 79866 79805 0.86 73796 73458 0.88 709989 84373 80863 80830 0.78 740340 737388 0.40 74976 73 0.38 846356 808 83060.40 75776 73954.40 78830 6947 687857 76448 744537. 7907 75053 0.44 7485 75730.6 69334 76304 7585.4 730 73458 0.06 7856 700.03 69654 76633 75583.37 73309 737388 0.59 736 740757.8 70035 77307 75606.0 734544 73954 0.68 74636 74655 0.04 77679 75730 0.7 775476 763070.60 74653 744537 0.39 74993 7505 0.6 78894 780864 773968 0.88 75030 7585 0.6 75345 753795 0.05 7503 700 0.56 78608 775775.3 75696 75583 0.4 75848 75670 0.30 7655 79956 7544 75606 0.4 7649 7605 0.3 7995 80886 7694 763070 0.5 76797 76505 0.37 73909 804653 769559 773968 0.57 76938 76665 0.35 73589 740757 0.66 86400 77555 775775 0.03 77309 77335 0.03 738595 74655.08 830536 805690 779656 74586 7505.0 6835 7439 5.74 6778 780705 788 0.4 74403 753795.3 73407 730345 0.8 699866 7439 3.3 78394 746840 75670.5 739736 73848 0.0 709556 730345.93 786437 786580 0.0 75933 7605. 746958 740348 0.88 74977 73848 3.5 78966 78845 0.9 75598 76505. 763685 75335.35 79344 740348.9 79063 760587 76665 0.79 769708 75960.36 73954 75335.9 79397 795995 0.6 764895 77335.0 7735 766947 0.57 7375 75960 3.00 798349 767897 788.8 77397 767700 0.7 744858 766947.97 804334 80905 0.8 77483 786580.5 78555 773968.47 75356 767700.8 80956 7838 78845 0.74 843 78909 3.0 76664 773968 0.96 8505 88875.68 7839 795995.63 83737 77564 78909.78 86889 808808 80905 0.73 730685 70049 98088 85096 88875.69 765587 79093 73800 70985.0 69444 768696 7378 7303 74855.09 69946 70985.4 77453 749845 73630 7585.43 704967 74855.40 7874 755705 743547 733755.3 764 7585.9 7 77500 73987 J = 0 J = J = J = 3 J = 3 J = 4 J = 5 J = 6 6

The resuts, presented n Tabe 7 above were cacuated n two steps. In the frst step, we obtaned the tota energes of the on by a pure anaytca method. Resuts of ths cacuaton are presented n coumn 3, and ther reatve accuracy compared to the measurements n Ref. 3 s shown n coumn 4. As expected, the pure theoretca method does not aow us to obtan an accuracy of % for the atomc eves of such a compcated on as Sn VIII. We have shaded those eves for whch the cacuatons are not accurate enough. In the second step, we used the 4 F dd = 0.8 F dd, F f ( dd ) = expermentay defned vaues and found that at f ( ) t ( ) 4 0.9 F t ( dd ), and ζ ( d ) =. 96 ζ ( d ) f 0 our vaues agree notceaby better, as shown by coumns t 5 and 6 n Tabe 7. In the formuas above, the subscrpt f means ft vaue, and the subscrpt t means theoretca vaue. At ast, we have cacuated the same eves by the Cowan code, whch despte beng purey theoretca, s st utzng scang factors near 0%, dependng upon the nucear charge of the eement [8]. Our cacuatons for the excted Sn VIII 4d 6 5p energy eves are presented n Appendx B. Expermenta vaues are taken from Ref. 3, and the ft vaues were 4 4 F dd 0.57 F dd F dd. F dd ζ d = 0. 44 ζ d, determned as f ( ) = t ( ), f ( ) = t ( ), f ( ) t ( ) 3 3 ζ ( p) = 0. 35 ζ ( p), F ( pd ) = 0.76 F ( pd ), and G ( pd ).4 G ( pd ) f t f t f =. As before, the theoretca vaues agreed wth the expermenta data wthn 5%, and once the expermenta vaues are taken nto account, the accuracy of our numerca smuaton becomes ess than % n most cases. In sum, our cacuated resuts compare favoraby wth avaabe expermenta data and a Cowan-average-term Hartree-Fock code, as we as other atomc codes wth reatvstc correctons. At the same tme, the accuracy of the purey theoretca methods s not wthn the requred %, so wthout the expermenta vaues, one cannot guarantee the cacuaton of accurate EUV transtons. We have attaned an accuracy wthn 3-7% of our purey theoretca cacuatons for Sn ons whch ack expermenta energy vaues. Once the energy eves are obtaned, we can sgnfcanty mprove our cacuatons to practcay as accurate as the measured data. t 7. DETAIED TIN EUV TRANSITIONS Emsson of ght takes pace when an atom changes from a state of hgher energy to a state of ower energy. Smary, absorpton resuts n an upward transton that s caused by the acton of the radaton fed on an atom. The energy of a transton from state to state j s formay defned as the dfference of tota energes of the atom or the on n states and j. As we known from the terature [6-8], the foowng are equvaent measures of tota strength of the spectrum ne: the probabty of a transton W j, the radaton ntensty of a spectrum ne I j, and the weghted spontaneous-emsson transton probabty g j Aj, expressed through the statstca weght of the eve g j = J + and the Ensten spontaneous emsson transton probabty rate A j. In our study, we normay use the non-dmensona (absorpton) oscator strength, whch s reated to the ne strength S by the foowng expresson: 7

8 ( ). 3 S J E E f j j + = (7.) As above, uness otherwse stated, we use atomc unts. The energy of the transton ( ) j E E s gven n Rydberg unts, + J s the degeneracy of the nta eve, and the ne strength s defned by the matrx eement of the wavefunctons of the nta and fna states wth the eectrc dpoe operator. Usng as before the ntermedate coupng approxmaton, one can determne the wavefunctons of both states n the form of a near combnaton of S-couped basc functons, where the coeffcents of the components are obtaned from the energy egenvectors for the correspondent states. Ths aows us to represent the ne strength n the form of the two dot products of the two component vectors to the dpoe-transton matrx. Such a technque reduces the probem of cacuatng the ne-strength matrx eements wth uncouped states to the matrx eements wth the S-couped bass functons. The mathematca expressons for the matrx eements are rather compex and engthy, so we drect the reader to the more compete theoretca pubcatons [7, 4, 5]. Usng tradtona atomc theory notaton for 6j- and 9j-symbos, fractona parentage coeffcents, and square brackets, we ony present the fna expressons for cacuaton of the ne strengths of typca transtons: ( ) ( ) [ ] ( ) ( )( ), { {,,,,,,, () 0.5 k k n n J S S k SS k n k n S P S S S S s S S S S S J J S J J S S nk D α α α α δ + + + + + + = (7.) ( ) ( ) [ ] ( ) ( ). {,,, () 0.5 n n J S S S n n S P S S J J S J J n D α α δ + + + = (7.3) In ths study, we do not have more than two open shes, the acute symbo s used for the core quantum numbers, the δ -functon expresses the seecton rue of nonchangeabe spn momentum, and the rada ntegra () P s defned by the matrx eement of the rada wavefunctons of partcpatng shes: ( ) ( ) ( ) + ± = 0, max, ()., max dr rp P P n n δ (7.4) Eectrc dpoe transtons can occur ony when the oscator strength (and, correspondenty, the ne strength) n (7.) s non-zero. From the propertes of the matrx eement n the defnton of the ne strength S, t foows that a transton can occur ony when the partcpatng states have opposte party, and. 0, ± = = J J J The transton 0 = J = J s not aowed. In cacuatng the EUV radaton from tn or xenon pasma, another dffcuty may take

pace, such as accountng for the transtons from the nner shes. As reported n recent studes [4], some nner eves of the ons wth ntermedate and heavy atomc numbers may have energes comparabe to the onzaton potenta, whch sgnfcanty decreases the accuracy of the HF method. For exampe, n xenon or tn pasma wth major radatng waveength regon of q q nterest near 3 nm, transtons from outer shes 4d 4d 5p must be consdered n q q combnaton wth the transtons 4d 4d 4 f and the transtons from the nner shes, such as 4 6 q 5 q+ p 4d 4 p 4d. The exstence of sem-emprca or expermenta nformaton on atomc energy eves may greaty hep n mprovng the accuracy of the ab nto smuaton. Unfortunatey, such nformaton for the ntermedate and hgh-z eements, especay ther hghy excted ons, s not fuy presented n the terature [6]. For exampe, no Sn expermenta data are avaabe for the EUV regon, except the work of Azarov and Josh [3], whch s not actuay dedcated to the 7 6 3.5 nm range and ony deas wth the on Sn VIII and ts transton 4d 4d 5p. Accordng to them, the transton has a very wde range of spttng (near 6 nm), wth the major nes beng concentrated around -3 nm (eft fgure n Fg. 5). As shown by the rght fgure n Fg. 5, resuts of our cacuaton show very good agreement n shape, wdth, and pace of the transton. Fgure 6 presents the resuts of our cacuaton of the EUV transtons for sx tn ons, from Sn VIII to Sn XIII. Detas of these cacuatons are presented n Appendx C, and the summary s shown n Tabe 0. V.I. Azarov and Y.N. Yosh [3] Fg. 5: Expermenta data vs. cacuaton for Sn VIII HEIGHTS-ATOM cacuaton 4 p 7 6 d 4d 5 transton We have determned that the 4d - 4f and 4d - 5p transtons ony party cover the EUV range of nterest. Among the fve ons startng from Sn IX, the hghest EUV emsson shoud be from the Sn XI 4d 4 3 4d 4 f transton and from the Sn XII 4d 3 4d 5 p transton. Smar resuts were obtaned by means of the Cowan code and are presented n Fg. 7. Despte the dfference n the hghest vaues of the oscator strength, the code aso predcts the parta coverage of the EUV range by severa transtons, when the hghest emsson corresponds to the Sn X 4d 5 4 4d 4 f and Sn XI 4d 4 3 4d 4 f transtons. 9

Fg. 6: Cacuaton of tn EUV transtons by the HEIGHTS-ATOM code 0

Fg. 7: Cacuaton of tn EUV transtons by the Cowan code

Tabe 0: Summary of the tn EUV transtons. Sn VIII Sn IX Sn X Sn XI Sn XII Sn XIII Ground Confguraton 4d 7 4d 6 4d 5 4d 4 4d 3 4d # eves n GC 9 34 37 34 9 9 Excted 4f Confguraton 4d 6 4f 4d 5 4f 4d 4 4f 4d 3 4f 4d 4f 4d 4f # eves n 4f EC 346 46 346 06 8 0 Transton array 4d 7-4d 6 4f 4d 6-4d 5 4f 4d 5-4d 4 4f 4d 4-4d 3 4f 4d 3-4d 4f 4d -4d 4f Tota oscator strength 6.89 5.89 4.80 3.7.69.7 Spttng, nm 4.04-.869.400-3.859.50-.036.963-0.05.770-8.38 3.09-6.6 Tota transtons 36 45 457 39 590 6 EUV transtons 0 5 63 0 4 5 Excted 5p Confguraton 4d 6 5p 4d 5 5p 4d 4 5p 4d 3 5p 4d 5p 4d 5p # eves n 5p EC 80 4 80 0 45 Transton array 4d 7-4d 6 5p 4d 6-4d 5 5p 4d 5-4d 4 5p 4d 4-4d 3 5p 4d 3-4d 5p 4d -4d 5p Tota oscator strength 0.80 0.69 0.58 0.46 0.35 0.3 Spttng, nm 5.580-7.5.506-46.903.66-.335.644-5.673.675-6.077.95-3.57 Tota transtons 456 97 759 57 388 46 EUV transtons 0 0 66 04 3 8. DETAIED TIN EUV OPACITIES As dscussed n Secton 4, the optca emsson and absorpton coeffcents are cacuated wthn the CRE mode wth ne spttng. Ths mode s based on the Hartree-Fock-Sater (HFS) method wth both eectrostatc and spn-orbt spttng of confguratons and spectra nes. Instead of the Racah theory of anguar moments for spttng the HFS confguraton average energes nto eves, we use the better approxmaton for the spt energy eves nstead, obtaned by the reaxed core Hartree-Fock method n ntermedate coupng. In ths way, the HFS data for the EUV ons, confguratons and transtons are removed, and substtuted the mproved data drecty n the CRE. Ths guarantees that the EUV range w be covered ony by those transtons that we have cacuated accuratey by the HF method, whe the other energy ranges w st be cacuated as prevousy. The sze of the generated EUV optca tabes s not very arge (around 0-40 MB of dsk space) due to the very narrow energy range and reatvey sma number of transtons wthn that range. Therefore, t s unnecessary to provde addtona modfcatons of the resut tabes, such as group averagng, as we dd prevousy to reduce the number of spectra ponts and the amount of requred dsk space. In Fg. 8 we present resuts of our computaton of the Sn emsson and absorpton coeffcents n the EUV range of 3.3-3.77 nm or 90.04-93.7 ev, for the typca EUV temperature and densty ranges, such as 0-30 ev and 0 6-0 9 cm -3. The coeffcents are very dense, when densty s around 0 7 cm -3, and become sparse when the densty s hgher (up to 0 9 cm -3 ) or ower (beow 0 6 cm -3 ). The resuts of ths cacuaton ndcate that, to maxmze the EUV output, the EUV source needs to operate wthn the named densty range.

CONCUSION Fg. 8: Tn pasma optca coeffcents n the EUV range The report revewed the major atomc and pasma methods we use wthn the comprehensve HEIGHTS-EUV package. The methods dffer n accuracy, competeness, and compcaton. The on states, popuatons of atomc eves, and optca coeffcents, such as emsson and absorpton coeffcents, are cacuated by means of the combnaton of the HFS atomc mode and CRE pasma mode wth spttng atomc eves. The accuracy of these methods s satsfactory to smuate the pasma magnetohydrodynamc behavor, but nsuffcent to generate the pasma spectroscopc characterstcs. Sgnfcant accuracy mprovement s acheved by usng the advanced SD HF n ntermedate coupng atomc method. However, emprca correcton of atomc energy eves s necessary to obtan the spectroscopc accuracy of the data and meet the % bandwdth requrement of eadng EUV source manufacturers wthn the 3.5 nm range. The appcabty of the methods s demonstrated for the tn ons and pasma. Deta atomc and pasma propertes, such as energy eves, oscator strengths, equaton of state, reatve on concentratons, wde range opactes, and EUV range opactes near 3.5 nm are cacuated, and compared wth the avaabe expermenta resuts. The presented resuts are used by the authors n smuaton the dynamcs and characterstcs of varous EUV sources. 3