An R-matrix approach for plasma modelling and the interpretation of astrophysical observation

Over the last decade an Auburn-Rollins-Strathclyde consortium has developed several suites of parallel R-matrix codes [1, 2, 3] that can meet the fundamental data needs required for the interpretation of astrophysical observation and/or plasma experiments. Traditionally our collisional work on light fusion-related atoms has been focused towards spectroscopy and impurity transport for magnetically confined fusion devices. Our approach has been to provide a comprehensive data set for the excitation/ionization for every ion stage of a particular element. As we progress towards a burning fusion plasma, there is a demand for the collisional processes involving tungsten, which has required a revitalization of the relativistic R-matrix approach. The implementation of these codes on massively parallel supercomputers has facilitated the progression to models involving thousands of levels in the close-coupling expansion required by the open d and f sub-shell systems of mid Z tungsten. This work also complements the electron-impact excitation of Fe-Peak elements required by astrophysics, in particular the near neutral species, which offer similar atomic structure challenges. Although electron-impact excitation work is our primary focus in terms of fusion application, the single photon photoionisation codes are also being developed in tandem, and benefit greatly from this ongoing work.

Electron-impact excitation of Ar 2+

Context: Emission from Ar III is seen in planetary nebulae, in H II regions, and from laboratory plasmas. The analysis of such spectra requires accurate electron impact excitation data. Aims: The aim of this work is to improve the electron impact excitation data available for Ar2+, for application in studies of planetary nebulae and laboratory plasma spectra. The effects of the new data on diagnostic line ratios are also studied. Methods: Electron-impact excitation collision strengths have been calculated using the R-Matrix Intermediate-Coupling Frame-Transformation method and the R-Matrix Breit-Pauli method. Excitation cross sections are calculated between all levels of the configurations 3s^23p^4, 3s3p^5, 3p^6, 3p^53d, and 3s^23p^3nl (3d ≤ nl ≤ 5s). Maxwellian effective collision strengths are generated from the collision strength data. Results: Good agreement is found in the collision strengths calculated using the two R-Matrix methods. The collision strengths are compared with literature values for transitions within the 3s^23p4 configuration. The new data has a small effect on Te values obtained from the I(λ7135 Å+ λ7751 Å)/ I(λ5192 Å) line ratio, and a larger effect on the Ne values obtained from the I(λ7135 Å)/I(λ9 μm) line ratio. The final effective collision strength data is archived online.

Atomic data for modelling fusion and astrophysical plasmas

Trends and focii of interest in atomic modelling and data are identified in connection with recent observations and experiments in fusion and astrophysics. In the fusion domain, spectral observations are included of core, beam penetrated and divertor plasma. The helium beam experiments at JET and the studies with very heavy species at ASDEX and JET are noted. In the astrophysics domain, illustrations are given from the SOHO and CHANDRA spacecraft which span from the solar upper atmosphere, through soft x-rays from comets to supernovae remnants. It is shown that non-Maxwellian, dynamic and possibly optically thick regimes must be considered. The generalized collisional-radiative model properly describes the collisional regime of most astrophysical and laboratory fusion plasmas and yields self-consistent derived data for spectral emission, power balance and ionization state studies. The tuning of this method to routine analysis of the spectral observations is described. A forward look is taken as to how such atomic modelling, and the atomic data which underpin it, ought to evolve to deal with the extended conditions and novel environments of the illustrations. It is noted that atomic physics influences most aspects of fusion and astrophysical plasma behaviour but the effectiveness of analysis depends on the quality of the bi-directional pathway from fundamental data production through atomic/plasma model development to the confrontation with experiment. The principal atomic data capability at JET, and other fusion and astrophysical laboratories, is supplied via the Atomic Data and Analysis Structure (ADAS) Project. The close ties between the various experiments and ADAS have helped in this path of communication.

The Rossiter-McLaughlin effect reloaded: Probing the 3D spin-orbit geometry, differential stellar rotation, and the spatially-resolved stellar spectrum of star-planet systems

When a planet transits its host star, it blocks regions of the stellar surface from view; this causes a distortion of the spectral lines and a change in the line-of-sight (LOS) velocities, known as the Rossiter-McLaughlin (RM) effect. Since the LOS velocities depend, in part, on the stellar rotation, the RM waveform is sensitive to the star-planet alignment (which provides information on the system’s dynamical history). We present a new RM modelling technique that directly measures the spatially-resolved stellar spectrum behind the planet. This is done by scaling the continuum flux of the (HARPS) spectra by the transit light curve, and then subtracting the infrom the out-of-transit spectra to isolate the starlight behind the planet. This technique does not assume any shape for the intrinsic local profiles. In it, we also allow for differential stellar rotation and centre-to-limb variations in the convective blueshift. We apply this technique to HD 189733 and compare to 3D magnetohydrodynamic (MHD) simulations. We reject rigid body rotation with high confidence (>99% probability), which allows us to determine the occulted stellar latitudes and measure the stellar inclination. In turn, we determine both the sky-projected (λ ≈ −0.4 ± 0.2◦) and true 3D obliquity (ψ ≈ 7+12 −4 ◦ ). We also find good agreement with the MHD simulations, with no significant centre-to-limb variations detectable in the local profiles. Hence, this technique provides a new powerful tool that can probe stellar photospheres, differential rotation, determine 3D obliquities, and remove sky-projection biases in planet migration theories. This technique can be implemented with existing instrumentation, but will become even more powerful with the next generation of high-precision radial velocity spectrographs.

Rates of E1, E2, M1, and M2 transitions in Ni II

Aims. We present rates for all E1, E2, M1, and M2 transitions among the 295 fine-structure levels of the configurations 3d⁹, 3d⁸4s, 3d⁷4s², 3d⁸4p, and 3d⁷4s4p, determined through an extensive configuration interaction calculation. 

Methods. The CIV3 code developed by Hibbert and coworkers is used to determine for these levels configuration interaction wave functions with relativistic effects introduced through the Breit-Pauli approximation. 

Results. Two different sets of calculations have been undertaken with different 3d and 4d functions to ascertain the effect of such variation. The main body of the text includes a representative selection of data, chosen so that key points can be discussed. Some analysis to assess the accuracy of the present data has been undertaken, including comparison with earlier calculations and the more limited range of experimental determinations. The full set of transition data is given in the supplementary material as it is very extensive. 

Conclusions. We believe that the present transition data are the best currently available.

Valence and L-shell photoionization of Cl-like argon using R-matrix techniques

Photoionization cross-sections are obtained using the relativistic DiracAtomic R-matrix Codes (DARC) for all valence and L-shell energy ranges between 27 and 270 eV. A total of 557 levels arising from the dominant configurations 3s²3p⁴, 3s³p5, 3p⁶, 3s²3p³[3d, 4s, 4p], 3p⁵3d, 3s²3p²3d², 3s³p4³d, 3s³p3³d2 and 2s²2p⁵3s²3p⁵ have been included in the targetwavefunction representation of the Ar III ion, including up to 4p in the orbital basis. We also performed a smaller Breit-Pauli (BP) calculation containing the lowest 124 levels. Direct comparisons are made with previous theoretical and experimental work for both valence shell and L-shell photoionization. Excellent agreement was found for transitions involving the ²P^o initial state to all allowed final states for both calculations across a range of photon energies. A number of resonant states have been identified to help analyse and explain the nature of the spectra at photon energies between 250 and 270 eV.

Ultraviolet emission lines of Si II in cool star and solar spectra

Recent atomic physics calculations for Si II are employed within the CLOUDY modelling code to analyse Hubble Space Telescope (HST) STIS ultraviolet spectra of three cool stars, β Geminorum, α Centauri A and B, as well as previously published HST/GHRS observations of α Tau, plus solar quiet Sun data from the High Resolution Telescope and Spectrograph. Discrepancies found previously between theory and observation for line intensity ratios involving the 3s²3p ²PJ-3s3p^{2 4}P^J' intercombination multiplet of Si II at ~ 2335 Å are significantly reduced, as are those for ratios containing the 3s²3p ²PJ-3s3p^{2 2}DJ ~ transitions at ~1816 Å. This is primarily due to the effect of the new Si II transition probabilities. However, these atomic data are not only very different from previous calculations, but also show large disagreements with measurements, specifically those of Calamai et al. for the intercombination lines. New measurements of transition probabilities for Si II are hence urgently required to confirm (or otherwise) the accuracy of the recently calculated values. If the new calculations are confirmed, then a long-standing discrepancy between theory and observation will have finally been resolved. However, if the older measurements are found to be correct, then the agreement between theory and observation is simply a coincidence and the existing discrepancies remain.

Electron-impact excitation of open d-shell ions

Astrophysics is driven by observations, and in the present era there are a wealth of state-of-the-art ground-based and satellite facilities. The astrophysical spectra emerging from these are of exceptional quality and quantity and cover a broad wavelength range. To meaningfully interpret these spectra, astronomers employ highly complex modelling codes to simulate the astrophysical observations. Important input to these codes include atomic data such as excitation rates, photoionization cross sections, oscillator strengths, transition probabilities and energy levels/line wavelengths. Due to the relatively low temperatures associated with many astrophysical plasmas, the accurate determination of electron-impact excitation rates in the low energy region is essential in generating a reliable spectral synthesis. Hence it is these atomic data, and the main computational methods used to evaluate them, which we focus on in this publication. We consider in particular the complicated open d- shell structures of the Fe-peak ions in low ionization stages. While some of these data can be obtained experimentally, they are usually of insufficient accuracy or limited to a small number of transitions.

A comparison of theoretical Mg VI emission line strengths with active-region observations from SERTS

R-matrix calculations of electron impact excitation rates in N-like Mg VI are used to derive theoretical electron-density-sensitive emission line ratios involving 2s²2p³ - 2s2p⁴ transitions in the 269-403 Å wavelength range. A comparison of these with observations of a solar active region, obtained during the 1989 flight of the Solar EUV Rocket Telescope and Spectrograph (SERTS), reveals good agreement between theory and observation for the 2s²2p^{3 4}S - 2s2p ^{4 4}p transitions at 399.28, 400.67, and 403.30 Å, and the 2s²2p^{3 2}p - 2s2p^{4 2}D lines at 387.77 and 387.97 Å. However, intensities for the other lines attributed to Mg VI in this spectrum by various authors do not match the present theoretical predictions. We argue that these discrepancies are not due to errors in the adopted atomic data, as previously suggested, but rather to observational uncertainties or mis-identifications. Some of the features previously identified as Mg VI lines in the SERTS spectrum, such as 291.36 and 293.15 Å, are judged to be noise, while others (including 349.16 Å) appear to be blended.

ALMA-resolved salt emission traces the chemical footprint and inner wind morphology of VY Canis Majoris

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

Localized structures in complex plasmas in the presence of a magnetic field

In this work, the general framework in which fits our investigation is that of modeling the dynamics of dust grains therein dusty plasma (complex plasma) in the presence of electromagnetic fields. The generalized discrete complex Ginzburg-Landau equation (DCGLE) is thus obtained to model discrete dynamical structure in dusty plasma with Epstein friction. In the collisionless limit, the equation reduces to the modified discrete nonlinear Schrödinger equation (MDNLSE). The modulational instability phenomenon is studied and we present the criterion of instability in both cases and it is shown that high values of damping extend the instability region. Equations thus obtained highlight the presence of soliton-like excitation in dusty plasma. We studied the generation of soliton in a dusty plasma taking in account the effects of interaction between dust grains and theirs neighbours. Numerical simulations are carried out to show the validity of analytical approach.