LOW Theoretical Background

The theoretical sensitivity is computed using direction, frequency, and LST-dependent SEFD values following the procedure described in Sokolowski et al. (2022). See the linked paper for more details.

Sensitivity of a beamformed observation

The LOW sensitivity calculator backend uses the modified radiometer equation to estimate the rms noise achievable in a coherently beamformed tied-array observation. Following Morello et al. (2020), the detection threshold of a coherently beamformed tied-array observation with \(\mathrm{N}_\mathrm{ant}\) is given by the equation

(1)\[\mathrm{S}_\mathrm{rms} = \frac{\beta}{\epsilon} \frac{\mathrm{SEFD}_\mathrm{I}}{\sqrt{\mathrm{t}_\mathrm{int} \Delta\nu \mathrm{N}_\mathrm{ant}^2 }} \sqrt{\frac{W}{P-W}}.\]

In the above equation,

  • \(\beta\) is the digitisation loss. Following Kouwenhoven & Voute (2001), we take \(\beta=1.05\) for 8-bit quantisation.

  • \(\epsilon\) is the search efficiency. For FFT-based searches, \(\epsilon=0.75\) Morello et al. (2020).

  • \(\mathrm{SEFD}_\mathrm{I}\) is the System Equivalent Flux Density for Stokes I in Jy.

  • \(\mathrm{t}_\mathrm{int}\) is the observation duration in seconds.

  • \(\Delta \nu\) is the observation bandwidth in Hz.

  • \(\mathrm{N}_\mathrm{ant}\) is the number of antennas.

  • \(\mathrm{W}\) is the observed pulse width in milliseconds.

  • \(\mathrm{P}\) is the pulse period in milliseconds.

A pulse observed at the telescope is broader than the intrinsic pulse width due to propagation effects in the interstellar medium (ISM) like dispersion and scattering.

Pulse broadening (in milliseconds) due to dispersion measure \(\mathrm{DM}\) (in \(\mathrm{pc}~\mathrm{cm}^{-3}\)) can be calculated using the relation (Lorimer & Kramer 2004)

\[\mathrm{W}_\mathrm{DM} = \mathrm{K}_\mathrm{DM} \times \mathrm{DM} \times \frac{\Delta \nu}{\nu^3},\]

where \(\mathrm{K}_\mathrm{DM}=8.3 \times 10^6~\mathrm{ms}\), \(\Delta \nu\) is the channel frequency resolution in MHz, and \(\nu\) is the representative frequency (usually taken to be the band centre) in MHz.

Pulse broadening due to scattering in the ISM (\(\mathrm{W}_\mathrm{scat}\)) is given by the relation (Bhat et al. 2004)

\[\log (\mathrm{W}_\mathrm{scat}) = -6.46 + 0.154 \log (\mathrm{DM}) + 1.07 \log (\mathrm{DM})^2 - 3.86 \log (\nu_\mathrm{GHz})\]

where \(\nu_\mathrm{GHz}\) is frequency in GHz.

The observed pulse width \(\mathrm{W}\) is related to the intrinsic pulse width \(\mathrm{W}_\mathrm{int}\) by the relation

\[\mathrm{W} = \sqrt{\mathrm{W}_\mathrm{int}^2 + \mathrm{W}_\mathrm{DM}^2 + \mathrm{W}_\mathrm{scat}^2}.\]

Note that the equation for beamformed sensitivity (1) can be expressed in terms of Stokes I continuum sensitivity. If \(\mathrm{S}_\mathrm{cont}\) is the continuum sensitivity, assuming that the observation is carried out using the same number of stations, bandwidth, and observation duration, the beamformed sensitivity can then be rewritten as

\[\mathrm{S}_\mathrm{rms} = \frac{\beta}{\epsilon} \mathrm{S}_\mathrm{cont} \sqrt{\frac{\mathrm{N}_\mathrm{ant}-1}{\mathrm{N}_\mathrm{ant}}} \sqrt{\frac{\mathrm{W}}{\mathrm{P}-\mathrm{W}}}.\]

Dynamic range warning

In addition to calculating the theoretical sensitivity, the LOW sensitivity calculator backend warns if it detects a bright (>10 Jy) off-axis source within the full width at half maximum (FWHM) of the LOW station beam calculated at the centre of the specified band. The backend uses a bright source catalogue derived from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey (Hurley-Walker et al. (2017)).