Issue 
A&A
Volume 638, June 2020



Article Number  A37  
Number of page(s)  14  
Section  Interstellar and circumstellar matter  
DOI  https://doi.org/10.1051/00046361/201833157  
Published online  09 June 2020 
Constraining the redshifts of unlocalised fast radio bursts^{★}
^{1}
Jodrell Bank Center for Astrophysics, The University of Manchester,
Alan Turing Building,
Manchester,
M13 9PL,
UK
email: charles.walker@postgrad.manchester.ac.uk; rene.breton@manchester.ac.uk
^{2}
School of Chemistry and Physics, University of KwaZuluNatal,
Westville Campus, Private Bag X54001,
Durban
4000, South Africa
email: ma@ukzn.ac.za
^{3}
NAOCUKZN Computational Astrophysics Center (NUCAC), University of KwaZuluNatal,
Durban
4000,
South Africa
Received:
4
April
2018
Accepted:
20
January
2020
Context. The relationship between the dispersion measures (DMs) and redshifts of fast radio bursts (FRBs) is of scientific interest. Upcoming commensal surveys may detect and localise many FRBs to the subarcsecond angular resolutions required for accurate redshift determination. Meanwhile, it is important to exploit sources accumulated with more limited localisation to their maximum scientific potential.
Aims. We present techniques for the DMredshift analysis of large numbers of unlocalised FRBs, accounting for uncertainties due to their extragalactic DM components, redshift dependences, and progenitor scenarios.
Methods. We reviewed the components comprising observed FRB DMs. We built redshiftscalable probability distribution functions for these components, which we combined in cases of multiple progenitor scenarios. Accounting for prior FRB redshift distributions we inverted these models, enabling FRB redshifts to be constrained.
Results. We illustrate the influence of FRB progenitors on their observed DMs, which may remain significant to redshift z ~ 3. We identify the FRB sample sizes required to distinguish between multiple progenitor scenarios. We place new, physically motivated redshift constraints on all catalogued FRBs to date and use these to reject potential host galaxies in the localisation area of an FRB according to various models. We identify further uses for DMredshift analysis using many FRBs. We provide our code so that these techniques may be employed using increasingly realistic models as our understanding of FRBs evolves.
Key words: techniques: miscellaneous / stars: neutron / pulsars: general / local insterstellar matter / galaxies: halos / intergalactic medium
Python code and examples are available at https://doi.org/10.5281/zenodo.1209920
© C. R. H. Walker et al. 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Many fast radio bursts (FRBs) have unconstrained redshifts. Since their discovery in 2007 (Lorimer et al. 2007), tens of these shortduration (~0.1–10 ms), dispersed radio bursts have been catalogued (Petroff et al. 2016)^{1}. Save for the repeating FRB 121102 (Spitler et al. 2014, 2016), their extragalactic origins have been inferred via their dispersion measures (DMs).
These DMs, related to the frequencydependent index of refraction and resulting arrival time delays experienced by radio waves which propagate through ionised electrons, are proportional to the integrated electron densities along their propagation paths. All FRBs have DMs measured in excess of that which may be afforded by electrons of the Milky Way alone.
Extragalactic electrons may lie in the intergalactic medium (IGM), the host galaxy of the FRB and the environment around the source (Macquart et al. 2015). According to relativistic theories, the observed time delay associated with these electrons is modified. This, along with sightline variance in IGM electrons within collapsed structures results in a redshiftdependent potential DM distribution for FRBs (Ioka 2003; McQuinn 2014; Macquart et al. 2015). This distribution is additionally intrinsically linked to the locations of FRBs within their host galaxies and thus to their progenitors.
Analysis of the FRB DMredshift relationship is of astronomical significance. Macquart et al. (2015) review the scientific potential of redshiftconstrained FRBs, including using 10^{4} as probes of the baryonic matter of the Universe and 10^{3} as cosmic rulers for measuring the geometry of the Universe beyond redshift z ~ 2. Statistical studies into the potential IGM (see Zhou et al. 2014; McQuinn 2014) and host galaxy (see Xu & Han 2015; Yang et al. 2017) contributions to FRB DMs have been undertaken to compare different models accounting for the structure of the IGM or the spacial distributions of FRBs within host galaxies, for example. In practice, many FRBs with welldetermined redshifts must be accumulated to exploit such studies fully.
Redshifts may be obtained spectroscopically via host galaxy association using interferometry (e.g. Tendulkar et al. 2017), however many FRBs have been (and will continue to be) detected by telescopes with large fields of view (FoVs) containing many galaxies, within which their locations may be unconstrained. The FoV of a Parkes Lband multibeam receiver (half power beam width; HPBW ~ 14.4′), for example, may contain 65 000 galaxies out to redshift 6 (Scoville et al. 2007)^{2}. Molonglo’sUTMOST may improve localisation to 15″ × 15″ (Caleb et al. 2017); ASKAP’s fly’seye mode has demonstrated localisation of an FRB to 8′ × 8′ (Bannister et al. 2017) and aims for <7″ localisation when fully operational (Macquart et al. 2010); and CHIME will improve localisation for potentially tens of FRBs per day to tens of arcminutes (Ng et al. 2017; CHIME/FRB Collaboration 2018). But the large FoVs and onsky time available to resolutionlimited, single dish and fly’seye surveys suggest that many FRBs will be discovered to worse than the subarcsecond accuracy necessary for unambiguous association with a single host (see e.g. Mahony et al. 2018).
Standard calculation of redshifts for such unlocalised FRBs involves estimation and removal of the sightlinedependent DM component attributed to Milky Way electrons (using Galactic electron models, e.g. NE2001, Cordes & Lazio 2002 or YMW, Yang et al. 2017). This may be followed by simplification of the theoretical relationship (see Ioka 2003) between redshift and the IGM DM component (see e.g. Petroff et al. 2016), assuming idealised hydrogen and helium ionisation fractions and an average line of sight. Finally, no host DM component (e.g. Petroff et al. 2016) or a fixed host component (e.g. ~ 100 pc cm^{−3}; Lorimer et al. 2007; Thornton et al. 2013) may be assumed. Associated uncertainties from, for example redshiftdependent variance in collapsed systems (e.g. galactic haloes) along IGM sightlines (potentially contributing ~ 180−400 pc cm^{−3} at redshift 1; McQuinn 2014) or source positions in host galaxies (potentially contributing thousands of pc cm^{−3} in the rest frame, Xu & Han 2015; Yang et al. 2017), may not be considered.
These uncertainties leave many catalogued FRBs with redshifts that are not ideally constrained, which hinders analysis of the DMredshift relationship. The purposes of this paper are as follows: Firstly, we aim to introduce a framework for performing physically motivated DMredshift analysis for large numbers of unlocalised FRBs, accounting for uncertainties such as (a) potential host contributions to DM (following e.g. Xu & Han 2015; Yang et al. 2017), (b) the cosmic variance of IGM structure (following e.g. McQuinn 2014), and (c) FRB progenitor evolution across cosmic time (following e.g. Zhou et al. 2014). Secondly, we apply physically motivated models to this framework, thus illustrating the effects of different progenitor scenarios on the probability distribution functions (PDFs) of observed DMs and redshifts of FRBs. Lastly, we discuss the implications of our models, to provide constrained redshift values for currently unlocalised FRBs and to identify further uses for FRB DMredshift analysis.
In Sect. 2 we summarise the extragalactic DM components which contribute to the observed arrival time delay between two FRB frequencies and review their respective redshift dependences. We discuss our choices for physically motivated PDFs of these DM components. In Sect. 3 we present a framework for the combination of these PDFs to retrieve the overall likelihood of an FRB of given redshift having a particular excess DM, or vice versa. We demonstrate this method using our models. Section 4 discusses the implications of our findings.
Unless otherwise stated, in this paper we adopt a spatially flat ΛCDM cosmology model with the cosmological parameters fixed to the bestfit values of Planck Collaboration VI (2020): Ω_{b} = 0.048, Ω_{m} = 0.309, Ω_{Λ} = 0.691, n_{s} = 0.9608, σ_{8} = 0.815, H_{0} = 67.3 km s^{−1} Mpc^{−1}.
2 Contributions to FRB dispersion measure
The observed time delay Δt_{obs} between twofrequencies ν_{1, obs} and ν_{2,obs} may be written as a function of dispersion measure DM_{obs} (Manchester & Taylor 1972). For an extragalactic FRB, DM_{obs} may be deconstructed into components (1)
representing contributions from electrons along the FRB light path in the MW, IGM, its host galaxy, and its local environment (dependent on its progenitor model), respectively (Macquart et al. 2015). The excess extragalactic contribution to DM_{obs} (i.e. DM_{IGM} + DM_{host} + DM_{local}) is referred to in this work as DM_{exc}.
2.1 Milky Way
For a given FRB, DM_{MW} is typically computed using Galactic electron distributions (e.g. NE2001 Cordes & Lazio 2002, YMW Yang et al. 2017) modelling the Milky Way interstellar medium (ISM). A contribution from the Galactic halo of the Milky Way (estimations in the range from 2–30 pc cm^3 as reviewed by Mahony et al. 2018) may also be considered. Removal of a completely accurate DM_{MW} leaves a remaining DM_{exc} component solely associated with extragalactic electron distributions. Uncertainties associated with Milky Way electron models may be considered (e.g. in Keane 2016). Through analysis of pulsars with parallaxdetermined distances, the potential of such models to misrepresent the DM_{MW} contribution to DM_{obs} has been demonstrated, the most extreme result of which could be the misclassification of an FRB as another transient phenomenon, or vice versa. In particular Keane (2016) highlight a source, J1354+24, as having a 20–40 % probability of being an FRB wrongly labelled as an RRAT. For the purposes of this paper we assume correct DM_{MW} values and discard these to focus on extragalactic, redshiftdependent contributions to uncertainty (see Sect. 4.1 for a discussion).
2.2 Intergalactic medium
While 40% of the baryonic matter budget of the Universe is readily accounted for (e.g. in galaxies, Xray corona, and warm intergalactic matter), the remaining portion is less easily detectable. Recent use of the thermal SunyaevZel’dovich effect (e.g. de Graaff et al. 2019; Tanimura et al. 2019) and quasar absorption line observations (e.g. Nicastro et al. 2018) reveal largescale filamentary structure in the IGM as a potential source of this remaining matter. As all IGM electrons along an FRB propagation path contribute to DM_{IGM}, FRBs may similarly serve as probes of elusive baryon distributions (McQuinn 2014).
The DM_{IGM} is source redshift (z_{s}) dependent. During standard DM_{IGM}redshift analysis, an ionisation fraction profile for the IGM is assumed, accounting for the reionisation history of both hydrogen and helium^{3}. A linear approximation to the resulting ⟨DM_{IGM}(z_{s})⟩ profile, which is the mean IGM contribution to FRB DMs as a function of redshift when averaged over all sightlines, is then used to estimate z_{s} for observed FRBs. These linear relationships take the form z_{s} = (DM_{obs} −DM_{MW})∕X pc cm^{−3}, where X is constant and dependent on the exact ionisation fractions assumed (see e.g. Petroff et al. 2016; Keane 2018; Zhang 2018). Additionally, sightline variability in the number of collapsed systems (e.g. galactic haloes) along FRB propagation paths results in a variance around ⟨DM_{IGM}(z_{s})⟩ for an FRB population (Ioka 2003; McQuinn 2014). To accurately contain FRB DMredshift constraints, this variance must be considered.
McQuinn (2014) simulate PDFs describing the likelihood of an FRB originating at source redshift z_{s} contributing an IGM component DM_{IGM} to the total observed DM_{obs}, according to various galactic halo electron distributions. Following this work, we consider a Gaussian approximation P(DM_{IGM}z_{s}) for this distribution with redshiftevolving mean DM_{IGM}(z_{s}) and variance . We make use of the matter power spectrum from public code CAMB^{4}, with Planck bestfitting cosmological parameters (Planck Collaboration VI 2020). More complex profiles may be substituted in the future. Our approximation, derived in full in Appendix A, is shown in Fig. 1.
Fig. 1 Probability distribution functions (P(DM_{IGM})) for FRBs of source redshift z_{s} = 1 accounting for DM_{IGM} sightline variance due to collapsed systems. The Gaussian approximation assumed in this work is compared to more complex McQuinn (2014) baryonic halo models (tracing top hat functions, NFW profiles, and cosmological simulations, respectively, obtained via. private communication). Profiles are normalised to ∫ P(DM z_{s} = 1) d DM = 1000. 
2.3 Host galaxy
As discussed by, for example Ioka (2003), Zhou et al. (2014) and Macquart et al. (2015), the arrival time delay Δt_{obs,host} between observing frequencies ν_{1,obs} and ν_{2,obs} for an FRB of source redshift z_{s} due to a restframe host DM contribution DM_{host, r} depends on the following:
 1.
The relationship between the frequency of the radio wave as emitted in the rest frame of a receding galaxy ν_{r} and the observed frequency ν_{obs} (2)
 2.
The relationship between a time delay measured in the rest frame of the galaxy Δt_{r,host} and the observed Δt_{obs,host} due to cosmological time dilation, (3)
The DM_{host} contributions to DM_{obs} have previously been investigated and found to be dependent on both FRB host galaxy types and the source locations within them. Ioka (2003) explored host galaxy components for GRBs, finding potential contributions of up to 10^{5} pc cm^{−3} for sources in starforming regions. Xu & Han (2015) evaluate FRB DM_{host} values from restframe spiral and elliptical galaxies by scaling NE2001, finding potential contributions of approximately thousands/tens of pc cm^{−3}, respectively. Yang et al. (2017) estimate the mean host contribution for 21 observed FRBs to be ~ 270 pc cm^{−3}. Accounting for DM_{host} variation during DMredshift analysis (particularly of low z_{s} events) is therefore necessary.
In this work we simulate restframe (z_{s} = z_{0} = 0) host galaxy electron distribution models, which we populate with FRBs according to various progenitor scenarios. By integrating electron densities along random sightlines we build restframe PDFs describing the likelihood of host galaxies contributing DM_{host} components to DM_{obs}. For demonstrative purposes, we model static galaxies which do not evolve with redshift^{5}. By letting P(U)dU = P(x)dx with DM_{host} = U(x) and DM_{host,r} = x the redshiftdependent PDF becomes (6)
Our models are described below.
2.3.1 Spiral galaxies
As simple spiral galaxy electron distribution models (see Schnitzeler 2012 for a review) have previously been proven sufficient to reproduce the DMs of Milky Way pulsars with known parallax distances to within a factor 1.5–2, rather than scaling NE2001 (see e.g. Xu & Han 2015) to demonstrate our DMredshift analysis framework, we opt for a simpler model comprised of two components. The first is the GBCa model (as denoted by Schnitzeler 2012), which was developed using 109 pulsars with DMindependently determined distances Gómez et al. (2001). This model is a sum of two ellipsoids made from decaying exponentials and describes the inner Milky Way Galactic disc, where the majority of the stellar mass of the Milky Way resides. It is well constrained within 4–12 kpc of the Galactic centre. We make the assumption that the model is acceptable within these central 4 kpc for the demonstrative purposes of this work. The second component is a model tracing neutral hydrogen (HI), which was developed by Kalberla & Dedes (2008). This model is a single decaying exponential and approximates the average midplane volume density of the Milky Way disc between 7 and 35 kpc. We combine these two components into a smoothly varying function as follows: (7)
where cylindrical coordinates r and h denote the radial distance from the Galactic centre and height above the Galactic plane, respectively. All relevant parameters are provided in Table 1. To reconcile the two components, the HI model is multiplied by an ionisation fraction of ~ 15%. At a radius ~50 kpc the model reaches an electron density consistent with the IGM at redshift zero as estimated by Ioka (2003). The model is shown in Fig. 2.
Fig. 2 Radial and height cross sections of our spiral galaxy electron density model n_{e,spiral} (r, h) (pink lines). Constituent components of the model tracing the inner MW Galactic disc (Gómez et al. 2001; Schnitzeler 2012) (blue lines) and an ionised outer HI disc model (Kalberla & Dedes 2008) are shown (brown lines). Shaded (blue/brown) regions highlight areas in which the inner/outer profiles are considered well constrained. The horizontal black dotted line indicates the approximate electron density for the IGM at z_{s} = 0 (Ioka 2003). 
Spiral galaxy electron distribution model parameters.
2.3.2 Elliptical galaxies
We use a β model to describe elliptical galaxy electron distributions, i.e. (8)
following Brown & Bregman (2001), Mamon & Łokas (2005). The model traces hot gas density within the virial radius as a function of galactrocentric radius r, central density ρ_{0}, β_{g} = 0.5, and r_{c} as follows: (9)
where R_{e} is the effective (halflight) galaxy radius and q = 10. Assuming hot gas traces free electrons, we modify the model to approach the local IGM electron density at the virial radius of a galaxy (~ 70 kpc) via multiplication with an ionisation fraction (F = 0.04). The relationships between total galaxy mass, central density, and effective radius for the β model are provided by Mamon & Łokas (2005). We model a galaxy of total virial mass 5 × 10^{10} M_{⊙} with ρ_{0} ≃ 6 × 10^{7} M_{⊙} kpc^{−3} and R_{e} = 3.2 kpc.
2.3.3 Distributions of FRBs within galaxies
Spacial distributions of FRBs within host galaxies rely on their progenitors. To model potential scenarios we populate our host galaxies with FRBs of different x,y,z coordinates and calculate DMs for randomly orientated observers at their boundaries. The resulting P(DM_{host}z_{s} = 0) PDFs are shown in Fig. 3. Our modelling choices are discussed below.
Stellar distributions in spiral galaxies
Stellar populations and phenomena potentially associated with FRBs include massive progenitor stars and superluminous supernovae (Tendulkar et al. 2017), magnetar outbursts (Pen & Connor 2015), pulsar supergiant pulses (Cordes & Wasserman 2016), collapsing (Falcke & Rezzolla 2014) or merging (Totani 2013) neutron stars, and pulsar companions (Mottez & Zarka 2014). Such progenitors should follow the spacial distributions of young stellar populations (e.g. OB stars) or neutron stars. We opted to model OB stars, pulsar populations of two ages (hereafter YPSRs, OPSRs) separated by a 8 Myr characteristic age, and millisecond pulsars (MSPs) separately to account for displacement due to, for example, supernova kicks.
Height distributions as a function of height h above the galactic plane (H_{*}(h), where * denotes the model type) for OB stars, YPSRs, and OPSRs are modelled following MaízApellániz (2001); Sun & Han (2004) and written as (10) (11) (12)
where the scale height h_{PSR} = h_{0} + σ_{PSR}t. The height distribution for MSPs is modelled after Lorimer (2013) as follows: (13)
Parameters are provided in Table 2.
Radial distributions as a function of galactocentric radius (R_{*} (r), where * denotes the model type) for YPSRs and OPSRs, are modelled with a gamma function following Yusifov & Küçük (2004) and with a Gaussian distribution for MSPs and OB stars following Lorimer (2013) and written as (14)
where X = r + r_{1} and X_{⊙} = r_{⊙} + r_{1}. (15)
Parameters are provided in Table 2.
Fig. 3 Restframe PDFs of relative likelihood against DM_{host} for multiple FRB progenitor scenarios. From top left to bottom right: sources follow stellar distributions (OB stars, YPSRs, OPSRs, MSPs) in spiral galaxies; homogeneous distributions (within spheres of radius H kpc) in spiral galaxies; elliptical distributions (confined to multiples of the galaxy effective radius, R_{e} = 3.2 kpc) in elliptical galaxies; and homogeneous distributions (within spheres of radius H kpc) in elliptical galaxies. 
Parameters used in pulsar distribution models.
Other distributions
We populate elliptical galaxies assuming their stellar mass distribution (and thus ionised electron distribution) to be a tracer for FRBs. We achieve this with a source distribution, following n_{e,elliptical}(r) normalised to the number simulated sources. We also generate homogeneously distributed sources within spheres of variable radius for both spiral and elliptical galaxies. Such distributions attempt to model central galactic halobased progenitors. At smaller limiting radii these distributions could serve as models of AGNrelated FRB mechanisms (e.g. relativistic jetplasma interactions Romero et al. 2016; Vieyro et al. 2017). Larger limiting radii could be consistent with progenitors located in more extended galactic haloes. We note that the most accurate DM distributions at the smallest limiting radii require electron distribution models with welldefined central profiles.
2.4 Burst environment
The DM_{local} contributions to DM_{obs} are also progenitordependent. Progenitor theories of FRB are diverse and may be separated into cataclysmic and noncataclysmic events. The repeating FRB 121102, by virtue of its repetition, must be noncataclysmic. Subarcsecond localisation has enabled study of its environment. This object is hosted by a lowmetallicity dwarf (Tendulkar et al. 2017) and is associated with a compact, persistent, offcentre (with respect to the optical centre of the galaxy) variable radio source within a starforming region, which may be a lowluminosity AGN that is powered supernova remnant or pulsar wind nebula (Chatterjee et al. 2017; Marcote et al. 2017; Bassa et al. 2017). The extreme and variable rotation measures of the burst indicate possible proximity of a neutron star and massive accreting black hole, and parallels to PSR J17452900 and Sagittarius A* have been drawn (Michilli et al. 2018)^{6}. Analysis suggests that FRB 121102’s DM_{local} could be small (Tendulkar et al. 2017). Such studies have potential implications for progenitor models, however FRB 121102 may not necessarily be representative of all FRBs. More accurate modelling of DM_{local} may arise as further FRBs are localised. In this work we assume DM_{local} = 0.
3 Excess electron model
In this section we present the framework for combining the PDFs of DM_{exc} components to retrieve a PDF, which describes the likelihood of an FRB at source redshift z_{s} to have an observed extragalactic DM component DM_{exc}, contributing to a measured time delay Δt_{obs,exc} between two observing frequencies ν_{1,obs} and ν_{2,obs} on Earth. We demonstrate this using our simulated IGM and progenitor models. Increasingly, accurate component models developed using FRB discoveries (or via other means) may allow these PDFs, which we refer to as excess electron models, to become increasingly useful tools for DMredshift analysis. In Sect. 3.1 we describe our method and in Sect. 3.2 we apply Bayes’ theorem to the excess electron model and discuss advantages of doing so.
3.1 P(DMz) methodology and results
Following Eq. (1), an FRB emitted from a source at redshift z_{s} has a DM_{exc} formed from the sum of its extragalactic components. Assuming these components to be random independent variables drawn from their respective underlying PDFs, P(DM_{IGM}z_{s}) and P(DM_{host}z_{s}), the convolution of these PDFs, (16)
describes the likelihood of an FRB of given source redshift z_{s} to have an observed extragalactic DM component DM_{exc}. This DM_{exc} results in a measured time delay Δt_{obs,exc} between two observing frequencies (ν_{1,obs} and ν_{2,obs}) on Earth.
Example excess electron models generated using our IGM (Sect. 2.2) and host galaxy (Sect. 2.3) PDFs are shown in Fig. 4. Comparisons between projections of these models for different values of z_{s} and P(DM_{IGM}z_{s}) illustrate the initially significant (and diminishing with redshift) influence of the host galaxy on P(DM_{exc}z_{s}). These comparisons show that the potential host galaxy contribution to DM_{exc} increases the likelihood of an FRB to have a DM_{exc} in excess of that predicted by conventional analysis. This excess may be particularly large for sources at low redshifts. In some cases (see e.g. subplots (a) and (c), which trace YPSR distributions in spiral galaxies and elliptical distributions in elliptical galaxies, respectively), the influence of the host diminishes significantly by z_{s} = 3, signified by the convergence of the P(DM_{exc}) curves with their IGM counterparts. In other cases (see e.g. subplots (c) and (f), where FRBs are confined to galactic regions of high electron density), the host’s potential influence is still significant at z_{s} = 3.
Differences between P(DM_{exc}z_{s}) profiles are also apparent between different progenitor models, particularly for sources with low redshifts. The DM_{exc} number density profiles acquired for large samples of redshiftconstrained FRBs could therefore prove valuable tools for FRB progenitor analysis (see Sect. 4.3). We note that each of our models considers a single progenitor scenario. In reality, a combination of these models may be more representative of the true FRB population.
3.2 P(zDM) methodology and results
Applying Bayes’ theorem to Eq. (16): (17)
where the PDF P(z_{s}DM_{exc}) describes the likelihood that an FRB with a measured time delay Δt_{obs,exc} between observing frequencies (ν_{1,obs} and ν_{2,obs}) on Earth due to a DM excess DM_{exc} will have a source redshift z_{s} according to the IGM and host models comprising P(DM_{exc} z_{s}). For practical application the model should be normalised to ∫ P(z_{s} DM_{exc}) dz = 1. In Eq. (17), P(DM_{exc}) is the probability of observing a DM excess and P(z_{s}) is the assumed prior FRB redshift distribution.
Figure 5 shows projections of P(z_{s}DM_{exc}) for FRBs of different DM_{exc} values assuming a progenitor population following a YPSR spacial distribution in spiral galaxies. The left subplot compares projections against P(z_{s}DM_{IGM}), illustrating the potential influence of the host galaxy on probable redshifts of an FRB. The comparisons illustrate the diminishing influence of the host with increasing DM_{exc} owing to the increased likelihood of a highredshift source. This result has implications for constraining redshifts using their DMs alone. Using such techniques, FRBs with larger DM_{exc} values are more accurately constrained. We note that the underlying progenitors of FRBs can significantly alter the redshift range available to an FRB (see Sect. 4.3), so the most accurate constraints result from use of realistic host and IGM models.
The assumed prior FRB redshift distribution P(z_{s}) may also influence the potential redshift ranges available to FRBs. The uniform prior (U(0,6)) considered across the generated redshift range in Fig. 5 (left) assumes that FRBs are equally distributed with redshift. Redshift distributions of FRBs rely on their progenitors, so this may not be the case. Zhou et al. (2014) previously considered the redshift distribution of FRBs to follow an Erlang distribution fit, i.e. (18)
where k = 2, γ = 1 which has been associated with GRBs by Shao et al. 2011 and accounts for a dearth of observed events in the nearby Universe (presumably owing to the relative scarcity of galaxies), and the decreased likelihood of detection of distant events by sensitivitylimited instruments. Figure 5 (right) compares P(z_{s}DM_{exc}) projections assuming uniform and Erlang priors for an FRB of DM_{exc} equal to that of FRB 121102. The Erlang prior suppresses the likelihood of its source originating at low redshifts. Both distributions are still consistent with spectroscopic measurements of FRB 121102’s host galaxy.
4 Discussion and conclusions
In previous sections we have reviewed the contributing components to observed FRB DMs, their redshift dependences, and have combined physically motivated PDFs for these components. Excess electron models resulting from such techniques allow more accurate constraints to be placed on unlocalised FRB redshifts than arise from standard practice. This is beneficial as an FRB’s observed DM will always be welldefined regardless of whether or not a spectroscopic redshift may be obtained. The threefold purposes of this work are as follows. Firstly, we seek to improve FRB DMredshift analysis at z_{s} > 2, by considering the redshiftevolving mean in DM_{IGM} and its associated variance due to galactic haloes along sightlines, the contribution of the host galaxy to DM_{exc} assuming different progenitor scenarios, and the effects of underlying prior FRB redshift distributions on their potential redshift ranges. Furthermore, we aim to illustrate the differences that physicallymotivated models have on FRB DMredshift probability distributions. Secondly, we place constraints on the redshifts of unlocalised FRBs using their DMs given the above considerations. Finally, we identify potential uses for DMredshift analysis in future FRB research.
Fig. 4 Excess electron models P(DM_{exc}z_{s}) illustrating the relationship between observed excess dispersion measure DM_{exc} and source redshift z_{s} for simulated FRB populations. Subplots consider progenitor scenarios following: (a) YPSR distributions in spiral galaxies, homogeneous distributions within galactocentric spheres of radius; (b) H = 0.2 kpc; (c) H = 0.004 kpc in spiral galaxies; elliptical distributions confined to regions (d) 1; (e) 0.1 times the effective radius (R_{e} = 3.2 kpc) of elliptical galaxies; and (f) homogeneous distributions within galactocentric spheres of radius H = 0.1 kpc in elliptical galaxies. Projections of each model for different z_{s} (solid lines) are compared with P(DM_{IGM}z_{s}) (dotted lines) to illustrate the influence of the host contribution. The twodimensional PDFs themselves (z_{s} against DM_{exc} with colour indicating relative likelihood) are inset into each subplot. Models are normalised to ensure ∫ P(DM_{exc}z_{s}) dDM = 1. 
4.1 Dispersion measure component analysis
For all FRB progenitor scenarios we consider, our P(DM_{host}z_{s} = 0) curves indicate potentially significant restframe dispersion (DM_{host} > 400 pc cm^{−3}, and in some cases >700 pc cm^{−3}) by the host galaxy (see Fig. 3). This result is consistent with previous spiral galaxy dispersion studies (e.g. Zhou et al. 2014; Xu & Han 2015). Our simulations also predict potentially nonnegligible dispersion by elliptical galaxies, which conflicts with past literature (e.g. Xu & Han 2015). This discrepancy may be a result of different modelling methods. Xu & Han (2015) scale a Milky Waysized elliptical galaxy (constructed by neglecting NE2001’s thin disc and spiral arm components) using H_{α} to trace ionised electrons, whereas we use a β model and confine sources to areas within the effective radius of the galaxy. We also note that our elliptical models are constructed to approach the local IGM electron density at their respective virial radii. However, excess matter past this point could in principle contribute to DM_{host}. Therefore future work could utilise an improved elliptical model obtained by theoretically deriving the electron number density profile of Eq. (8) and deducing the true radius at which it equals the cosmic mean.
Figure 4 indicates that the relative host galaxy contribution to DM_{exc} is most influential at low z_{s} and diminishes as z_{s} increases. This effect, due to cosmological time dilation and host galaxy recession (Ioka 2003), is independent of progenitor/host model and has previously been noted to inadequately suppress host galaxy contributions > 700 pc cm^{−3} for FRBs of z_{s} < 2 (Zhou et al. 2014; Macquart et al. 2015). Our models show that in some cases progenitor/host scenarios may significantly influence an FRB’s potential DM_{exc} range out toz_{s} > 3 (e.g. when sources are confined to regions of high electron density).
Our chosen models may not accurately represent the origins of individual FRBs. FRB 121102, the only spectroscopically constrained source to date, originates in a dwarf galaxy (Tendulkar et al. 2017), which is a scenario we have not considered. During the analysis we hold our galaxy masses fixed, which could be improved by considering realistic host galaxy mass distributions (see e.g. Luo et al. 2018). As our assumed masses lie close to the top of the mass range, such considerations may further suppress the overall host contribution to P(DM_{exc}z_{s}) relative to that of the IGM. Future DMredshift analysis would also ideally incorporate a DM_{local} component. Finally, we note that we idealistically assume perfectly removed DM_{MW} components in this paper. Uncertainties in Milky Way disc and halo electron density distribution models must be considered for the most accurate DMredshift analysis of FRBs.
Fig. 5 P(z_{s}DM_{exc}) projections after application of Bayes’ theorem, illustrating potential source redshift z_{s} ranges for FRBs with observed excess electron components DM_{exc}. The model chosen assumes FRB progenitors to follow YPSR distributions in spiral galaxies. Left subplot: DM_{exc} projections (solid lines) against P(z_{s}DM_{IGM}) (dashed lines), illustrating the influence of the host contribution. A uniform (U(0,6) prior redshift distribution is assumed. Right subplot: projections for an FRB of DM_{exc} equal to thatof FRB 121102. Uniform (pink) and Erlang (P(z_{s}) = P(z_{s};k, γ)) prior redshift distributions are compared. Vertical lines indicating FRB 121102’s redshift as predicted conventionally (dashdotted) and as observed spectroscopically (dashed) are included. Models are normalised to ensure ∫ P(z_{s} DM) dz = 1. 
4.2 Redshift constraints for unlocalised FRBs
Table B.1 presents revised source redshifts for all catalogued FRBs to date, constrained to bounds drawn from their 95% confidence intervals. The table includes progenitor scenarios chosen to highlight the impact of the host/progenitor on an FRB’s potential redshift range. Our chosen progenitor scenarios are as follows:
 a)
Distributions neglecting a DM_{host} component. For FRBs of z_{s} < 2, these distributions essentially recover conventionally calculated redshifts.
 b)
Sources tracing young pulsar populations in Milky Waylike spiral galaxies.
 c)
Sources confined to galactocentric spheres of radius H = 0.2 kpc in Milky Waylike spiral galaxies.
 d)
Sources tracing the gas profiles of 5 × 10^{10} M_{⊙} elliptical galaxies, confined to areas within their R_{e} = 3.2 kpc effective radii.
 e)
Sources confined to galactocentric spheres of radius H = 0.1 kpc in 5 × 10^{10} M_{⊙} elliptical galaxies.
In all cases a U(0,6) prior FRB redshift distribution has been assumed. The true cosmological FRB redshift distribution depends on FRB phenomenology.
We provide our models and a Jupyter notebook with this work so that our results may be reproduced, and so that others may obtain redshifts for FRBs with our models or substitute their own P(DM_{IGM}z_{s}), P(DM_{host}z_{s} = z_{0}), and P(z_{s}) distributions. We propose that such techniques be considered for general use during future FRB analysis to better constrain the redshifts of unlocalised FRBs.
4.3 Potential uses for DMredshift analysis
In this section, we identify potential implications of the DMredshift relationship and uses for the framework presented in this paper.
Progenitor distinction
Figure 4 illustrates the potentially significant impact of FRB host galaxies on their possible DM_{exc} values. Therefore observed DM_{exc} distributions for FRBs constrained to z_{s} ranges may be useful for ruling out progenitor scenarios. We illustrate the approach for the set of models described in Sect. 4.2. We consecutively draw DM_{exc} samples of increasing size from a model in the set held to be true and perform 10 000 Bayesian model selection tests against the other models. We define the threshold for unambiguous distinction of a “true” model from all others to be when 99.7% of tests result in “decisive” evidence (i.e. log_{10}(k) > 2, where k is the Bayes factor) in favour of the true model. We find that unambiguous distinction of any true model from the remainder of the set requires between 200−300, 300−400, or 700−800 FRBs with accurate DM_{exc} values and source redshifts between 0.45 < z_{s} < 0.65, 0.95 < z_{s} < 1.05, or 1.45 < z_{s} < 1.65, respectively. To distinguish between models of sources tracing stellar distributions (e.g. model b) and sources confined to regions of high electron density (e.g. models c, e), fewer than 50 FRBs are required for any of these redshift ranges.
Fig. 6 P(z_{s}DM_{exc}) projections after application of Bayes’ theorem for an FRB of DM_{exc} equal to that of FRB 171020. Subplots consider progenitor scenarios following: (a) distributions neglecting a host galaxy DM component, (b) YPSR distributions in spiral galaxies, (c) homogeneous distributions within galactocentric spheres of radius H = 0.2 kpc in spiral galaxies, (d) elliptical distributions confined to the effective radii (R_{e} = 3.2 kpc) of elliptical galaxies, and (e) homogeneous distributions within galactocentric spheres of radius H = 0.1 kpc in elliptical galaxies. Uniform (pink lines) and Erlang (grey lines) prior redshift distributions are compared. Vertical lines indicate the conventionally calculated redshift of the FRB. Shaded regions indicate bounds drawn from 68 to 95% confidence intervals. Models are normalised to unity. 
Host galaxy rejection
An observed FRB of unconstrained z_{s} always has a welldefined DM_{obs}. After subtracting an accurate DM_{MW}, progenitor model and prior redshift distributiondependent P(z_{s}DM_{exc}) curves may be used to constrain z_{s} to a statistical significance. Galaxies within the FRB localisation region, but outside the selected redshift range may then be discounted as potential hosts. We implement this technique for the lowDM_{obs} (= 114 pc cm^{−3}) FRB 171020, motivated by the work of Mahony et al. (2018).
Figure 6 shows P(z_{s}DM_{exc}) projections for the burst assuming the progenitor scenarios in Sect. 4.2. Table B.1 contains numerical constraints for their 95% confidence regions. The potential available z_{s} range is shown to be significantly influenced by both progenitor model and P(z_{s}). Models (a)−(d) are broadly consistent with Mahony et al. (2018) calculations^{7} but with larger upper limits (0.098 < z_{s, max} < 0.143). Model e differs significantly from these ranges (1.325 < z_{s} < 1.432). This result may be understood by considering the large z_{s} necessary for relativistic effects to suppress sufficiently the large DM_{host} acquired from an electrondense environment. A search of the NASA Extragalactic Database (NED) within the localisation region of the FRB using our revised z_{s} ranges yields 2148 objects, 13 of which have spectroscopically constrained redshifts. Of these, 11 fall outside our ranges for models (a)−(d). The remaining two galaxies are discussed in Mahony et al. (2018). If FRB 171020 were to have originated from a model (e)like progenitor, these two galaxies could also be discounted as potential hosts.
Milky Way uncertainty mitigation
By observing small areas of sky, surveys may mitigate uncertainty as a consequence of DM_{MW} sightline variations. Over a small enough area, the Milky Way contribution to DM_{obs} could be considered uniform, resulting in a constant systematic offset in the DM_{exc} distribution for all FRBs observed in any given direction. Current and future FRB searches capable of piggybacking deep extragalactic surveys (e.g. VFASTR, MEETRAP, LOFTe; Wayth et al. 2011; BurkeSpolaor et al. 2016; Stappers 2016; Walker et al. 2018)could potentially provide such sets of FRBs. Assuming FRBs have no preferential latitude dependence (Bhandari et al. 2018), this approach would yield similar numbers of events in any pointing direction.
DM_{MW} analysis
By considering the above approach, sufficient numbers of pointings by deep extragalactic FRB searches over large sky areas collecting sufficiently large numbers of FRBs would allow the disentanglement of DM_{exc} and DM_{MW}. Assuming the true underlying P(DM_{exc}z_{s}) to be directionindependent, the DM_{obs} distribution for FRBs of given z_{s} observed in different directions could be compared and searched for systematic offsets.
We demonstrate this technique by analysing how well an underlying DM_{MW} distribution may be recovered for a simulated FRB population. To do so, we generated an allsky FRB population, with DM_{MW} and DM_{exc} components drawn from the NE2001 model and a generic normal distribution, respectively. The latter distribution was chosen for simplicity, as the true distribution is irrelevant in this situation (see below). We then divided the sky into regions small enough to have constant DM_{MW} using a HEALPix tessellation; we experimented with NSide = 4, 8, and 16. Assuming the P(DM_{obs}z_{s}) distribution to be the same in all directions except for positiondependent offset due to the DM_{MW} contribution, we can recover the relative distribution of this Milky Way component across the sky by calculating the average DM within each tile. An intrinsic uncertainty in the overall zero point of DM_{MW} exists owing to the unknown mean DM of P(DM_{exc}z_{s}). However, the relative reconstruction of DM_{MW} is independent of the extragalactic probability distribution, thus justifying our earlier choice of a simple normal distribution. Our simulations show that we can recover DM_{MW} to ~37 (nFRB/tile)^{−1/2} pc cm^{−3} for any skytile tile ≲200 sq. deg., provided that it is far enough from the Galactic plane (i.e. ~ 5°) for the DM gradient within it to be fairly small.
The same reconstruction technique can also be applied to fluxlimited surveys (e.g. CRAFT, CHIME/FRB; Macquart et al. 2010; CHIME/FRB Collaboration 2018) since P(DM_{exc}F > F_{lim}) is also positionindependent. For a representative allsky survey containing 10 000 FRBs, it would be possible to reconstruct the sky DM_{MW} at a precision ~4.6 pc cm^{−3}, i.e. ~10% relative uncertainty, at a 215 sq. deg. resolution. Interestingly, it would then be possible to remove the MW contribution from all FRBs and combine these in order to obtain a proper P(DM_{exc}F > F_{lim}) that is accurate at the ~2% level as far as recovering the second moment is concerned.
As cosmological analysis of future redshiftconstrained FRB populations relies on precise removal of DM_{MW}, this application is a potentially crucial case for obtaining large numbers of FRBs with a limited localisation accuracy.
DM_{IGM}–DM_{host} disentanglement
As DM_{IGM}(z_{s}) is a function integrated over redshift while DM_{host}(z_{s}) is evaluated at a specific source redshift, the two components may theoretically be disentangled using redshiftconstrained FRBs. Techniques adapted by analogy to frequencydependent component separation of Galactic foreground and cosmic microwave background (CMB) emission (see, e.g. Planck Collaboration IX 2016) may be of use. However, even though the redshift dependence of these components is well known, their true underlying PDFs is not readily extractable. Thus component separation may require forward modelling using realistic template P(DM_{IGM}z_{s}) and P(DM_{host}z_{s}) distributions or a decomposition into basis functions such as a Gaussian mixture.
4.4 Concluding remarks
We present in this paper a statistical framework for the exploration of the FRB DMredshift relationship, which provides the basis for several analyses. Firstly, it facilitates the assessment of host galaxy contributions to FRB DMs using physically motivated models. We find that all of our host models may potentially contribute restframe DM_{host} values > 400 pc cm^{−3}. By applying relativistic effects in an expanding Universe, we demonstrate the diminishing influence of DM_{host} with increasing redshift. But for extreme progenitor scenarios, in which FRBs originate in electrondense regions (e.g. near to galactic centres), the DM_{host} component may still significantly influence P(DM_{exc}z_{s}) profiles out to z_{s} > 3. We place limits on the number of semiredshiftconstrained FRBs required to distinguish between our models, finding that < 50 events may be required to identify galactocentric from stellar progenitor populations. Secondly, the framework may be used for constraining the redshifts of unlocalised FRBs more rigorously than by using conventional techniques. By consulting P(z_{s}DM) probability distributions generated using physically motivated progenitor models, we provide constraints for all FRBs catalogued to date. We use these constraints to reject host galaxy candidates within the localisation region of the lowDM_{obs} FRB 171020. A repository containing Python code and examples of the DMredshift analysis demonstrated in this work is provided online^{8}. Finally, such techniques may aid DM_{MW} uncertainty mitigation and DM_{obs} component disentanglement. Assuming a directionindependent P(DM_{exc}z_{s}), the DM_{MW} contribution to DM_{obs} could be identified with a search for sightlinevariable offsets in histogrammed DM_{exc} distributions obtained by fluxlimited surveys (e.g. CHIME/FRB). Surveys yielding many FRBs within small enough redshift ranges could potentially put FRBs without spectroscopic redshift measurements to use with this technique. It may also be possible to separate DM_{IGM} and DM_{host} using their respective redshift dependences.
Acknowledgements
We thank Francesco Pace for vital conversation and instruction during this project, Matthew McQuinn for providing useful IGM comparison models, and Joseph Callingham, Jason Hessels, Cherry Ng and the anonymous referee for their suggestions for improving this paper. C.R.H.W. acknowledges support from a UK Science and Technology Facilities Council studentship. R.P.B. acknowledges support from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 715051; Spiders). Y.Z.M. is supported by the National Research Foundation of South Africa with Grant no. 105925 and no. 104800.
Appendix A Variance in DM_{IGM}
For our PDF of likelihood of a FRB of source redshift z_{s} to contribute DM_{IGM} to the overall observed dispersion measure DM_{obs}, we assume a Gaussian (A.1)
where the n_{e}(z) is the threedimensional electron number density at redshift z, where the variance is (A.3)
following McQuinn (2014). In this equation, χ is the comoving distance factoring out Universal expansion, d χ = c dz∕H(z), where c is the speed of light and H(z) is the Hubble factor at z; is the mean electron density at z = 0; the matter power spectrum at z for wavenumber k is P_{e}(k, z); and k_{⊥} is the perpendicular component (Dodelson 2003).
Following Eqs. (3) and (4), and Appendix (A.1)–(A.7) in Ma & Zhao (2014), n_{e}(z) may be written as (A.4)
where Ω_{b} = 0.048 is the fractional baryon density (Planck Collaboration VI 2020), m_{p} is the proton mass, μ_{e} ≃ 1.14 is the mean mass per electron, is the gas density at redshift z, where ρ_{cr} = 1.879 h^{2} × 10^{−29} g cm^{−3} is the Universe’s present critical density (Dodelson 2003), and (A.5)
where χ_{e} is the number ratio of ionised to total electrons, Y_{p} ≃ 0.24 is the helium mass fraction, and is the number of ionised electrons per helium atom. The quantity may range between 0–2 therefore χ_{e} may range from 0.86–1. In this work we assume χ_{e} = 1.
Eq. (A.2) may be written more explicitly as (A.7)
Substituting relevant quantities (A.8)
where (Mo & White 2002).
Likewise, following McQuinn (2014) the sightline variance in DM_{IGM} is (A.9)
where by Eq. (A.4) (A.10)
Appendix B Redshift constraints for existing FRBs
DMderived source redshifts z_{s} for catalogued FRBs (Col. 1) calculated via conventional techniques (z_{s} ~DM_{exc}∕(1200 pc cm^{−3}) Petroff et al. 2016; Ioka 2003); (Col. 2), and from excess electron models described in Sect. 4.2 with bounds drawn from their 95% confidence intervals (Cols. 3−7).
References
 Bannister, K. W., Shannon, R. M., Macquart, J.P., et al. 2017, ApJ, 841, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Bassa, C. G., Tendulkar, S. P., Adams, E. A. K., et al. 2017, ApJ, 843, L8 [NASA ADS] [CrossRef] [Google Scholar]
 Bhandari, S., Keane, E. F., Barr, E. D., et al. 2018, MNRAS, 475, 1427 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, B. A., & Bregman, J. N. 2001, ApJ, 547, 154 [NASA ADS] [CrossRef] [Google Scholar]
 BurkeSpolaor, S., Trott, C. M., Brisken, W. F., et al. 2016, ApJ, 826, 223 [NASA ADS] [CrossRef] [Google Scholar]
 Caleb, M., Flynn, C., Bailes, M., et al. 2017, MNRAS, 468, 3746 [NASA ADS] [CrossRef] [Google Scholar]
 Chatterjee, S., Law, C. J., Wharton, R. S., et al. 2017, Nature, 541, 58 [NASA ADS] [CrossRef] [Google Scholar]
 CHIME/FRB Collaboration (Amiri, M., et al.) 2018, ApJ, 863, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Cordes, J. M., & Lazio, T. J. W. 2002, ArXiv eprints [astroph/0207156] [Google Scholar]
 Cordes, J. M., & Wasserman, I. 2016, MNRAS, 457, 232 [NASA ADS] [CrossRef] [Google Scholar]
 de Graaff, A., Cai, Y.C., Heymans, C., & Peacock, J. A. 2019, A&A, 624, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dodelson, S. 2003, Modern Cosmology (Cambridge, US: Academic Press) [Google Scholar]
 Falcke, H., & Rezzolla, L. 2014, A&A, 562, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gómez, G. C., Benjamin, R. A., & Cox, D. P. 2001, AJ, 122, 908 [NASA ADS] [Google Scholar]
 Ioka, K. 2003, ApJ, 598, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Kalberla, P. M. W., & Dedes, L. 2008, A&A, 487, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Keane, E. F. 2016, MNRAS, 459, 1360 [NASA ADS] [CrossRef] [Google Scholar]
 Keane, E. F. 2018, Nat. Astron., 2, 865 [NASA ADS] [CrossRef] [Google Scholar]
 Lorimer, D. R. 2013, IAU Symp., 291, 237 [NASA ADS] [Google Scholar]
 Lorimer, D. R., Bailes, M., McLaughlin, M. A., Narkevic, D. J., & Crawford, F. 2007, Science, 318, 777 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Luo, R., Lee, K., Lorimer, D. R., & Zhang, B. 2018, MNRAS, 481, 2320 [NASA ADS] [CrossRef] [Google Scholar]
 Ma, Y.Z., & Zhao, G.B. 2014, Phys. Lett. B, 735, 402 [NASA ADS] [CrossRef] [Google Scholar]
 Macquart, J.P., Bailes, M., Bhat, N. D. R., et al. 2010, PASA, 27, 272 [NASA ADS] [CrossRef] [Google Scholar]
 Macquart, J. P., Keane, E., Grainge, K., et al. 2015, Advancing Astrophysics with the Square Kilometre Array (Berlin: Springer), 55 [CrossRef] [Google Scholar]
 Mahony, E. K., Ekers, R. D., Macquart, J.P., et al. 2018, ApJ, 867, L10 [NASA ADS] [CrossRef] [Google Scholar]
 MaízApellániz, J. 2001, AJ, 121, 2737 [Google Scholar]
 Mamon, G. A., & Łokas, E. L. 2005, MNRAS, 363, 705 [NASA ADS] [CrossRef] [Google Scholar]
 Manchester, R. N., & Taylor, J. H. 1972, Astrophys. Lett., 10, 67 [NASA ADS] [Google Scholar]
 Marcote, B., Paragi, Z., Hessels, J. W. T., et al. 2017, ApJ, 834, L8 [NASA ADS] [CrossRef] [Google Scholar]
 Masui, K., Lin, H.H., Sievers, J., et al. 2015, Nature, 528, 523 [NASA ADS] [CrossRef] [Google Scholar]
 McQuinn, M. 2014, ApJ, 780, L33 [NASA ADS] [CrossRef] [Google Scholar]
 Michilli, D., Seymour, A., Hessels, J. W. T., et al. 2018, Nature, 553, 182 [NASA ADS] [CrossRef] [Google Scholar]
 Mo, H. J., & White, S. D. M. 2002, MNRAS, 336, 112 [NASA ADS] [CrossRef] [Google Scholar]
 Mottez, F., & Zarka, P. 2014, A&A, 569, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ng, C., Vanderlinde, K., Paradise, A., et al. 2017, ArXiv eprints, [arXiv:1702.04728] [Google Scholar]
 Nicastro, F., Kaastra, J., Krongold, Y., et al. 2018, Nature, 558, 406 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Pen, U.L., & Connor, L. 2015, ApJ, 807, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Petroff, E., Barr, E. D., Jameson, A., et al. 2016, PASA, 33, e045 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration IX 2016, A&A, 594, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VI 2020, A&A, in press, https://doi.org/10.1051/00046361/201833910 [Google Scholar]
 Romero,G. E., del Valle, M. V., & Vieyro, F. L. 2016, Phys. Rev. D, 93, 023001 [NASA ADS] [CrossRef] [Google Scholar]
 Schnitzeler, D. H. F. M. 2012, MNRAS, 427, 664 [NASA ADS] [CrossRef] [Google Scholar]
 Scoville, N., Aussel, H., Brusa, M., et al. 2007, ApJS, 172, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Shao, L., Dai, Z.G., Fan, Y.Z., et al. 2011, ApJ, 738, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Spitler, L. G., Cordes, J. M., Hessels, J. W. T., et al. 2014, ApJ, 790, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Spitler, L. G., Scholz, P., Hessels, J. W. T., et al. 2016, Nature, 531, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Stappers, B. 2016, in Proceedings of MeerKAT Science: On the Pathway to the SKA. 2527 May, 2016 Stellenbosch, South Africa (MeerKAT2016), Online at https://pos.sissa.it/cgibin/reader/conf.cgi?confid=277, id.10, 10 [Google Scholar]
 Sun, X. H., & Han, J. L. 2004, MNRAS, 350, 232 [NASA ADS] [CrossRef] [Google Scholar]
 Tanimura, H., Hinshaw, G., McCarthy, I. G., et al. 2019, MNRAS, 483, 223 [NASA ADS] [CrossRef] [Google Scholar]
 Tendulkar, S. P., Bassa, C. G., Cordes, J. M., et al. 2017, ApJ, 834, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Thornton, D., Stappers, B., Bailes, M., et al. 2013, Science, 341, 53 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Totani, T. 2013, PASJ, 65, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Vieyro, F. L., Romero, G. E., BoschRamon, V., Marcote, B., & del Valle, M. V. 2017, A&A, 602, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Walker, C. R. H., Breton, R. P., Harrison, P. A., et al. 2018, IAU Symp., 337, 422 [NASA ADS] [Google Scholar]
 Wayth, R. B., Brisken, W. F., Deller, A. T., et al. 2011, ApJ, 735, 97 [NASA ADS] [CrossRef] [Google Scholar]
 Xu, J., & Han, J. L. 2015, Res. Astron. Astrophys., 15, 1629 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, Y.P., Luo, R., Li, Z., & Zhang, B. 2017, ApJ, 839, L25 [NASA ADS] [CrossRef] [Google Scholar]
 Yusifov, I., & Küçük, I. 2004, A&A, 422, 545 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zhang, B. 2018, ApJ, 867, L21 [NASA ADS] [CrossRef] [Google Scholar]
 Zhou, B., Li, X., Wang, T., Fan, Y.Z., & Wei, D.M. 2014, Phys. Rev. D, 89, 107303 [NASA ADS] [CrossRef] [Google Scholar]
The true ratio of ionised to total electrons may be between 0.86 and 1 (see Appendix A).
See for example Luo et al. (2018) for an approach to galaxy simulation accounting for a mass ranges and evolution with redshift.
FRB 110523’s observed magnetisation and scattering properties may also indicate a similar environment (Masui et al. 2015).
All Tables
DMderived source redshifts z_{s} for catalogued FRBs (Col. 1) calculated via conventional techniques (z_{s} ~DM_{exc}∕(1200 pc cm^{−3}) Petroff et al. 2016; Ioka 2003); (Col. 2), and from excess electron models described in Sect. 4.2 with bounds drawn from their 95% confidence intervals (Cols. 3−7).
All Figures
Fig. 1 Probability distribution functions (P(DM_{IGM})) for FRBs of source redshift z_{s} = 1 accounting for DM_{IGM} sightline variance due to collapsed systems. The Gaussian approximation assumed in this work is compared to more complex McQuinn (2014) baryonic halo models (tracing top hat functions, NFW profiles, and cosmological simulations, respectively, obtained via. private communication). Profiles are normalised to ∫ P(DM z_{s} = 1) d DM = 1000. 

In the text 
Fig. 2 Radial and height cross sections of our spiral galaxy electron density model n_{e,spiral} (r, h) (pink lines). Constituent components of the model tracing the inner MW Galactic disc (Gómez et al. 2001; Schnitzeler 2012) (blue lines) and an ionised outer HI disc model (Kalberla & Dedes 2008) are shown (brown lines). Shaded (blue/brown) regions highlight areas in which the inner/outer profiles are considered well constrained. The horizontal black dotted line indicates the approximate electron density for the IGM at z_{s} = 0 (Ioka 2003). 

In the text 
Fig. 3 Restframe PDFs of relative likelihood against DM_{host} for multiple FRB progenitor scenarios. From top left to bottom right: sources follow stellar distributions (OB stars, YPSRs, OPSRs, MSPs) in spiral galaxies; homogeneous distributions (within spheres of radius H kpc) in spiral galaxies; elliptical distributions (confined to multiples of the galaxy effective radius, R_{e} = 3.2 kpc) in elliptical galaxies; and homogeneous distributions (within spheres of radius H kpc) in elliptical galaxies. 

In the text 
Fig. 4 Excess electron models P(DM_{exc}z_{s}) illustrating the relationship between observed excess dispersion measure DM_{exc} and source redshift z_{s} for simulated FRB populations. Subplots consider progenitor scenarios following: (a) YPSR distributions in spiral galaxies, homogeneous distributions within galactocentric spheres of radius; (b) H = 0.2 kpc; (c) H = 0.004 kpc in spiral galaxies; elliptical distributions confined to regions (d) 1; (e) 0.1 times the effective radius (R_{e} = 3.2 kpc) of elliptical galaxies; and (f) homogeneous distributions within galactocentric spheres of radius H = 0.1 kpc in elliptical galaxies. Projections of each model for different z_{s} (solid lines) are compared with P(DM_{IGM}z_{s}) (dotted lines) to illustrate the influence of the host contribution. The twodimensional PDFs themselves (z_{s} against DM_{exc} with colour indicating relative likelihood) are inset into each subplot. Models are normalised to ensure ∫ P(DM_{exc}z_{s}) dDM = 1. 

In the text 
Fig. 5 P(z_{s}DM_{exc}) projections after application of Bayes’ theorem, illustrating potential source redshift z_{s} ranges for FRBs with observed excess electron components DM_{exc}. The model chosen assumes FRB progenitors to follow YPSR distributions in spiral galaxies. Left subplot: DM_{exc} projections (solid lines) against P(z_{s}DM_{IGM}) (dashed lines), illustrating the influence of the host contribution. A uniform (U(0,6) prior redshift distribution is assumed. Right subplot: projections for an FRB of DM_{exc} equal to thatof FRB 121102. Uniform (pink) and Erlang (P(z_{s}) = P(z_{s};k, γ)) prior redshift distributions are compared. Vertical lines indicating FRB 121102’s redshift as predicted conventionally (dashdotted) and as observed spectroscopically (dashed) are included. Models are normalised to ensure ∫ P(z_{s} DM) dz = 1. 

In the text 
Fig. 6 P(z_{s}DM_{exc}) projections after application of Bayes’ theorem for an FRB of DM_{exc} equal to that of FRB 171020. Subplots consider progenitor scenarios following: (a) distributions neglecting a host galaxy DM component, (b) YPSR distributions in spiral galaxies, (c) homogeneous distributions within galactocentric spheres of radius H = 0.2 kpc in spiral galaxies, (d) elliptical distributions confined to the effective radii (R_{e} = 3.2 kpc) of elliptical galaxies, and (e) homogeneous distributions within galactocentric spheres of radius H = 0.1 kpc in elliptical galaxies. Uniform (pink lines) and Erlang (grey lines) prior redshift distributions are compared. Vertical lines indicate the conventionally calculated redshift of the FRB. Shaded regions indicate bounds drawn from 68 to 95% confidence intervals. Models are normalised to unity. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.