An efficient technique for pre-selecting low-redshift damped Lyα systems

The number of z ∼ 1 damped Lyα systems (DLAs, log  N (H  i ) ≥ 20.3) per unit redshift is approximately 0.1, making them relatively rare objects. Large, blind QSO surveys for low-redshift DLAs are therefore an expensive prospect for space-borne ultraviolet telescopes. Increasing the efficiency of these surveys by pre-selecting DLA candidates based on the equivalent widths (EWs) of metal absorption lines has previously been a successful strategy. However, the success rate of DLA identification is still only ∼35 per cent when simple EW cut-offs are applied, the majority of systems having 19.0 < log  N (H  i ) < 20.3. Here, we propose a new way to pre-select DLA candidates. Our technique requires high-to-moderate-resolution spectroscopy of the Mg  ii λ2796 transition, which is easily accessible from the ground for 0.2 ≲ z ≲ 2.4. We define the D -index, the ratio of the line equivalent width to velocity spread, and measure this quantity for 19 DLAs and eight sub-DLAs in archival spectra obtained with echelle spectrographs. For the majority of absorbers, there is a clear distinction between the D -index of DLAs compared with sub-DLAs (Kolmogorov–Smirnov probability = 0.8 per cent). Based on this pilot data sample, we find that the D -index can select DLAs with a success rate of up to 90 per cent, an increase in selection efficiency by a factor of 2.5 compared with a simple EW cut. We test the applicability of the D -index at lower resolution and find that it remains a good discriminant of DLAs for full width at half-maximum (FWHM) ≲ 1.5 Å. However, the recommended D -index cut-off between DLAs and sub-DLAs decreases with poorer resolution and we tabulate the appropriate D -index values that should be used with spectra of different resolutions.     

Non-Sobolev modelling of radiation-pressure-driven flows in active galactic nuclei

We present a new general scheme for calculating the structure and dynamics of radiation-pressure-driven photoionized flows. The new method goes one step beyond the Sobolev approximation. It involves a numerical solution of the radiative transfer in absorption lines, including the effects of differential expansion and line interactions such as line locking and blanketing. We also present a new scheme for calculating the radiation pressure due to trapped line photons in finite, differentially expanding flows. We compare our results for the radiation pressure force with those obtained using the Sobolev approximation and show the limitations of the latter. In particular, we demonstrate that the Sobolev method gives a poor approximation near discontinuity surfaces and its neglect of line blanketing can lead to erroneous results in high-velocity flows. We combine the newly calculated radiation pressure force with self-consistent photoionization and thermal calculations to study the dynamics and spectral features of broad absorption-line flows and highly ionized gas flows in active galactic nuclei (AGN). A comparison with Sobolev-type calculations shows that the latter overestimates the terminal velocity of the flow and, conversely, underestimates its opacity. We also show that line locking on broad emission lines can have a significant effect on the dynamics and spectral features of AGN flows.     

An HLLC Riemann solver for relativistic flows – II. Magnetohydrodynamics

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The Harten–Lax–van Leer contact wave (HLLC) solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave.
Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the constrained transport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one- and two-dimensional numerical test problems.