Application of factor analysis and geostatistics to identify processes controlling the distribution and bioaccessibility of trace elements in Northern Ireland soils

Geologic and environmental factors acting over varying spatial scales can control
trace element distribution and mobility in soils. In turn, the mobility of an element in soil will affect its oral bioaccessibility. Geostatistics, kriging and principal component analysis (PCA) were used to explore factors and spatial ranges of influence over a suite of 8 element oxides, soil organic carbon (SOC), pH, and the trace elements nickel (Ni), vanadium (V) and zinc (Zn). Bioaccessibility testing was carried out previously using the Unified BARGE Method on a sub-set of 91 soil samples from the Northern Ireland Tellus1 soil archive. Initial spatial mapping of total Ni, V and Zn concentrations shows their distributions are correlated spatially with local geologic formations, and prior correlation analyses showed that statistically significant controls were exerted over trace element bioaccessibility by the 8 oxides, SOC and pH. PCA applied to the geochemistry parameters of the bioaccessibility sample set yielded three principal components accounting for 77% of cumulative variance in the data
set. Geostatistical analysis of oxide, trace element, SOC and pH distributions using 6862 sample locations also identified distinct spatial ranges of influence for these variables, concluded to arise from geologic forming processes, weathering processes, and localised soil chemistry factors. Kriging was used to conduct a spatial PCA of Ni, V and Zn distributions which identified two factors comprising the majority of distribution variance. This was spatially accounted for firstly by basalt rock types, with the second component associated with sandstone and limestone in the region. The results suggest trace element bioaccessibility and distribution is controlled by chemical and geologic processes which occur over variable spatial ranges of influence.

Modification of the formation of high-Mach number electrostatic shock-like structures by the ion acoustic instability

The formation of unmagnetized electrostatic shock-like structures with a high Mach number is examined with one- and two-dimensional particle-in-cell (PIC) simulations. The structures are generated through the collision of two identical plasma clouds, which consist of equally hot electrons and ions with a mass ratio of 250. The Mach number of the collision speed with respect to the initial ion acoustic speed of the plasma is set to 4.6. This high Mach number delays the formation of such structures by tens of inverse ion plasma frequencies. A pair of stable shock-like structures is observed after this time in the 1D simulation, which gradually evolve into electrostatic shocks. The ion acoustic instability, which can develop in the 2D simulation but not in the 1D one, competes with the nonlinear process that gives rise to these structures. The oblique ion acoustic waves fragment their electric field. The transition layer, across which the bulk of the ions change their speed, widens and their speed change is reduced. Double layer-shock hybrid structures develop.

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.

Electron impact excitation of Be-like ions: a comparison of DARC and ICFT results

Emission lines of Be-like ions are frequently observed in astrophysical plasmas, and many are useful for density and temperature diagnostics. However, accurate atomic data for energy levels, radiative rates (A-values) and effective electron excitation collision strengths ($\Upsilon$) are required for reliable plasma modelling. In general it is reasonably straightforward to calculate energy levels and A- values to a high level of accuracy. By contrast, considerable effort is required to calculate $\Upsilon$, and hence it is not always possible to assess the accuracy of available data. Recently, two independent calculations (adopting the $R$-matrix method) but with different approaches (DARC and ICFT) have appeared for a range of Be-like ions. Therefore, in this work we compare the two sets of $\Upsilon$, highlight the large discrepancies for a significant number of transitions and suggest possible reasons for these.

Temporal structure of attosecond pulses from laser-driven coherent synchrotron emission

The microscopic dynamics of laser-driven coherent synchrotron emission transmitted through thin foils are investigated using particle-in-cell simulations. For normal incidence interactions, we identify the formation of two distinct electron nanobunches from which emission takes place each half-cycle of the driving laser pulse. These emissions are separated temporally by 130 attoseconds and are dominant in different frequency ranges, which is a direct consequence of the distinct characteristics of each electron nanobunch. This may be exploited through spectral filtering to isolate these emissions, generating electromagnetic pulses of duration ~70 as.

Photoionization of Co+ and electron-impact excitation of Co2 + using the Dirac R-matrix method

Modelling of massive stars and supernovae (SNe) plays a crucial role in understanding galaxies. From this modelling we can derive fundamental constraints on stellar evolution, mass-loss processes, mixing, and the products of nucleosynthesis. Proper account must be taken of all important processes that populate and depopulate the levels (collisional excitation, de-excitation, ionization, recombination, photoionization, bound–bound processes). For the analysis of Type Ia SNe and core collapse SNe (Types Ib, Ic and II) Fe group elements are particularly important. Unfortunately little data is currently available and most noticeably absent are the photoionization cross-sections for the Fe-peaks which have high abundances in SNe. Important interactions for both photoionization and electron-impact excitation are calculated using the relativistic Dirac atomic R-matrix codes (DARC) for low-ionization stages of Cobalt. All results are calculated up to photon energies of 45 eV and electron energies up to 20 eV. The wavefunction representation of Co III has been generated using GRASP0 by including the dominant 3d⁷, 3d⁶[4s, 4p], 3p⁴3d⁹ and 3p⁶3d⁹ configurations, resulting in 292 fine structure levels. Electron-impact collision strengths and Maxwellian averaged effective collision strengths across a wide range of astrophysically relevant temperatures are computed for Co III. In addition, statistically weighted level-resolved ground and metastable photoionization cross-sections are presented for Co II and compared directly with existing work.