Time-dependent R -matrix theory for ultrafast atomic processes

We describe an ab initio nonperturbative time-dependent R-matrix theory for ultrafast atomic processes. This theory enables investigations of the interaction of few-femtosecond and -attosecond pulse lasers with complex multielectron atoms and atomic ions. A derivation and analysis of the basic equations are given, which propagate the atomic wave function in the presence of the laser field forward in time in the internal and external R-matrix regions. To verify the accuracy of the approach, we investigate two-photon ionization of Ne irradiated by an intense laser pulse and compare current results with those obtained using the R-matrix Floquet method and an alternative time-dependent method. We also verify the capability of the current approach by applying it to the study of two-dimensional momentum distributions of electrons ejected from Ne due to irradiation by a sequence of 2 as light pulses in the presence of a 780 nm laser field.

Coherent and stochastic contributions of compound resonances in atomic processes: Electron recombination, photoionization, and scattering

In open-shell atoms and ions, processes such as photoionization, combination (Raman) scattering, electron scattering, and recombination are often mediated by many-electron compound resonances. We show that their interference (neglected in the independent-resonance approximation) leads to a coherent contribution, which determines the energy-averaged total cross sections of electron- and photon-induced reactions obtained using the optical theorem. In contrast, the partial cross sections (e.g., electron recombination or photon Raman scattering) are dominated by the stochastic contributions. Thus, the optical theorem provides a link between the stochastic and coherent contributions of the compound resonances. Similar conclusions are valid for reactions via compound states in molecules and nuclei.

Many-electron relativistic calculation and interpretation of atomic processes in time dependent heavy-ion scattering

The time dependence of a heavy-ion-atom collision system is solved via a set of coupled channel equations using energy eigenvalues and matrix elements from a self-consistent field relativistic molecular many-electron Dirac-Fock-Slater calculation. Within this independent particle model we give a full many-particle interpretation by performing a small number of single-particle calculations. First results for the P(b) curves for the Ne K-hole excitation for the systems F{^8+} - Ne and F{^6+} - Ne as examples are discussed.

Entanglement analysis of atomic processes and quantum registers

During recent years, quantum information processing and the study of N−qubit quantum systems have attracted a lot of interest, both in theory and experiment. Apart from the promise of performing efficient quantum information protocols, such as quantum key distribution, teleportation or quantum computation, however, these investigations also revealed a great deal of difficulties which still need to be resolved in practise. Quantum information protocols rely on the application of unitary and non–unitary quantum operations that act on a given set of quantum mechanical two-state systems (qubits) to form (entangled) states, in which the information is encoded. The overall system of qubits is often referred to as a quantum register. Today the entanglement in a quantum register is known as the key resource for many protocols of quantum computation and quantum information theory. However, despite the successful demonstration of several protocols, such as teleportation or quantum key distribution, there are still many open questions of how entanglement affects the efficiency of quantum algorithms or how it can be protected against noisy environments. To facilitate the simulation of such N−qubit quantum systems and the analysis of their entanglement properties, we have developed the Feynman program. The program package provides all necessary tools in order to define and to deal with quantum registers, quantum gates and quantum operations. Using an interactive and easily extendible design within the framework of the computer algebra system Maple, the Feynman program is a powerful toolbox not only for teaching the basic and more advanced concepts of quantum information but also for studying their physical realization in the future. To this end, the Feynman program implements a selection of algebraic separability criteria for bipartite and multipartite mixed states as well as the most frequently used entanglement measures from the literature. Additionally, the program supports the work with quantum operations and their associated (Jamiolkowski) dual states. Based on the implementation of several popular decoherence models, we provide tools especially for the quantitative analysis of quantum operations. As an application of the developed tools we further present two case studies in which the entanglement of two atomic processes is investigated. In particular, we have studied the change of the electron-ion spin entanglement in atomic photoionization and the photon-photon polarization entanglement in the two-photon decay of hydrogen. The results show that both processes are, in principle, suitable for the creation and control of entanglement. Apart from process-specific parameters like initial atom polarization, it is mainly the process geometry which offers a simple and effective instrument to adjust the final state entanglement. Finally, for the case of the two-photon decay of hydrogenlike systems, we study the difference between nonlocal quantum correlations, as given by the violation of the Bell inequality and the concurrence as a true entanglement measure.

Atomic scale simulation of deformation in magnesium single crystals

The deformation behaviour of magnesium single crystals under plane strain conditions has been examined using molecular dynamics modelling. The simulations were based on an existing atomic potential for magnesium taken from the literature. A strain of 10% was applied at rates of 3x10⁹s^-1 and 3x10⁷s^-1. The simulations predicted the formation of mechanical twins that accommodated extension in the c-axis direction of the hexagonal unit cell. However, the predicted twin is not of the same kind found in magnesium, but is that commonly observed in titanium. It is believed that further analysis of the physical properties predicted by this interatomic potential will shed more light on the atomic processes controlling twinning in Magnesium alloys. It also highlights the need for improvements to the interatomic potential such that more accurate deformation behaviour can be attained.

An investigation of Fe XV emission lines in solar flare spectra

Previously, large discrepancies have been found between theory and observation for Fe XV emission line ratios in solar flare spectra covering the 224-327 angstrom wavelength range, obtained by the Naval Research Laboratory's S082A instrument on board Skylab. These discrepancies have been attributed to either errors in the adopted atomic data or the presence of additional atomic processes not included in the modelling, such as fluorescence. However our analysis of these plus other S082A flare observations (the latter containing Fe XV transitions between 321-482 angstrom), performed using the most recent Fe XV atomic physics calculations in conjunction with a chianti synthetic flare spectrum, indicate that blending of the lines is primarily responsible for the discrepancies. As a result, most Fe XV lines cannot be employed as electron density diagnostics for solar flares, at least at the spectral resolution of S082A and similar instruments (i.e.similar to 0.1 angstrom). An exception is the intensity ratio I(3s3p P-3(2)-3p(2) P-3(1))/I(3s3p P-3(2)-3p(2) D-1(2))=I(321.8 angstrom)/I(327.0 angstrom), which appears to provide good estimates of the electron density at this spectral resolution.

Numerical simulation of the charge balance in an EBIT

The electron beam ions traps (EBITs) are widely used to study highly charged ions (HCIs). In an EBIT, a high energy electron beam collides with atoms and ions to generate HCIs in the trap region. It is important to study the physics in the trap. The atomic processes, such as electron impact ionisation (EI), radiative recombination (RR), dielectronic recombination (DR) and charge exchange (CX), occur in the trap and numerical simulation can give some parameters for design, predict the composition and describe charge state evolution in an EBIT [Phys. Rev. A 43 (199 1) 4861]. We are presently developing a new code, which additionally includes a description of the overlaps between the ion clouds of the various charge-states. It has been written so that it can simulate experiments where various machine parameters (e.g. beam energy and current) can vary throughout the simulation and will be able to use cross- sections either based on scaling laws or derived from atomic structure calculations. An object-oriented method is used in developing the new software, which is an efficient way to organize and write code. (C) 2003 Elsevier Science B.V. All rights reserved.

Electron-impact excitation of Cr II A theoretical calculation of collision and effective collision strengths for forbidden transitions

Context. Absorption or emission lines of Cr II are observed in a wide variety of astrophysical spectra and accurate atomic data are urgently needed to interpret these lines. Many of these data are impossible to measure experimentally and a full theoretical treatment is the only means by which these data can be obtained.

Aims. In this paper, we present collision strengths and effective collision strengths for electron-impact excitation of Cr II for forbidden transitions among the lowest-lying 74 fine-structure levels. Effective collision strengths have been computed for 18 individual electron temperatures of astrophysical importance, ranging from 2000-100 000 K.

Methods. The parallel suite of R-matrix packages, RMATRX II, which has recently been extended to allow for the inclusion of relativistic effects, has been used in the present work to compute the collision strengths and effective collision strengths for electron-impact excitation of Cr II. We concentrate in this publication on low-lying forbidden lines among the lowest 74 jj fine-structure levels with configurations 3d(5) and 3d(4)4s, although atomic data has been evaluated for all 39 060 transitions among the 280 jj levels of configurations 3d(5), 3d(4)4s and 3d(4)4p. This work constitutes the largest evaluation ever performed for this ion involving 1932 coupled channels.

Results. Collision and effective collision strengths are presented for all transitions among the lowest 74 J pi states of Cr II and comparisons made with the work of Bautista et al. (2009). While the effective collision strengths agree well for some transitions, significant discrepancies exist for others. We believe that the present atomic data represents the most accurate, most sophisticated and most complete data set for electron-impact excitation of Cr II and we would recommend them to astrophysicists and plasma physicists in their application work. We would expect that the effective collision strengths presented for the important low-lying forbidden lines are accurate to within 15%.

Electron-impact excitation of Sc II: collision strengths and effective collision strengths for fine-structure transitions.

Accurate fine-structure atomic data for the Fe-peak elements are essential for interpreting astronomical spectra. There is a severe paucity of data available for Sc II, highlighted by the fact that no collision strengths are readily available for this ion. We present electron-impact excitation collision strengths and Maxwellian averaged effective collision strengths for Sc II. The collision strengths were calculated for all 3916 transitions amongst 89 jj levels (arising from the 3d4s, 3d2, 4s2, 3d4p, 4s4p, 3d5s, 3d4d, 3d5p, 4p2 and 3d4f configurations), resulting in a 944 coupled channel problem. The R-matrix package RMATRXII was utilized, along with the transformation code FINE and the external region code PSTGF, to calculate the collision strengths for a range of incident electron energies in the 0 to 8.3 Rydberg region. Maxwellian averaged effective collision strengths were then produced for 27 temperatures lying within the astrophysically significant range of 30 to 105 K.
The collision strengths and effective collision strengths were produced for two different target models. The purpose was to systematically examine the effect of including open 3p correlation terms into the configuration interaction expansion for the wavefunction. The first model consisted of all 36 CI terms that could be generated with the 3p core closed. The second model incorporated an additional six configurations which allowed for single-electron excitations from within the 3p core. Comparisons are made between the two models and the results of Bautista et al., obtained by private communication. It is concluded that the first model produced the most reliable set of collision and effective collision strengths for use in astrophysical and plasma applications.

nero- a post-maximum supernova radiation transport code

The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a non-local thermodynamic equilibrium (NLTE) treatment of radiation transport has rarely been studied. In this paper, we present a new SN radiation transport code, nero, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare nero to the radiation transport code of Jerkstrand, Fransson & Kozma and to the nebular code of Mazzali et al. All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, nero is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' Type Ia supernova (SN Ia) 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post-maximum epochs. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

Comparing Optical Thomson Scattering With Emission Spectroscopy In He Plasmas Induced By Laser Breakdown At 1 Atmosphere

We have carried out optical Thomson scattering measurements from a laser induced breakdown in He at 1 atmosphere. The breakdown was created with a Nd:YAG laser with 9ns pulse duration and 400mJ pulse energy focused into a chamber filled with He. A second harmonic Nd: YAG laser with 9ns pulses and up to 80mJ energy was used to obtain temporally and spatially resolved data on the electron density and temperature. In parallel experiments, we measured the emission of the 447.1nm line from He I. Initial results suggest good agreement between densities inferred but full Abel inversion is needed for conclusive results.

Electron impact excitation of Mg VIII: Collision strengths, transition probabilities and theoretical EUV and soft X-ray line intensities for Mg VIII

Context: Mg VIII emission lines are observed in a range of astronomical objects such as the Sun, other cool stars and in the coronal line region of Seyfert galaxies. Under coronal conditions Mg VIII emits strongly in the extreme ultraviolet (EUV) and soft X-ray spectral regions which makes it an ideal ion for plasma diagnostics.

Aims. Two theoretical atomic models, consisting of 125 fine structure levels, are developed for the Mg VIII ion. The 125 levels arise from the 2s(2)2p, 2s(2)p2, 2p(3), 2s(2)3s, 2s(2)3p, 2s(2)3d, 2s2p3s, 2s2p3p, 2s2p3d, 2p(2)3s, 2p(2)3p and 2p(2)3d configurations. Electron impact excitation collision strengths and radiative transition probabilities are calculated for both Mg VIII models, compared with existing data, and the best model selected to generate a set of theoretical emission line intensities. The EUV lines, covering 312-790 angstrom, are compared with existing solar spectra (SERTS-89 and SUMER), while the soft X-ray transitions (69-97 angstrom) are examined for potential density diagnostic line ratios and also compared with the limited available solar and stellar observational data.

Methods. The R-matrix codes Breit-Pauli RMATRXI and RMATRXII are utilised, along with the PSTGF code, to calculate the collision strengths for two Mg VIII models. Collision strengths are averaged over a Maxwellian distribution to produce the corresponding effective collision strengths for use in astrophysical applications. Transition probabilities are also calculated using the CIV3 atomic structure code. The best data are then incorporated into the modelling code CLOUDY and line intensities generated for a range of electron temperatures and densities appropriate to solar and stellar coronal plasmas.

Results. The present effective collision strengths are compared with two previous calculations. Good levels of agreement are found with the most recent, but there are large differences with the other for forbidden transitions. The resulting line intensities compare favourably with the observed values from the SERTS-89 and SUMER spectra. Theoretical soft X-ray emission lines are presented and several density diagnostic line ratios examined, which are in reasonable agreement with the limited observational data available.

Electron-impact excitation of FeII - Collision strengths and effective collision strengths for low-lying fine-structure forbidden transitions

Context. Absorption or emission lines of Fe II are observed in many astrophysical spectra and accurate atomic data are required to interpret these lines. The calculation of electron-impact excitation rates for transitions among even the lowest lying levels of Fe II is a formidable task for theoreticians.

Aims. In this paper, we present collision strengths and effective collision strengths for electron-impact excitation of Fe II for low-lying forbidden transitions among the lowest 16 fine-structure levels arising from the four LS states 3d(6)4s D-6(e), 3d(7) F-4(e), 3d(6)4s D-4(e), and 3d(7) P-4(e). The effective collision strengths are calculated for a wide range of electron temperatures of astrophysical importance from 30-100 000 K.

Methods. The parallel suite of Breit-Pauli codes are utilised to compute the collision cross sections for electron-impact excitation of Fe II and relativistic terms are included explicitly in both the target and the scattering approximation. 100 LS or 262-jj levels formed from the basis configurations 3d(6)4s, 3d(7), and 3d(6)4p were included in the wavefunction representation of the target, including all doublet, quartet, and sextet terms. Collision strengths for a total of 34191 individual transitions were computed.

Results. A detailed comparison is made with previous theoretical works and significant differences were found to occur in the effective collision strengths, particularly at low temperatures.