Atomic and Molecular Aspects of Astronomical Spectra Arxiv:1211.1236V1 [Physics.Atom-Ph] 3 Nov 2012
Total Page:16
File Type:pdf, Size:1020Kb
Atomic and Molecular Aspects of Astronomical Spectra Taha Sochi A dissertation submitted to the Department of Physics & Astronomy University College London in fulfilment of the requirements for the degree of Doctor of Philosophy arXiv:1211.1236v1 [physics.atom-ph] 3 Nov 2012 September 2012 Declaration The work presented in this thesis is the result of my own effort, except where otherwise stated. All sources are acknowledged by explicit reference. Taha Sochi ......................... 1 Abstract In the first section of this thesis, we present the atomic part of our investigation. AC2+ atomic target was prepared and used to generate theoretical data required in the investigation of recombination lines that arise from collisions between elec- trons and ions in thin plasma found in planetary nebulae and other astrophysical objects. The R-matrix method of electron scattering theory was used to describe the C2+ plus electron system. Theoretical data concerning bound and autoioniz- ing states were generated in the intermediate-coupling approximation by R-matrix and Autostructure codes and compared to experimental data. The comparison revealed very good agreement. These theoretical data were then used to generate dielectronic recombination data for C+ which include transition lines, oscillator strengths, radiative transition probabilities, as well as emissivity and dielectronic recombination coefficients. The data were cast in the form of a line list, called SS1, containing 6187 optically-allowed transitions which include many C ii lines observed in astronomical spectra. The data were validated by comparison to C+ recombination data obtained from a number of sources in the literature. This line list was used to analyze the spectra from a number of astronomical objects, mainly planetary nebulae, and identify their electron temperature where the observational data were obtained from the literature. The electron temperature investigation was also extended to include free electron energy analysis which uses observational and theoretical data of FF and FB transitions to investigate the long-standing problem of discrepancy between the results of recombination and forbidden lines analysis and its possible connection to the electron distribution (Maxwellian or non-Maxwellian). In the course of this investigation two elaborate methods, one 2 for finding and analyzing resonances (K-matrix method) and the other for analyz- ing and identifying electron temperature from astronomical spectra (least squares minimization), were employed. A computer program for atomic transition analy- sis was also developed and used as the main tool for obtaining the line list and analyzing the observational spectra. In the second section of the thesis we present the results of our molecular in- vestigation; the generation of a comprehensive, calculated line list of frequencies + + and transition probabilities for the singly deuterated isotopologue of H3 ,H2D . The line list, which is the most comprehensive one of its kind in existence, contains over 22 million rotational-vibrational transitions occurring between more than 33 thousand energy levels and covers frequencies up to 18500 cm−1. All energy levels with rotational quantum number, J, up to 20 are considered, making the line list useful for temperatures up to at least 3000 K. About 15% of these levels are fully assigned with approximate rotational and vibrational quantum numbers. The list is calculated using a previously proposed, highly accurate, ab initio model imple- mented in a high-accuracy computer code based on a two-stage discrete variable representation (DVR) approach. Various consistency checks were carried out to test and validate the results. All these checks confirmed the accuracy of the list. A temperature-dependent partition function, valid over a more extended temper- ature range than those previously published, and cooling function are presented. Temperature-dependent synthetic spectra for the temperatures T = 100; 500; 1000 and 2000 K in the frequency range 0-10000 cm−1 were also generated and presented graphically. 3 Acknowledgements and Statements I would like to thank my supervisors Professor Peter Storey and Professor Jonathan Tennyson for their kindness, guidance and advice. It is a great privi- lege and gratifying experience to work with them over several years and learn from their deep knowledge and vast expertise in atomic, molecular and astro physics. I also would like to thank the external examiner Professor Nigel Badnell from the University of Strathclyde, and the internal examiner Dr Sergey Yurchenko from the University College London for their constructive remarks and corrections and valuable recommendations for possible future work. Apart from RMATRX-I and Autostructure atomic codes and DVR3D molecular code which were used as the main computing tools, I also used RESMESH and another piece of code from Pete to calculate f-values from photoionization cross sections. The rest of tools and techniques belong to me. Anyway, I assume full responsibility for any error or failure. The emphasis in this thesis is on the practical aspects of my work in an attempt to provide a full and credible account of how this investigation was conducted and how the conclusions were reached. Bare minimum of theory will be given when necessary for motivation and background as the theory is readily available in the literature and mainly belongs to other contributors. Finally, I should acknowledge the encouragement and motivation which I re- ceived from family and friends while doing this PhD with special tribute to my wife. 4 The following articles are based on results and activities related directly to the research presented in this thesis: • Taha Sochi. Emissivity: A program for atomic transition calculations. Com- munications in Computational Physics, 7(5):1118-1130, 2010. + • Taha Sochi; Jonathan Tennyson. A computed line list for the H2D molecular ion. Monthly Notices of the Royal Astronomical Society, 405(4):2345-2350, 2010. • Taha Sochi; Peter J. Storey. Dielectronic Recombination Lines of C+. Atomic Data and Nuclear Data Tables (Accepted), 2012. • Peter J. Storey; Taha Sochi. Electron Temperature of Planetary Nebulae from C ii Recombination Lines. Submitted, 2012. 5 Contents Declaration . .1 Abstract . .2 Acknowledgements and Statements . .4 Contents . .6 List of Figures . 10 List of Tables . 13 Nomenclature . 18 Part I: Atomic Physics 24 1 Introduction 25 1.1 Recombination Processes . 25 1.1.1 Radiative Recombination . 26 1.1.2 Dielectronic Recombination . 27 1.2 Literature Review . 29 1.2.1 Theory of Recombination Lines . 29 1.2.2 Recombination Lines of Carbon . 33 1.3 Recombination versus Forbidden Lines Analysis . 38 2 Physics of Scattering and Recombination 41 2.0.1 Close Coupling Approximation . 42 2.0.2 R-matrix Method . 45 2.0.2.1 Inner and Outer Regions . 47 2.0.2.2 Matching of Solutions . 48 2.0.2.3 Bound States . 48 2.0.2.4 Resonances . 48 2.0.2.5 Oscillator Strengths . 49 2.0.2.6 Photoionization Cross Sections . 50 2.0.2.7 Relativistic Effects . 50 6 2.0.3 RMATRX-I Code . 51 3 Dielectronic Recombination Lines of C ii 55 3.1 Theoretical Background . 56 3.2 Autoionizing and Bound States Calculations . 62 3.2.1 Methods for Investigating Resonances . 63 3.2.1.1 QB Method . 64 3.2.1.2 Time-Delay Method . 66 3.2.1.3 K-Matrix Method . 67 3.2.2 Preparing Target . 68 3.2.3 Developing and Implementing K-matrix Method . 69 3.2.4 Calculations . 72 3.3 f-Values Calculations . 75 3.3.1 f-Values for Free-Free Transitions . 76 3.3.2 f-Values for Free-Bound Transitions . 77 3.3.3 f-Values for Bound-Bound Transitions . 82 3.4 Recombination and Emissivity Calculations . 82 3.4.1 Practical Aspects . 83 3.5 Comparison to Previous Work . 84 4 Electron Temperature of Astronomical Objects 89 4.1 Method . 91 4.1.1 Least Squares Minimization . 95 4.2 Astronomical Objects . 98 4.2.1 NGC 7009 . 99 4.2.2 NGC 5315 . 105 4.2.3 NGC 7027 . 109 4.2.3.1 Baluteau et al ..................... 109 4.2.3.2 Zhang et al ...................... 110 4.2.4 IC 418 . 115 4.2.5 NGC 2867 . 120 4.2.6 DQ Herculis 1934 . 123 4.2.7 CPD - 56◦8032 . 126 4.2.8 He 2-113 . 129 4.3 Electron Distribution . 132 4.4 Epilogue . 136 7 Part II: Molecular Physics 140 + 5 Computed Line List for H2D 141 5.1 Introduction . 141 5.1.1 Literature Review . 144 5.1.2 Line Lists . 146 5.2 Method . 147 5.3 Validation and Primary Results . 154 5.3.1 Comparison to Experimental and Theoretical Data . 154 5.3.2 Partition Function . 155 5.3.3 Cooling Function . 158 5.3.4 Synthetic Spectra . 161 6 Conclusions and Future Work 164 6.1 Conclusions . 164 6.2 Recommendations for Future Work . 167 Bibliography 169 A Using Lifetime Matrix to Investigate Single-Channel Resonances190 B Error Analysis 195 B.1 Observed Error . 196 B.2 Fixed Percentage Error . 197 B.3 Poisson Error . 197 C Emissivity Program Documentation 199 C.1 Theoretical background . 200 C.2 Input and Output Files . 202 C.2.1 `Input' Data File . 202 C.2.2 `ResEmis' Input Data File . 209 C.2.3 `ELEVEmis' Input Data File . 210 C.2.4 `Transitions' Output Data File . 211 C.3 Tests . 211 D Input Data for R-matrix and Autostructure 213 D.1 R-matrix Input Data . 213 D.1.1 STG1 . 213 D.1.2 STG2 . 214 8 D.1.3 STGJK . 215 D.1.4 STG3 . 216 D.1.5 STGF . 216 D.1.6 STGB . 217 D.1.7 STGBF . 217 D.1.8 STGBB . 217 D.1.9 STGQB . 217 D.2 Autostructure Input Data . 218 D.2.1 Polarizability . 218 D.2.2 Transition Probabilities . 218 D.2.3 f-values . 218 E Tables 220 Bound States .