X-RAY SOURCES A N D SHOCK COMPRESSION SCHEMES FOR PHOTOABSORPTION EDGE SPECTROSCOPY. By S T E V E N R. D Y K E B.Sc.(Hon.), Physics, University of V i c t o r i a , 1993 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF T H E REQUIREMENTS FOR T H E D E G R E E OF M A S T E R OF SCIENCE in T H E FACULTY OF GRADUATE STUDIES DEPARTMENT OF PHYSICS We accept this thesis as conforming to the required standard T H E UNIVERSITY OF BRITISH COLUMBIA September 1995 © S T E V E N R. D Y K E , 1995 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of B r i t i s h C o l u m b i a , I agree that the L i b r a r y shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. D E P A R T M E N T OF PHYSICS T h e University of B r i t i s h C o l u m b i a 2075 Wesbrook Place Vancouver, C a n a d a V 6 T 1W5 Date: Abstract A scheme to probe the temporal evolution of the K-shell photoabsorption edge i n shock compressed aluminum is presented in two components. In the first component, a one dimensional hydrodynamic code coupled to a non-local thermodynamic equilibrium calculation is developed and used to model shock propagation and target preheat in laser irradiated C H - A l - S i multilayer targets. Single and multiple shocks are used to create conditions in the a l u m i n u m by which density, temperature and ionization effects on K-shell photoabsorption can be isolated. T h e final result consists of the o p t i m a l laser parameters and target designs for creating well characterized plasma states, necessary for accurate atomic modeling of the K-edge. T h e second component describes the results of x-ray spectra observed from 12 elements across the periodic table, in search of an emission source to backlight the K-edge measurement. Spectral flatness and high emission levels in the region of 1520-1600 e V around the K-edge of aluminum are the criteria for a suitable backlighter. Lead, gold and y t t r i u m seem promising as backlighter sources, and the emission lines observed in silver provide a convenient spectral calibration source. ii Table of Contents Abstract ii Table of Contents iii List of Tables vii List of Figures viii Acknowledgement xv 1 Introduction 1.1 Laser-Matter Interactions 1 1.2 Photoabsorption Edge Research 2 1.3 Present Work 6 1.3.1 Strategy to Unfold K-Edge Dependencies . . 7 1.3.2 Backlighter Sources 7 1.4 2 1 : Thesis Organization 9 Physical Processes in Laser-Matter Interactions 10 2.1 Laser Driven Shock Waves 10 2.1.1 Shock Wave Formation 10 2.1.2 Shock Compression of Solids and the Rankine-Hugoniot Relations 13 2.1.3 Impedance-Mismatch Technique 18 2.2 Atomic Physics . . . : . iii 21 2.3 2.2.1 Introduction 21 2.2.2 A t o m i c Processes 22 2.2.3 Non-Local Thermodynamic Equilibrium Model 25 R a d i a t i o n Emission, Absorption and Transport 2 3.1 Bremsstrahlung Emission 2.3.2 Recombination Emission 2.3.3 Line Emission : 2.4 3 . 28 28 '. . . . 29 30 R a d i a t i o n Transport 31 N u m e r i c a l Simulations 35 3.1 Introduction . . 35 Laser-Target Code ( L T C ) 36 3.2.1 P h y s i c a l Content i n L T C 36 3.2.2 Numerical Methods in L T C 41 3.2.3 Comment on L T C 42 .3.2 3.3 Physics Content in the N L T E M o d e l 42 3.3.1 Calculation of State Population • • 42 3.3.2 Rate Coefficients 3.3.3 Photoexcitation 47 3.3.4 Line Transport 48 3.3.5 C o n t i n u u m Transport . . 43 . 50 3.4 Modifications for Present Work 50 3.5 Comparison of N L T E Calculation w i t h R A T I O N Code 53 3.5.1 Energy Levels and Level C o u p l i n g 53 3.5.2 Rate Coefficients 54 3.5.3 Results of Comparison . iv 55 3.6 4 Survey of Simulation Parameters . 61 Results of R a d i a t i o n - H y d r o d y n a m i c Simulations 64 4.1 64 4.2 General Results . .'. .' 4.1.1 Hydrodynamics . 64 4.1.2 Radiation Transport 64 4.1.3 Laser Absorption 67 4.1.4 Flux Limiter 4.1.5 Final Choice of Absorption Routine and Flux Limiter . . . . . . . : . 71 Optimized Shock Compression Schemes 72 4.2.1 Laser Pulse and Target Design .72 4.2.2 Density and Temperature Conditions 74 4.2.3 Discussion of Individual Compression Schemes 74 5 Experimental Facility 6 69 5.1 Laser Facility 5.2 Focal Conditions 5.3 X-Ray Spectrometers 5.4 Targets 5.5 Processing of Data 95 95 . . . . . . 97 99 . 102 102 Results and Discussion of X - R a y Spectra 110 6.1 X-Ray Spectra 110 6.1.1 Bismuth, Lead and Gold 114 6.1.2 Tungsten 118 6.1.3 Tantalum 6.1.4 Samarium . . . 118 122 v 6.2 7 6.1.5 Tin 122 6.1.6 Silver 122 6.1.7 Niobium 125 6.1.8 Yttrium 127 6.1.9 Germanium 127 6.1.10 Magnesium 131 Discussion , 131 S u m m a r y and Conclusions 136 7.1 Summary of Present Work 136 7.2 Conclusions 137 7.3 Future W o r k 137 Bibliography 139 List of Tables 3.1 ' 3.2 Photoabsorption edge energies for the different ionization stages of carbon: 51 Ionization potentials for the ground and excited states included in the N L T E model. . 52 3.3 Transitions and oscillator strengths used in the N L T E calculation. . . . . 53 4.1 Pressures and speeds characteristic of single shock waves in the C H layer. 66 4.2 Pressures and speeds characteristic of single shock waves i n the A l layer. 66 4.3 Pressures in the C H and a l u m i n u m layers when different laser absorption routines are used 4.4 68 Dependence of peak pressure w i t h i n and shock arrival time at the alum i n u m layer on the value of the flux limiter used. Simulation parameters are described i n the text. . 4.5 71 Compression and temperature conditions obtained in the shock compression schemes, along w i t h final target and laser pulse parameters. 5.1 76 Summary of sample elements used i n this work in order of decreasing atomic number. . . . 5.2 A l u m i n u m wavelengths used in deriving the spectral dispersion relation. . 6.1 E x p e r i m e n t a l and theoretical energies of lines of an a l u m i n u m spectrum when the spectral calibration equation is applied vii 105 108 115 List of Figures 1.1 Qualitative picture of a laser-target interaction. 3 1.2 , A general compression (density )-temperature space used to illustrate the method of decoupling density, temperature, and ionization effects on K edge position. 2.1 8 Schematic of the hydrodynamics in a laser irradiated target, p is the target density, and T and Tj are the electron and ion temperatures. The critical e density layer is denoted by n . . . cr 11 2.2 Illustration of the shock build up process in a laser driven shock wave. . . 2.3 Schematic of the shock discontinuity. P is pressure, p is density, E is 14 total internal energy, D is shock speed and u is particle speed. Variables subscripted with a 0 represent the undisturbed solid. Particle speed ahead of the shock front is zero. . . . 2.4 Hugoniot (solid), isotherm (dot-dash), and isentrope (dash) 2.5 Schematic of the impedance mismatch technique, pi, Pi, Ei,Usi, 17 • 19 and U Pi are the density, pressure, total internal energy, shock speed, and particle speed in the two regions. At the interface, P 4- P — P and U i = U 2- • 0 2.6 x 2 p P 20 Shaded region represents the parameter space accessible in laser-produced plasma experiments. Applicability of certain models is shown. For example, the line labeled L T E Z=5 corresponds to the range of validity of the 2.7 local thermodynamic equilibrium model for Boron 26 Radiation transport in planar geometry. 32 viii 3.1 Level coupling scheme used i n this work. S is ionization, a is radiative and three-body recombination, aoi X~ l is dielectronic recombination, X and are collisional excitation and de-excitation, respectively, and A is the spontaneous decay 3.2 44 Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma w i t h an electron density of 1 0 19 cm - 3 . T h e dashed lines w i t h circles correspond to R A T I O N calculations, whereas the solid lines w i t h squares correspond to N L T E calculations 3.3 56 Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma w i t h an electron density of 1 0 20 cm . T h e dashed lines - 3 w i t h circles correspond to R A T I O N calculations, whereas the solid lines w i t h squares correspond to N L T E calculations. 3.4 57 Ionizations produced by the present N L T E calculation and R A T I O N for . a carbon plasma w i t h an electron density of 1 0 21 cm . The. dashed lines - 3 w i t h circles correspond to R A T I O N calculations, whereas the solid lines w i t h squares correspond to N L T E calculations 3.5 58 Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma w i t h an electron density of 1 0 22 c m ~ . T h e dashed lines 3 w i t h circles correspond to R A T I O N calculations, whereas the solid lines w i t h squares correspond to N L T E calculations 3.6 59 Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma w i t h an electron density of 1 0 23 cm - 3 . T h e dashed lines w i t h circles correspond to R A T I O N calculations, whereas the solid lines 3.7 w i t h squares correspond to N L T E calculations. 60 Five-zone scheme used i n target design. . . 62 ix 4.1 Snapshot of hydrodynamic profiles w i t h i n the target at 200 picoseconds of simulation time. T h e shock wave is shown propagating through the C H layer, imminent upon the C H - A 1 boundary. T h i s particular target consists of 27 fim C H - 3 pm A l - 6 //m S i . T h e laser irradiance is 3.75 x 1 0 14 W/cm . 65 2 4.2 P l o t of reflectivity as a function of time for both the electromagnetic wave solver and inverse bremsstrahlung absorption routines. 4.3 . 70 V i s u a l aid to the important parameters in laser pulse design. T h e rise and fall times between levels is fixed at 100 ps 4.4 73 Illustration of the densities and temperatures achieved in the shock compression schemes 4.5 .. Compression (i.e. . 75 density) as a function of time w i t h i n the a l u m i n u m layer for the reference state. T h e decompression prior to shock arrival is due to radiative preheat. T h e time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. 4.6 . 78 Temperature as a function of time w i t h i n the a l u m i n u m layer for the reference state. Note the temperature increase prior to shock arrival due to radiation preheat. T h e time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. 4.7 . . . . . . . . . . . 79 Pressure as a function of time w i t h i n the a l u m i n u m layer for the reference state. T h e time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. .' x ; 80 4.8 Compression (i.e. density) as a function of time w i t h i n the aluminum layer for the high average ionization state. Note the much higher level of decompression caused by radiative preheat in this case as opposed to other simulation results. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line 4.9 83 Temperature as a function of time w i t h i n the aluminum layer for the high average ionization state. R a d i a t i o n preheat is clearly an important effect for this laser irradiance. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line 84 4.10 Pressure as a function of time w i t h i n the aluminum layer for the high average ionization state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line 85 4.11 Average ionization as a function of shock pressure for a l u m i n u m along the single shock Hugoniot. These values are obtained from the Q E O S model. 4.12 Compression (i.e. 86 density) as a function of time w i t h i n the a l u m i n u m layer for the isochoric state: T h e time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line 88 4.13 Temperature as a function of time w i t h i n the a l u m i n u m layer for the isochoric state. T h e low laser irradiance used for the first pulse leads to a situation where radiation preheat is not a significant factor. The time i n terval over which quasi-steady state conditions are achieved is indicated by the horizontal line 89 4.14 Pressure as a function of time w i t h i n the aluminum layer for the isochoric state. T h e time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line xi 90 4.15 Compression (i.e. density) as a function of time w i t h i n the aluminum layer for the isothermal state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line 4.16 Temperature as a function of time w i t h i n the aluminum layer for the isothermal state. The low laser irradiance used for the first pulse leads to a situation where radiation preheat is not a significant factor. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line 4.17 Pressure as a function of time w i t h i n the aluminum layer for the isothermal state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line 5.1 Schematic of the laser facility. O S C is the oscillator, P A is the preamplifier and A l - 4 are the amplifiers. I R M is an infrared mirror and SF1-6 are the vacuum spatial filters. A P is the hard aperture and S A is the saturable absorber. P l - 8 are polarizers, B S R is a beam splitter-reflector, N D is a neutral density optical absorber, I F is an interference filter, and G E N T E K is the piezoelectric beam energy monitor. S H G is the second harmonic generation crystal 5.2 Sample laser pulse. The horizontal scale is 1 nanosecond per division. . . 5.3 E x p e r i m e n t a l arrangement for scanning the focal distribution. N D is a neutral density optical absorber, and I F is the interference filter described in the text. xii 5.4 Images of the laser spot at best focus. Attenuation of the top picture is . 100 times greater than the attenuation of the bottom picture, and 10 times greater than the attenuation of the center picture. T h e magnification of these images is 165 times 5.5 . Target chamber and P E T crystal spectrometer. T h e diameter of the target chamber is approximately 90 cm. , 5.6 F i l m cassette. T h e film plane is indicated by the dotted line 5.7 Arrangement of the x-ray spectrometer. 103 104 T h e path length R = a + b of the longest wavelength reference line is shown. 6.1 107 A print of a sample a l u m i n u m and tantalum exposure. Some of the known lines are labeled (refer to Figures 6.2, 6.3, and 6.8) Ill 6.2 Sample aluminum reference spectrum. Some common lines are labeled. 6.3 Enlarged aluminum reference spectrum showing the clearly resolved satellite lines qr, a-d, and klj of A l X I 6.4 116 Sample lead spectrum. T h e position of the a l u m i n u m K-edge is shown by 117 Sample gold spectrum. T h e position of the a l u m i n u m K-edge is shown by the dotted line 6.7 ' 119 Sample tungsten spectrum. T h e position of the aluminum K-edge is shown by the dotted line 6.8 112 Sample bismuth spectrum. T h e position of the aluminum K-edge is shown the dotted line 6.6 . 113 by the dotted line 6.5 101 120 Sample tantalum spectrum. T h e position of the a l u m i n u m K-edge is shown by the dotted line 121 xiii 6.9 Sample samarium spectrum. T h e position of the a l u m i n u m K-edge is shown by the dotted line 123 6.10 Sample t i n spectrum. T h e position of the a l u m i n u m K-edge is shown by the dotted line 124 6.11 Sample silver spectrum. The position of the a l u m i n u m K-edge is shown by the dotted line. . 126 6.12 Sample niobium spectrum. The position of the aluminum K-edge is shown by the dotted line. ..." 128 6.13 Sample y t t r i u m spectrum. The position of the a l u m i n u m K-edge is shown by the dotted line 6.14 Sample germanium spectrum. 129 The position of the a l u m i n u m K-edge is shown by the dotted line 6.15 Sample magnesium spectrum. 130 The position of the a l u m i n u m K-edge is shown by the dotted line 132 6.16 Moseley plot used to aid in line classification xiv 133 Acknowledgement I would like to thank D r . Andrew N g for his unending enthusiasm and support for this project. His mentorship has taught me many valuable lessons about how science should be done. I would also like to thank Andrew Forsman, Kristin Smith, and George Pinh'o for all of their expert assistance in performing experiments and debugging computer code. T h e unbounded love and support of my family and friends has carried me through many rough spots over the last two years.- For this I a m truly grateful and appreciative. xv Chapter 1 Introduction 1.1 Laser-Matter Interactions T h e study of laser-produced plasmas and laser-matter interactions is closely connected to the development of laser technology [1]. It was the advent of the Q-switching technique [2] that set the stage for creating conditions previously inaccessible w i t h i n matter. Pressures of many tens of megabars, accompanied by temperatures of millions of degrees and compressions of several times solid density, are achievable when a high-power laser beam is focussed onto a target. Traditionally, only much lower pressure regimes had been accessible, and these were obtained through the use of high velocity impacts or explosive detonations. Underground nuclear explosions [3] have created pressures greater than those obtained with any other means, but experiments of this type are hardly practical. Now that extreme pressures can be routinely generated in the laboratory w i t h commercially available laser sources, the opportunity exists to map out the behaviour of matter to conditions approaching those inside stellar objects. In addition to providing an understanding of fundamental plasma physics, research involving laser-produced plasmas has found application i n many other fields, including atomic physics, astrophysics, inertial confinement fusion, x-ray lasers and x-ray lithography. A detailed account of the use of laser-produced plasmas can be found in the work of H o r a [4]. A brief overview of the processes involved in laser-target interactions at this point is 1 Chapter 1. Introduction 2 instructive. Figure 1.1 shows a schematic of the laser-target interaction. Light incident from a high-power laser is focussed onto a solid target, and initially absorbed w i t h i n the skin depth of the front surface. A s the intensity of the laser radiation increases, material from the target surface evaporates and expands outward (ablation), reaching speeds in excess of 1 0 c m / s . A t laser irradiances greater than 10 W / c m , the ablated material 7 9 2 is ionized [5], and forms a low density, high temperature plasma called the corona. T h e laser light penetrates the coronal plasma only up to the critical density. A t this point, the electron density yields a plasma frequency equal to the laser frequency and the incident laser light is reflected. Electron thermal conduction, x-ray radiation, and suprathermal electrons are possible mechanisms for carrying absorbed laser energy from the critical density layer to the ablation front, enhancing the ablation process. T h e momentum of the expanding plasma is balanced by the formation of a shock wave which propagates into the target and produces a dense, strongly coupled plasma. These plasmas typically have electron densities of ~ 1 0 23 cm - 3 , temperatures of 1-40 e V , and ion-ion coupling constants much larger than one [6]. T h e range of plasma densities and temperatures attained i n the different regions of a laser-heated solid can easily span several orders of magnitude. 1.2 Photoabsorption Edge Research T h e shock waves generated by laser-matter interactions provide a test-bed for high pressure equation of state research. In the past, however, this research has focussed on shock speed [7, 8, 9, 10] and temperature [10] measurements, which yield little information on the electronic structure, ionization state, ionization potential, ion-ion, or ion-electron correlation strength [11, 12] of dense plasmas. These properties are mutually dependent, and are important for opacity and equation of state calculations. Chapter 1. Introduction Figure 1.1: Qualitative picture of a laser-target interaction. Chapter 1. Introduction 4 One of the first indications that K-shell photoabsorption edge effects could play a role in x-ray transmission was reported i n the experiment of N g et al. [13]. In this work, a significant increase i n x-ray transmission through shock compressed a l u m i n u m targets (as compared to cold aluminum targets) was attributed to a blue shift of the K - s h e l l photoabsorption edge of the a l u m i n u m . Such a shift would allow a larger set of photon energies to pass through the target. However, the x-ray transmission measurement was critically dependent upon knowledge of the shock breakout time, and it was eventually realized that the shock breakout time had not been correctly determined. Thus, the increased x-ray transmission was not due to a blue shift of the K - s h e l l photoabsorption edge as had been thought. T h e first experiment to probe the atomic physics of a dense plasma was that reported by B r a d l e y et al. [14]. In this pioneering work, the structure and position of the K - shell photoabsorption edge of chlorine was measured as multilayer B i - K C l targets were radiatively heated and subsequently shock compressed by laser ablation. T h e observed red shift of 8 ± 4 e V during the compression phase suggested pressure ionization of the 3p level, also evidenced by the appearance of the ls-3p absorption line. T h i s experiment was also the first demonstration of the depression of ionization potential [14] as the plasma density increased past solid density. T h e edge shift was modeled as the competition between an energy increase i n the edge from ionization effects and an energy decrease due to continuum lowering. T h e edge broadening, however, was attributed to Stark broadening and gradients i n temperature and density w i t h i n the shock compressed region. W h i l e this experiment represents the first successful effort in unraveling the atomic physics of strongly coupled plasmas, it is not without problems. T h e use of emission from the front side plasma for the absorption measurement precludes a determination of the true position of the unshifted photoabsorption edge. In the experiment, the best Chapter 1. Introduction 5 measurement of the unshifted edge position was still red-shifted by 3 e V from high resolution measurements on cold material [15], possibly a result of higher levels of x-ray heating than had been anticipated. More importantly, however, is the appearance of the ls-3p absorption line on top of the K-edge. A s this line appears, it becomes difficult to determine whether the absorption is due to the K-shell or due to the Sp vacancy, thus adding another level of uncertainty to the measurement. Finally, in an effort to characterize the shock, a one dimensional hydrocode which did not include any type of radiation transport physics was used. Moreover, questions regarding the accuracy of the high pressure equation of state of chlorine could expose problems w i t h interpreting the edge shift as due to x-ray or shock heating processes. Shock induced red shifts of 7 ± 2 e V i n the K-shell photoabsorption edge of a l u m i n u m were later reported by D a S i l v a et al. [16] and interpreted w i t h the solid state model of G o d w a l et al. [17] and the density functional theory approach of Perrot [18]. T h e use of a relatively thick, uniform target i n the experiment by D a S i l v a et al. [16] served two purposes. Target uniformity allowed plasma states to be characterized by shock speed measurements, while target thickness was chosen to reduce the level of radiative preheat. Emission from the front side plasma passing through the target provided the absorption measurement. Unfortunately, the limited amount of emission at the start of the laser pulse left the initial unshifted edge position unobserved. Moreover, the limited level of transmitted x-rays forced the data to be recorded only on Polaroid film, leading to only a semi-quantitative assessment of intensity levels and edge position. More definitive measurements of the K - s h e l l photoabsorption edge profile are needed. A l u m i n u m represents the ideal material for these experiments, since its high pressure equation of state is well established [19] and thus the shock compressed system can be accurately modeled. W e l l characterized plasma states are necessary to test atomic physics models. Chapter 1. Introduction 1.3 6 Present W o r k T h e structure and position of \ h e K-edge in a l u m i n u m (1560 e V in undisturbed material [20]) depends upon density, temperature, and ionization state [21], as well as effects such as continuum lowering [22] and electron degeneracy. T h e dense plasma behind the shock front is not hot enough to emit an appreciable flux of x-rays, and hence emission spectroscopy is not a viable technique for probing this region. A promising alternative is to use absorption spectroscopy and study the temporal evolution of the spectral profile of an inner shell photoabsorption edge. Ideally, the x-rays for the absorption measurement would be provided by a backlighter, a laser-produced plasma situated behind the sample which is to be probed. W e l l characterized plasma states, necessary for accurate atomic modeling, are produced by coupling tailored laser pulses to specific targets. T h e a i m of this thesis research is to provide the necessary background to carry out such an experiment. Target and laser pulse design is done with a one-dimensional radiation-hydrodynamics computer code. A non-local thermodynamic equilibrium ( N L T E ) model is employed to calculate ion level populations and emitted radiation. C o n t i n u u m radiation is transported through the target using multigroup diffusion [23], whereas line radiation is transported using the concept of photon escape probabilities [23, 24]. T h e piecewise parabolic method [25, 26] is used in conjunction w i t h equations representing the conservation of mass, momentum, and energy, to model shock propagation and hydrodynamics w i t h i n the target. T h e targets consist of C H - A l - S i sandwiches, where the C H and S i layer thicknesses are varied but the A l thickness is fixed at 3 microns. T h e C H layer on the laser deposition side provides pressure enhancement of the shock wave as it passes from the C H to the a l u m i n u m layer, as well as preventing the aluminum from releasing out the front side Chapter 1. Introduction 7 prematurely. T h e silicon layer allows quasi-steady state conditions to be achieved in the aluminum before the shock emerges from the rear surface. T h i s quasi-steady state is required for ~100 picoseconds, which, given the temporal resolution of modern streak cameras, should provide sufficient time for measurements. 1.3.1 Strategy to U n f o l d K - E d g e Dependencies T h e routes to decoupling the effects of density, temperature and ionization state on K edge energy are represented as three branches originating from a reference state in a density-temperature space (refer to Figure 1.2). T h e term 'reference state' is merely a label for the state w i t h which conditions on the three branches are compared in order to deduce the effects of different hydrodynamic conditions on the K-edge position. T h e branch for studying density effects includes the reference state plus one point at higher density (but the same temperature), created by launching multiple shocks into the target. Likewise, multiple shock compression is used to create a point of the same density as the reference point, but lower temperature. A greater average ionization resultant.from the higher temperatures created by large amplitude single shock waves allows the study of ionization state effects on K - s h e l l photoabsorption. Hence, the simulations have been undertaken to provide the.optimal pulse and target characteristics to achieve the above four states necessary for unfolding the individual dependencies. 1.3.2 Backlighter Sources T h e ideal backlighter would be spectrally flat (i.e. no spectral lines or emission edges), and as bright as possible in the region surrounding the K-edge, 1520-1600 eV. W i t h these characteristics in m i n d , a search for a suitable backlighter has been carried out. Spectral emission from 12 elements ( M g , Ge, Y , N b , A g , Sn, S m , T a , W , A u , P b , and B i ) has been recorded w i t h a flat pentaerythritol ( P E T ) crystal spectrometer. A n aluminum Chapter 1. Introduction 1 0.0 1 1 1 8 I 1 1 1 | 1 • 1 1.0 1 2.0 1 — i i 3.0 i — | — i i I i 4.0 1 — i — i — | — i — i — i — i — | i i i I i i 5.0 i i ' 1 — i I i 6.0 i 1 i 1 — i 1 — | i I . 7.0 . r . -. I 8.0 Compression Figure 1.2: A general compression(density)-temperature space used to illustrate the method of decoupling density, temperature, and ionization effects on K-edge position. Chapter 1. Introduction 9 spectrum superimposed upon the sample spectrum was used for energy calibration with an accuracy of ± 1 e V in the region around 1560 e V . A t o m i c number differences between the samples was less than ten, providing adequate coverage in the region of the periodic table studied. 1.4 Thesis Organization Chapter 2 reviews the important processes in laser-matter interactions, notably shock propagation and the emission, absorption and transport of radiation. A description of the non-local thermodynamic equilibrium model is also given. Chapter 3 describes the computer codes used to model the ablation process and optimize the laser pulse and target designs. Chapter 4 summarizes the results of the numerical simulations. Details of the experiment are given in Chapter 5, and the x-ray spectra are discussed i n Chapter 6. C o n c l u d i n g remarks are made in Chapter 7. Chapter 2 P h y s i c a l Processes in L a s e r - M a t t e r Interactions T h i s chapter presents a brief review of the physics involved in the formation of shock waves when laser light is focussed onto a solid target. Shock propagation in multilayered targets is discussed in terms of the Rankine-Hugoniot equations and the impedancemismatch technique. As well, emphasis is placed on describing the physics involved in the emission, absorption and transport of radiation through the plasma. 2.1 Laser D r i v e n Shock Waves 2.1.1 Shock Wave F o r m a t i o n Features of the ablation process along w i t h associated density and temperature profiles appear in Figure 2.1. For short wavelength (A < 1/um), medium intensity (I ~ few x 1 0 15 W/cm 2 ) laser radiation , the dominant absorption mechanism is inverse bremsstrahlung, or free-free absorption [28]. In this process, electrons oscillating in the electric field of the incident light wave collide w i t h ions, and thus convert their directed energy of motion into random thermal energy. T h e atoms of the target material are quickly heated and ionized, forming a dense plasma on the target surface. D u r i n g this initial plasma formation, there exists a dense plasma bounded by solid target on one side and vacuum on the other. This step profile in density quickly decays into a rarefaction wave, and ions are accelerated into the vacuum at the ion acoustic speed [29]. T h e ejected material forms a low density 10 Chapter 2. Physical Processes in Laser-Matter Interactions 11 POSITION Figure 2.1: Schematic of the hydrodynamics i n a laser irradiated target, p is the target density, and T and are the electron and ion temperatures. T h e critical density layer is denoted by n . e cr Chapter 2. Physical Processes in Laser-Matter Interactions (n < 1 0 e 23 12 c m ) , high temperature ( T > 100 eV) region called the coronal plasma (in - 3 e analogy w i t h the coronal layer of stellar objects), labeled as region 1 of Figure 2.1. T h e laser light penetrates the coronal plasma only to the critical density layer , the point at which the electron density yields a plasma frequency resonant with the laser frequency and the incident light is reflected. Analytically, the critical density is given by where eo is the p e r m i t t i v i t y of free space, m respectively, and u t is the laser frequency. e and e are the electron mass and charge, T h e reflection process can take place over spatial scales as short as a few laser wavelengths [30]. Note that the critical density is proportional to the inverse square of the laser wavelength, and hence shorter wavelength radiation propagates further into the target before being reflected. Energy transport beyond the critical density layer is accomplished by electron thermal conduction, x-ray radiation emitted by the hot corona, and suprathermal electrons[27] generated at the critical density layer. T h e steep temperature gradient which exists between the cold solid and the hot coronal plasma gives rise to a large thermal flux towards the interior of the target. T h i s occurs in region 2 of Figure 2.1, which is known as the ablation zone joins the cold solid at the conduction or ablation zone. T h e ablation front. Energy carried to the ablation front is balanced by a return flux of accelerated and heated target materials into the corona. T h e ablation front moves through the target at a constant rate provided that the mass ablation rate remains constant. Once the energy flux between the critical density layer and the solid target is established, shock formation arises as a consequence of conservation of mass, momentum and energy. T h e momentum of the ablated material must be balanced, and i n this way a shock wave is produced (region 3 of Figure 2.1). T h e amplitude of the shock wave w i l l depend Chapter 2. Physical Processes in Laser-Matter Interactions 13 upon the laser irradiance, and as the irradiance increases (e.g. the rising edge of a laser pulse) a series of shock waves of increasing amplitude are launched into the target. T h e shock wave is a compressional disturbance, and the first of this series w i l l travel at slightly over the local sound speed, compressing the target. T h e next shock wave i n the series will travel at the sound speed of the compressed material, and hence propagates faster than the previous. (Recall that the sound speed is given by c s = )• T h i s effect leads to a pile up zone, where multiple shocks of increasing amplitudes and velocities are propagating into the target. T h e latter waves w i l l catch up w i t h the earlier ones, leading to a steepening of the shock profile, as illustrated in Figure 2.2. T h e many small shock waves have thus coalesced into a s i n g l e well-defined shock front, characterized as a propagating discontinuity i n density, temperature and pressure. Constant ablation pressure sustains the shock, so that further acceleration of the shock front requires an increase in the ablation pressure. 2.1.2 Shock Compression of Solids and the Rankine-Hugoniot Relations Shock compression of solids is sufficiently different from the case of liquids and gases to warrant a brief review. A lucid account may be found in Zel'dovich and Raizer [31]. T h e essential difference between a solid and a liquid (or gas) is the strong interparticle interaction present i n a solid. Atoms in condensed matter are bound together by an elastic potential which resists expansion and contraction, while at the same time m i n i m i z i n g the potential energy of the system. T h i s interparticle. force is effective over distance scales of only Angstroms, and yet accounts for the ability of a solid to support shear stresses, whereas a liquid or gas cannot. T h e interparticle forces present in a solid are usually accounted for by the inclusion of pressure and energy terms arising solely from this elastic force. In this way, the .total Chapter 2. Physical Processes in Laser-Matter Interactions Figure-2:2: Illustration of the shock build up process in a laser driven shock Chapter 2. Physical Processes in Laser-Matter Interactions 15 pressure and internal energy is given by (2.2) P = P + P c E = E C t + (2.3) E t where P and E are the elastic force terms and Pt and Et are the thermal terms. Recall c c that pressure arises from the transfer of momentum of particles participating i n thermal motion. T h e non-thermal pressure and energy terms arise from the interparticle binding forces alone, and remain finite even as the temperature approaches absolute zero. A n . i m p o r t a n t aspect i n the shock compression of solids is the distinction between strong and weak shock waves. In the weak shock l i m i t , the principal contribution to the internal energy is E , and the shock speed is related to the compressibility of the c material at standard conditions. T h i s approaches the l i m i t of an acoustic wave, and the compression of the material is of the order of a few percent. T h e pressure range for which shocks in solids are weak is up to ~0.1 M B a r . In the strong shock l i m i t , however, the energy imparted to the material by the shock wave far exceeds the material's elastic potential. T h e resulting shock heated material behaves as an ideal gas, and the ideal gas equation of state may be applied. In this case, the thermal term E t dominates. Pressures of the order of ~100 M B a r are needed to generate strong shock waves in solids. For pressures intermediate between these limits, the thermal and elastic terms contribute approximately equally. T h e locus of thermodynamic states which can be reached in the shock compression process is determined by the equation of state ( E O S ) of the material. T h e E O S is generally represented as a two-dimensional slice of the three-dimensional surface of pressure, density and temperature. Reversible thermodynamic processes are represented as continuous paths between initial arid final states on the E O S surface. T h e irreversible process Chapter 2. Physical Processes in Laser-Matter Interactions 16 of shock compression is represented as a j u m p discontinuity from initial to final states on the E O S surface. For a given initial state, the properties of the final state are constrained by the conservation of mass, momentum and energy across the shock front (refer to Figure 2.3). Mass and momentum conservation are expressed as V D-u V D 0 ^ ^ Du P - P = — 0 where D is the shock speed and u is the particle speed behind the shock front. (2.4) , „ (2.5) s The equation for the conservation of energy across the shock front relates the internal energies of the initial and final states Eo, E , to the pressure and volume jumps across the discontinuity E - E 0 = ±x(P-P )(V -V). 0 0 (2.6) Equations 2.4- 2.6 form the Rankine-Hugoniot relations, and can be derived from completely general arguments. T h e y are independent of the aggregate state of the material through which the shock propagates, and are valid for both spherical and planar shock waves [32]. N o assumptions of continuity of flow variables is made i n their derivation. These relations can be obtained in a mathematically more formal manner by considering the shock front discontinuity as the l i m i t i n g case of very large but finite gradients. One then integrates the hydrodynamic equations of motion across this layer as the thickness tends to zero. T h e Rankine-Hugoniot relations constitute a system of three equations w i t h five unknowns. If the equation of state, P = P(p, E) is known for the material, it may be added as a fourth equation. T h i s allows the system to be reduced to one equation w i t h one free parameter, which yields the locus of thermodynamic states (i.e. the Hugoniot Chapter 2. Physical Processes in Laser-Matter Interactions 17 A Pi P 1 E 1 - Po D P o E 0 Figure 2.3: Schematic of the shock discontinuity. P is pressure, p is density, E is total internal energy, D is shock speed and u is particle speed. Variables subscripted w i t h a 0 represent the undisturbed solid. Particle speed ahead of the shock front is zero. Chapter 2. Physical Processes in Laser-Matter Interactions 18 curve) reached by a single shock for specific initial conditions. In Figure 2.4, a qualitative example of the Hugoniot, isotherm and isentrope are shown for comparison. 2.1.3 I m p e d a n c e - M i s m a t c h Technique T h e impedance-mismatch technique was first demonstrated for laser-target interactions by Vesser et al. [33]. Subsequent studies by Holmes et al. [34] and Cottet et al. [35] have verified the usefulness of this technique for pressure enhancement of shock waves at the interface of dissimilar materials. T h e important parameter characterizing the process is the material-dependent shock impedance, pU . s W h e n the shock wave reaches the interface of the two media, reflected and transmitted waves are generated. Whether these waves are reflections or rarefactions depends upon the relative shock impedances of the two materials. Specifically, i f P of the incident pressure pulse, and Pi, P 2 0 is the amplitude those of the reflected and transmitted pulses (refer to Figure 2.5) then [36] • p _ 1 ~ {pUsh ~ (pUs)l 7~TT\ 1 I TT\ p X r ° , \ x Z - ' ) {pU )2 + {pU )i s s 2{pU ) {pU )2 + [pU )i s s 2 s where the subscripts 1 and 2 denote the regions on either side of the interface. Con- servation of momentum and energy require that the pressure and particle velocity be continuous across the interface. Pressure enhancement can be achieved by propagating a shock into a material which has a higher shock impedance than the material from which the shock originated. C o n versely, impedance mismatch explains the low pressure release wave associated w i t h shock breakout. W h i l e technically the impedance mismatch technique refers to shock propagation Chapter 2. Physical Processes in Laser-Matter Interactions A Figure 2.4: Hugoniot (solid), isotherm (dot-dash), and isentrope (dash). Chapter 2. Physical Processes in Laser-Matter Interactions 20 INTERFACE Pi P 1 E P2 1 P E 2 2 Figure 2.5: Schematic of the impedance mismatch technique. p Pi, E U i, and U i are the density, pressure, total internal energy, shock speed, and particle speed in the two regions. At the interface, P + P = P and U = U . h Q x 2 pl p2 h S P Chapter 2. Physical Processes in Laser-Matter Interactions 21 through dissimilar materials, the idea can also be applied to uniform targets. For example, consider the case of two shock waves propagating through a uniform target. T h e second shock travels in the previously compressed and heated material, and is just about to overtake the first. Just after coalescence, a shock wave of a new pressure and speed is formed, which then traverses the undisturbed material. propagating rarefaction is generated. A t this point a counter- T h e rarefaction is due simply to the fact that there now exists a shock wave of a new pressure, and hence speed, traveling in a material w i t h different density. Arguments using this idea of shock impedance can be used to explain the 'ringing' phenomena associated w i t h shock compression of multi-material targets, and w i l l be evident when inspecting the hydrodynamic profiles inside targets presented in Chapter 4. 2.2 Atomic Physics 2.2.1 Introduction R a d i a t i o n from laser produced plasmas can be identified as originating from essentially three different regions w i t h i n the target [37]. Emission from the high temperature, low density coronal plasma consists mainly of faint continuum emission. T h i s is due p r i m a r i l y to the fact that the most abundant ionic species are either fully ionized or are ground state H-like or He-like ions. T h e greatest part of the radiation originates from the region between the critical density layer and the ablation front. Here, the plasma density and temperature are high enough to cause intense x-ray emission. T h e remaining target material, including the plasma state behind the shock front, is not hot enough to emit appreciable amounts of x-ray radiation, but rather, emits as a blackbody source [37]. For medium to high atomic number elements, the emission is dominated by boundbound transitions [38, 39], whereas for low atomic numbers, the emission tends to be Chapter 2. Physical Processes in Laser-Matter Interactions 22 composed p r i m a r i l y of recombination continuum [40]. T h e relatively low temperatures needed to fully ionize low atomic number materials leads to fewer bound electrons from which line emission can originate. A s well, in low atomic number plasmas, bound-free and free-free emission dominate over bound-bound line radiation due to electron collisional quenching of the line radiation. T h e population density contributing to continuum radiation lies below.the collisional quenching l i m i t [40], and therefore continuum radiation is unaffected. Furthermore, in high atomic number elements as much as 50% of the total incident laser radiation can be converted to x-rays [42, 72], making the plasma an intense source of radiation. T h i s radiation represents a significant energy transport mechanism and can lead to several effects. X - r a y s can be deposited near the ablation front, enhancing the ablation process. X - r a y radiation also represents an energy loss mechanism, which can lead to a reduction i n ablation pressure and hence shock speed. More importantly, however, x-fays propagating into the target can preheat the material before the shock arrives, complicating the characterization of the shocked state. T h e first step to quantitatively describe the radiation emission is the development of a model to calculate the level populations of the ions i n the plasma. T h i s requires knowledge of all the atomic processes and associated rate coefficients important in the regime of interest. T h e processes considered important in this work are described below, followed by a discussion of the non-local thermodynamic equilibrium model used. 2.2.2 A t o m i c Processes Collisional Ionization Collisional ionization N (j) z + e=* N (k) z+1 + e+ e (2.9) Chapter 2. Physical Processes in Laser-Matter Interactions 23 occurs when an ion of charge Z, state j, becomes ionized through a collision with a free electron. State k is some excited state or ground state different from state j. Three-Body Recombination As the inverse process of collisional ionization, three-body recombination consists of an electron being captured by an ion. N (j) z + e + e=> N -\k) z + e. (2.10) The non-recombining electron satisfies momentum conservation. Three-body recombination rates are calculated using the method of detailed balance with collisional ionization rates and the resulting rate coefficients are proportional to the square of the electron density. This process can generally be neglected for plasma densities below 10 16 cm - 3 [43]. Radiative Recombination Radiative recombination is the emission of a photon upon the capture of a free electron by an ion. N \j) +. => N ~\k) + hu z z e (2lll) This process represents the major source of x-ray radiation emitted over a large range of plasma conditions. Radiative recombination rates are generally obtained through detailed balance with photo-ionization rates. i Chapter 2. Physical Processes in Laser-Matter Interactions 24 Photoexcitation and Photoionization A s the plasma density increases to the point where reabsorption of radiation becomes optically thick. In this scenario, population distribu- important, the plasma is known as tions in the various excited levels of the different ionization stages are altered by photoexcitation and photoionization, which i n turn affects the spectral emission [96]. T h i s introduces a coupling between the radiation field and the population distribution which must be considered to obtain accurate spectral information. Physically, the processes of photoexcitation and photoionization are represented as N (j) z + hu =» N (k) z (2.12) and N (J) + hi/=> N {k) + e , z respectively. z+l A n y model which ignores these effects is only treating the (2.13) optically thin case, and may not provide an accurate calculation of radiation emission or transport levels. Collisional Excitation and De-excitation Collisional excitation and its inverse process of de-excitation, N (j) z + e & N \k) z + e (2.14) occur when electron-ion collisions induce a transition in the state of the bound electrons. Detailed balance may be used to calculate the de-excitation rate when the excitation rate is known. Chapter 2. Physical Processes in Laser-Matter Interactions 25 Radiative Transistions A l s o termed bound-bound transitions, radiative transitions N (j)=>N {k) z + hv z (2.15) are responsible for the line emission of atoms. B o u n d electrons make transitions to lower energy levels, releasing the excess energy i n the form of radiation. Dielectronic Recombination T h e three-step process of dielectronic recombination starts w i t h the capture of an electron by an ion into a level of large principal quantum number. T h e excess energy excites a bound electron N (j)+e^N - (f,k') z z . l (2.16) The state consisting of the two excited electrons is unstable, and can undergo radiative decay N ~ (f, z 1 k') N ~ (k) z 1 + hu . (2.17) Generally, the stabilizing radiative process is the decay of level j' into an excited state which lies below the ionization l i m i t . T h e outer electron causes the radiated photon wavelength to be longer than that of the unperturbed resonance, giving rise to satellite lines in atomic spectra. 2.2.3 Non-Local Thermodynamic Equilibrium Model The shaded region of Figure 2.6 shows the extent of parameter space accessible in laserproduced plasma experiments. It is immediately apparent that only a very general model may be used to adequately describe the atomic processes present over the entire range Chapter 2. Physical Processes in Laser-Matter Interactions 26 TEMPERATURE [«V] Figure 2.6: Shaded region represents the parameter space accessible in laser-produced plasma experiments. A p p l i c a b i l i t y of certain models is shown. For example, the line labeled L T E Z=5 corresponds to the range of validity of the local thermodynamic equil i b r i u m model for B o r o n . Chapter 2. Physical Processes in Laser-Matter Interactions 27 of density and temperature. T h e coronal model, which assumes that ionization is determined by a balance between collisional ionization and radiative recombination, is only applicable where the radiative decay rate is much greater than the collisional decay rate 52A ».n 'E X .. jk e f k<j (2.18) jk :k<j M a n y atomic processes are neglected with this model, making it a crude approximation except for very low density plasmas. T h e local thermodynamic equilibrium ( L T E ) model, on the other hand, is valid only at sufficiently high densities [44]. In this model, the charge state is determined from the Saha-Boltzmann equation [45] from which the ion level populations are calculated using the Boltzmann relation Here, wf and xf a r e the statistical weight and ionization energy, respectively, of level i. T h e high electron densities needed to make the assumptions of the L T E model valid preclude its use for modeling the physics of most plasmas. . Instead, a collisional-radiative equilibrium ( C R E ) or non-local.thermodynamic equilibrium ( N L T E ) model must be employed. T h i s model is applicable over all density and temperature regimes, and gives equivalent results to the coronal and local thermodynamic equilibrium models in the limit of low and high electron densities, respectively. In this 'non-local thermodynamic equilibrium model, an atomic rate equation of the form ^ L U = Y,W N -Y:W. }N kj l k k l j . is written for each ground and excited state included in the model. 'W~ k (2.20) k Here, W j and are the rates for populating and depopulating, respectively, the level Nj. k These coefficients include all of the physical processes described in the previous section. If the Chapter 2. Physical Processes in Laser-Matter Interactions 28 hydrodynamic characteristics of the plasma change on a time scale greater than that of the important atomic processes, E q u a t i o n 2.20 reduces to the steady state form o =E^*-Ew7* ^ •• (- ) 1 k 2 21 k A n estimate of the time required to reach steady state has been suggested [46], t where the laser wavelength X ss , (2.22) is in microns. For 0.53 pm light, t L T h i s estimate for t = 8.98 x 1 0 - A £ 1 0 ss ss ~ 200 picoseconds. allows the rate equation 2.20 to be used in for few-nanosecond time scale simulations, but could pose problems for the shorter time scale simulations (~ 1 nanosecond) presented in Chapter 4. In principle, the non-local thermodynamic equilibrium model can be used for any plasma conditions provided accurate atomic rates for all processes are known. are cases, however, where accurate atomic rates are unavailable. There In this model, the processes of collisional excitation and collisional ionization are included despite possibly large uncertainties in their atomic rates. Fortunately, Salzmann [47] has concluded that parameters sensitive to the average ionization , such as internal energy, free-free and free-bound radiation, are generally insensitive to the rate coefficients. 2.3 2.3.1 R a d i a t i o n Emission, A b s o r p t i o n and Transport Bremsstrahlung Emission T h e r m a l bremsstrahlung emission arises from the acceleration of an electron in the electrostatic field of an ion, w i t h the electron as the primary emitter. T h e power radiated per unit volume per unit frequency interval is given by /•oo I {y) fJ = hvY,n N z e ry I J 0 v f(E)do dE e ff (2.23) Chapter 2. Physical Processes in Laser-Matter Interactions 29 where v and f(E) are the electron velocity and distribution functions, respectively. T h e e differential cross-section for bremsstrahlung emission, doff, is usually expressed as the classical cross-section times a correction factor, gff, called the Gaunt factor. Integration of Equation 2.23 over a Maxwellian electron distribution of temperature T yields e 6 8 x 10~ 7 //W = — ^ 1 / 2 47 e£ %f n le e*p(-hv/T ) z2N e (2.24) Z where / / / has the unit of [W cm"" H z ] . 3 - 1 T h e inverse bremsstrahlung absorption coefficient, Kff(u), can be evaluated by app l y i n g KirchhofT's law [48] Zff(v)=*ff(y)B(v) (2.25) where the bremsstrahlung emission per unit solid angle, Cff(v), is calculated from E q u a - tion 2.23 and B(u) is the Planck blackbody function „, 2/uA N 1 - 1 ' = — eMhv/kT) 2.3.2 ' ( 2 2 6 ) Recombination Emission In the process of recombination, a free electron of energy E is captured by an ion into a bound level of ionization energy x f emitted. - 1 , and a photon of energy hv = (E + Xn~ ) * l s T h e energy of the emitted photon is thus a function of the electron energy E, and only photons with energy E > Xn~ l a r e generated. T h i s is the source of the recombination edges observed i n free-bound continuum spectra [49, 50]. U s i n g the hydrogenic ion cross-section [51], the power radiated per unit volume per unit frequency interval by free-bound radiation is W = 6 8x1 0 n J2Z>N eM~hv/T ) z le f -v-Z-i - 4 7 e e Z k^-g^Mx^/Te) £ n>n ee n m i n I (2.27) Chapter 2. Physical Processes in Laser-Matter Interactions in units of [W c m [52] and £ n - 3 H z ] . In E q u a t i o n 2.27, g 30 is the average free-bound Gaunt factor - 1 fb is the number of available states for recombination. T h e sum n > n in is m taken over the different shells for which hv > Xn' • 1 The photoionization absorption coefficient Hf (u) can be calculated using detailed b balance in conjunction w i t h Kirchhoff's law and the Planck function. T h i s leads to = 1.3 x where U f Z-l\2 v u C - v,-3. 1 — exp Z-1, ? f (hv (2.28) is the partition function and /?/(, is i n c m . z 2.3.3 TJZ i oTKTZ« 1^ - >\X ^ ^) -,40 - 1 Line Emission Line emission results from the transition of an electron between bound energy levels of the ion. Specifically, for upper level u and lower level I, the radiated power per unit volume per unit frequency interval is simply Ibbi") = N A hu (j) (u) , u also in units of [W c m - 3 ul Q (2.29) e H z ] . A \ is the transition rate (Einstein-A coefficient) between - 1 u levels u and I, 4>e(v) is the normalized emission line profile function, and Nu is the population of the upper level. U s i n g the relationships between the Einstein A and B coefficients [48], the line absorption coefficient can be written •' where (j> { ) v a ^ - '(j,- k) '' KM= A N ! M ) - (2 30) is the normalized absorption line profile function, and gu and gi are the, statistical weights for the two levels. Chapter 2. Physical Processes in Laser-Matter Interactions 2.4 31 R a d i a t i o n Transport T h e simplest radiation transport problem to solve is that of the optically thin plasma, where reabsorption of radiation is neglected. In high-density laser-produced plasmas, however, this approximation is often invalid. L i n e intensities and line shapes can be significantly modified by absorption w i t h i n the plasma. More importantly, though, atomic level populations are modified by photoexcitation and photoionization, leading to a coupling between the radiation field and the atomic level populations. If the radiation field is assumed to adjust instantaneously to changes in the plasma temperature and density, then the radiative transfer equation [53] takes the particularly simple form dl {s) u *«/( ) S ds where t v K„(S)I (S) U (2.31) and K are the total emission and absorption coefficients, and include a l l of V the processes described in the preceding section. E q u a t i o n 2.31 is simply the statement that the intensity increment dl over the path ds is equal to the increase e„(s)ds v emission minus the decrease K (s)I (s)ds u l/ due to due to absorption. Consider the case of radiation transport in planar geometry as depicted in Figure 2.7. U s i n g [i = cos 6 (9 being the angle between the path of an emitted photon and the x axis) and the optical depth r v (2.32) along w i t h the source function S u (2.33) allows the radiative transfer equation to be written as dI {T ,n) u v dr v = S -IV V l (2.34) Chapter 2. Physical Processes in Laser-Matter Interactions Figure 2.7: R a d i a t i o n transport in planar geometry. 32 Chapter 2. Physical Processes in Laser-Matter Interactions 33 Equation 2.34 can be integrated to obtain the formal solution of the transfer equation + jT exp { ^ ~ ) IU(T*, A*) = MO", p) exp M ^ d < . . ( 2 3 5 ) The special case of a constant source function S yields the solution v J „ ( T „ , p) = S + exp x [7,(0, p) - S„] v . (2.36) T w o l i m i t i n g cases exist for this solution. In the optically thick l i m i t (T„ > > 1), the emitted intensity reduces to I„{T ,p) = S v v , (2.37) while the optically t h i n limit (T„ < < 1) yields I (T ,p) = (\-T /p)I {Q,p) V v v v + is v . (2.38) In the numerical calculations presented in Chapter 4, the optically thick l i m i t described above is obtained by the inclusion of photoexcitation and photoionization processes. In calculating the energy transfer i n planar geometry, it is useful to consider I as v consisting of forward (/+) and backward (/~) components, subject to the boundary conditions /+(0,AX) = 0 7;(TJ,//) = 0 . , (2.39) The energy flow per unit time per unit area perpendicular to the r axis is then given by F {T ) V v = F+- F~ (2.40) where the radiative fluxes are defined as F* = 2TT [ ±l Jo I±pdp . (2.41.) T h e absorbed power per unit volume is then calculated w i t h E (x) abs HF =^ . (2.42) Chapter 2. Physical Processes in Laser-Matter Interactions 34 For numerical calculations it is convenient to consider energy absorption as occurring in one of n finite w i d t h cells using the expression Eg" = F„{T?) - FrW ) . 1 • (2.43) T h e energy absorbed in cell n due to emission from cell 0 (including absorption i n cell 0) is approximated with E a bs n = fAx [E (T:-T»/2)-E (T^-T°j2)} T 0 2 where I is the radiated power per unit volume and E 2 0 integral. 2 (2.44) is the second order exponential T h e factor of 1/2 accounts for the assumed isotropic forward and backward emission. E q u a t i o n 2.44 is used extensively in the numerical calculations discussed in the following chapters. Chapter 3 N u m e r i c a l Simulations T h i s chapter presents a description of the numerical codes used to simulate both the hydrodynamics of the laser ablation process and the emission, absorption and transport of radiation through the target. A survey of the relevant input parameters used for the simulations is presented at the end of the chapter. 3.1 Introduction T h e hydrodynamics of the laser ablation process are modeled w i t h the Laser-Target Code ( L T C ) , a substantially modified version of the M E D U S A inertial confinement fusion code [54]. L T C has been described in considerable detail elsewhere [55] and therefore only a brief overview of the important physics is presented here. T h e emphasis in L T C is placed on the accurate modeling of the shock compression process only, and so important physical processes such as radiation emission and transport have been neglected. These more complex processes are treated in H Y R A D , a code which couples the accurate hydrodynamic modeling of L T C w i t h a detailed calculation of the emission, absorption and transport of radiation through the target. H Y R A D uses a non-local thermodynamic equilibrium ( N L T E ) model to determine the ion state populations and emitted radiation. Line and continuum radiation can be transported as in an optically t h i n plasma, or can be coupled to the population distribution to produce a simulation corresponding to an optically thick plasma. Detailed tests have been performed for the N L T E model [56], and thus only the important aspects of the code are reviewed 35 Chapter 3. Numerical Simulations 36 here. Computer codes are necessary in modeling laser ablation, as the many simultaneous physical processes preclude the use of simpler analytic formulae [28, 57]. Simulations, however, are only approximate, and involve assumptions and simplifications. Numerical results should not be taken as the final word, as only comparison w i t h experiment can validate the level of complexity used in a simulation. 3.2 Laser-Target C o d e ( L T C ) 3.2.1 Physical Content in L T C T h e one-dimensional hydrocode L T C solves the plane-parallel fluid equations for a singlefluid, two-temperature plasma. Electrons and ions are assumed to have the same velocity, but each species maintains a characteristic temperature. T h e most complete description involves a radiation field at a t h i r d temperature in equilibrium w i t h the electrons and ions. T h e fluid equations are solved in the Lagrangean formalism, where the calculation follows the time evolution of the individual fluid elements. T h e time (t) and Lagrangean coordinate (m) are regarded as the independent variables. Specifically, the Lagrangean coordinate m(r, t) (in k g / m ) is defined in terms of the laboratory coordinates (r,t) 2 and the density profile p(r, t) by [25] (3.1) where Ri(t) is the position of the free (rear) surface and r is the position of the cell under consideration. E q u a t i o n 3.1 is simply an expression for the total mass of the Lagrangean cells situated between Ri(t) and r . T h e physics associated with the laser ablation process and subsequent hydrodynamic Chapter 3. Numerical Simulations 37 motion of the material is contained w i t h i n a set of coupled differential equations representing conservation of mass, momentum and energy. In terms of m and t , these equations are dv du dt dm du ( d(P + Pi) e dt dt + _ Q (3-3) dt + 1 dt dV + J (3.4) dt U 1 dF _ c (3.2) e P 2 Q dm d(U/ ) dt _ dt ~ d(fU) dt dm (3.5) ^ _dF + dm + C K p pF X c T h e principal dependent variables which are solved for are specific volume V = 1/p, velocity w, material internal energies E ei (3.6) (3.7) fluid Ei (electron and ion components), radiation energy density U, and radiation flux F. Equations 3.2 and 3.3 represent fluid continuity and momentum conservation, respectively. T h e momentum equation considers the total hydrodynamic pressure, as indicated by the sum of the electron and ion components. R a d i a t i o n pressure is negligible com- pared to typical shock pressures, and is thus excluded from the momentum equation. Equations 3.4 and 3.5 express energy conservation, for electrons and ions, i n terms of the first law of thermodynamics. T h e energy source term Q is composed of the various energy transport and deposition processes included in the model. Equations 3.6 and 3.7 describe the interaction of the radiation field w i t h the target material. T h e Eddington factor / in E q u a t i o n 3.6 is restricted to the range | < / < 1. T h e detailed derivations of Equations 3.6 and 3.7 can be found in the work of M i h a l a s and M i h a l a s [53]. T h e radiation energy density and flux equations are included here Chapter 3. Numerical Simulations 38 only for completeness as radiation emission and transport w i t h i n the target is modeled in substantially greater detail w i t h the N L T E model described in the following section. A number of secondary quantities are required to complete the description, and these can be obtained once the principal dependent variables are defined. Most important in the prescription is the equation of state. Generally, the E O S is expressed in the form T = T (V,E ), P = P (V,E ) (3.8) T = T (V,E ), P = P (V,E ) (3.9) e e i e i e i e i i e i where V is the specific volume and T , T j and P ,Pi are the electron and ion temperae e ture and pressure components, respectively. Other important quantities include the ion density, P Ui — (3.10) Am P where A is the average atomic mass number of the target material and m p is the proton mass. T h e electron density is given by n =< Z > x m , e (3-11) and the C o u l o m b logarithm [54] is given by I n A = 1 6 . 3 4 - l n ( T / n J / < Z >~ ) 3 2 1 2 e l . (3.12) T h e average ionization < Z > can be obtained by either interpolation from an appropriate table (such as the S E S A M E [89] database) or by calculation w i t h an atomic model. The ionization for certain regions w i t h i n the target is calculated w i t h the N L T E model. Values for the remainder of the regions are obtained from a table. T h e laser wavelength A is needed to define the critical density 1.1 x 1 0 A 2 15 (3.13) Chapter 3. Numerical Simulations where n CT has the unit of m 39 and the wavelength is in meters. - 3 Other necessary data such as electrical and thermal conductivities may be obtained from tabular formats so that the hydrocode can be applied equally well to a variety of target materials. T h e source terms which enter the conservation of energy equations are given by Qi = Hi- Q = (3.14) ei +X (3.15) where H is the heat flow due to thermal conduction, K is the rate of energy transfer e H E +K K ei +A las rad EI Ai between the ions and electrons, and term as is the rate of absorption of laser energy. T h e X d accounts for radiation transport. ra Electron and ion heat conduction, H and H , are calculated w i t h the classical heat E { conduction model H = -V-KVT P (3.16) where K is the appropriate thermal conductivity and / is the flux limiter. For electronic heat conduction, K can be interpolated from a table (such as S E S A M E [89]) or calculated w i t h the Spitzer expression [60] T K = 1.995 x I O 0.095 ' °- ' (1 + 0.24Z) - 9 Z E + 2 4 5/2 e (3.17) ZlnA T h e ion thermal conductivity is calculated w i t h the Spitzer equation [60] j,5/2 m = 4.3 x 1 0 - 1 -fZlnA . 2 v (3.18) ' Electron-ion energy exchange is governed by the term [54] K = 0.59 x 1 0 ~ n ( T j - T )T^ M~ Z 8 EI 2 e E V 2 In A (3.19) where M and Z are the atomic mass number and charge state, respectively, of the ions. Chapter 3. Numerical Simulations 40 T w o methods can be employed to calculate the absorption of the incident laser radiation. T h e first method is adequate for ~ 5 0 0 picosecond or longer time-scale simulations, whereas the second method is better suited for few-hundred picosecond or shorter timescale calculations. Results from these two methods are discussed in greater detail i n the next Chapter. In the first method, laser energy is absorbed v i a inverse bremsstrahlung up to the critical density layer. T h e absorbed laser power at position (r, t) is given by A (r,t) =--^ (3.20) las where $ t ( r , t) is the laser power reaching point r at time t. U s i n g the absorption coefficient of Stallcop and B i l l m a n [61] p 2 13.51 a A A = Z\nA (1 ( 1-- /?)V2 3 / 2 22 T ^ ( 3 2 1 ) where /3 = n /n , $ ( r , t) is calculated w i t h e cr L rRo ®L(r,t) = <l> {R ,t) 1 — exp J —a(r',t)dr' L 0 (3.22) Laser energy which is not absorbed upon reaching the critical density layer is reflected back out of the plasma, where it can be re-absorbed on the way out. Laser energy escaping the plasma is lost. In the second method, known as the electromagnetic wave solver, absorption of laser energy is accomplished by solving the Helmholtz equations [62] for electromagnetic waves in a plasma [63]. T h e initial cold solid and subsequent hot plasma are characterized by a temporally and spatially dependent complex dielectric function 247T e(z, i ) = 1 + : — x . a ( c j ) where o(u>) is the electrical conductivity. (3.23) A t each time step of the calculation, the Helmholtz equations are solved numerically to satisfy boundary conditions corresponding Chapter 3. Numerical Simulations 41 to incident and reflected waves i n vacuum, and an evanescent wave inside the target. T h e solution yields the complex electric field amplitude, from which the reflectivity and energy deposition rate < E • J > = Re(o)\E\ /2 2 (3.24) are calculated. T h e energy deposition rate is then incoporated into the hydrodynamic equations as the heat source term T h e radiation source term X (t) 3 Tad A[ . as for cell j is calculated w i t h the N L T E segment of H Y R A D . T h i s is discussed in greater detail in the following sections, but briefly, T h e terms E £ and E™ account for line and continuum radiation deposition respectively. l 3.2.2 Numerical Methods in L T C T h e zoning scheme used i n L T C is similar to that of M E D U S A , consisting of a mesh of cells of mass A m . Associated w i t h each cell is a cell-centered average value for each principal dependent variable. Solution of the fluid equations is then obtained in two phases for each time-step. In the first phase, the piecewise parabolic method [25, 26] is used to solve a R i e m a n n shock tube problem at the interface of each cell, thus advancing the hydrodynamic motion of the mesh. T h i s defines the quantities r and u, and hence V and p , at the advanced time level. T h e Crank-Nicholson [54] time-centered differencing scheme is then used to solve the energy equations i m p l i c i t l y for E e and E^ Solution of these equations is done iteratively to account for the nonlinear dependencies of the source terms and equation of state on the thermodynamic variables. One to three iterations are sufficient to achieve a convergence accuracy on the order of one percent. Chapter 3. Numerical Simulations 3.2.3 42 C o m m e n t on L T C T h e principle emphasis of L T C is to model the basic laser ablation process. Accordingly, several mechanisms have not been treated. H o t electron generation and transport has been ignored since for irradiances < 1 0 1 5 W / c m , this process may not account for a 2 significant fraction of energy transport [58, 59]. T h e various parametric processes (e.g. stimulated B r i l l o u i n scattering, stimulated R a m a n scattering, two plasmon decay, etc.) which can develop i n the coronal plasma are also not considered. Most importantly, L T C does not contain the level of sophistication necessary to model the atomic physics of the coronal plasma. For this, the N L T E package described i n the previous chapter is coupled to L T C . T h e next section discusses the details of the N L T E model developed. 3.3 3.3.1 Physics Content i n the N L T E M o d e l C a l c u l a t i o n of State P o p u l a t i o n T h e collisional-radiative equilibrium model described i n the previous chapter is used to calculate the i o n state populations. T h e population of a given ion state N (j) z is calculated w i t h a rate equation of the form + neN ^(J> 4i,r + = [neN - (j")S& » z z l d £ N ^A z + n ] T N (t)X^ + n £ z tj e i>j -N {j)[n S f> ' z z e J z e i>j + n af +£ A +n £ e " .(3-26) JV (i)X ] z y i<j tj +n £ X] e i<j where Z e i<j tj i>j S^_ y, is the collisional ionization rate from N ~ {j") to N (j), a is the sum of z x z X radiative, dielectronic and three-body recombination rates, a =a R +a D I + na , e 5B (3.27) Chapter 3. Numerical Simulations 43 Aij is the spontaneous decay rate from level i to j, and X - , X ^ y 1 are the collisional excitation and de-excitation rates from level i to j, respectively. C o u p l i n g between all states is permitted w i t h i n the model, but for the sake of computational speed, only coupling between an excited state and the ground state of the next ionization stage is considered, as illustrated in Figure 3.1. T h i s approximation is valid for electron densities less than 1 0 22 cm - 3 [64]. For the near solid-density plasmas encountered i n the following simulations, the coupling should also include non-adjacent ionization stages. However, a coupling scheme this complete would be too computationally intensive to be feasible, and so the above approximation must be employed. The average ionizations predicted by this N L T E calculation and a code which couples select non-adjacent ionization stages are examined in Section 3.5. T h e solution of the rate equations can be further simplified if the characteristic time scale for changes in the plasma temperature and density is assumed to be long compared to the time scale of the atomic processes. T h i s allows the steady state approximation dN (j)/dt — 0 to be used. z Once the coefficients for the rate equation have been calculated, the population dist r i b u t i o n of the ions is determined. T h e electron density dependence of the equations requires iteration of the solution until convergence to w i t h i n a specified amount is achieved." 3.3.2 Rate Coefficients Ionization T h e ionization rate coefficient is calculated from the result of Landshoff and Perez [65] (3.28) 0.915 . L(l + 0 . 0 6 4 * T / x f )* e where S l j ' 1 1 0.42 (1 + 0MT /xf) \ 2 e has units of [cm /s] and f is the number of outer shell electrons. 3 Figure 3.1: Level coupling scheme used in this work. S is ionization, a is radiative and three-body recombination, a is dielectronic recombination, X and X~ are collisional excitation and de-excitation, respectively, and A is the spontaneous decay. l D I Chapter 3. Numerical Simulations 45 Three-Body Recombination B e i n g the inverse process of collisional ionization means that the three-body recombination rate can be calculated by the method of detailed balance. Thus, where g(Z, 0) is the multiplicity of the ground state and g(Z — the state formed by the recombination process. (OSZB)Z O 1 j n a s is the multiplicity of units of [cm /s]. 6 Radiative Recombination T h e radiative recombination rate, in [cm /s], is calculated using the method of Seaton 3 [66], (<XR)zfl = d 5 2 x l(T Z(£ 1 4 1 / 2 ( 0 . 4 3 + 0.5 In 0 + O . 4 7 0 _ 1 / 3 ) , (3.30) where <f> = xf /(kT ) . e Dielectronic Recombination T h e dielectronic recombination rate, in [cm /s], is calculated w i t h the Burgess and Merts 3 [67] prescription 2.4 x 1 0 - (« /)K''°= ' ^,/2 C 9 _ ^ , (-E* WEWJl-P^) (3.31) where = fiff? ^ B{z) z E* P- ) 20 32 Xij 1 + 0 . 0 1 5 Z / ( Z + 1)^ (3.33) 3 and A(Z, j) = x l l{\ + 0.105a; + 0.015a: ) An = 0 = 0.5z / ( l + 0.21a: + 0.30a; ) An ^ 0 x 2 1 / 2 2 2 (3.34) Chapter 3. Numerical Simulations 46 where the definition x = (1 + Z)*f- has been used. In these equations, Xij a n d hj a r e the excitation energy and oscillator strength respectively of the i —> j transition of the recombining ion of charge Z. Spontaneous E m i s s i o n T h e spontaneous emission rate is derived from the oscillator strength. ^% 1.5 xV' For transition ' (3 35) where Qj and <?, are the statistical weights of the upper and lower levels, respectively, is the oscillator strength of the transition, and A is the wavelength of the emitted photon. Collisional Excitation and De-excitation T h e collisional excitation rate is calculated w i t h the semi-classical method of impact parameters [68]. Specifically, Xn = 1.578 x 10~* FII < y XijTe' where is the oscillator strength. exp ( ^ ) v 7 e (3.36) y < ^ • > is the four-parameter thermally averaged Gaunt factor taken from the work of Mewe [69]. Detailed balance is used to calculate the collisional de-excitation rate ^"=IHi)^ • B o t h X^ and X^ 1 -> 3(3 7 have units of [cm /s]. 3 Continuum Lowering In addition to the atomic processes listed above, the N L T E model also accounts for continuum lowering. T h i s effect reduces the ionization potential of each state xf by an Chapter 3. Numerical Simulations 47 (3.38) <z > <z> 2 Z' = K (3.39) z'(z' + i)jv; = 1.937 x 10~ T ± 3.3.3 -,1/2 (3.40) e Photoexcitation O p t i c a l p u m p i n g of ion populations i n the plasma is incorporated into the model through the concept of photon escape probability. T h e region under consideration is divided into n cells. Photons emitted in cell j have a finite probability Cji of being absorbed in cell i. If Nf is the total upper level population of cell i, and A \ is the spontaneous transition rate between the levels of interest, then u the upper level population of cell % is governed by a rate equation of the form (3.41) T h e first two terms on the right hand side describe mechanisms responsible for populating and depopulating the level N™, and hence include a l l of the atomic processes previously described. T h e t h i r d term accounts for the photoexcitation of the level by line radiation emitted from a l l cells j into cell i. T h e occurrence of a photon being emitted and reabsorbed in the same region is implicit in this final term. A s before, the steady state approximation dN^/dt = 0 is used to simplify calculations. A line photon emitted i n cell j has a probability ji P = e(Tji) P ~ P (Tji e + AT;) (3.42) Chapter 3. Numerical Simulations 48 of being absorbed in cell i. T h e angle-averaged probability P {r) that a photon traverses e an optical depth r without being absorbed or scattered is defined by [70] 1 r°° = 2J ^MrtJdy "( ) P T (3.43) 0 where ip and </>„ are the normalized emission and absorption line profiles and E is the u 2 second order exponential integral. T h e total line center optical depth A T , is related to the line center absorption coefficient KQ and the cell w i d t h Axi by A r , = KoAo-i . . ; (3.44) The following approximate expressions for P , adopted for the model, e Fe = (1 + 1.861607r + 0 . 8 1 7 3 9 3 r ) 2 0 -1 0 = 0.286ro (ln(1.95r ))/ V r < 3 ' 0 1 0 (3.45) r > 3 0 have been calculated assuming a Doppler line profile. T h e probability Cji is then calculated with where the factor of 1/2 accounts for the assumed isotropic forward and backward emission of radiation. 3.3.4 Line Transport Non-resonant absorption of line photons is accounted for i n the following manner. T h e optical depth of cell i at photon energy h u • A ^ W 0 is calculated w i t h = [«//("<>)+ «f>o.)] AXi (3.47) where the bremsstrahlung absorption coefficient /«// is given by 3 4 x 10 n 6 «// = ' ^1/2 , I ~ exp(-/w/ /T )] x Y.9f Z N £ 1 2 0 e f Z (3.48) Chapter 3. Numerical Simulations 49 and.K^(^o) is the total inner shell photoabsorption coefficient as described in [56]. Neglecting photoexcitation, the number of photons N? absorbed in cell j due to line emission from cell i is evaluated w i t h N[ =lNrA [E (T^ i UL + Ar^/2)-E (T^ 2 + ^/2 2 + Arr)}. (3.49) B y defining the continuum escape probability Pj% (r) as l P^(T) = E (T™ + ATZ/2 + T), 2 (3.50) the coupling coefficient Cji used to couple resonant absorption of line radiation to the ion population distribution can be approximated w i t h A ppe C» = NH— —— (3.51) which includes all of the absorption processes. Equation 3.51 employs the definitions AP** = P% (0) - Pg ( A r . ) (3.52) and AP = P j f ( 0 ) - P^(Ari) cn . (3.53) Consequently the number of photons absorbed due to all other processes is AP ^pcn cn ji N b - ji N P ^p pe + • (3-54) which represents an absorbed energy of E% = hv N$' 0 in the calculation of X . rad (3.55) Chapter 3. Numerical Simulations 3.3-5 50 Continuum Transport The method of multi-group diffusion [23] is used to transport bremsstrahlung and recombination radiation. T h i s allows the radiation w i t h i n each energy b i n to be treated in a manner similar to that used for line transport. Specifically, for bin k, extending from u k to v +u k E%{k) = \mAxi [E (T™ + A r f / 2 ) - E (r% + 2 + Ar™)] 2 (3.56) where Ii(k) is given by r \l {v) + Ti(k)= The ff + I {v))dv fb . (3.57) spectral intensity per unit frequency interval per unit volume for bremsstrahlung and recombination radiation is calculated w i t h IfM = J/2 "« E Z N g 2 z ff exp(-hu/T ) (3.58) e and 6 8 x 1 W = ' fl~ 47 "« E Z N 2 Z exp(-^/T ) £ e P v- 2 - 1 ^ r - 9 expfe*"/^) 1 f b (3.59) respectively. 3.4 Modifications for Present Work The original implementation of the N L T E calculation was specific to homogeneous alum i n u m targets. A number of changes were necessary to enable H Y R A D to handle the multi-material targets used for the present work. The C H layer on the laser deposition side is expected to dissociate completely into carbon and hydrogen ions. F o r this reason, the photoabsorption edge energies required Chapter 3. Numerical Simulations 51 Ionization Stage K - E d g e (eV) Reference C+° C+ 283.84 310.00 [20] [71] 320.13 354.31 390.67 489.99 [71] [72] [71] [71] 1 C+ • C+ 2 C +4 G +5 3 Table 3.1: Photoabsorption edge energies for the different ionization stages of carbon. are those of the various ionization stages of carbon atoms alone. Hydrogen ions do not contain any photoabsorption edges. Table 3.1 summarizes the K-edge energies used i n calculating the inner shell optical depth in the C H layer. Opacities needed for the radiation transport are calculated by fitting fourth order polynomials to existing opacity data [20] in energy regions between absorption edges, as suggested by Biggs and L i g h t h i l l [73]. In the case of C H (Z=3.5), the opacity is weighted by atomic mass, i.e. 86% for carbon and 14% for hydrogen. T h e N L T E calculation described previously made no reference to the specifics of the target material. M a t e r i a l dependence is introduced through the atomic input data. A t o m i c d a t a for only one species could be included i n the C R E calculation, and therefore the hydrogen ions in the C H plasma have been ignored. T h e carbon atomic data includes all ionization states plus four excited states i n each of the C + 4 and C + excited states. 5 Table 3.2 summarizes the electronic configurations considered by the present calculation. Electronic configurations i n the lower ionization stages of carbon have been ignored since the relatively small ionization potentials coupled w i t h moderate-to-high laser irradiances means that the majority of carbon atoms w i l l be fully ionized. Table 3.3 lists the chosen line transitions w i t h corresponding oscillator strengths. Oscillator strengths have been obtained from the work of Wiese [74]. Chapter 3. Numerical Simulations Ionization Stage +o c C+ 1 c+ 3 +2 C C +4 4 C + C+ 4 4 4 C+ 5 C + C + C +5 Ionization Potential (eV) 2p 2p 11.260 24.384 2s 2s Is 2p 2 47.888 64.494 2 392.09 84.128 37.502 2 3p 4p 5p ls 2p 5 +6 —' +5 C + C Outer Electron 3p 4p 5p C + C 5 . 52 21.088 13.433 489.99 122.36 54.255 30.425 19.468 0.00 Table 3.2: Ionization potentials for the ground and excited states included i n the N L T E model. Chapter 3. Numerical Simulations Transition ls ls ls ls 2 2 2 2 53 Oscillator Strength -ls5p -ls4p -ls3p -ls2p ls-5p ls-4p ls-3p ls-2p 0.024 0.051 0.141 0.647 0.014 0.029 0.079 0.415 Table 3.3: Transitions and oscillator strengths used in the N L T E calculation. 3.5 C o m p a r i s o n of N L T E Calculation w i t h R A T I O N Code To test the validity of the present N L T E calculation, average ionizations produced by our calculation for a carbon plasma were compared to those produced w i t h the widely used [75] atomic physics code R A T I O N [76]. A brief description of the physics included in R A T I O N is presented here, w i t h an emphasis on how R A T I O N differs from our N L T E calculation. A more detailed description of the level structure and coupling considered by the R A T I O N code can be found i n References [76] and [106]. 3.5.1 Energy Levels and Level Coupling In the calculation of ion level populations by R A T I O N , a distinct emphasis is placed on the Li-like through fully-stripped ionization stages. Ion species of neutral through Be-like are represented simply as a single ground state and its associated ionization potential. There are a total of 88 energy levels included in the R A T I O N calculation for carbon, and the ionization stages of Li-like through fully stripped comprise 84 of these levels. Li-like levels include detailed nl configurations for all / levels up to n = 4, w i t h level energies calculated from a method due to Edlen [77]. Energies for the n — 5 level are Chapter 3. Numerical Simulations 54 calculated w i t h the formula of Vainstein and Safronova [78]. Energies of levels with principal quantum number, n = 6-10 are calculated w i t h another method of E d l e n [79] He-like energy levels in states w i t h principal quantum number n = 2-7 are taken from the d a t a compilation of Scofield [80]. For levels between n — 7 and n = 10, the hydrogenic formula ' E(n,Z) = I (Z) p is used, where I (Z) P (3.60) n . is the ionization potential of the He-like species, also obtained from Scofield [80]. E q u a t i o n 3.60 is used for H-like energy levels w i t h n = 2-10, where I (Z) p is again taken from the work of Scofield [80]. T h e energy levels used i n the N L T E calculation of this work were taken from the d a t a tables of K e l l y [81]. Similar to the R A T I O N code, our calculation contains an expanded level structure in higher ionization stages, specifically 4 excited levels in each of the C and C + 5 + 4 ionization stages (refer to Table 3.2). However, the N L T E calculation contains a total of only 15 states, far fewer t h a n considered by R A T I O N . T h e 88 energy levels considered by R A T I O N represent specific atomic transitions. M a n y of these transitions occur between excited levels of adjacent and non-adjacent ionization states, in contrast to our N L T E calculation which couples excited states to ground states of the next ionization stage (refer to Figure 3.1). T h e coupling scheme in R A T I O N thus represents a more complete description of the atomic processes in the plasma. 3.5.2 Rate Coefficients A wide variety of rate coefficients are employed by R A T I O N , typically w i t h the choice of equation w i t h i n a particular process dependent upon the ion species. This represents a level of complexity not included in our calculation. Chapter 3. Numerical Simulations 55 T h e collisional ionization rate is calculated w i t h the prescription of L o t z [82], and detailed balance w i t h this rate is used to determine the three-body recombination rate. T h e spontaneous radiative recombination rate for all ion stages except H-like is calculated with the formula of Spitzer [83]. For H-like ions, this rate is determined by the method of Seaton [66]. T h e stimulated radiative recombination, photoionization, and photoexcitation rates, for all ion species, are calculated according to the methods outlined in the work of M i h a l a s [84]. T h e collisional excitation rate for H-like species employs the formula found in the work of V a n Regemorter [85], whereas the methods of Cochrane and M c W h i r t e r [86] are used for ionization stages between neutral and He-like. A u t o i o n ization rates for He-like ions are calculated from the work of Vainstein and Safronova [87], whereas this rate for Li-like ions is taken from the work of Gabriel [88]. Inverse autoionization rates for the two ion species are determined by the principal of detailed balance. T h e R A T I O N code contains all of the atomic processes considered by our N L T E calculations, plus many that are not, namely, autoionization, inverse autoionization and stimulated radiative recombination. T h e only rate coefficient formula common to R A T I O N and our N L T E calculation is the spontaneous radiative, recombination rate of Seaton [66]. 3.5.3 Results of Comparison Figures 3.2 to 3.6 illustrate the difference in average ionizations produced by the two calculations for a carbon plasma over the temperature range 50 to 550 e V and electron density range 1 0 19 to 1 0 23 cm - 3 . T h e greatest difference between the two calculations occurs for high density and low temperature plasma conditions. electron density of 1 0 23 cm - 3 For example, at an and plasma temperature of 50 eV, R A T I O N predicts an average ionization of 3.60, whereas the N L T E calculation predicts a value of 3.31. T h e Chapter 3. Numerical Simulations 56 6.00 5.95 h •I CO N 5.90 'c o CD O) CO S 5.85 5.80 - 5.75 100.0 200.0 300.0 400.0 500.0 600.0 Temperature (eV) Figure 3.2: Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma with an electron density of 1 0 c m . The dashed lines with circles correspond to R A T I O N calculations, whereas the solid lines w i t h squares correspond to N L T E calculations. 19 - 3 Chapter 3. Numerical Simulations 6.00 |—-n 1 • 57 r < 5.80 - 5 75 I 0.0 , I 100.0 lL I 200.0 i _l i I * 300.0 400.0 , L 500.0 i I . 600.0 Temperature (eV) Figure 3.3: Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma with an electron density of 10 c m . T h e dashed lines w i t h circles correspond to R A T I O N calculations, whereas the solid lines with squares correspond to N L T E calculations. 20 - 3 Chapter 3. Numerical Simulations 58 6.00 Temperature (eV) Figure 3.4: Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma w i t h an electron density of 10 c m . T h e dashed lines w i t h circles correspond to R A T I O N calculations, whereas the solid lines w i t h squares correspond to N L T E calculations. 21 - 3 Chapter 3. Numerical Simulations 59 6.5 -.---.-0-T---0 5.5 c o CO _N 'c o CD U) C i _O CD > < 4.5 3.5 1 0.0 1 • 100.0 1 1 200.0 • • 300.0 1 1 400.0 i • 500.0 I 600.0 Temperature (eV) Figure 3.5: Ionizations produced by the present NLTE calculation and RATION for a carbon plasma with an electron density of 10 cm" . The dashed lines with circles correspond to RATION calculations, whereas the solid lines with squares correspond to NLTE calculations. 22 3 Chapter 3. Numerical Simulations 2 0 I ' 0.0 i i 100.0 i J 200.0 i 60 i i i 300.0 400.0 J i 500.0 i 600.0 Temperature (eV) Figure 3.6: Ionizations produced by the present N L T E calculation and R A T I O N for a carbon plasma with an electron density of 1 0 cm"" . The dashed lines w i t h circles correspond to R A T I O N calculations, whereas the solid lines w i t h squares correspond to N L T E calculations. 23 3 Chapter 3. Numerical Simulations 61 difference between the two values is always less than 10%, and decreases as electron density and temperature increase. A t near-solid density and high temperature conditions the two calculations give very similar results (refer to Figures 3.5 and 3.6). Thus, even though the level coupling scheme used in our N L T E calculation is theoretically only valid for n < 1 0 e 22 cm - 3 [64], the fact that it predicts ionizations very close to those predicted by R A T I O N (which considers a much broader coupling scheme) allows the N L T E model to be used for the near soliddensity plasmas encountered in these simulations. 3.6 Survey of Simulation Parameters In all simulations, the numerical mesh is divided into five regions consisting of 150-300 cells (refer to Figure 3.7). T h e number of cells used varies w i t h target thickness, and is chosen to simultaneously maximize both numerical resolution and computational speed. T h e first three regions on the laser deposition side consist of the C H compound. T h e first of these three regions is uniformly zoned w i t h a thickness per cell of ~ . 0.05 A. Cell w i d t h in this region must be less than 10% of the laser wavelength to achieve sufficient resolution in the ablation process. T h e second region is geometrically zoned to provide a smooth transition between the finely zoned front side region and the more coarsely zoned t h i r d region. T h e bulk of the C H lies in the t h i r d region, constructed w i t h a w i d t h per cell of 0.2 u.m. Simulations w i t h a more finely zoned target did not produce quantitatively different results, but required more computing time. T h e fourth region contains the required 3.0 pm of a l u m i n u m . T h e code more easily handles constant volume per cell across dissimilar material interfaces (as opposed to constant mass per cell), and hence the aluminum layer spans 15 cells. Following the convention used at the Chapter 3. Numerical Simulations LASER CH CH CH Al Figure 3.7: Five-zone scheme used in target design. Si Chapter 3. Numerical Simulations 63 C H - A 1 interface, the cell w i d t h in the silicon layer is also set to 0.2 //m per cell. The thickness and hence number of cells in this fifth region depends upon the specific target, and therefore no single set of values was used for all of the simulations. T h e S E S A M E [89] equation of state table 7590 describes the thermodynamics of the C H layer. Electrical and thermal conductivities are calculated w i t h the Spitzer formulae as tabular d a t a is currently unavailable in our laboratory. T h e ionization state and opacity of each C H cell is calculated w i t h the N L T E model. T h e Q E O S [90] equation of state is used to generate E O S data along w i t h electrical and thermal conductivities for the aluminum. Ionization state was also obtained w i t h Q E O S , but opacity values were calculated using the N L T E model. E q u a t i o n of state data, electrical and thermal conductivities and ionization state for the silicon layer is interpolated from the Q E O S tables. Opacity calculations follow from the Biggs and L i g h t h i l l approach used i n the N L T E model. A l l simulations considered the following processes i n the C H layer : collisional ionization, radiative, three-body, and dielectronic recombination, collisional excitation , collisional de-excitation, spontaneous de-excitation and finally, continuum lowering. T h e full radiation transport package, consisting of continuum transport, line transport, and photoexcitation, was used for the C H , aluminum, and silicon layers. Results of the simulations are discussed in the following chapter. Chapter 4 Results of Radiation-Hydrodynamic Simulations T h i s chapter presents a survey of the results of the radiation-hydrodynamic computer code simulations. F i n a l parameters of target and laser pulse design are summarized, along w i t h density and temperature profiles corresponding to each of the shock compression processes. 4.1 4.1.1 General Results Hydrodynamics Figure 4.1 shows a snapshot of the hydrodynamic profiles w i t h i n the target during a simulation of the laser ablation process. Tables 4.1 and 4.2 present a summary of the pressures and shock speeds produced i n the C H and aluminum layers from single shock wave simulations. M a n y features of the hydrodynamic evolution of the shock wave agree w i t h previous work [38],[91]-[97]. 4.1.2 Radiation Transport T h e most obvious artifact of radiation transport is the radiatively heated zone ( R H Z ) . T h i s is the optically thick region of nearly constant density and temperature occurring immediately behind the shock compressed material. T h e explanation for the formation of this region involves analysis of the mean free path of x-rays emitted from the front side plasma. A s the shock propagates through the target, the w i d t h of the compressed 64 Chapter 4. Results of Radiation-Hydrodynamic Simulations 1000.00 71—i—l—i 1—l—l—r - 65 - | — l — l — l — l — I — l — i — i — i — i — i — i — r - Pressure (MBar) Density (x100kg/m3j| Temperature (eV) CORONA \ 100.00 SHOCK Al Si / RHZ 10.00 1.00 0.10 CH 0.01 0.0 10.0 20.0 30.0 40.0 Position (microns) Figure 4.1: Snapshot of hydrodynamic profiles w i t h i n the target at 200 picoseconds of simulation time. T h e shock wave is shown .propagating through the C H layer, imminent upon the C H - A 1 boundary. T h i s particular target consists of 27 C H - 3 pmA\ - 6 /j,m S i . T h e laser irradiance is 3.75 x 1 0 W / c m . 14 2 Chapter 4. Results of Radiation-Hydrodynamic Simulations Irradiance(10 13 W/cm ) 2 1.25 2.50 Pressure(MBar) 2.36 3.93 9.29 36.2 6.25 37.5 Speed(10 4 66 m/s) 1.75 2.43 3.02 6.75 Table 4.1: Pressures and speeds characteristic of single shock waves in the C H layer. Irradiance(10 13 1.25 2.50 6.25 37.5 W/cm ) 2 Pressure(MBar) 3.90 6.49 15.8 61.5 Speed(10 4 m/s) 1.50 1.88 2.72 6.02 Table 4.2: Pressures and speeds characteristic of single shock waves in the A l layer. region becomes larger than the mean free path for some subset of wavelengths emitted from the front side plasma. Hence, a larger portion of x-ray deposition occurs behind the shock front, rather than ahead of it. In uniform a l u m i n u m targets, this effect can actually eliminate a significant fraction of preheat [38]. This is clearly not the case for C H , arid so x-ray preheat is an important factor for a l l times during the hydrodynamic evolution. T h e w i d t h of the R H Z grows linearly in time, similar to the results obtained by Salzmann et al. [38] in homogeneous a l u m i n u m targets. T h e radiation front associated with the R H Z , propagating behind the shock front, is responsible for the second, delayed heating process observed in rear surface emission studies [96, 98, 99]. For moderate irradiances, the shock wave outruns the radiation front. However, for irradiances of ~ 1 0 16 W/cm 2 and higher, the radiation front w i l l run ahead of the shock front [97], greatly heating and decompressing the target material so that the shock process is relatively ineffective. T h e levels of radiation emission and transport are consistent w i t h previous studies Chapter 4. Results of Radiation-Hydrodynamic Simulations 67 using laser heated carbon foils [40]. For laser irradiances of 1.25, 2.50, and 6 . 2 5 x l 0 1 3 W / c m , the percentages emitted by the various processes breaks down approximately 2 as : 6% free-free, 87% bound-free, and 7% bound-bound. F o r the . 3 . 7 5 x l O 14 W/cm case, the percentages are : 7% free-free, 89% bound-free, and 4% bound-bound. 2 The fraction of bremsstrahlung and recombination (recombination is p r i m a r i l y into K shell levels [40]) increases w i t h increasing laser irradiance, and the level of bound-bound radiation decreases accordingly. T h i s can be attributed to the higher average ionization and hence fewer available electrons to make bound-bound transitions w i t h i n the ions. T h e bound-free intensity greatly exceeds the bound-bound intensity due to,electron collisional quenching of the line radiation [40]. T h e population density contributing to continuum radiation lies below the collisional quenching density, and hence is unaffected. T h e small difference i n the amount of line radiation between this work and that of Ref. [40] can be attributed to the number of bound transitions included in the model. However, since continuum radiation accounts for over 90% of the emitted radiation, a difference of a few percent in the level of line radiation cannot be expected to have any significant effect on target hydrodynamics. Finally, the ionization burn wave [96] which induces transparency to certain emission lines in a l u m i n u m targets is not expected to have a significant effect on radiation transport through the C H layer. Most of the emission lies above the 490 e V K-shell energy of C 6 + [40], and hence w i l l be attenuated by inner shell photoionization processes regardless of the ionization state. 4.1.3 Laser Absorption For certain sets of simulation parameters, the two methods of calculating the deposition of laser energy w i t h i n the target (refer to Section 3.2.1) predict different hydrodynamic conditions. T h e question then arises as to the validity of each of the absorption routines. Chapter 4. Results of Radiation-Hydrodynamic Simulations Routine % D u m p at n E M Wave Solver 68 C H . Pressure(MBar) A l Pressure(MBar) — 15.3 26.5 c r Inv. Brem. 100 26.7 45.5 Inv. Inv. Inv. Brem. Brem. Brem. 10 1 0.1 19.3 15.6 9.6 32.7 27.1 16.0 Table 4.3: Pressures in the C H and a l u m i n u m layers when different laser absorption routines are used. The inverse bremsstrahlung absorption routine contains a free parameter which dictates the amount of laser energy deposited at the critical density surface. T h i s energy dump is crucial i n the initial formation of the plasma, but can be expected to become less important for laser pulses w i t h nanosecond risetimes since a sizable fraction of laser energy w i l l reach the critical density layer only at very early times during the pulse. O n the other hand, for laser pulses w i t h sub-nanosecond risetimes, the absorption process will show a significant dependence on the value of this parameter. T h e calculation which solves the Helmholtz equations, on the other hand, contains no such free parameter. To study the effect of the energy dump parameter contained in the inverse bremsstrahlung absorption routine, a series of simulations were performed. of intensity 2 x l 0 Si target. 1 4 W/cm 2 A single trapezoidal pulse was assumed'incident on a 10.9/xm C H - 2/xm A l - 10pm T h e two absorption calculations were employed, and the switch to control the energy dump in the inverse bremsstrahlung routine was varied from 100% to 0.1%. N o r m a l incidence w i t h s-polarized light (i.e. E parallel to the plane of incidence) was used in the electromagnetic wave solver calculation. T h e flux limiter, another free parameter, was fixed at 0.1, a value suggested by detailed Fokker-Plank calculations [100]. Table 4.3 presents a summary of the shock wave pressures i n the C H and a l u m i n u m layers corresponding to the different simulation parameters. Chapter 4. Results of Racliation-Hydrodynamic Simulations 69 T h e inverse bremsstrahlung absorption routine with 1% dump at the critical density layer predicts roughly the same hydrodynamic conditions as the electromagnetic wave solver. T h e reason for this is evident in Figure 4.2, where the reflectivity of the absorbing plasma is plotted against time. D u r i n g the first ~ 1 0 0 picoseconds, the two calculations predict, very different absorptions. However, at ~200 picoseconds, the absorption levels are roughly the same, and by 300 picoseconds, they are nearly identical. T h e spikes occuring in the reflectivity level of the inverse bremmstrahlung absorption routine are a result of the finite number of cells used to model the plasma density gradient in the absorbing region. Hence, for few-hundred picosecond time-scale and longer simulations, the two absorption routines will lead to the same hydrodynamic conditions. However, for simulations on time-scales shorter than ~300 picoseconds, care must be taken to choose the correct absorption mechanism. 4.1.4 Flux Limiter One of the largest uncertainties in computer simulations of this type concerns value of the flux limiter, / , to be used. A value of / = 0.6, corresponding to a free-streaming flux for a Maxwellian electron velocity distribution, is almost certainly too large for the present laser irradiances. Lower values of / limit the amount of electronic heat conduction from the critical density surface to the cold solid, thus producing a much hotter coronal plasma. A s less heat is transported to the ablation front, the the shock pressure w i t h i n the target is reduced and the shock arrival at the a l u m i n u m layer is delayed. To obtain the general trend of shock characteristics on flux limiter value, three simulations have been performed w i t h the following parameters : 2.0 x 1 0 14 W/cm 2 on 10.9/xm C H - 2yum A l - 10/xm S i . Absorption of laser energy is calculated w i t h the electromagnetic wave solver routine. T h e value for / is varied, w i t h a l l other simulation parameters oo 1 — — —•— — —•— — — — — —•—>— 1 0.0 1 1 1 1 100.0 1 J 1 1 1 1 200.0 300.0 Time (ps) • • • 1 • • •—• I 400.0 500.0 Figure 4.2: P l o t of reflectivity as a function of time for both the electromagnetic wave solver and inverse bremsstrahlung absorption routines. Chapter 4'. Results of Radiation-Hydrodynamic Simulations Flux Limiter Peak Pressure ( M B a r ) Shock A r r i v a l T i m e (ps) 0.60 0.10 0.03 36.9 26.5 18.5 214 235 247 Table 4.4: Dependence of peak pressure w i t h i n and shock arrival time at the aluminum layer on the value of the flux limiter used. Simulation parameters are described in the text. held constant. T h e results for conditions w i t h i n the aluminum layer are summarized in Table 4.4. T h e different values of the flux limiter have significant effects on target hydrodynamics, w i t h shock speed and pressure showing a larger variation than compression or temperature. T h e variation in shock speed w i l l shift the temporal position of the steady state conditions and so this factor must be carefully considered during experimental design. Furthermore, the flux limiter will also play a critical role i n schemes where the t i m i n g between multiple shock waves is important. B y experimentally characterizing shock transit through aluminum and C H targets, and comparing w i t h numerical simulations, the correct value for the flux limiter can be deduced. 4.1.5 F i n a l Choice of A b s o r p t i o n R o u t i n e and F l u x L i m i t e r T h e absorption routine which solves the Helmholtz equations, assuming plane-polarized light, was chosen for the final calculations. T h i s routine, coupled w i t h a flux limiter of 0.1, predicts hydrodynamic conditions similar to those of more well established hydrocodes [101]. 71 Chapter 4. Results of Radiation-Hydrodynamic Simulations 4.2 72 O p t i m i z e d Shock Compression Schemes 4.2.1 Laser Pulse and Target Design T h e laser pulse consists of one or two intensity plateaus, w i t h linear rise and fall times between each level. T h e rise and fall times are fixed at 100 picoseconds, in accordance w i t h current shaped laser pulses. Moreover, this short, linear rise time is capable of generating shock waves at intensities below which shock wave formation is not possible w i t h Gaussian laser pulses. T h e individual intensities 7; and plateau lengths Tj of each laser pulse are varied to achieve the appropriate hydrodynamic conditions. Figure 4.3 presents a visual aid to the relevant laser parameters. A l t h o u g h the essential target characteristics were described in the introduction, a more detailed description would prove useful. T h e C H layer on the laser deposition side serves three purposes. Firstly, it provides an absorption medium to minimize the effects of radiation preheat on the a l u m i n u m layer. Secondly, by acting as a material of low shock impedance, the impedance mismatch between this layer and the a l u m i n u m layer allows enhancement of the shock pressure as the shock wave passes through the C H - a l u m i n u m boundary. Finally, the C H provides support to the shock compressed a l u m i n u m layer, preventing release out the front side before steady state is obtained for the necessary time. B y varying the C H thickness, the proper delay between multiple shocks, rarefactions and release waves can be achieved. W h i l e increasing the thickness of this layer would certainly reduce the level of radiative preheat, excessive thickness would degrade the absorption measurement. T h e C H thicknesses presented as final target designs represent a trade off between preheat and absorption measurement factors. T h e sole restriction in target design concerns the thickness of the a l u m i n u m layer, which is preset to 3.0 yum. T h i s value represents a compromise between transmission level of the absorption signal (from the backlighter) and observable effect. If the layer is too Chapter 4. Results of Radiation-Hydrodynamic Simulations 100 X1 ioo x 1 0 0 2 Time (ps) F i gure 4.3: Visual aid to the important parameters in laser pulse design. The rise and fall times between levels is fixed at 100 ps. 73 Chapter 4. Results of Radiation-Hydrodynamic Simulations 74 thick, the intensity of the transmitted x-rays may be too faint to detect w i t h the desired accuracy. Conversely, if the a l u m i n u m is too thin, the amount of material displaying the photoabsorption edge shift may be too little for a definitive measurement. T h e silicon layer serves only to prevent premature shock breakout and subsequent release of the aluminum layer. T h i s layer must be made as t h i n as possible to minimize the effects of absorption and hence signal loss, and yet must maintain steady state in the aluminum layer for the required time. T h e similar atomic properties between this layer and the aluminum make silicon the natural choice for target design. Moreover, the use of silicon could aid in target fabrication ; bulk silicon could easily be etched to the desired thickness, and the a l u m i n u m and polystyrene layers deposited on top. 4.2.2 Density and Temperature Conditions Table 4.5 and Figure 4.4 summarize the density and temperature conditions achieved in the aluminum layer w i t h the different shock compression processes. T h e quoted values represent the average value of five of the Lagrangean mesh points (regularly spaced across the a l u m i n u m layer) at a time half way through the steady state plateau. Steady state is considered to have been reached after the rear aluminum mesh has been shocked up and the rarefaction generated at the aluminum-silicon interface (described later) has traversed the entire a l u m i n u m layer. T h e uncertainties are the difference between the highest and lowest values achieved by any one of the Lagrangean mesh points in the layer (during the steady state time) and the average value described above. 4.2.3 Discussion of Individual Compression Schemes To aid in future target design, a brief discussion of the important processes i n each of the shock compressions of Figure 4.4 is given. Emphasis is placed on the description of the shock and rarefaction waves propagating inside the target. Chapter 4. Results of Radiation-Hydrodynamic Simulations 100 ~1 1 1 r- ~| 1 T 1 75 f~ H 1 I l - High Z* 5 > CD zs "CcO o •I—• I— Reference Isotherm 10 CD CL E CD dx^ Isochore t x ^ J 1.5 2.0 2.5 1 i i i I 3.0 i i i i I 3.5 _l I L_ 4.0 Compression Figure 4.4: Illustration of the densities and temperatures achieved in the shock compression schemes. Chapter 4. Results of Radiation-Hydrodynamic Simulations 76 State P/PAI T(eV) Reference. 2.4±0.2 8.5±0.5 18-3-4 700 , — 6.25 , — Isochore 2.4±0.3 3.3±0.3 2.5±0.1 900 , 600 900 , 600 1.25 , 2.25 2.0±0.1 22-3-5 22-3-5 Isotherm 3.7±0.4 Pdx > Tfa 2.2±0.1 8.5±0.8 4.9±0.3 27-3-7 27-3-7 900 , 600 900 , 600 2.5 , 10.0 2.5 , 10.0 H i g h Z* 3.1±0.3 29.0±2.0 27-3-6 700 , — 37.5 , — Ptx, Ttx 1 C H - A l - S i nm r l , 2 (P'S) Ii,2 1 0 13 W/cm 2 1.25 , 2.25 Table 4.5: Compression and temperature conditions obtained in the shock compression schemes, along w i t h final target and laser pulse parameters! Chapter 4. Results of Radiation-Hydrodynamic Simulations 77 T h e most important factor to consider when choosing the compression and temperature conditions of the reference state, and hence the isothermal and isochoric paths, is radiative preheat. T h e relatively low compressions and temperatures chosen to map out the compression-temperature space can be obtained w i t h sufficiently low laser irradiances that radiative preheat does not seriously degrade steady-state conditions. O n the other hand, the plasma conditions generated should still produce a measurable K-edge shift [16, 17, 56]. T h e graphs in Figures 4.5 to 4.17 represent the time evolution of 5 Lagrangean mesh cells w i t h i n the aluminum layer. T h e first and last cells of the a l u m i n u m layer are included, as well as three others which are equally spaced among the remaining cells. T i m e zero on the graph corresponds to the start of the laser pulse. T h e time interval over which quasi-steady state conditions are assumed to he present is indicated by a horizontal line above the plot of the Lagrangean mesh. In a l l of the figures, the cells are labeled as follows : front side (first cell) — , fourth cell • - -, eigth cell , twelfth cell - — , and rear side (fifteenth cell) - • - . Reference State Plots of compression, temperature and pressure w i t h i n the a l u m i n u m layer during the shock compression process associated w i t h the reference state appear in Figures 4.5 to 4.7. A s stated previously, the term 'reference state' is merely a label for the state w i t h which conditions on the three branches are compared in order to deduce the effects of varying hydrodynamic conditions on the K-edge energy. M u c h of the analysis' of this state is common to the other shock produced states, and is only described once. R a d i a t i o n transport is a noticeable factor in this compression scheme, manifesting itself in a variety of ways. F i r s t , there are the 0.8 e V and ~ 1 0 % heated and decompressed conditions i n the a l u m i n u m prior to shock arrival. However, these conditions are not Figure 4.5: Compression (i.e. density) as a function of time w i t h i n the aluminum layer for the reference state. The decompression prior to shock arrival is due to radiative preheat. T h e time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. Figure 4.6: Temperature as a function of time w i t h i n the aluminum layer for the reference state. Note the temperature increase prior to shock arrival due to radiation preheat. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. Chapter 4. Results of Radiation-Hydrodynamic Simulations Figure 4.7: Pressure as a function of time within the aluminum layer for the reference state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. 80 Chapter 4. Results of Radiation-Hydrodynamic Simulations uniform across the layer, and as expected, the higher temperatures and decompressions experienced by the cells closer to the laser deposition side of the target result i n lower peak shock compression. Connected to the non-uniformity in radiation preheat is the spread of densities and temperatures observed across the a l u m i n u m layer at any time during the steady state plateau. T h i s range arises as a consequence of the cells starting at different i n i t i a l conditions in the shock compression process. F i n a l l y , radiation emitted from the front side plasma causes continual decompression of the a l u m i n u m layer after the shock compression process is complete. T h e heating due to radiation absorption is offset by the cooling associated w i t h the expansion process, and so even though the density drops, the temperature remains nearly constant. Single shock processes are the easiest to describe simply because there are fewer counter- and co-propagating waves w i t h i n the target. T h e decompression (from rear to front) visible immediately after the rear-most a l u m i n u m cell is shock compressed can be attributed to a rarefaction wave generated at the aluminum-silicon interface. T h i s effect is apparent in all of the simulations, resulting from an impedance mismatch between layers common to a l l the target designs considered in this work. After this rarefaction has traversed the entire a l u m i n u m layer, the reference state (po,T ) forms, existing for ~150 0 picoseconds. T h e lifetime of this state is limited by shock breakout. T h e rarefaction wave associated w i t h release from the target rear surface propagates back into the aluminum, illustrated as the rear mesh being the first to decompress. Theoretically, a thinner silicon layer could be used to force the steady state to a 100 picosecond lifetime, and improve the absorption measurement by l i m i t i n g signal loss through the silicon. However, simulations with 3//m of silicon show that not enough support is provided for the shock compressed aluminum, and it begins to decompress almost immediately. Four microns of silicon appears to be the m i n i m u m useful thickness. 81 Chapter 4. Results of Radiation-Hydrodynamic Simulations 82 High Average Ionization Density, temperature and pressure conditions necessary to achieve a high average ionization (relative to those achieved in the other shock compression schemes) are presented in Figures 4.8 to 4.10. Radiation transport is clearly an important factor for this irradiance, as the aluminum layer is heated to nearly 4 eV and decompressed by ~20% prior to shock arrival. The purpose of producing such a high shock pressure is to raise the average ionization of the aluminum layer to ~3.5. Although Figure 4.11 suggests that the average ionization should be greater than 4 (given the shock pressure of this process), QEOS is known to overestimate the average ionization, and a value of ~3.5 is more likely [41]. Figure 4.11 serves to demonstrate the large shock pressure increases necessary to generate small increases in the average ionization. This shock compression process allows a study of ionization effects on the K-edge position. Since single shocks produce higher temperatures for a given compression than multiple shocks, the single shock used to form this state will produce a density not too different from that of the Reference state, aiding in isolating the ionization effects. Once the temperature effects have been studied (from the isochoric branch), these can be unfolded from the observed K-edge dependence in this scheme to arrive at the true ionization effect. The irradiation conditions used to generate this state represent the limit of practicality for obtaining any form of steady state conditions in experiments of this type. The 10% variation between highest and lowest densities achieved over the 100 picosecond plateau pushes the limits of acceptability. Originally, a state containing mostly A l be produced, but the irradiance needed was almost 6 x l 0 1 4 + 4 ions was to W / c m , and the final shock 2 compressed state displayed 30% variation over a 100 picosecond plateau. Chapter 4. Results of Radiation-Hydrodynamic Simulations 100.0 200.0 300.0 400.0 500.0 600.0 700.0 83 800.0 Time (ps) Figure 4.8: Compression (i.e. density) as a function of time within the a l u m i n u m layer for the high average ionization state. Note the much higher level of decompression caused by radiative preheat in this case as opposed to other simulation results. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. 0.0 100.0 I I I I I 200.0 I II I II I 300.0 400.0 I I I I I I I I I I II I 500.0 600.0 700.0 I I I I I l l 800.0 Time (ps) Figure 4.9: Temperature as a function of time within the aluminum layer for the high average ionization state. Radiation preheat is clearly an important effect for this laser irradiance. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. " Chapter 4. Results of Radiation-Hydrodynamic Simulations Q l~i i i "100.0 : I I 200.0 i i i I i i 300.0 | ' i i i i i 400.0 ii i i i i i 500.0 i i 85 ,,i 600.0 ; I , , , 700.0 800.0 Time (ps) Figure 4.10: Pressure as a function of time within the aluminum layer for the high average ionization state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. Chapter 4. Results of Radiation-Hydrodynamic Simulations 14.0 Trr] 1 1—i—r~i-ri-T | 1 86 1—i i i i i ij 12.0 10.0 c o "CO N 8.0 'c o CD O) CO \ CD 6.0 > < 4.0 2.0 0.0 _i 0.1 i i _ J 1.0 1 LII I I I I I 10.0 - 1 u_d_ 100.0 _l 1000.0 I I I I I I I 10000.0 Shock Pressure (MBar) Figure 4.11: Average ionization as a function of shock pressure for aluminum along the single shock Hugoniot. These values are obtained from the QEOS model. Chapter 4. Results of Radiation-Hydrodynamic Simulations 87 Isochoric State Figures 4.12 to 4.14 present hydrodynamic conditions w i t h i n the a l u m i n u m layer during the double shock process used to obtain the point on the isochore of Figure 4.4. T h e compression associated w i t h each plateau of the laser pulse is easily recognizable, and remains uncomplicated by the presence of radiation effects. One obvious benefit of using a low laser irradiance is the lack of significant radiation preheat. T h i s is clearly illustrated in Figure 4.13. Neither of the two shock compressions are complicated by the decompressed and heated initial conditions, and hence these results are closely approximated by simulations which do not include any form of radiation transport. T h e relative t i m i n g of the two shocks used to create the hydrodynamic conditions (p , T ) leads to the formation of an 80 picosecond steady state (pt ,T ) t t x tx prior to the arrival of the second shock. T h i s state is formed just after the rarefaction wave generated by the first shock at the aluminum-silicon interface has propagated back through the aluminum, decompressing it. A l t h o u g h the density and temperature conditions created at this first plateau do not lie on either the isotherm or isochore of the reference state, they do provide an opportunity to obtain K-edge data without using a separate pulse-target scheme. T h e 200 picosecond wide plateau at (pt,T ) provided by the second shock is ultimately t limited by the thickness of the silicon layer. T h e rarefaction wave associated w i t h shock breakout decompresses the a l u m i n u m , demonstrated in Figure 4.12 as the density drops in order of rear to front. Isothermal State Compression, temperature and pressure conditions associated w i t h the state isothermal to the reference state appear i n Figures 4.15 to 4.17. T h e analysis of the shock and Chapter 4. Results of Radiation-Hydrodynamic Simulations 3.0 88 I | 1 I T'l" I II I I j I I I I I I I I [ [• I I 1 I | | | | | | | | | I I I l 1 2.0 c o CO CO CD Q. E o O 1.0 0.0 1000.0 1 I I 1 1 1 1200.0 1400.0 I 1 i l l''I II I I 1600.0 1800.0 Time (ps) 1 I 2000.0 2200.0 Figure 4.12: Compression (i.e. density) as a function of time within the aluminum layer for the isochoric state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. Chapter 4. Results of Radiation-Hydrodynamic 0.0 1000.0 1 1 1 1 1 1 1 1 1 1 1200.0 1 , 1 1 1400.0 Simulations 1 1 1 1 1600.0 89 1 1800.0 2000.0 ' 2200.0 Time (ps) Figure 4.13: Temperature as a function of time w i t h i n the a l u m i n u m layer for the isochoric state. T h e low laser irradiance used for the first pulse leads to a situation where radiation preheat is not a significant factor. T h e time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. Chapter 4. Results of Radiation-Hydrodynamic Simulations 90 10.0 8.0 - CO CD ZJ CO CO CD Ql 6.0 - 4.0 h 2.0 \- 0.0 1000.0 1200.0 1400.0 1600.0 1800.0 2000.0 2200.0 Time (ps) Figure 4.14: Pressure as a function of time within the aluminum layer for the isochoric state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. Chapter 4. Results of Radiation-Hydrodynamic Simulations 91 rarefaction waves generated i n this shock process is identical to that of the isochoric state. A g a i n , radiation effects do not play a significant role. T w o quasi-steady state plateaus are produced which can generate K-edge data. T h e first plateau, occurring between 1300 and 1400 picoseconds, contains hydrodynamic conditions which lie on the isochoric branch, and provide a convenient median point to the reference state and the previous shock compression scheme. T h e second plateau, occurring between 1490 and 1600 picoseconds, is limited by the thickness of the aluminum layer. Final Comment T h e six distinct hydrodynamic conditions generated by the four shock compression schemes described here should provide enough data to unravel the density, tempera- ture, and ionization state dependencies of the K-edge position. Three of the main points, namely the reference, isochoric, and isothermal, lie along nearly exact isochoric and isothermal branches from the reference point. T h e point labeled (pdxiTdx) in Figure.4.4 is a rather fortuitous by-product of the isothermal state, and serves to fill out the isochoric branch. A n o t h e r point along the isothermal branch would be ideal, but one could not be found which did not suffer greatly from radiation preheat. Chapter 4. Results of Radiation-Hydrodynamic 900.0 1100.0 1300.0 92 Simulations 1500.0 1700.0 1900.0 Time (ps) Figure 4.15: Compression (i.e. density) as a function of time w i t h i n the a l u m i n u m layer for the isothermal state. The time interval over which quasi-steady state conditions are achieved is.indicated by the horizontal line. Chapter 4. Q Results of Radiation-Hydrodynamic I ii 900.0 i i i i i i i I i 'i i 1100.0 iii 93 Simulations i 1300.0 1500.0 Time (ps) iiiii 1700.0 , 1 1900.0 Figure 4.16: Temperature as a function of time within the a l u m i n u m layer for the isotherm a l state. T h e low laser irradiance used for the first pulse leads to a situation where radiation preheat is not a significant factor. The time interval over which quasi-steady " state conditions are achieved is indicated by the horizontal line. Chapter 4. Results of Radiation-Hydrodynamic Simulations 900.0 1100.0 1300.0 1500.0 1700.0 94 1900.0 Time (ps) Figure 4.17: Pressure as a function of time w i t h i n the aluminum layer for the isothermal state. The time interval over which quasi-steady state conditions are achieved is indicated by the horizontal line. Chapter 5 Experimental Facility T h i s chapter presents a description of the laser facility and analysis techniques used in obtaining x-ray spectra from various elements. 5.1 Laser Facility A schematic of the laser facility appears in Figure 5.1. T h e laser chain consists of a Quantel N d : Y A G oscillator and pre-amplifier followed by four Quantel Nd:Glass amplifiers. T h e oscillator is passively Q-switched w i t h a saturable absorber to produce a single 2.3 ns (full w i d t h at half maximum) temporally and spatially Gaussian pulse at 1.064 /xm. A second saturable absorber is placed between the oscillator and the rest of the chain to prevent premature seeding of the oscillator by spontaneous emission from the amplifiers. T o reduce beam nonuniformities incurred i n the amplification process, a vacuum spatial filter follows each amplifier. T h e horizontal polarization of the beam is preserved by pairs of thin film polarizers situated between each amplifier and the following spatial filter. T h e amplified 1.064 p:m pulse then passes through a deuterated potassium dihydrogen phosphate ( K D P type II) crystal where approximately 50% of the beam energy is converted to 0.532 / i m light. B e a m energies of up to 15 J at 0.532 u,m can be obtained from this system. T h e remaining infrared beam is dumped, and the 0.532 //m light is steered to the target chamber and beam diagnostic equipment by a series of dichroic mirrors. T h e pulse shape and beam energy of every laser shot is monitored by measuring the 95 5. Experimental Chapter BSR/ Facility LENS h ' G E N T E K _ LENS BSR 96 ND I -| | PHOTODIODE BSR BSR BSR\ T A4 SF 5 IRM P7, P8 SHG SF 6 F3SR SR / I R DUMP A3 SF 4 IRM P4'P3 A2 SF 3 P2 P1 A 1 SF 2 7 ^ X- PA H = h \ v SF 1 AP I R M \ - 0 — I I OSC SA Figure 5.1: Schematic of the laser facility. O S C is the oscillator, P A is the preamplifier and A l - 4 are the amplifiers. I R M is an infrared mirror and SF1-6 are the vacuum spatial filters. A P is the hard aperture and S A is the saturable absorber. P l - 8 are polarizers, B S R is a beam splitter-reflector, N D is a neutral density optical absorber, I F is an interference filter, and G E N T E K is the piezoelectric beam energy monitor. S H G is the second harmonic generation crystal. Chapter 5. Experimental Facility 97 leakage through a beam splitter-reflector. A Gentek [102] piezoelectric detector is used to obtain the energy measurement. T h e absolute energy calibration of the Gentek is accomplished by firing the full beam into a Scientech [103] calorimeter and comparing the calorimeter reading w i t h the piezoelectric reading. T h e pulse shape is obtained w i t h a H a m m a m a t s u [104] photodiode coupled to a Tektronix 7104 oscilloscope. A sample laser pulse appears in Figure 5.2. Shot to shot variations in laser power are less than ten percent. 5.2 Focal Conditions In order to achieve known irradiation conditions, it is first necessary to map out the focal distribution of the lens. T h e following procedure is used. A surface deposition graticule is placed in the target plane, and using an incandescent lamp in conjunction w i t h an f/1.4 camera lens, an image of the graticule is obtained on P o l a r o i d film. A 5300A interference filter (100A bandwidth) is placed in front of the Polaroid camera to help avoid chromatic aberrations (see Figure 5.3). T h e image of the graticule is used to determine the magnification of the image relay system as well as focus the camera lens on the target plane. T h e graticule is removed and a high reflectance mirror is placed at the entrance of the target chamber to divert nearly all of the laser light into the Scientech calorimeter. T h i s is to avoid damaging the image relay system, as the full laser chain must be fired in order to obtain the most accurate focal distribution. Effects such as thermal lensing, which can alter the focal position by several millimeters, appear only when the full chain is fired. A n image of the focussed beam at the target plane is then obtained as the^ focussing lens position is varied. Different neutral density filters placed in front of the Polaroid camera allow images of the focal spot corresponding to different intensity levels on film, thus expanding the dynamic range of the measurement. Chapter 5. Experimental Facility Figure 5.2: Sample laser pulse. The horizontal scale is 1 nanosecond per division. Chapter 5. Experimental Facility 99 Figure 5.4 shows the smallest focal spot obtained (best focus) at different intensity levels. T h e calculation of the irradiance on target is somewhat arbitrary, as there is considerable freedom in choosing the size of the spot radius to use. T h e present work uses the spot radius corresponding to the 60% level, as has been used in previous work [56, 55]. W i t h this definition the irradiance of the beam on target is calculated w i t h 0.6 1 = xE — , (5.1) where E is the beam energy, r is the full w i d t h half m a x i m u m duration of the laser pulse, and r 6 0 is the radius of the spot which contains 60% of the beam energy, in this case 30/im. T h i s yields an irradiance on target of 1.5 x 1 0 14 W/cm . 2 Originally, three different irradiances were to be used to investigate how plasma conditions affected spectral emission. However, the lower irradiances achieved by attenuating the beam w i t h neutral density filters, and enlarging the spot size on target did not produce sufficient x-ray emission levels for analysis. T h e only irradiation condition which could produce sufficient emission was that corresponding to best focus. 5.3 X-Ray Spectrometers Figure 5.5 shows the layout of the target chamber. A single flat pentaerythritol ( P E T ) B r a g g diffraction crystal is used to record emission between 7.3A and 8.3k . B-10 foil [105].. is used to protect the fragile crystal from plasma debris, and provides an absorption edge on the exposure which is used to verify the wavelength calibration. Spurious film exposure due to secondary fluorescence is prevented by i) inserting a plastic shield between the plasma and the spectrometer crystal, ii) fabricating the equipment from a low atomic number material such as plastic, and iii) covering aluminum parts w i t h t h i n layers of plastic. T h e emission spectra are recorded on K o d a k P F 4 6 9 film, which consists of K o d a k emulsion number 2530 on an acetate base. T h e film is cut into strips approximately 3 Chapter 5. Experimental Facility 100 Figure 5.3: Experimental arrangement for scanning the focal distribution. N D is a neutral density optical absorber, and I F is the interference filter described in the text. Chapter 5. Experimental Facility 101 Figure 5.4: Images of the laser spot at best focus. Attenuation of the top picture is 100 times greater than the attenuation of the bottom picture, and 10 times greater than the attenuation of the center picture. T h e magnification of these images is 165 times. Chapter 5. Experimental Facility 102 cm long by 1 cm wide and housed inside a portable, light tight film cassette. T h e film cassette is machined from a solid block of aluminum and the front surface (excluding the entrance window) covered w i t h plexiglass (refer to Figure 5.6). T h e entrance window of the cassette is 0.6 cm tall by 2.0 cm wide and is covered w i t h a strip of beryllium, allowing the sealed unit to be light tight for normal room light operation, while still allowing x-ray emission to reach the emulsion. A removable mask allows the a l u m i n u m reference spectrum to be recorded on one half of the emulsion, and then by flipping the mask, the sample spectrum can be recorded on the other half of the emulsion. The film cassette is not moved between reference shot and sample shot, and so w i t h b o t h exposures on the same piece of film, an accurate wavelength calibration can be obtained. E a c h sample element spectrum is obtained w i t h two laser shots i n an effort to improve spectral quality. 5.4 Targets Target thickness is compensated for by translating the target mount axially along the beam line. T h i s ensures that the different spectra can be overlaid on a single exposure while still obtaining an accurate wavelength calibration. Table 5.1 presents a summary of the different targets used i n this work. 5.5 Processing of Data E a c h piece of exposed film is carefully developed and dried, and the resulting negative is mounted between glass slides to avoid damage. T h e spectra are then measured using a microdensitometer. T h e microdensitometer uses a photomultiplier tube to measure the film exposure of a 1 m m t a l l by 8 pm wide region, and the output, i n the form of a current, is recorded Figure 5.5: Target chamber and P E T crystal spectrometer. T h e diameter of the target chamber is approximately 90 cm. Chapter 5. Experimental Facility CASSETTE BODY 104 ACCESS CAP PLASTIC COVERING BERYLLIUM WINDOW /////. / / / / / FILM PLANE BERYLLIUM WINDOW Figure 5.6: F i l m cassette. T h e film plane is indicated by the dotted line. Chapter 5. Experimental Facility 105 Element A t o m i c Number Sample T y p e Bismuth Lead Gold Tungsten Tantalum Samarium Tin Silver Niobium Yttrium Germanium Aluminum Magnesium 83 82 50/im on 125/xm plastic 79 74 73 62 50 47 41 39 32 13 12 100/xm foil 25u,m foil 25/xm foil 25/xm foil 30/im foil 50/xm foil 50//m foil 100/xm foil 100/xm foil 500//m wafer 25/jm foil 25/xm foil Table 5.1: Summary of sample elements used in this work in order of decreasing atomic number. Chapter 5. Experimental Facility 106 for every 5 pm step of the scanning stage: T h e dynamic range of the microdensitometer has been determined to be approximately 100 [106] whereas the dynamic range of the film is only 30, and hence there is no loss of information w i t h this technique. T h e data from the microdensitometer scan represents the film exposure as a function of distance along the scan. T h e first step i n the calibration procedure consists of converting this distance measurement into a wavelength scale. T h i s spectral calibration is derived using three lines whose wavelengths and spatial relationship on the exposure are accurately known. A list of the calibration lines used i n this work appears in Table 5.2. T h e calibration lines should span as large a portion of the recorded spectrum as possible. Figure 5.7 shows the spectrometer arrangement used to derive the spectral dispersion at the film plane. For photons of wavelength A, the B r a g g diffraction angle 0 is given by A = 2dsin0 (5.2) where d is the lattice spacing of the crystal (4.37.1 A for P E T ) . U s i n g this equation, and three calibration lines, one can derive a and R and the following relation to convert the spatial coordinate on film to a wavelength : (5.3) where X is the spatial film coordinate. T h e angle QQ'IS the B r a g g diffraction angle of the longest wavelength reference line. T h e variable R is the path length from the plasma to the film plane for this reference line, and a is the angle of intersection this line makes w i t h the film plane. T h e choice of calibration lines is somewhat arbitrary, although more consistent results are obtained when satellite lines are used instead of resonance lines.. Resonance lines have high opacity and tend to overlap lines from other, more highly excited atomic states. Thus, while the resonance lines are dominant, it is more difficult to accurately locate Chapter 5. Experimental Facility 107 PLASMA o CRYSTAL Figure 5.7: Arrangement of the x-ray spectrometer. longest wavelength reference line is shown. The path length R = a + b of the Chapter 5. Experimental Facility 108 Transition Wavelength (A) Reference #12 7.2710 7.3635 7.7571 7.8574 7.8683 7.8717 [107] [107] [108] [107], [107] [107] 8.0713 [107] #11 He a abed k lj op Table 5.2: A l u m i n u m wavelengths used i n deriving the spectral dispersion relation. their center position. T h i s inaccuracy can mean the difference between values for R , a and # which are consistent w i t h the experimental set up, a n d those which are not. 0 T h i s sensitivity is due almost entirely to the small overall size of the spectrometer. The longest path length is only a few times the entire exposure w i d t h on film. A s well, the sphericity of the emitting plasma introduces curvature to the spectral lines on film. T h i s curvature translates into extra line w i d t h when the exposure is scanned, hence leading to uncertainty i n determining the true line position. T h e second step i n the analysis is to convert the film exposure obtained from the microdensitometer scan into relative intensities. T h i s is done using a best-fit p o l y n o m i a l obtained from previous work [106], and is expected to be accurate to w i t h i n 20% i n the relative intensity range of interest. Specifically, the relative intensity is obtained from the film exposure using 1 = 0.02819/ - 0 . 2 5 6 2 / + 0 . 8 3 6 9 / - 0.2298/ - 3.455 4 3 2 (5.4) where X is the relative intensity and / is the film exposure. T h e final step is to correct the relative intensity for the transmission through the b e r y l l i u m window on the film cassette and through the B-1'0 foil used to protect the crystal. T h e beryllium window correction follows exponential attenuation, where the opacity Chapter 5. Experimental Facility 109 is calculated from a best fit polynomial i n the energy region of interest. Unfortunately, the exact composition of the B-10 was unavailable, and so the attenuation properties had to be determined experimentally. Relative intensities of lead spectra taken w i t h and without the B-10 foil are compared to derive an effective attenuation as a function of wavelength, which is applied to the data. T h e following chapter presents the x-ray spectra obtained by the methods outline above. Chapter 6 Results and Discussion of X - R a y Spectra T h i s chapter presents the results of the x-ray spectra obtained from various elements by the methods outlined in the previous chapter. Comparison w i t h previously published spectra .is used to determine select transitions, but many other transitions remain unclassified. Suitability of elements as backlighters is also discussed. 6.1 X - R a y Spectra A print of a sample a l u m i n u m and tantalum exposure appears in Figure 6.1. Some of the known lines for both spectra are labeled (Refer to Figures 6.2, 6.3, and 6.8). The most important point to note about the following spectra is that they are neither time nor space resolved. Rather, they represent the average spectral emission of the plasma over the entire emission volume. For example, any density or temperature diagnostics made from the a l u m i n u m satellite line ratios represent only average conditions. T y p i c a l a l u m i n u m reference spectra appear in Figures 6.2 and 6.3. Figure 6.3 shows that the spectrometer was able to resolve the k from the l,j satellite lines of Li-like a l u m i n u m . T h i s suggests a resolving power A / A A of 2300. Higher resolving power could not be demonstrated as no other closely spaced lines appeared i n this region. U s i n g the wavelength calibration equation (5.3) and the range of values obtained for a and R during processing, the derived uncertainty o f the spectrometer is ± 1 e V in the region of the K-edge. T h i s estimate is supported by two pieces of evidence. The wavelength calibration w i t h typical values of a and R is applied to a random aluminum 110 Chapter 6. Results and Discussion of X-Ray Spectra 111 Aluminum Tantalum Figure 6.1: A print of a sample a l u m i n u m and tantalum exposure. lines are labeled (refer to Figures 6.2, 6.3, and 6.8). Some of the known Chapter 6. Results and Discussion of X-Ray Spectra 112 Figure 6.2: Sample aluminum reference spectrum. Some common lines are labeled. Chapter 6. Results and Discussion of X-Ray Spectra 113 Figure 6.3: Enlarged aluminum reference spectrum showing the clearly resolved satellite lines qr, a-d, and klj of A l XI. Chapter 6. Results and Discussion of X-Ray Spectra 114 spectrum and agreement is seen to be excellent (see Table 6.1). T h e second piece of evidence is the consistent appearance of the K-edge on the sample spectra between 1559.2 and 1600.8 e V , visible prior to correction for the B-10 attenuation. E a c h microdensitometer scan has had the same film calibration curve applied in converting film exposure to relative intensity. A l t h o u g h the accuracy of this curve is only 20% for the intensities seen here, comparisons can still be made between the different spectra. T h e spectra are discussed i n order of decreasing atomic number. 6.1.1 B i s m u t h , L e a d and G o l d T h e spectra of bismuth, lead and gold appear i n Figures 6.4 to 6.6. T h e observed emissions are remarkably similar given the range of atomic numbers that these elements span, and therefore the spectra are discussed together. T h e most prominent feature of the three spectra is the appearance of transition arrays between 1545 and 1635 e V . It is interesting to note that these arrays appear at nearly the same center energies and w i t h almost identical widths. T h i s would suggest that the arrays come from transitions which are relatively insensitive to the nuclear structure and number of shielding electrons, and hence are derived from inner subshell transitions. Previous work [109, 110, 111, 120] in analysing the spectra from these elements i n different wavelength regions allows certain transitions to be ruled out. Candidates for the observed arrays are the N i I like ion transitions [111] of 3 d - 3 d 1 0 9 / 2 nf , 3d -3d| nf 1 0 5 / 2 / 2 , and 3 p 3 d - 3 p i j 3 d n d 6 7 / 2 1 0 1 0 / 2 5 / 2 with principal quantum number n equal to 4 or 5. T h e average intensity of the bismuth emission i n the ± 3 0 e V region around the K edge is nearly 1.5 times that of lead and gold, although a l l three display roughly the same amount of fractional variation. T h e prominent line structure of bismuth in the region centered at 1555 e V , which is not seen in lead or gold, make the latter two more attractive as backlighters. In fact, bismuth has been used previously as a backlighter [11, 12] i n the Chapter 6. Results and Discussion of X-Ray Spectra Transition Experiment (eV) Theory (eV) #12 #11 1706.2 1684.1 1705.5 1684.2 He 1598.3 1589.3 1588.5 1580.1 1578.1 1575.5 1574.5 1536.3 1598.7 1588.5 1588.0 a y s,t q,r a-d k 1J o,P 1580.0 1578.3 1576.1 1575.3 1536.4 115 • Table 6.1: E x p e r i m e n t a l and theoretical energies of lines of an a l u m i n u m spectrum when the spectral calibration equation is applied. Chapter 6. Results and Discussion of X-Ray Spectra 116 Bismuth 5.0 ~i—i—i—F—i—i—i—i—I—r T 1111I111 J - 4.0 CO -4—' c Z) -e 3.0 < 0 CD > 2.0 •4—' _CJJ CD DC 1.0 Q Q I—i—"—i—i—i—i—i—i—i—l—i—i—i—J i i i i i l t i i i i i i i -t i i i I i i i i i i i i i i i 725.0 1675.0 1625.0 1575.0 1525.0 Energy (eV) Figure 6.4: Sample bismuth spectrum. The position of the aluminum K-edge is shown by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 117 Lead 5.0 ~i—i—i—i—i—i—i—i—I—i—r 4.0 CO -e 3.0 < CO c CD CD > 2.0 CD DC 1.0 0.0 1725.0 J—i—i—I—I—I—i—i—i 1675.0 I I I i j 1625.0 1575.0 Energy (eV) 1525.0 Figure 6.5: Sample lead spectrum. The position of the aluminum K-edge is shown by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 118 K-edge region, but the sample bismuth spectrum obtained by these researchers had less spectral resolution. It is unknown how the spectral structure seen in Figure 6.4 would affect the results of this previous work. On the basis of observed line structure, lead and gold provide more suitable backlighter characteristics than bismuth. This structure, however, will be enhanced or degraded depending on irradiance, and therefore none or all of these elements may turn out to be useful at other irradiances. 6.1.2 Tungsten The observed x-ray emission spectrum of tungsten appears in Figure 6.7. Analyses exist for this element in the ranges of 2.6-3.3 keV [109] and 1.6-1.9 keV [118]. Four of the spectral lines seen in this work have been identified previously, and have been labeled according to Kef. [118]. The g,h and i transitions originate from the Cu I like 3d 4s-3d 4s4p array. The (3 line group belongs to the 3d 4s -3d 4s 4p Zn I like 10 9 10 2 9 2 transition. The lower energy lines remain unclassified, but most likely belong to inner subshell transition arrays. As with previous elements, the line structure around the K edge makes tungsten unsuitable as a backlighter, even though it is one of the more intense emitters studied in this work. 6.1.3 Tantalum Tantalum, whose x-ray spectrum appears in Figure 6.8, is one of the more actively studied elements in the classification of unresolved transition arrays [111, 117, 118, 119]. Again, the lines are labeled in accordance with Ref. [118]. The Cu I like transition, 3d 4s-3d 4s4p, labeled with lowercase letters, gives rise to 10 9 the majority of lines in this energy region. The two N i I like resonance lines are the most prominent emission lines, in agreement with previous work [109]—[121]. A few lines escape Chapter 6. Results and Discussion of X-Ray Spectra Figure 6.6: Sample gold spectrum. the dotted line. 119 T h e position of the a l u m i n u m K-edge is shown by Chapter 6. Results and Discussion of X-Ray Spectra Figure 6.7: Sample tungsten spectrum. The position of the aluminum K-edge is by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 121 Tantalum 16.0 ~i—r—\—r i—i—i—r- 14.0 12.0 Z) JQ < CO c 10.0 8.0 CD 0 > > 05 CD 6.0 rr 4.0 2.0 h J—i—i—i—i—i—i—I—i—i—i—i—i—i—i i i I i i i 0.0 1725.0 1675.0 1625.0 i i i i t_ i I i i I 1575.0 1525.0 Energy (eV) Figure 6.8: Sample tantalum spectrum. The position of the aluminum K-edge is shown by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 122 classification, but like tungsten, probably originate from inner subshell transition arrays. Unfortunately, tantalum displays more line structure near 1560 e V than is acceptable for a backlighter element. 6.1.4 Samarium A sample samarium spectrum appears in Figure 6.9. B u r k h a l t e r et al. [116] have obtained spectra of the rare earth elements S m through Y b , but have published results only for G d . Comparison of the S m spectrum w i t h the published G d spectrum does not lead to any line identifications. T h e intense transition array structure over the observed energy range makes the S m spectrum ideal as d a t a for transition array modeling, but useless as a backlighter. 6.1.5 Tin Spectra of laser produced t i n plasmas below 1240 e V have been published previously [113, 114, 115] but Figure 6.10 shows the first spectrum obtained above 1500 e V . Transition arrays appear superimposed on a sloped background of bremsstrahlung radiation. Reference [113] allows certain transitions to be ruled out, so that candidates for the observed emission are transition arrays of the types 3 d - 3 d fe fe_1 5 f and 3 d - 3 d fe fc_1 5p in Sn X X X I I I to Sn X X I I I ions. W h i l e the extensive transition array structures around the K-edge render Sn unsuitable as a backlighter, they do provide a testbed for further transition array modeling. 6.1.6 Silver Previous spectral analyses for this element have concentrated on a much lower energy regime [123, 126, 124], whereas Figure 6.IT shows the spectrum in the region 1 5 0 0 - 1725 Chapter 6. Results and Discussion of X-Ray Spectra 123 Figure 6.9: Sample samarium spectrum. The position of the aluminum K-edge is shown by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra Tin I Q I i ' ' ' i i i' i i i i 1725.0 1675.0 iii iii i , 1625.0 1575.0 Energy (eV) i; i • • • i 1525.0 Figure 6.10: Sample t i n spectrum. The position of the a l u m i n u m K-edge is shown the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 125 eV. Reference [126] suggests that the observed spectrum may originate from 2s-3p, 2p-3d, and 2p-3s transitions in A g + 3 9 . However, detailed line classifications do not exist, and none are offered here. Theoretically, the temporal evolution of the intensity of the 1560 eV line of silver could be used to monitor the a l u m i n u m K-edge shift. A s the K-edge shifts towards the red, the lower energy components experience less attenuation, and the line would appear more broad and w i t h increased intensity. Conversely, as the K-edge shifts towards the blue, the line would suffer greater attenuation, and would appear more narrow and with less intensity. T h i s method, however, suffers from several serious drawbacks. F i r s t of all, the w i d t h of the emission line is much less than the magnitude of the expected edge shift, making this technique useful only during the initial moments of the experiment, and not over the entire measurement period. Furthermore, the K-edge measurement is of the time-resolved form, and thus the temporal dependence of the emission of the 1560 e V line of silver would be required, adding yet another level of complexity to the measurement. Finally, the attenuation of this line only provides K-edge information right at 1560 e V , whereas a spectrally flat backlighter would provide a K-edge measurement over a much broader spectral range. However, silver may find a use as a spectral calibration source. M a n y intense, well defined emission lines surround the K-edge, and the resonance at 1560 e V may serve as a convenient reference point by which to calibrate spectral measurements. 6.1.7 Niobium Spectra of highly ionized N b have been obtained from vacuum spark measurements, [123] and like those of silver, these measurements have focussed on a much lower energy region of the spectrum. Figure 6.12 presents the N b spectrum in the region of 1500 - Chapter 6. Results and Discussion of X-Ray Spectra 126 Silver 5.0 ~~| i i i—i—i—i—i—i—i—]—i—i—i—i—i—i—i—i—i—|—r~ -i—i—i—r~ 4.0 cn 'E -Q < 3.0 CO c CD CD 2.0 as CD rr 1.0 0.0 1725.0 1675.0 1625.0 1575.0 Energy (eV) 1525.0 Figure 6.11: Sample silver spectrum. The position of the a l u m i n u m K-edge is shown by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 127 1725 e V . N o line identifications can be made, although the existence of transition arrays seems likely, as is the case w i t h medium to high atomic number elements. A s atomic number decreases, however, there should be a trend toward non-overlapping resonance line structures and transition array structures should become less frequent. 6.1.8 Yttrium Previous analyses of Y t from vacuum spark and laser produced plasma spectra [123, 126] have been done in the 120 - 1200 e V range. Figure 6.13 shows the spectrum in a much higher energy region. Comparison w i t h previous work does not lead directly to any line classifications. T h e relatively small amount of line structure coupled w i t h small fractional intensity variations make Y t a good candidate for a backlighter if the low emission levels seen here can be increased by increasing the laser irradiance. 6.1.9 Germanium Figure 6.14 presents the spectrum of Ge in the range 1500 - 1725 e V . N o previous analyses exist for the spectrum of this element in this energy range, and no line classifications are offered here. T h e spectrum most likely consists of unresolved transition arrays superimposed upon a sloped background of bremsstrahlung radiation. T h e intense spectral emission over the entire range makes G e unusable as a backlighter source, which is unfortunate given its unparalleled level of brightness. Interestingly, G e has been used as a backlighter for shock waves studies where spectral structure is not i n h i b i t i n g [127]. T h i s spectrum provides a rich source of d a t a to be modeled for further transition array studies. Chapter 6. Results and Discussion of X-Ray Spectra 128 Niobium 5.0 - i — i — i — i — i — r - i—i—i—i—i—i—i—i—i—j—i—i—i—i—i—i—i—i—i—j—i—i—i—i— 4.0 'c 3 -E 3.0 < CO c CD CD > 2.0 03 CD DC 1.0 | Q I—I 1 1 — — II 1 — — I I 1— — II—I—I—I I I I I I I I I ' ' 1725.0 1675.0 1625.0 I i i ii i i i i i I i i i. 1575.0 1525.0 Energy (eV) Figure 6.12: Sample n i o b i u m spectrum. T h e position of the a l u m i n u m K-edge is shown by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 129 Yttrium 5.0 4.0 CO 'c -a 3.0 < CO c CD CD 2.0 03 CD rr 1.0 0.0 1725.0 • • • . i 1675.0 J—i—i—i—I—i—i—i—i—i i i i i i i i } i • i i i i i 1625.0 1575.0 Energy (eV) 1525.0 Figure 6.13: Sample y t t r i u m spectrum. T h e position of the a l u m i n u m K-edge is shown by the dotted line. Chapter 6. Results and Discussion of X-Ray Spectra 130 Germanium 30.0 25.0 CO E 20.0 JD < % c 15.0 o c: 0 > 03 •g CD io.o DC 5.0 0.0 1725.0 —I—I—I—I—I—I—I—I 1 I 1675.0 I 1625.0 —1—I—I—I } I 1 I I l_ 1575.0 1525.0 Energy (eV) Figure 6.14: Sample germanium spectrum. shown by the dotted line. The position of the aluminum K-edge is Chapter 6. Results and Discussion of X-Ray Spectra 6.1.10 131 Magnesium A typical M g spectrum appears in F i g . 6.15. T h e emission lines at 1700, 1660 and 1580 eV are the He^, H e , and He^ resonance lines of M g X I , respectively. T h e structure at 7 the low energy edge of the recorded spectrum consists of the 1-6 satellite lines of M g X . Aside from the line structure w i t h i n 30 e V on either side of the K-edge, the lack of intensity of the observed spectrum makes magnesium unusable as a backlighter. 6.2 Discussion T h e aim of this work was to find a suitably flat and bright spectral source for use as a backlighter. Classification of spectral features for each element was never a priority. T h e computational resources needed to model the unresolved transition arrays present in these spectra were not available, and hence the primary means by which lines are identified for high-Z elements could not be used. Moreover, much of the previous work done in classifying these spectra has been concentrated i n different spectral regions than the one studied in this work. T h e multitude of possible transitions originating from the many coexisting ionization stages in the plasma makes line identifications a labourious task. Certain trends can be employed to aid in the identification process, but these are often of limited use. One method employs the systematic Z 2 dependence of atomic features on atomic number by means of a Moseley plot [116], as illustrated in F i g . 6.16 (taken from Reference [116]). T h e square root of the transition energy for a select transition is plotted against atomic number. Interpolation from the resulting nearly linear fit is used to determine the energy of the same transition for different elements. O f course, this method is only approximate, and relies heavily on previous calculations. Another technique used to unravel atomic spectra involves the systematic addition Chapter 6. Results and Discussion of X-Ray Spectra Magnesium 14.0 ~i——i—i—i—i—i—r i i i i i I—i—i—i—i—i—i—i—i—i—i—i—i—m—i—i—i—r~ 12.0 — CO 10.0 "c Z) -Q < 8 0 'co c CD I CD 6.0 .> " + •— < _C0 CD DC 4.0 2.0 0.0 1725.0 1675.0 1625.0 1575.0/ Energy (eV) Figure 6.15: Sample magnesium spectrum. shown by the dotted line. i i i i ii 1525.0 The position of the aluminum K-edge Chapter 6. Results and Discussion of X-Ray Spectra 133 45 -I 1 28 30 : I 40 I 50 i 60 i 70 ATOMIC NUMBER Figure 6.16: Moseley plot used to aid in line classification. Chapter 6. Results and Discussion of X-Ray Spectra 134 and subtraction of spectator electrons [113, 116, 121]. T h e energy of a specific transition for a single ionization state is calculated. Spectator electrons added to or subtracted from various shells w i l l shift the transition energy by constant amount, allowing identification of the transition in different ionization states. For example, removal of successive 3d electrons i n G d + 3 6 shifts the 3 d - 4 f transition up by +39 e V . T h i s has been used to identify 3d-4f lines in G d + 3 T and G d + 3 8 [116]. T h i s method, however, suffers from a form of de- generacy. Different transitions originating in different ionization stages can lead to nearly identical energies [109]. A s well, line shifting due to asymmetric Doppler broadening [113] can shift spectral lines, m a k i n g the spectator electron technique inaccurate. It should be noted that the spectral classifications which are offered in the previous sections are by no means final. Instead, they have been determined by extrapolation and comparison w i t h previous analyses, and simply serve as a starting point for future work. One of the most important aspects of these spectra, as was alluded to earlier, was the dependence of line structure and continuum emission strength on laser irradiance. Increasing the irradiance should lead to greater ion populations in the higher ionization states. Unfortunately, much of the spectra seen in this work remains unclassified as to specific line transitions and even from which ionization species they originate. Hence, increasing the laser power could have one of two effects. T h e line structure already observed could become enhanced, rendering a l l of the elements studied here unusable, or, the relative populations of the different ionization stages could shift so that the line emission would appear in a different spectral region. T h e spectra observed here would then seem to flatten out, as the plasma becomes more like a continuum emitter. Likely candidates to be used as a backlighter source are narrowed to three. Lead and gold show reasonable levels of brightness, w i t h little observed line structure. However, the intensity variation displayed by these elements around 1560 e V could prove to be a problem. Y t t r i u m , on the other hand, shows excellent spectral flatness, and would be Chapter 6. Results and Discussion of X-Ray Spectra 135 suitable provided it is sufficiently bright. Finally, although not usable as a backlighter, silver seems promising as an emission source by which to calibrate spectrograph^ instruments in an experimental situation. Chapter 7 Summary and Conclusions T h i s final chapter presents a summary of the work undertaken i n this thesis project and the conclusions which can be drawn from it. Suggestions for future work are also offered. 7.1 Summary of Present Work T h e temporal evolution of the K photoabsorption edge of aluminum under shock compressed conditions allows the study of the atomic physics of hot, dense matter. This thesis work has concentrated on providing the necessary background information for the design of an efficient experiment to probe this physics regime. In order to decouple to effects of density, temperature and ionization state on the position of the K-edge, four separate shock compression schemes are used. Each scheme couples a tailored laser pulse to a specific C H - A l - S i target i n order to produce a well characterized plasma state w i t h i n the aluminum layer. B y comparing three of these states back to a reference state, density, temperature and ionization state effects can be individually studied. T h e shock compression schemes are designed w i t h a one dimensional hydrodynamic code incorporating non-local thermodynamic equilibrium physics. T h e N L T E model considers the atomic physics of the C H layer only, but line and continuum radiation, as well as inner-shell photoionization processes are considered for a l l cells i n the simulation. T h e target designs represent the most efficient thicknesses in terms of absorption measurement signal, radiation preheat, and laser pulse length. T h e laser pulse is modeled following 136 Chapter 7. Summary and Conclusions 137 presently available shaped pulses. T h e second aspect of a K-edge experiment involves the x-ray emission used for the photoabsorption measurement. In search of a suitable source, the x-ray emission between 1525 and 1725 e V of 12 elements between M g and B i has been obtained w i t h a P E T crystal spectrometer. Accuracy of the spectroscopic measurement is ± 1 e V . Suitable backlighter sources are identified, and a wealth of new data for unresolved transition array studies is produced. 7.2 Conclusions T h e pulse-target schemes summarized in Chapter 4 present a viable scheme for study of density, temperature and ionization effects on K-edge position. A l l steady state plasma conditions used to generate K-edge data have sufficiently long lifetimes and small variations i n hydrodynamic conditions to generate well characterized plasma states. Candidates for suitable x-ray emission sources are narrowed to three. Lead and gold show reasonable levels of brightness, w i t h little observed line structure. Y t t r i u m , on the other hand, shows excellent spectral flatness, but displays rather low emission levels. Finally, although not usable as a backlighter, silver seems promising as an emission source by which to calibrate spectrographic instruments i n an experimental situation. 7.3 Future Work Ideally, the N L T E calculation would include a l l of the layers of the target, as opposed to just the C H layer as considered here. T h i s would result in a far more accurate calculation of the ionization in the a l u m i n u m layer than is currently provided by the Q E O S or S E S A M E tables. A s well, some improvement would be expected on the absolute levels of radiation emission, as now line radiation could originate from the a l u m i n u m and silicon Chapter 7. Summary and Conclusions 138 layers. To gain further understanding of target preheat, fast electron generation could also be incorporated into the code. T h i s effect may prove significant for higher irradiances, and eventually needs to be considered. T h e most important future work to be carried out in terms of the x-ray spectra would involve a study of dependence of emission on laser irradiance. Different irradiances can be expected to shift relative ion abundances in the different ionization states, thus enhancing or degrading line emission and/or intensity levels in the energy region surrounding the K-edge of a l u m i n u m . Finally, a different target configuration than considered here, consisting of C H - S i - A l Si, may prove more effective in reducing levels of radiative preheat. Pressure enhancement could still be achieved at the CH-silicon interface, but the increased radiation absorption by the silicon layer may allow a much cleaner shock compression of the aluminum layer. Bibliography [1] V . P. Silin, Sov. Phys. J . E . T . P . 21, 1127 (1965) ; M . V . G o l d m a n , Annals of Physics 38, 117 (1966) ; P. K . K a w and J . M . Dawson, Phys. F l u i d s 12, 2586 (1969) [2] F . J . M c C l u n g and R . W . Hellwarth, J . A p p l . Phys. 33, 828 (1962) [3] C . E . Ragan, Phys. Rev. A 21, 458 (1980) [4] H . Hora, Plasmas at High Temperatures and Densities Applications and Implication of Laser-Plasma Interaction, Springer-Verlag, B e r l i n (1991) [5] P. Mulser, R . Sigel, and S. W i t k o w s k i , Phys. Rep. 6, 187 (1973) [6] P. Celliers, A . N g , G . X u , and A . Forsman, Phys. R e v . Lett. 68, 2305 (1992) [7] L . Vesser and J . Solem, Phys. R e v . Lett. 40, 1291 (1978) [8] R . J . Trainor, J . W . Shaner, J . W . Aubebach, and N . C . Holmes, Phys. Rev. Lett. 42,1154 (1979) [9] F . Gottet, J . P. R o m a i n , R . Fabbro, and B . Faral, Phys. Rev. Lett. 52, 1884 (1984) [10] A . N g , D . Parfeniuk, and L . DaSilva," Phys. R e v . Lett. 54, 2604 (1985) [11] T . A . H a l l , A . Djaoui, R . W . Eason, C . L . Jackson, and B . Shiwar, Phys. R e v . Lett. 60, 2034 (1988) [12] R . W . Eason, D . K . Bradley, J . D . Kilkenny, and G . N . Greaves, J . Phys. C 17, 5067 (1984) [13] A . N g , Y . G a z i t , F . A d a m s , D . Parfeniuk, P. Celliers, and L . DaSilva, P l a s m a Phys. Cont. N u c l . Fusion 3, 187 (1986) [14] D . K . Bradley, J . Kilkenny, S. J . Rose, and J . D ; Hares, Phys. R e v . Lett. 59, 2995 (1987) [15] L . G . P a r a t t and E . L . Jossem, Phys. R e v . 97, 916 (1955) [16] L . DaSilva, A . N g , B . K . G o d w a l , G . C h i u , F . Cottet, M . Richardson, P. Jaanimagi, and Y . T . Lee, Phys. Rev. Lett. 62, 1623 (1989) 139 Bibliography 140 [17] B . K . G o d w a l , A . N g , L . D a S i l v a , Y . T . Lee, and D . A . L i b e r m a n , Phys. Rev. A , 40, 4521 (1989) [18] F . Perrot and M . W . C . Dharma-Wardana, Phys. Rev. Lett. 71, 797 (1993) [19] W . J . Nellis, J . A . Moriarty, N . W . Ashcroft, N . C Holmes, and G . R . Gathers, Phys. R e v . L e t t . 6 0 , 1414 (1988) [20] B . L . Henke, P. Lee, T . J . Tanaka, R . L . Shimabukuro, and B . J . Fujikawa, A t . D a t a N u c l . D a t a Tables, 27, 1 (1982) [21] J . C . Slater, Quantum V o l . 1. Theory of Atomic Structure, [22] J . Stewart and K . P y a t t , Astrophys. J . 144, M c G r a w - H i l l , New Y o r k , 1960, 1203 (1966) [23] D . Duston, R . W . C l a r k , J . Davis, and J . P. Apruzese, Phys. Rev. A , 27, 1441 (1983) [24] J . P. Apruzese, J . Davis, D . Duston, and K . G . Whitney, J . Quant. Spectros. R a d i a t . Transfer 23, 479 (1980) [25] P. Collela and P. R . Woodward, J . C o m p . Phys. 54, 175 (1984) [26] P. Collela and H . G l a z , J . C o m p . Phys. 59, 264 (1985) [27] W . K r u e r , The Physics (1991) of Laser-Plasma Interactions, Addison-Wesely, New York [28] P. M o r a , Phys. F l u i d s 25, 1051 (1982) [29] P , M o r a and R . Pellat, Phys. F l u i d s 22, 2300 (1979) [3.0] C . V a n Kessel, Z. Naturforsch 30a, 1581 (1975) [31] Y . B . Zel'dovich and Y . P. Raizer, Physics of Shockwaves and High Temperature Hydrodynamic Phenomena, W . D . Hayes and R . F . Probstein eds., Academic Press, New Y o r k (1966) [32] J . S. Rinehart, Stress Transients in Solids, Addison-Weseley, New Y o r k (1965) [33] L . R . Vesser, J . Solem, and A . Lieber, A p p l . Phys. Lett. 35, 761 (1979) [34] N . C . Holmes, R . J . Trainor, R . A . Anderson, L . R . Vesser, and G . A . Reeves, Shock Waves in Condensed Matter, W . J . Nellis, L . Seaman and R . D . G r a h a m eds., A I P Press, New Y o r k (1982), p. 160 Bibliography 141 [35] F . Cottet, M . Hallouin, J . P. R o m a i n , R . Fabbro, B . Faral, and H . Pepin, A p p l . Phys. Lett. 47, 678 (1985) [36] N . L . Hickerson, High-Velocity New Y o r k (1970), p. 23 Impact Phenomena, R . K i n s l o w ed., Academic Press, [37] S. Eliezer, A . K r u m b e i n , and D . Salzmann, J . Phys. D 11, 1963 (1978) [38] D . Salzmann, H . Szichman, and A . K r u m b e i n , Phys. F l u i d s 30, 515 (1987) [39] D . Salzmann and A . K r u m b e i n , J . A p p l . Phys. 49, 3229 (1978) [40] D . Duston, R . Clark, and J . Davis, Phys. Rev. A 31, 3220 (1985) [41] M . W . C . Dharma-Wardana, Private Communication [42] T . M o c h i z u k i , T . Yabe, K . Okada, M . Hamada, N . Ikeda, S. Kiyokawa, and C . Yamanaka, Phys. Rev. A 33, 525 (1986) [43] C . De Michelis and M . M a t t i o l i , Nuclear Fusion 21, 677 (1981) [44] R . M c W h i r t e r , Plasma Diagnostic Academic Press, New York (1966) [45] H . G r i e m , Techniques, R . Huddlestone and S. Leonard eds., Plasma Spectroscopy, M c G r a w - H i l l , New York (1964) [46] R . M c W h i r t e r and H . P. Summers, Applied Atomic Collision Physics Huddlestone and S. Leonard eds., Academic Press, New Y o r k (1984) Vol. 2, R . [47] D . Salzmann, J . Quant. Spectros. Radiat. Transfer 27, 359 (1981) [48] G . B . R y b i c k i and A . P. L i g h t m a n , and Sons, New York (1979) Radiative Processes in Astrophysics, John W i l e y [49] P. G . Burkhalter, F . C . Y o u n g , B . H . R i p i n , S. E . Bodner, and D . J . Nagel, Phys. Rev. A 15, 1191 (1977) [50] B . F . Young, W . H . Goldstein, A . L . Osterheld, R . E . Stewart, and G . E . Busch, J . Phys. B 22, L533 (1989) [51] M . H . K e y and R . J . Hutcheon, Advances in Atomic and Molecular Bates arid B . Bederson eds., V o l . 16, Academic Press, New York (1980) Physics, D . [52] W . J . Karzas and R . Latter, Astrophys. J . Suppl. Ser. 6, 167 (1961) [53] D . Mihalas and B . W . Mihalas, Foundations University Press, New York (1984) of Radiation Hydrodynamics, Oxford Bibliography 142 [54] J . P. Christiansen, D . E . T . F . Ashby, and K . V . Roberts, C o m p u t . Phys. C o m m . 7, 271 (1974) 1 [55] P. Celliers, P h . D . Thesis, U . B . C , 1987 (unpublished) [56] L . D a S i l v a , P h . D . Thesis, U . B . C , 1988 (unpublished) [57] W . M . Manheimer and D . G . Colombant, Phys. F l u i d s 25, 1644 (1982) [58] J . H . Nuckolls, L . W o o d , A . Thiessen and G . Zimmerman, Nature 239, 139 (1972) [59] C . E . M a x and K . Estabrook, C o m m . P l a s m a Phys. Controlled Fusion 5, 239 (198G) ; W . K r u e r , C o m m . P l a s m a Phys. Controlled Fusion 6, 167 (1981) [60] L . Spitzer and R . H a r m , Phys. Rev. 89, 977 (1952) [61] J . R . Stallcop and K . W . B i l l m a n , P l a s m a Phys. 16, 722 (1973) [62] M . B o r n and E . Wolf, Principles of Optics, 6th ed., Pergamon Press, Oxford (1980) [63] A . N g , P. Celliers, A . Forsman, R . More, Y . Lee, F . Perrot, M . Dharma-Wardana, and G . Rinker, Phys. Rev. Lett 72, 3351 (1994) [64] Y . T . Lee, J . Quant. Spectros. R a d i a t . Transfer 38, 131 (1987) [65] R . K . Landshoff and J . D . Perez, Phys. Rev. A 20, 1619 (1976) [66] M . J . Seaton, M o n . Not. R . A s t r o n . Soc. 119, 81 (1959) [67] A . Burgess, Astrophys. J . 141, 1588 (1965) [68] E . O r a n and J . Davis, J . A p p l . Phys. 45, 2480 (1974) [69] R . Mewe, A s t r o n . Astrophys. 20, 215 (1980) [70] R . C . Canfield, A . N . M c C l y m o n t , and R . C . Puetter, Methods in Radiative Transfer, Cambridge University Press, New Y o r k (1984) [71] R . F . R e i l m a n and S. T . Manson, Astrophys. J . Suppl. Ser. 40, 815 (1979) [72] P. Celliers and K . E i d m a n n , Phys. Rev. A 41, 3270 (1990) [73] F . Biggs and R . L i g h t h i l l , Sandia Laboratories Report N o . S C - R R - 7 1 0507, 1971 (unpublished) [74] W . L . Wiese, Atomic Transition U . S . G o v t . P r i n t . Off. (1966) Probabilities - a Critical Data Compilation, Wash. Bibliography 143 [75] B . F . Y o u n g , R . Stewart, and G . Busch, Phys. Rev. Lett. 61, 2851 (1988) ; J . Bailey, R . Stewart, J . Kilkenney, T . P h i l i p s , and R . Lee, J . Phys. B 19, 2639 (1986) ; R . Kauffman, R . Lee, and K . Estabrook, Phys. Rev. A 35, 4286 (1986) [76] R . W . Lee, B . L . W h i t t e n , and R . E . Stout, J . Quant. Spectroscopty R a d i a t . Transfer 32, 91 (1984) [77] B . E d l e n , Phys. Scr. 19, 255 (1976) [78] L . Vainstein and U . Safronova, Phys. Scr. 31, 519 (1985) [79] B . Edlen, Phys. Scr. 17, 565 (1978) [80] J . Scofield, Lawrence Livermore Laboratory Internal Document U C D - 1 6 8 4 8 (1975) [81] R . Kelly, J . Phys. C h e m . Ref. D a t a , V o l . 16, Suppl. 1, 1967 [82] W . L o t z , Z . Phys. 206, 205 (1967); 216, 241 (1968); 220, 486 (1969) [83] L . Spitzer, A p . J . 107, 6 (1948) [84] D . Mihalas,Stellar Atmospheres, 2nd ed., W . Freeman, San Francisco (1978) [85] H . V a n Regemorter, A p . J . 136, 906 (1962) [86] D . Cochrane and R . M c W h i r t e r , Phys. Scr. 28, 25 (1983) [87] L . Vainstein and U . Safronova, A t . D a t a N u c l . D a t a Tables 21, 49 (1978) [88] A . G a b r i e l , M o n . Not. Royal Astro. Soc. 160, 99 (1972) [89] S E S A M E Library, Los Alamos N a t i o n a l Laboratory, Los Alamos [90] R . M . More, K . H . Warren, D . A . Y o u n g , and G . B . Z i m m e r m a n , Phys. F l u i d s 31, 3059 (1988) [91] S. Jackel, D . Salzmann, A . K r u m b e i n , and S. Eliezer, Phys. F l u i d s 26, 3138 (1983) [92] S. Eliezer, L . Gitter, A . K r u m b e i n , and D . Salzmann, Phys. Lett. 86A, 464 (1981) [93] D . Salzmann, S. Eliezer, A . K r u m b e i n , and L . Gitter, Phys. Rev. A 28, 1738 (1983) [94] R . Marchand, C . Capjack, and Y . T . Lee, Phys. F l u i d s 31, 949 (1988) [95] R . J . Trainor and Y . T . Lee, Phys. F l u i d s 25, 1898 (1982) [96] D . Duston, R . C l a r k , J . Davis, and J . Apruzese, Phys. Rev. A 27, 1441 (1983) Bibliography 144 [97] R. Schmalz, J. Meyer-ter-Vehn, and R. Ramis, Phys. Rev. A 34, 2177 (1986) [98] E . McLean, S. Gold, J. Stamper, R. Whitlock, H . Griem, S. Obenschain, S. Gitomer, and M . Matzen, Phys. Rev. Lett. 45, 1246 (1980) [99] A. Ng, D. Parfeniuk, L . DaSilva, D . Pasini, Phys. Fluids 28, 2915 (1985) [100] A . R. Bell, R. G . Evans, and D. J. Nicholas, Phys. Rev. Lett 46, 243 (1981); G . J. Rickard, A . R. Bell, and E . M . Epperlein, Phys. Rev. Lett. 62, 2687 (1989) [101] L A S N E X Computer Code, property of Lawrence Livermore National Laboratory. [102] Gentek Inc., 2625 Dalton St., Ste-Foy Quebec, G1P 3S9, Canada [103] Scientech Inc., 5649 Arapahoe Ave., Boulder C O 80303, U S A [104] Hammamatsu photodiode model R1193U [105] B-10 Foil, Alexander Vacuum Research, Greenfield Mass., U S A [106] A . Forsman, M.Sc. Thesis, U . B . C . 1992 (unpublished) [107] D . Duston, J.E. Rogerson, J. Davis, and M . Blaha Phys. Rev. A 28, 2968, 1983 [108] C. Chenais-Popovics, C. Fievet, J.P. Geindre and J.C. Gauthier, Phys. Rev. A 40, 3194 1989 [109] A . Zigler, M . Klapisch, and P. Mandelbaum, Phys. Lett. A 117, 31 (1986) 110] S. Kiyokawa, T. Yabe, and K . Okada, Phys. Rev. Lett. 54, 1999 (1985) 111] N . Tragin, J. Geindre, P. Monier, J. Gauthier, C. Chenais-Popovics, J. Wyart, and C Bauche-Arnoult, Phys. Scr. 37, 72 (1988) r 112] M . Klapisch, H . Mora, A . Bar-Shalom, P. Mandelbaum, J. Schwob, A . Zigler, and S. Jackel, Phys. Lett. 79A, 67 (1980) 113] P. Burkhalter, U . Feldman, and.R. Cowan, J . Opt. Soc. A m . 64, 1058 (1974) 114] B . Edlen, Atomic Spectra, Springer-Verlag, Berlin (1964) 115] E . Alexander, M . Even-Zohar, B . Fraenkel, and S. Goldsmith, J. Opt. Soc. A m . 61, 508 (1971) 116] P. Burkhalter, D. Nagel, R. Whitlock, Phys. Rev. A 9, 2331 (1974) 117] P. Audebert, J. Gauthier, J. Giendre, C. Chenais-Popovics, J. Wyart, Phys. Rev. A 32, 409 (1985) Bibliography 145 [118] P. Mandelbaum, M . Klapisch, A . Bar-Shalom, J. Schwob, and A. Zigler, Phys. Scr. 27, 39 (1982) [119] M . Klapisch, P. Mandelbaum, A . Zigler, Phys. Scr. 34, 51 (1985) [120] M . Busquet, D . Pain, J. Bauche, and E . Luc-Koenig, Phys. Scr. 31, 137 (1985) [121] C. Bauche-Arnoult, E . Luc-Koenig, J. Wyart, J. Giendre, P. Audebert, and C. Chenais-Popovics, Phys. Rev. A 33, 791 (1986) [122] M . Klapisch, P. Mandelbaum, A . Bar-Shalom, A. Zigler, and S. Jackel, J. Opt. Soc. A m . 71, 1276 (1981) [123] N . Schweitzer, M . Klapisch, J. Schwob, M . Finkenthal, A . Bar-Shalom, P. Mandelbaum, and B . Fraenkel, J. Opt. Soc. A m . 71, 219 (1981) [124] M . Even-Zohar, Some euv Spectra from Laser-Produced Plasmas, Rep. C E A - E U R A T O M FC-768 [125] M . Klapisch, E . Meroz, P. Mandelbaum, A . Zigler, C. Bauche-Arnoult, and J. Bauche, Phys. Rev. A 25, 2391 (1982) [126] P. Burkhalter, C. Dozier, D. Nagel, Phys. Rev. A 15, 700 (1977) [127] T. Endo, K . Shigemori, H . Azechi, K . Mima, M . Sato, Phys. Rev. Lett. 74, 3608 (1995)
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- X-ray sources and shock compression schemes for photoabsorption...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
X-ray sources and shock compression schemes for photoabsorption edge spectroscopy Dyke, Steven R. 1995
pdf
Page Metadata
Item Metadata
Title | X-ray sources and shock compression schemes for photoabsorption edge spectroscopy |
Creator |
Dyke, Steven R. |
Date Issued | 1995 |
Description | A scheme to probe the temporal evolution of the K-shell photoabsorption edge in shock compressed aluminum is presented in two components. In the first component, a one dimensional hydrodynamic code coupled to a non-local thermodynamic equilibrium calculation is developed and used to model shock propagation and target preheat in laser irradiated CH-Al-Si multilayer targets. Single and multiple shocks are used to create conditions in the aluminum by which density, temperature and ionization effects on K-shell photoabsorption can be isolated. The final result consists of the optimal laser parameters and target designs for creating well characterized plasma states, necessary for accurate atomic modeling of the K-edge. The second component describes the results of x-ray spectra observed from 12 elements across the periodic table, in search of an emission source to backlight the K-edge measurement. Spectral flatness and high emission levels in the region of 1520-1600 eV around the K-edge of aluminum are the criteria for a suitable backlighter. Lead, gold and yttrium seem promising as backlighter sources, and the emission lines observed in silver provide a convenient spectral calibration source. |
Extent | 6390484 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-02-10 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0085443 |
URI | http://hdl.handle.net/2429/4367 |
Degree |
Master of Science - MSc |
Program |
Physics |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 1995-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_1995-0539.pdf [ 6.09MB ]
- Metadata
- JSON: 831-1.0085443.json
- JSON-LD: 831-1.0085443-ld.json
- RDF/XML (Pretty): 831-1.0085443-rdf.xml
- RDF/JSON: 831-1.0085443-rdf.json
- Turtle: 831-1.0085443-turtle.txt
- N-Triples: 831-1.0085443-rdf-ntriples.txt
- Original Record: 831-1.0085443-source.json
- Full Text
- 831-1.0085443-fulltext.txt
- Citation
- 831-1.0085443.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0085443/manifest