Resonant States in Negative Ions

Resonant states are multiply excited states in atoms and ions that have enough energy to decay by emitting an electron. The ability to emit an electron and the strong electron correlation (which is extra strong in negative ions) makes these states both interesting and challenging from a theoretical point of view. The main contribution in this thesis is a method, which combines the use of B splines and complex rotation, to solve the three-electron Schrödinger equation treating all three electrons equally. It is used to calculate doubly excited and triply excited states of 4S symmetry with even parity in He-. For the doubly excited states there are experimental and theoretical data to compare with. For the triply excited states there is only theoretical data available and only for one of the resonances. The agreement is in general good. For the triply excited state there is a significant and interesting difference in the width between our calculation and another method. A cause for this deviation is suggested. The method is also used to find a resonant state of 4S symmetry with odd parity in H2-. This state, in this extremely negative system, has been predicted by two earlier calculations but is highly controversial. Several other studies presented here focus on two-electron systems. In one, the effect of the splitting of the degenerate H(n=2) thresholds in H-, on the resonant states converging to this threshold, is studied. If a completely degenerate threshold is assumed an infinite series of states is expected to converge to the threshold. Here states of 1P symmetry and odd parity are examined, and it is found that the relativistic and radiative splitting of the threshold causes the series to end after only three resonant states. Since the independent particle model completely fails for doubly excited states, several schemes of alternative quantum numbers have been suggested. We investigate the so called DESB (Doubly Excited Symmetry Basis) quantum numbers in several calculations. For the doubly excited states of He- mentioned above we investigate one resonance and find that it cannot be assigned DESB quantum numbers unambiguously. We also investigate these quantum numbers for states of 1S even parity in He. We find two types of mixing of DESB states in the doubly excited states calculated. We also show that the amount of mixing of DESB quantum numbers can be inferred from the value of the cosine of the inter-electronic angle. In a study on Li- the calculated cosine values are used to identify doubly excited states measured in a photodetachment experiment. In particular a resonant state that violates a propensity rule is found.

Heating the Early Universe : Numerical Methods and Their Analysis

During the epoch when the first collapsed structures formed (6<z<50) our Universe went through an extended period of changes. Some of the radiation from the first stars and accreting black holes in those structures escaped and changed the state of the Intergalactic Medium (IGM). The era of this global phase change in which the state of the IGM was transformed from cold and neutral to warm and ionized, is called the Epoch of Reionization.In this thesis we focus on numerical methods to calculate the effects of this escaping radiation. We start by considering the performance of the cosmological radiative transfer code C2-Ray. We find that although this code efficiently and accurately solves for the changes in the ionized fractions, it can yield inaccurate results for the temperature changes. We introduce two new elements to improve the code. The first element, an adaptive time step algorithm, quickly determines an optimal time step by only considering the computational cells relevant for this determination. The second element, asynchronous evolution, allows different cells to evolve with different time steps. An important constituent of methods to calculate the effects of ionizing radiation is the transport of photons through the computational domain or ``ray-tracing''. We devise a novel ray tracing method called PYRAMID which uses a new geometry - the pyramidal geometry. This geometry shares properties with both the standard Cartesian and spherical geometries. This makes it on the one hand easy to use in conjunction with a Cartesian grid and on the other hand ideally suited to trace radiation from a radially emitting source. A time-dependent photoionization calculation not only requires tracing the path of photons but also solving the coupled set of photoionization and thermal equations. Several different solvers for these equations are in use in cosmological radiative transfer codes. We conduct a detailed and quantitative comparison of four different standard solvers in which we evaluate how their accuracy depends on the choice of the time step. This comparison shows that their performance can be characterized by two simple parameters and that the C2-Ray generally performs best.