Fine-structure photoionization cross sections of Fe II

We present the first calculation of fine-structure photoionization cross sections for the ground state of singly ionized Fe. These large-scale ab initio calculations, limited to the near-threshold region, were performed in the close-coupling approximation using a Dirac–Coulomb R -matrix method implemented within a modified version of the DARC package. Our calculated cross sections reproduce in detail the resonance structures observed in previous experimental determinations.

Photoionization cross-sections for Fe II

The LS R-matrix method was used to compute new photoionization cross sections for Fe II. Results are compared with available experimental data and with previous calculations of the cross section. We also present the first fine-structure photoionization data for this ion obtained with the fully-relativistic DARC codes.

Electron-impact excitation and ionization of W 3+ for the determination of tungsten influx in a fusion plasma

Tungsten will be employed as a plasma facing material in the ITER fusion reactor under construction in Cadarache, France; therefore, there is a significant need for accurate electron-impact excitation and ionization data for the ions of tungsten. We report on the results of extensive calculations of ionization and excitation for W 3+ that are intended to provide the atomic data needed for the determination of impurity influx diagnostics of tungsten in several existing tokamak reactors. The electron-impact excitation rate coefficients for this study were determined using the relativistic R -matrix method. The contribution to direct electron-impact ionization was determined using the distorted-wave approximation, the accuracy of which was verified by an R -matrix with pseudo states calculation. Contributions to total ionization from excitation autoionization were also generated from the relativistic R -matrix method. These results were then employed to calculate values of ionization per emitted photon, or SXB ratios, for four carefully selected spectral lines; these data will allow the determination of impurity influx from tungsten facing surfaces. For the range of densities of importance in the edge region of a tokamak reactor, these SXB ratios are found to be nearly independent of electron density but vary significantly with electron temperature.

Collisional-radiative calculations of He line emission in low-temperature plasmas

We present spectral modeling results for neutral helium. Our underlying atomic data contains radiative transition rates that are generated from atomic structure calculations and electron-impact excitation rates, that are determined from both the standard R-matrix method and the R-matrix with pseudostates RMPS method. In this paper, we focus on transitions of particular importance to diagnostic line ratios. For example, our calculated rate coefficient for the electron-impact transition 1s3s 1S→1s3p 1P, which has a pronounced effect on the 728.1 nm diagnostic spectral line, is found to be in good agreement with previous experimental mea- surements. We also consider transitions from the 1s2 1S ground and 1s2s 3S terms to terms of the n=4 shell. They are found to be affected significantly by coupling of the bound states to the target continuum continuum coupling, which is included in our RMPS calculation, but not in our standard R-matrix calculation. We perform collisional-radiative calculations to determine spectral line intensity ratios for three ratios of particular interest, namely the 504.8 nm/471.3 nm, 492.2 nm/471.3 nm, and 492.2 nm/504.8 nm line ratios. Comparing our results determined from the RMPS excitation rates with those from the standard R-matrix excitation rates, we find that continuum coupling affects the rate coefficients significantly, leading to different values for all three line ratios. We also compare our modeling results with spectral measurements taken recently on the Auburn Helicon plasma device, finding that the ground and metastable populations are not in equilibrium, and that the experimental measurements are more consistent with the 1s2s 3S metastable term populations being short lived in the plasma.

The effects of radiative cascades on the x-ray diagnostic lines of Fe16+

We present complete collisional-radiative modelling results for the soft x-ray emission lines of Fe16+ in the 15 Å–17 Å range. These lines have been the subject of much controversy in the astrophysical and laboratory plasma community. Radiative transition rates are generated from fully relativistic atomic structure calculations. Electron-impact excitation cross sections are determined using a fully relativistic R-matrix method employing 139 coupled atomic levels through n = 5. We find that, in all cases, using a simple ratio of the collisional rate coefficient times a radiative branching factor is not sufficient to model the widely used diagnostic line ratios. One has to include the effects of collisional-radiative cascades in a population model to achieve accurate line ratios. Our line ratio results agree well with several previous calculations and reasonably well with tokamak experimental measurements, assuming a Maxwellian electron-energy distribution. Our modelling results for four EBIT line ratios, assuming a narrow Gaussian electron-energy distribution, are in generally poor agreement with all four NIST measurements but are in better agreement with the two LLNL measurements. These results suggest the need for an investigation of the theoretical polarization calculations that are required to interpret the EBIT line ratio measurements.

Effective collision strengths for forbidden transitions among the 3s23p3 fine-structure levels of Cl III

Effective collision strengths for the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 3s²3p³ configuration of Cl III are presented. The calculation employs the multichannel R-matrix method to compute the electron-impact excitation collision strengths in a close-coupling expansion, which incorporates the lowest 23 LS target eigenstates of Cl III. These states are formed from the 3s²3p³, 3s3p⁴, 3s²3p²3d and 3s²3p²4s configurations. The Maxwellian-averaged effective collision strengths are presented graphically for all 10 fine-structure transitions over a wide range of electron temperatures appropriate for astrophysical applications [log T(K) = 3.3 - log T(K) = 5.9]. Comparisons are made with the earlier seven-state close-coupling calculation of Butler & Zeippen, and in general excellent agreement is found in the low-temperature region where a comparison is possible [log T(K) = 3.3 - log T(K) = 4.7]. However, discrepancies of up to 30 per cent are found to occur for the forbidden transitions which involve the ⁴S^o ground state level, particularly for the lowest temperatures considered. At the higher temperatures, the present data are the only reliable results currently available.

Effective collision strengths for fine-structure forbidden transitions among the 2s22p3 levels of S x

Effective collision strengths for electron-impact excitation of the N-like ion S x are calculated in the close-coupling approximation using the multichannel R-matrix method. Specific attention is given to the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^O, ²D^o and ²P^o levels in the 2s²2p³ ground configuration. The total (e^- + ion) wavefunction is expanded in terms of the 11 lowest LS eigenstates of S x, and each eigenstate is represented by extensive configuration-interaction wavefunctions. The collision strengths obtained are thermally averaged over a Maxwellian distribution of velocities, for all 10 fine-structure transitions, over the range of electron temperatures log T(K) = 4.6-6.7 (the range appropriate for astrophysical applications). The present effective collision strengths are the only results currently available for these fine-structure transition rates.

Effective collision strengths for fine-structure forbidden transitions among the 2s22p3 levels of Ne IV

Effective collision strengths for electron-impact excitation of the N-like ion NeIV are calculated in the close-coupling approximation using the multichannel R-matrix method. Specific attention is given to the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 2s²2p³ ground-state configuration. The expansion of the total wavefunction incorporates the lowest 11 LS eigenstates of NeIV, consisting of eight n = 2 terms with configurations 2s²2p³, 2s2p⁴ and 2p⁵, together with three n = 3 states of configuration 2s²2p²3s. We present in graphical form the effective collision strengths obtained by thermally averaging the collision strengths over a Maxwellian distribution of velocities, for all 10 fine-structure transitions, over the range of electron temperatures log T(K) = 3.6 to log T(K) = 6.1 (the range appropriate for astrophysical applications). Comparisons are made with the earlier, less sophisticated close-coupling calculation of Giles, and excellent agreement is found in the limited temperature region where a comparison is possible [log T(K) = 3.7 to log 7(K) = 4.3]. At higher temperatures the present data are the only reliable results currently available.

Effective collision strengths for fine-structure forbidden transitions among the 3s 23p 3 levels of Ar IV

The multichannel R-matrix method is used to compute electron impact excitation collision strengths in Ar IV for all fine-structure transitions among the ⁴S°, ²D° and ²P° levels in the 3s ²3p ³ ground configuration. Included in the expansion of the total wavefunction are the lowest 13 LS target eigenstates of Ar iv formed from the 3s ²3p ³, 3s3p ⁴ and 3s ²3p ²3d configurations. The effective collision strengths, obtained by averaging the electron collision strengths over a Maxwellian distribution of electron velocities, are presented for all 10 fine-structure transitions over a wide range of electron temperatures of astrophysical interest (T _e = 2000-100 000 K). Comparisons are made with an earlier 7-state close-coupling calculation by Zeippen, Butler & Le Bourlot, and significant differences are found to occur for many of the forbidden transitions considered, in particular those involving the ⁴S° ground state, where discrepancies of up to a factor of 3 are found in the low-temperature region. © 1997 RAS.

Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

R-matrix-Floquet theory of multiphoton processes:VIII. A linear equations method

A new linear equations method for calculating the R-matrix, which arises in the R-matrix-Floquet theory of multiphoton processes, is introduced. This method replaces the diagonalization of the Floquet Hamiltonian matrix by the solution of a set of linear simultaneous equations which are solved, in the present work, by the conjugate gradient method. This approach uses considerably less computer memory and can be readily ported onto parallel computers. It will thus enable much larger problems of current interest to be treated. This new method is tested by applying it to three-photon ionization of helium at frequencies where double resonances with a bound state and autoionizing states are important. Finally, an alternative linear equations method, which avoids the explicit calculation of the R-matrix by incorporating the boundary conditions directly, is described in an appendix.

Stiffness analysis of polymeric composites using the finite element method

The ability to predict the mechanical behavior of polymer composites is crucial for their design and manufacture. Extensive studies based on both macro- and micromechanical analyses are used to develop new insights into the behavior of composites. In this respect, finite element modeling has proved to be a particularly powerful tool. In this article, we present a Galerkin scheme in conjunction with the penalty method for elasticity analyses of different types of polymer composites. In this scheme, the application of Green's theorem to the model equation results in the appearance of interfacial flux terms along the boundary between the filler and polymer matrix. It is shown that for some types of composites these terms significantly affect the stress transfer between polymer and fillers. Thus, inclusion of these terms in the working equations of the scheme preserves the accuracy of the model predictions. The model is used to predict the most important bulk property of different types of composites. Composites filled with rigid or soft particles, and composites reinforced with short or continuous fibers are investigated. For each case, the results are compared with the available experimental results and data obtained from other models reported in the literature. Effects of assumptions made in the development of the model and the selection of the prescribed boundary conditions are discussed.