Energy levels, radiative rates and electron impact excitation rates for transitions in C III

We report energy levels, radiative rates (A-values) and lifetimes for the astrophysically important Be-like ion C III. For the calculations, 166 levels belonging to the n ≤ 5 configurations are considered and the GRASP (General-purpose Relativistic Atomic Structure Package) is adopted. Einstein A-coefficients are provided for all E1, E2, M1 and M2 transitions, while lifetimes are compared with available measurements as well as theoretical results, and no large discrepancies noted. Our energy levels are assessed to be accurate to better than 1 per cent for a majority of levels, and A-values to better than 20 per cent for most transitions. Collision strengths are also calculated, for which the Dirac Atomic R-matrix Code (DARC) is used. A wide energy range, up to 21 Ryd, is considered and resonances resolved in a fine energy mesh in the thresholds region. The collision strengths are subsequently averaged over a Maxwellian velocity distribution to determine effective collision strengths up to a temperature of 8.0 ×10[5]K, sufficient for most astrophysical applications. Our data are compared with the recent R-matrix calculations of Fernández-Menchero et al., and significant differences (up to over an order of magnitude) are noted for several transitions over the complete temperature range of the results.

Energy levels, radiative rates and electron impact excitation rates for transitions in Fe XIV

Energies and lifetimes are reported for the lowest 136 levels of Fe XIV, belonging to the (1s(2)2s(2)2p(6)) 3s(2)3p, 3s(3)p(2), 3s(2)3d, 3p(3), 3s(3)p(3)d, 3p(2)3d, 3s(3)d(2), 3p(3)d(2) and 3s(2)4l configurations. Additionally, radiative rates, oscillator strengths and line strengths are calculated for all electric dipole (E1), magnetic dipole (M1), electric quadrupole (E2) and magnetic quadrupole (M2) transitions. Theoretical lifetimes determined from these radiative rates for most levels show satisfactory agreement with earlier calculations, a swell as with measurements. Electron impact excitation collision strengths are also calculated with the Dirac atomic R-matrix code (DARC) over a wide energy range up to 260 Ryd. Furthermore, resonances have been resolved in a fine energy mesh to determine effective collision strengths, obtained after integrating the collision strengths over a Maxwellian distribution of electron velocities. Results are listed for all 9180 transitions among the 136 levels over a wide range of electron temperatures, up to 10(7.1) K. Comparisons are made with available results in the literature, and the accuracy of the present data is assessed.

Electron-phonon thermalization in a scalable method for real-time quantum dynamics

We present a quantum simulation method that follows the dynamics of out-of-equilibrium many-body systems of electrons and oscillators in real time. Its cost is linear in the number of oscillators and it can probe time scales from attoseconds to hundreds of picoseconds. Contrary to Ehrenfest dynamics, it can thermalize starting from a variety of initial conditions, including electronic population inversion. While an electronic temperature can be defined in terms of a nonequilibrium entropy, a Fermi-Dirac distribution in general emerges only after thermalization. These results can be used to construct a kinetic model of electron-phonon equilibration based on the explicit quantum dynamics.

Intermediate-coupling R -matrix calculations of electron-impact excitation of Fe 5+

For applications to laboratory and astrophysical plasmas, there is a great need for accurate electron-impact excitation data between individual levels in the lower charge-state ions of iron. Recently, we have reported on the first intermediate-coupling R -matrix calculation of electron-impact excitation in Fe 4+ , in which the close-coupling expansion of the target included levels from both ground and excited configurations (Ballance et al 2007 J. Phys. B: At. Mol. Opt. Phys. [/0953-4075/40/23/f01] 40 F327 , 2008 Europhys. News 39 14). In this paper, we present the results of two large intermediate-coupling Dirac R -matrix calculations of electron-impact excitation of Fe 5+ . The results from the two calculations, which differ only in the configuration–interaction expansions of the target, are compared. These comparisons provide some indication of the accuracy of the calculations and the resulting data should be useful in modelling plasmas containing iron.

Photoionization cross section calculations for the halogen-like ions Kr + and Xe +

Photoionization cross section calculations on the halogen-like ions; Kr + and Xe + have been performed for a photon energy range from each ion threshold to 15 eV, using large-scale close-coupling calculations within the Dirac--Coulomb R -matrix approximation. The results from our theoretical work are compared with recent measurements made at the ASTRID merged-beam set-up at the University of Aarhus in Denmark and from the Fourier transform ion cyclotron resonance trap method at the SOLEIL synchrotron radiation facility in Saint-Aubin, France Bizau et al (2011 J. Phys. B: At. Mol. Opt. Phys. 44 055205) and the advanced light source M{ü}ller (2012 private communication), Aguliar et al (2012 J. Phys.: Conf. Ser . at press). For each of these complex ions our theoretical cross section results over the photon energy range investigated are seen to be in excellent agreement with experiment. Resonance energy positions and quantum defects of the prominent Rydberg resonances series identified in the spectra are compared with experiment for these complex halogen-like ions.

Photoionization cross sections for the trans-iron element Se + from 18 to 31 eV

Absolute photoionization cross-section calculations are presented for Se + using large-scale close-coupling calculations within the Breit--Pauli and Dirac--Coulomb R -matrix approximations. The results from our theoretical work are compared with recent measurements (Esteves 2010 PhD Thesis publication number AAI3404727, University of Reno, NV, USA; Sterling et al 2011 J. Phys. B: At. Mol. Opt. Phys. 44 025701; Esteves et al 2011 Phys. Rev. A 84 013406) made at the advanced light source (ALS) radiation facility in Berkeley, CA, USA. We report on results for the photon energy range 18.0--31.0 eV, which spans the ionization thresholds of the 4 S o 3/2 ground state and the low-lying 2 D o 5/2,3/2 and 2 P o 3/2,1/2 metastable states. Metastable fractions are inferred from our present work. Resonance energies and quantum defects of the prominent Rydberg resonances series identified in the spectra are compared for the 4p → n d transitions with the recent ALS experimental measurements made on this complex trans-iron element.

PAMOP: Petascale atomic, molecular and optical collision calculations

Petaflop architectures are currently being utilized efficiently to perform large scale computations in Atomic, Molecular and Optical Collisions. We solve the Schr\"odinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach. We briefly outline the parallel methodology used and implemented for the current suite of Breit-Pauli and DARC codes. In this report, various examples are shown from our theoretical results compared with experimental results obtained from Synchrotron Radiation facilities where the Cray architecture at HLRS is playing an integral part in our computational projects.

Dielectronic recombination of W 35 +

Accurate data for dielectronic recombination (DR) of the ions of tungsten are of significant interest in the modelling of tungsten impurity transport and radiative power loss in current tokamaks and in ITER. However, the complexity of the atomic structure for many of these ions makes level-resolved DR calculations untenable on currently available computers, especially for open d- and f-subshell ions. The majority of DR data presently available for ITER modelling are based on an average-atom approximation. To improve upon these baseline calculations, we investigate the use of the configuration-average distorted-wave (CADW) method to calculate DR rate coefficients for complex open d-shell systems. The aim is to produce rate coefficients that are sufficiently accurate in terms of modelling, yet greatly reduced in term of computational complexity compared to level-resolved calculations. In this paper, we consider the DR of W 35 + . Initially, we carry out several large-scale level-resolved calculations for the DR associated with the 4d → 4f and 4p → 4d excitations in this ion, using both the level-resolved distorted-wave and Dirac R -matrix methods. These calculations allow us to test the validity of the CADW approach on these same excitations by comparing cross sections and rate coefficients. These comparisons demonstrate that the CADW method is relatively accurate in relation to these level-resolved methods for the temperature range for which W 35 + should exist in a collisionally ionized plasma. We then present results for CADW rate coefficients for both Δ n = 0 and Δ n = 1 excitations for this ion. This study indicates that it is now feasible to generate a much improved comprehensive set of DR data for the entire tungsten isonuclear sequence.

Electron-impact excitation of Xe 26+ and its resultant spectral signature

We have carried out a 129 close-coupling level Dirac-Coulomb R-matrix calculation for the electron-impact excitation of Ni-like Xe. We have utilized this data to generate the spectral signature of Xe26+ in terms of feature photon-emissivity coefficients (F-PεCs). We have compared these F-PεCs with those generated using semi-relativistic plane-wave Born excitation data, which forms the heavy species baseline for the Atomic Data and Analysis Structure (ADAS), We find that the Born-based F-PεCs give a reasonable qualitative description of the spectral signature but that, quantitatively, the R-matrix-based F-PεCs differ by up to a factor of 2. The spectral signature of heavy species is key to diagnosing hot plasmas such as will be found in the International Thermonuclear Experimental Reactor.

Experimental and theoretical photoionization cross sections for Se + from 18 eV to 31 eV

Absolute Se photoionization cross-section measurements and Dirac-Coulomb R -matrix calculations are reported for the photon energy range 18.0 eV – 31.0 eV, which spans the ionization thresholds of the 4 S 0 3/2 ground state and the low-lying 2 D 0 3/2,5/2 and 2 P 0 1/2,3/2 metastable states. The determination of the photoionization and recombination properties of n -capture element ions is motivated by their astrophysical detection and the importance of their elemental abundances in testing theories of nucleosynthesis and stellar structure.

High-resolution photoionization of Xe + ions: Experiment and theory

High-resolution photoionization measurements of Xe + ions have been performed at the Advanced Light Source in Berkeley, California, USA. The experimental cross sections are compared with results from Dirac-Coulomb R-matrix calculations.

Dynamic analysis of earthquake phenomena by means of pseudo phase plane

This paper analyses earthquake data in the perspective of dynamical systems and its Pseudo Phase Plane representation. The seismic data is collected from the Bulletin of the International Seismological Centre. The geological events are characterised by their magnitude and geographical location and described by means of time series of sequences of Dirac impulses. Fifty groups of data series are considered, according to the Flinn-Engdahl seismic regions of Earth. For each region, Pearson’s correlation coefficient is used to find the optimal time delay for reconstructing the Pseudo Phase Plane. The Pseudo Phase Plane plots are then analysed and characterised.

Analysis of Forest Fires by means of Pseudo Phase Plane and Multidimensional Scaling Methods

Forest fires dynamics is often characterized by the absence of a characteristic length-scale, long range correlations in space and time, and long memory, which are features also associated with fractional order systems. In this paper a public domain forest fires catalogue, containing information of events for Portugal, covering the period from 1980 up to 2012, is tackled. The events are modelled as time series of Dirac impulses with amplitude proportional to the burnt area. The time series are viewed as the system output and are interpreted as a manifestation of the system dynamics. In the first phase we use the pseudo phase plane (PPP) technique to describe forest fires dynamics. In the second phase we use multidimensional scaling (MDS) visualization tools. The PPP allows the representation of forest fires dynamics in two-dimensional space, by taking time series representative of the phenomena. The MDS approach generates maps where objects that are perceived to be similar to each other are placed on the map forming clusters. The results are analysed in order to extract relationships among the data and to better understand forest fires behaviour.

Dynamic Analysis and Pattern Visualization of Forest Fires

This paper analyses forest fires in the perspective of dynamical systems. Forest fires exhibit complex correlations in size, space and time, revealing features often present in complex systems, such as the absence of a characteristic length-scale, or the emergence of long range correlations and persistent memory. This study addresses a public domain forest fires catalogue, containing information of events for Portugal, during the period from 1980 up to 2012. The data is analysed in an annual basis, modelling the occurrences as sequences of Dirac impulses with amplitude proportional to the burnt area. First, we consider mutual information to correlate annual patterns. We use visualization trees, generated by hierarchical clustering algorithms, in order to compare and to extract relationships among the data. Second, we adopt the Multidimensional Scaling (MDS) visualization tool. MDS generates maps where each object corresponds to a point. Objects that are perceived to be similar to each other are placed on the map forming clusters. The results are analysed in order to extract relationships among the data and to identify forest fire patterns.

Dynamic Analysis and Pattern Visualization of Forest Fires

This paper analyses forest fires in the perspective of dynamical systems. Forest fires exhibit complex correlations in size, space and time, revealing features often present in complex systems, such as the absence of a characteristic length-scale, or the emergence of long range correlations and persistent memory. This study addresses a public domain forest fires catalogue, containing information of events for Portugal, during the period from 1980 up to 2012. The data is analysed in an annual basis, modelling the occurrences as sequences of Dirac impulses with amplitude proportional to the burnt area. First, we consider mutual information to correlate annual patterns. We use visualization trees, generated by hierarchical clustering algorithms, in order to compare and to extract relationships among the data. Second, we adopt the Multidimensional Scaling (MDS) visualization tool. MDS generates maps where each object corresponds to a point. Objects that are perceived to be similar to each other are placed on the map forming clusters. The results are analysed in order to extract relationships among the data and to identify forest fire patterns.