CRAFTS is a multi-purpose drift scan survey conducted with FAST using a 19-beam receiver operating at 1.05–1.45 GHz, deployed in May 2018 and conducting blind FRB searches using multiple pipelines34. FRB 20190520B was discovered on 16 November 2019 in archived CRAFTS data, which are in 1-bit filterbank format with 196-μs time resolution and 0.122-MHz frequency resolution. In this first discovery observation, 3 bursts were detected in 10 s, and another burst was detected 20 s later. These 4 bursts from the drift scan survey gave a preliminary location for the source within an approximately 5-arcmin-diameter region. Taking the pointing location from FAST, 2 follow-up observations were performed with FAST on 25 April and 22 May in 2020 with 19-beam mode in which 15 bursts were detected. A monthly observation campaign was then conducted by FAST using the approximately 100-mas localization from the VLA (next section). After some regular telescope maintenance, 10 observations were performed spanning from 30 July 2020 to 19 September 2020 in which 60 more bursts were detected.

Bursts were detected from FRB 20190520B in each FAST monitoring observation. We list the properties of those bursts detected by FAST in Supplementary Table 1. The burst arrival time is in modified Julian date (MJD) format and has been transformed to the arrival time at the Solar System barycentre at 1.5 GHz with the DM values from Supplementary Table 1. The observed DMobs is measured using the method from ref. 9 and the code from ref. 35. We use the DM value that maximizes the structure of the highest signal-to-noise ratio (S/N) burst on a given date, and report this same DM for all other bursts detected on that date. Apparent epoch-to-epoch DM variations may simply be the result of a combination of variable effects, including intrinsic pulse structure, time–frequency drift, flux distribution in frequency and scattering (which affects bursts chromatically).

The FAST data were searched using a Heimdall-based pipeline34,36. For the FRB blind search in the 19-beam drift scan survey, the polarizations were added and only the total intensity (Stokes I) was recorded. The trial DM range was from 100 pc cm−3 to 5,000 pc cm−3 and we matched the pulse width using a boxcar search up to 200 ms. The candidates with S/N > 7 and present in fewer than 4 adjacent beams were manually examined for further inspection. After we identified FRB 20190520B as a new FRB source, the follow-up burst search was done with a narrower DM range (100–2,000 pc cm−3). Following the VLA localization, the tracking observations only recorded data from the central beam but with all the Stokes parameters and higher time resolution (about 50 μs). Candidates with S/N > 7 correspond to a fluence threshold of 9 mJy ms for a burst width of 1 ms.

The pulse widths were estimated by a Gaussian profile fit if the burst showed no evidence of scattering (see section on scattering below). A subpulse is recognized if the bridge between two closely spaced peaks drops more than 5σ below the higher peak. If the burst shows a scattering tail, the pulse width is derived from a combined fit for the Gaussian width and scattering time. We roughly estimate the bandwidth of each burst by dividing the whole bandpass into 50-MHz subbands and identifying the subbands containing burst emission. The burst fluences were determined using the bandwidth, temporal width and amplitude of each burst.


Following the FAST detection of bursts from FRB 20190520B, VLA observations were performed from July–November 2020 with director’s discretionary time (DDT) project 20A-557 at the most reliable position determined using FAST (good to 5 arcmin). Most of the observations were in the B array configuration (with a maximum baseline of 11.1 km), with the exceptions that the array configuration was BnA on MJD 59,161 and BnA->A on MJD 59,167 and MJD 59,169. In total, 11.4 h were spent on the source at multiple bands. The total bandwidths for the L (1.5 GHz), S (3 GHz) and C (5.5 GHz) bands were 1,024 MHz, 2,022 MHz and 2,022 MHz, respectively, with 1,024 channels, corresponding to channel bandwidths of approximately 1 MHz, 2 MHz and 2 MHz, respectively.

The details of the observations are given in Extended Data Table 1. The telescopes were pointed at the field centred at (RA, dec.)(J2000) = (16 h 02 min 01 s, −11° 17′ 28″). We used the ‘realfast’ search system at the VLA to search for bursts from FRB 20190520B in our VLA observations. Observations on MJD 59,169 contributed to the imaging at S band, but the ‘realfast’ system was not run on this day owing to a system error.

The ‘realfast’ search system has been described in detail in refs. 12,37, but here we discuss it briefly. Using a commensal correlator mode, the visibilites are sampled with 10-ms resolution and distributed to the ‘realfast’ graphics-processing-unit cluster to search for transients. The search pipeline rfpipe38 then applies the online calibration, de-disperses the visibilities and forms images at a set of different temporal widths. The 8σ fluence limit of a 10-ms image is 0.29 Jy ms, 0.18 Jy ms and 0.13 Jy ms for the L, S and C bands, respectively. Candidates with an image S/N greater than the threshold triggered the recording of fast sampled visibilities in windows ranging from 2 s to 5 s centred on each candidate. For each candidate, time–frequency data corresponding to the maximum pixel in the image are then extracted and post-processed. The candidate is then classified by fetch, a graphics-processing-unit-based convolutional neural network39.

Promising candidates selected visually by the ‘realfast’ team go through an offline analysis to refine the candidate parameters and to improve its detection significance, if possible. Several methods are tried within this analysis: offline search for the transient using a finer DM grid, varying radio frequency interference (RFI) flagging thresholds, changing image size (as the real-time search is done on non-optimal image sizes owing to computational limitations), subband search and so on (see section 2.4 of ref. 37 for more details). The pixel sizes of images formed by the real-time system for the L, S and C bands were 0.9″, 0.48″ and 0.27″, respectively, with an image size of 2,048 × 2,048 pixels corresponding to field sizes of 0.5°, 0.27° and 0.15°, respectively. During refinement, we searched the data with smaller pixel sizes: 0.5″ × 0.75″ for the L band, 0.38″ × 0.28″ for the S band and 0.27″ × 0.18″ for the C band. The significance of an astrophysical transient improves when the data is de-dispersed at a DM closer to the true DM of the candidate (see section 2 of ref. 40). Therefore, we re-ran the search with a finer DM grid at 0.1% fractional sensitivity loss, compared with 5% used by the real-time system. Noise-like events or RFI are sensitive to RFI flagging and image gridding parameters, and so they cannot be reproduced on refinement and are discarded. Sometimes, the transient signal is present only in a fraction of the whole frequency band, so we also re-run our search only using the relevant frequencies, which further improves the detection significance. We applied these techniques on all the candidates selected from these observations to refine their parameters.

Localization of bursts


The ‘realfast’ system makes several assumptions during calibration and imaging to improve the computational efficiency for the real-time imaging and search. Moreover, the point spread function (PSF) of the interferometer is not deconvolved from the image, so the real-time system forms ‘dirty’ images. The PSF shape makes the images more difficult to visually interpret and model.

To address these issues, we used a top-hat function to select the raw, de-dispersed, fast-sampled visibilities containing only the burst signal to re-image the burst data using the Common Astronomy Software Application package (CASA v5.6.2-3)41. Observations of 3C 286 (acquired before the FRB observations) were used to calibrate the flux density scale, bandpass and delays. Complex gain fluctuations over time were calibrated with observations of the calibrator J1558−1409. We performed phase-reference switching at intervals of 16 min, 13 min and 12 min for the L, S and C bands, respectively, consistent with the nominal phase-coherence timescales for the VLA, so that after calibration any short-timescale phase variations are negligible. Therefore, the systematic errors for the burst positions on short timescales are of the same magnitude as those for the deep imaging, as discussed below.

Determining properties of individual bursts

After calibration, the CASA task tclean was used to generate an image for each burst and estimate the S/N. Most of the bursts are spectrally confined, so for each burst we select the spectral window range that produces the highest image S/N. We then use the CASA task imfit around the FRB position to fit an elliptical Gaussian to the source in the radio image and measure the centroid location, peak flux density and 1σ image-plane uncertainties. Hereafter, we refer to these image-plane uncertainties as statistical errors on burst positions.

The average burst positions are calculated for each frequency band separately, by weighting each detection by the inverse of the position fit errors (statistical errors) reported by CASA. Statistical errors are inversely proportional to the S/N, and therefore high significance detection is expected to have smaller fit errors. The total positional error at each frequency band is obtained by adding this statistical error in quadrature with the systematic error in that band (determined using the deep radio images described below). The total positional error (RA error, dec. error) is (0.25″, 0.32″) for the L band, (0.28″, 0.17″) for the S band and (0.12″, 0.09″) for the C band. In the final analysis, the errors are dominated by systematic uncertainties at each frequency band (Extended Data Tables 2, 3). The final burst position was obtained by taking a weighted average of the burst positions at each of the three frequency bands, using the inverse of the total error at each frequency band as the weight. The best estimate of the burst position is RA = 16 h 02 min 04.272 s, dec. = −11° 17′ 17.32″ (J2000). We estimate the error on this position to be 0.10″ and 0.08″. We calculated the reduced chi-square value of individual burst positions with respect to this best estimate of the FRB position. The reduced chi-square was given by ({chi }^{2}=(1/8){sum }_{i}{({theta }_{i}-{theta }_{{rm{best}}})}^{2}/{sigma }_{i}^{2}), where ({theta }_{i}) refers to the coordinates of a burst, ({theta }_{{rm{best}}}) refers to the best estimate of the FRB position, ({sigma }_{i}) is the total error on a burst (statistical and systematic errors added in quadrature) and the sum is over all the 9 localized bursts. We obtained a reduced chi-square value of 0.53 and 0.77 for RA and dec. respectively.

We confirmed that there were no short-timescale phase errors by performing intermediate-timescale imaging of the continuum data. We imaged segments of 5–30 s, such that at least 1–3 sources could be detected at each frequency, and inspected the position variations of those sources over time. At all frequencies, the positions were stable and consistent with the statistical uncertainty (radiometer noise). Offsets of the sources were within the range of the systematics quoted for the PRS in the next section.

Owing to the low time resolution of the VLA data (10 ms), we do not perform any sophisticated modelling of the burst properties. The properties of the VLA bursts, including the flux densities, were obtained using CASA and the ‘realfast’ system, and are given in Extended Data Table 4. Here the reported widths should be considered as upper limits. The frequency extent of the burst signal was visually determined and reported in the last column. We also could not estimate a structuremaximizing DM for the VLA bursts because of the coarse time resolution. The DMs reported in this table are the values that maximize the S/N. This also explains the apparent variability seen in the DM values.

Persistent radio source

Deep radio images and PRS

The VLA campaign obtained two epochs at 1.5 GHz and six epochs at 3 GHz and 5.5 GHz, which resulted in an on-source time of about 3 h for 1.5 GHz, and about 4 h for both 3 GHz and 5.5 GHz.

Along with the ‘realfast’ output, the VLA visibilities with 3 s or 5 s sampling times were saved and analysed to search for persistent radio emission. This is done in parallel with saving the high-time-resolution (10 ms) data around the burst that was used in the previous section. We used the same data reduction, flagging and calibration approach that was used for burst localization. We then performed further flagging on the target and then subsequently imaged its Stokes I data using the CASA deconvolution algorithm tclean. To balance sensitivity while reducing sidelobes from a nearby bright source, we imaged with a Briggs weighing scheme (robust=0). In addition, self-calibrations were performed for all observations to correct considerable artefacts from the close-by bright sources in the field. We made use of the CASA task imfit to measure source flux densities by fitting an elliptical Gaussian model in the image plane.

We stacked observations at each central frequency in the uv plane and then imaged the Stokes I intensity, resulting in deep images at 1.5 GHz, 3 GHz and 5.5 GHz with root-mean-square (r.m.s.) noise of 9.0 μJy per beam, 4.5 μJy per beam and 3.0 μJy per beam, respectively. At 1.5 GHz, 3 GHz and 5.5 GHz, the upper limits of the deconvolved sizes of the PRS are as large as 1.4″ × 0.89″, 0.51″ × 0.14″ and 0.36″ × 0.1″, respectively. Thus, a conservative upper limit of the size of the PRS from VLA observations is 0.36″, which corresponds to 1.4 kpc at the angular diameter distance of 809 Mpc. The obtained positions for the PRS are shown in Extended Data Table 2 and the positions of the bursts and PRS are shown in Extended Data Table 1. The systematic offsets on these positions are estimated in the next section.

Systematic offsets

To determine the systematic errors on the coordinates of the PRS that we determined from the deep images, we ran the PyBDSF (https://www.astron.nl/citt/pybdsf/index.html) package to extract radio sources from the deep images. We then cross-matched the detected point sources in the deep images with the sources listed in the optical PanSTARRS survey DR142. We identified the radio point sources using the following criteria: (1) the peak intensity (Jy per beam) of a source should be 0.7, 0.5 and 0.5 times higher than its integrated flux (Jy) for the 1.5-GHz, 3-GHz and 5.5-GHz images, respectively; (2) the S/N (peak intensity/local r.m.s. noise) of a source should be greater than 5.

In total, we detected 375, 113 and 43 sources in the 1.5-GHz, 3-GHz and 5.5-GHz deep images, respectively. The selected sources were also visually checked to make sure that they are ‘point like’ sources in the deep images. We searched for matching optical sources within radii of 0.5″, 1.0″ and 2.0″. Going from 1.0″ to 2.0″, the additional cross-matched sources at each band are consistent with chance coincidence. Therefore, we adopt a 1.0″ cross-match radius, finding 136, 31 and 9 sources with optical counterparts in the PanSTARRS DR1 catalogue at 1.5 GHz, 3 GHz and 5.5 GHz, respectively. These cross-matched sources were used to determine the astrometry for FRB 20190520B and its PRS. By subtracting PanSTARRS coordinates from VLA coordinates, we estimate the systematic offsets in the 1.5-GHz, 3-GHz and 5.5-GHz positions, as listed in Extended Data Table 2. The systematic offsets are consistent with zero mean, and their uncertainties dominate the uncertainties of the PRS position.

Determining the position of the PRS

To determine the position of the PRS in the three frequency bands, we followed a procedure similar to what was used to determine the burst positions (see ‘Localization of bursts’). The best estimate for the PRS position is (RA, dec.)(J2000) = (16 h 02 min 04.261 s, −11° 17′ 17.35″). The error on this position is estimated to be 0.10″ and 0.05″ for RA and dec., respectively.

Variability and spectrum of the PRS

The flux density of the source measured at each epoch is shown in Extended Data Table 1. The measured flux densities show variations that are mostly consistent with measurement errors, but there are roughly 2σ variations in the 2.5-GHz subband that, if real, could be refractive interstellar scintillation or intrinsic variations, or both. To study the spectrum of the PRS, we split each of the observations into two 0.5 GHz/1 GHz subbands. Then we measured the average flux density at each of the subbands over the campaign. The multiband data were fit with a power-law model (({S}_{nu }propto {nu }^{alpha }), where Sv is the observed flux density at the frequency v, and α is the spectral index), yielding an average spectral index for the PRS of −0.41 ± 0.04 (Extended Data Fig. 2).

Chance coincidence association of the PRS

In the 1.5-GHz deep image, we detected 8 ‘point like’ sources, including the PRS, within 5 arcmin of the phase centre and with a flux density no less than that of the PRS (at 260 μJy) based on our point-source selection criteria described earlier. There is an additional bright source with a flux density of a few tens of millijansky in the region, but it was not classified as a ‘point like’ source based on our criteria.

To estimate the chance coincidence probability of the PRS with the bursts, we compared the solid angle corresponding to the uncertainty of burst localization and the average solid angle occupied by each of the eight ‘point like’ sources in the FRB field of view. The solid angle corresponding to each of the 8 ‘point like’ sources is roughly estimated as Ssource = π(5/60)2/8 steradians (Sr). The offset between the average position of the nine bursts and the position of the PRS at 5.5 GHz, which is best constrained when taking both statistical and systematic errors into account, is about 0.06 arcsec. This, along with a statistical error of 0.01 arcsec and a conservative estimate of the systematic error of 0.12 arcsec, can be used to estimate the offset between the PRS and the FRB position. We conservatively estimate this offset to be 0.19 arcsec. The solid angle corresponding to the offset therefore can be estimated as Soffset= π(0.19/60/60)2/8 Sr. The ratio between Soffset and Ssource gives the chance of coincident association of the PRS with the FRB position to be about 3 × 10−6.

Galaxy photometry and redshift determination

The deep R′-band (5,427–7,041 Å) images obtained by CFHT/MegaCam were stacked from archival observation data taken in 2014–2015 by the Canada France Hawaii Telescope (CFHT) archival pipeline MEGAPIPE, with a total of about 3.6 h on the field. The level 3 (flux calibrated) images were retrieved for our analysis.

NIR J-band images of the FRB–20190520B field were taken under a relatively poor seeing condition (about 1.3″) through a Subaru target-of-opportunity observation on 5 August 2020. A total of 1.4 h of observations were used for the final combined J-band image shown in Fig. 2. A J-band source of 22.07 ± 0.14 mag (AB) was detected at about 1″ southeast of the burst location, and is possibly stellar emission from the host galaxy. A faint source 2.5″ north has 22.87 ± 0.26 mag in the J band. Neither of these two sources were detected in the Ks-band image, with a 5σ limit of 21.74 mag (1.1 h).

An optical spectrum was obtained with the Double Spectrograph (DBSP) on the Palomar 200-inch telescope on 24 July 2020 using a 1″ slit width. This observation was executed before the CFHT/MegaCam archival data on FRB 20190520B field were found, and only PanSTARRS images were used for observation planning. The slit of DBSP was set to cover the PRS emission at RA = 16 h 02 min 04.27 s, dec. = −11° 17′ 17.5″ detected by the VLA in L-band on 22 July 2020. No clear optical counterpart was detected in any of the five band images of PanSTARRS from DR142. The slit was guided by the nearby M-star at RA = 16 h 02 min 04.48 s, dec. = −11° 17′ 19.1″ as reported in the PanSTARRS DR1 catalogue, with i = 20.4 mag, and 3.4″ due east of the PRS coordinates. The slit was set to a position angle of 108.5°, ensuring that both the PRS and the M-star fell within the slit. The observations with 2 × 900 s exposures were carried out under photometric sky conditions and subarcsecond seeing. The two-dimensional spectrum was generated using the image reduction and analysis facility (IRAF), including bias removal, flat-fielding and reduction of other instrumental effects. The one-dimensional spectrum was extracted from a 1.5″ window. The standard star BD+284211 was used for telluric correction and flux calibration. The DBSP one-dimensional spectrum is shown in Extended Data Fig. 3. The flux scale of the spectrum does not include the slit loss and registration error of PanSTARRS coordinates of the M-type star. The [O III] 5,007 Å line and the Hα line are both well detected (>5σ). The two emission lines are narrow, each with a full-width at half-maximum of about 10 Å. The redshift derived from these two spectral lines is z = 0.241 ± 0.001.

A follow-up Keck LRIS spectroscopic observation was carried out on 25 August 2020 under reasonable weather and seeing conditions (1.1″). The 1.5″ slit was set at a position angle of 160° to the extended optical emission seen in the MegaCam R′-band image around the FRB 20190520B location. A total exposure of 3,600 s was obtained. The emission lines Hα, Hβ, [O III] 4,859 Å, and [O III] 5,007 Å are well detected, indicating that the extended R′-band structure has the same redshift of z = 0.241. After the spectrum was corrected for Galactic foreground extinction, the line flux ratio between Hα and Hβ was used to estimate an extinction of AV = 0.80 mag, assuming a case-B line ratio of 2.86, yielding an Hα flux ({F}_{{rm{H}}{rm{alpha }}}=(42.0pm 0.5)times {10}^{-17}{{rm{ergs}}}^{-1}{{rm{cm}}}^{-2}.) The extinction-corrected Hα luminosity is ({L}_{{rm{H}}{rm{alpha }}}=7.4pm 0.2times {10}^{40}{{rm{ergs}}}^{-1}).

Chance association probability of the host galaxy

We use the approach of ref. 43 to estimate the chance coincidence probability of the association between the galaxy (J160204.31−111718.5) and FRB 20190520B. Assuming a uniform surface distribution of galaxies, the probability of chance coincidence follows a Poisson distribution; that is, (P=1-exp (-{n}_{i})), where ni is the mean number density of galaxies brighter than a specified R′-band magnitude in a circle of given radius, determined by the half light radius of the galaxy and the burst localization error region. This number density is estimated using the results from ref. 44. From previous discussion, we assume the localization error on the FRB position to be about 0.1″, and the size of the host galaxy to be 0.5″. The R′-band magnitude of the possible host is difficult to estimate as it is significantly affected by the emission lines. Thus, we conservatively estimate the R′-band magnitude to be (lesssim )23.3″ mag assuming a flat spectral energy distribution in vLv between the R′ band and the J band, where v and Lv are frequency and luminosity, respectively. The probability of a chance association between FRB 20190520B and J160204.31−111718.5 is estimated to be about 0.8%, independent of constraints from the association with a rare PRS or the observed FRB scattering. The next nearest galaxy is 6.5″ away from the FRB and has a chance coincidence probability of >20%.

FAST burst sample analysis

Repetition rate

FRB 20190520B was found to be regularly active with two or more bursts detected in each monitoring observation with FAST. Owing to the possible clustering behaviour of FRB emission, a Weibull distribution was used for the waiting time δ (the separation between adjacent detected bursts)45

$${mathscr{W}}=frac{k}{lambda }{left(frac{delta }{lambda }right)}^{k-1}{{rm{e}}}^{-{(delta /lambda )}^{k}},$$


where k is the shape parameter and λ is the scale parameter of the distribution. The Weibull distribution reduces to a Poisson distribution when k = 1. The mean waiting time, (E(delta )=lambda {Gamma }(1+1/k)), where Γ is the gamma function,  implies a burst rate (r=1/E(delta )). Extended Data Fig. 4 shows the inferred parameter distributions obtained with a Markov chain Monte Carlo. We find the burst rate of FRB 20190520B is (r={4.5}_{-1.5}^{+1.9}{{rm{h}}}^{-1}) with shape parameters (k={0.37}_{-0.04}^{+0.04}) when all 79 bursts above 9 mJy ms are used (left panel in Extended Data Fig. 4). Waiting times shorter than 1 s include substructure in individual bursts, so we also analyse only longer wait times δ ≥ 1 s, yielding (r={5.3}_{-1.0}^{+1.1},{{rm{h}}}^{-1}) with shape parameters (k={0.76}_{-0.08}^{+0.09}) (right panel in Extended Data Fig. 4).

Short- and long-timescale periodicity search

For a range of trial periods (P) and period derivatives ((dot{P})), the arrival times of the 79 FAST bursts were folded to calculate the phase of each burst for a given trial period. For the short-timescale periodicity search, the arrival times were folded at P between 1 ms and 1,000 s and (dot{P}) between 10−12 s s−1 and 1 s s−1. For the long-timescale periodicity search, the trial range of period P was from 2 d to 365 d and the range for period derivative (to account for spindown) was (dot{P}) from 10−12 d d−1 to 1 d d−1. The longest contiguous phase segment without bursts was calculated for each fold trial as a signature for possible periodicity.

For the short-timescale periodicity search, the bursts spread to a phase window larger than 60% of one period, which indicates no periodicity pattern in the trial range. For the long-timescale periodicity search, the longest phase segment without bursts showed two maxima about 0.6 near 67 d and 169 d. These two maxima are reproduced by folding the MJDs of all observation sessions, which indicates that the maxima reflect the sampling window rather than any burst periodicity. Thus, no long or short period of FRB 20190520B was detected.

Energy distribution

A 1-K equivalent noise calibration signal was injected before each observation session to obtain a high-quality flux density and energy calibration measurement for each detected burst. The injected noise calibration signal was used to scale the data into temperature units, yielding a nearly constant r.m.s. radiometer noise to within 6%. We then converted from units of kelvin to jansky using the zenith-angle-dependent gain curve, provided by the observatory through quasar measurements46. For most observations, the zenith angle is less than 20°, which corresponds to a stable gain of about 16 K Jy−1.

Assuming bursts have spectral index of about 0, we calculated their isotropic equivalent burst energies, E, following equation (9) of ref. 17:

$$E={10}^{39}{rm{erg}}times frac{4{rm{pi }}}{1+z}{left(frac{{D}_{{rm{L}}}}{{10}^{28}{rm{cm}}}right)}^{2}left(frac{{F}_{nu }}{{rm{Jy}},{rm{ms}}}right)left(frac{{nu }_{{rm{c}}}}{{rm{GHz}}}right),$$


where Fv = Sv × Weq is the specific fluence in units of erg cm−2 Hz−1 or Jy ms, Sv is the peak flux density, which has been calibrated with the noise level of the baseline, giving the flux measurement for each pulse at a central frequency of vc = 1.25 GHz, Weq is the equivalent burst duration, and the luminosity distance DL = 1,218 Mpc corresponds to a redshift z = 0.241 for the source of FRB 20190520B, using standard cosmological parameters47. The fluence-width distribution at 1.25 GHz for FRB 20190520B bursts can be seen in Extended Data Fig. 5. The histogram of burst energies (Extended Data Fig. 6) exhibits a bump that we fit with a log-normal function.

DM inventory analysis

The observed DM can be separated into four primary components (all in the observer’s frame):



where DMMW is the contribution from the Milky Way’s interstellar medium, DMhalo is the contribution from the Milky Way halo, DMhost the contribution from the host galaxy including its halo and any gas local to the FRB source, and DMIGM is the contribution from the IGM.

We use the NE2001 model16 to evaluate DMMW = 60 pc cm−3 (compared with 50 pc cm−3 from the YMW16 model48) as the mean of a uniform distribution with a generous ±40% width to conservatively represent the uncertainty in DMMW for the direction of FRB 20190520B. For the MW halo contribution, we use a uniform distribution from 25 pc cm−3 to 80 pc cm−3. Together, the Milky Way disk and halo components yield a total range from 61 pc cm−3 to 164 pc cm−3 with a mean of 113 pc cm−3 and r.m.s. uncertainty of 17 pc cm−3.

For the IGM, we use the Λ cold dark matter model (ΛCDM, where Λ is the cosmological constant) to calculate the mean DM contribution49

$${{rm{DM}}}_{{rm{IGM}}}(z)=frac{3c{H}_{0}{{Omega }}_{{rm{b}}}{f}_{{rm{IGM}}}}{8{rm{pi }}G{m}_{{rm{p}}}}{int }_{0}^{z}frac{chi (z)(1+z){rm{d}}z}{{[{{Omega }}_{{rm{m}}}{(1+z)}^{3}+{{Omega }}_{{Lambda }}]}^{frac{1}{2}}},$$


where the free electron number per baryon in the Universe is (chi (z)approx 7/8) (which is assumed to be constant for FRB redshifts), the normalized matter density is ({{Omega }}_{{rm{m}}}=0.315pm 0.007), the dark-energy fraction is ({{Omega }}_{{Lambda }}simeq 1-{{Omega }}_{{rm{m}}}), the baryonic fraction is ({{Omega }}_{{rm{b}}}={h}^{-2}times (0.02237pm 0.00015),,)the fraction of baryons in the IGM is fIGM, and the Hubble constant is ({H}_{0}equiv htimes 100,{rm{km}},{{rm{s}}}^{-1}=67.36pm 0.54,{rm{km}},{{rm{s}}}^{-1}{{rm{Mpc}}}^{-1}) (ref. 47). Using values for the speed of light c, the gravitational constant G and the proton mass mp, the resulting expression

$${{rm{DM}}}_{{rm{IGM}}}(z)approx 978,{rm{pc}},{{rm{cm}}}^{-3},{f}_{{rm{IGM}}}{int }_{0}^{z}frac{(1+z){rm{d}}z}{{[{{Omega }}_{{rm{m}}}{(1+z)}^{3}+{{Omega }}_{{Lambda }}]}^{frac{1}{2}}},$$


yields ({{rm{DM}}}_{{rm{IGM}}}(0.241)=248{f}_{{rm{IGM}}}{{rm{pccm}}}^{-3}).

Using ({{rm{DM}}}_{{rm{IGM}}}(z)) as the mean, we calculate the range of values for a given redshift using a log-normal distribution with variance ({sigma }_{{{rm{DM}}}_{{rm{IGM}}}}^{2}(z)={{rm{DM}}}_{{rm{IGM}}}(z){{rm{DM}}}_{{rm{c}}}), where ({{rm{DM}}}_{{rm{c}}}={{rm{50pccm}}}^{-3}) is chosen to provide cosmic variance consistent with published simulations. This gives ({sigma }_{{{rm{DM}}}_{{rm{IGM}}}}(0.241)=111sqrt{{f}_{{rm{IGM}}}}{{rm{pccm}}}^{-3}). The parameters for the log-normal distribution are then ({sigma }_{{{rm{lnDM}}}_{{rm{IGM}}}}={{mathrm{ln}[1+{({sigma }_{{{rm{DM}}}_{{rm{IGM}}}}(z)/{{rm{DM}}}_{{rm{IGM}}}(z))}^{2}]}}^{1/2}) and ({mu }_{{{rm{lnDM}}}_{{rm{IGM}}}}=,mathrm{ln},[{{rm{DM}}}_{{rm{IGM}}}(z)]-{sigma }_{{{rm{lnDM}}}_{{rm{IGM}}}}^{2},/2). It is noted that the log-normal distribution is asymmetric. Equation (5) and the log-normal distribution are also used to estimate the median DMIGM and its inner 68% uncertainty range shown in Fig. 3.

To constrain the host-galaxy DM, we combine the MW and IGM estimates and their uncertainties with the measured DM averaged over all bursts, which is ({{rm{DM}}}_{{rm{obs}}}=1,,204.7pm 4.{{rm{0pccm}}}^{-3}). The cumulative distribution function (CDF) of the log-normal distribution for the IGM contribution yields the range ({{rm{DM}}}_{{rm{IGM}}}={195}_{-70}^{+110}{{rm{pccm}}}^{-3}). We then calculate the posterior distribution of DMhost after marginalizing over the IGM and Milky Way distributions and using a flat, uninformative prior for DMhost. The median and probable ranges are calculated from the corresponding CDF. Using fIGM = 0.85 for the baryon fraction in the IGM17,50, we obtain a median value for DMhost and 68% probable interval ({{rm{DM}}}_{{rm{host}}}={903}_{-111}^{+72}{{rm{pccm}}}^{-3}). When we vary fIGM from 0.6 to 1, the total range of values for DMIGM is about 80 pc cm−1 to 350 pc cm−1 and for DMhost is about 1,020 pc cm−1 to 745 pc cm−1, where we have taken the median values resulting from each value of fIGM and extended the ranges using (half of the) corresponding 68% probable ranges for the smallest and largest values of fIGM.

For comparison, we also give redshift estimates for different values of fIGM if a fixed value of ({{rm{DM}}}_{{rm{host}}}=50{{rm{pccm}}}^{-3}) is used (as is often found in the literature) along with the above quoted mean value for the MW contribution. For fIGM ranging from 0.6 to 1, we obtain a range of redshift values from 2.2 to 0.9 (including the 68% redshift interval derived from the redshift CDF for each value of fIGM), much larger than for the redshift of 0.241 for the identified host galaxy.

Constraints on DMhost from Hα measurements

The DM of the host galaxy is independently estimated from its Hα emission by converting the extinction-corrected Hα flux, ({F}_{{rm{H}}{rm{alpha }}}=(42.0pm 0.5)times {10}^{-17}{{rm{ergs}}}^{-1}{{rm{cm}}}^{-2}), to a Hα surface density of 224 ± 3 Rayleighs in the source frame at z = 0.241, assuming host-galaxy dimensions of 0.5 × 0.5 arcsec (ref. 51). These host-galaxy dimensions are only a rough estimate based on the size of the Hα-emitting region in the Keck image, but because the image is seeing limited, the assumed dimensions may be even smaller, which would only serve to increase the Hα surface density in the following calculations. The Hα flux and estimated surface density are similar to those found for the host galaxy of FRB 20121102A51. The Hα surface density in the source frame (S{({rm{H}}{rm{alpha }})}_{{rm{s}}}) is related to the Hα surface density in the observer frame as (S{({rm{H}}{rm{alpha }})}_{{rm{s}}}={(1+z)}^{4}S{({rm{H}}{rm{alpha }})}_{{rm{o}}}), and is related to the emission measure (EM) in the source frame by

$${{rm{EM}}}_{{rm{s}}}=2.75,{rm{pc}},{{rm{cm}}}^{-6}{T}_{4}^{0.9}S{({rm{H}}{rm{alpha }})}_{{rm{s}}}approx 616pm 7,{rm{pc}},{{rm{cm}}}^{-6}times {T}_{4}^{0.9},left[frac{S{({rm{H}}{rm{alpha }})}_{{rm{s}}}}{224pm 3,{rm{R}}}right],$$


where we have used the redshift z = 0.241 and T4 is the temperature in units of 104 K. The EM is related to the DM in the ionized cloudlet model51 by ({rm{EM}}=[zeta (1+{{epsilon }}^{2})/f]{{rm{DM}}}^{2}/L), where (zeta ) represents cloud–cloud variations in the mean density, ε2 represents the variance of density fluctuations in a cloud, f is the filling factor and L is the path length through the gas sampled by the FRB. Using this relation, we obtain the corresponding source frame DM

$${{rm{DM}}}_{{rm{s}}}approx 392pm 3,{rm{pc}},{{rm{cm}}}^{-3}times {T}_{4}^{0.45},{left(frac{L}{5{rm{kpc}}}right)}^{1/2},{left[frac{f/0.1}{zeta (1+{{epsilon }}^{2})/2}right]}^{1/2},{left[frac{S{({rm{H}}{rm{alpha }})}_{{rm{s}}}}{224pm 3{rm{R}}}right]}^{1/2},$$


where we have adopted fiducial values of f = 0.1 and (zeta ={{epsilon }}^{2}=1), which are typical for the warm ionized medium in the Milky Way. For the path length through the gas, we adopt a fiducial value of L = 5 kpc, which is based on the apparent size of the Hα-emitting region in the host galaxy. However, given that the optical image is seeing limited, and that the orientation of the galaxy relative to the FRB line of sight is not known, this path length could be as large as 10 kpc or as small as 0.1 kpc. For a range of L from 0.1 kpc to 10 kpc, we find that DMs could be as small as 55 pc cm−3 or as large as 560 pc cm−3, for the same fiducial values of T4, f, (zeta ) and ϵ2. The estimated DM is also highly sensitive to the temperature: for a range of T4 from 0.5 to 5, DMs could range from 290 pc cm−3 to 810 pc cm−3, keeping all other parameters fixed at their fiducial values. In the observer’s frame, the measured DM contribution of the host galaxy is smaller by a factor of (1/(1+z)), yielding a nominal value of ({{rm{DM}}}_{{rm{host}},{rm{coeff}}}approx 316pm 2{{rm{pccm}}}^{-3}) for the coefficient in equation (7). The quoted errors account for only measurement errors in the Hα flux. To match the inferred value of DMhost ≈ 900 pc cm−3 requires that the three factors in equation (7) involving T4, L and f combine to a factor of about 3, which may easily be explained by a higher temperature or the broad range of possible values for L. Regardless, the large Hα EM affirms that the FRB DM receives a significant contribution from ionized gas in the host galaxy, but it is unclear whether the diffuse, Hα-emitting gas can account for the entire host-galaxy contribution or whether the FRB’s local environment contains significant ionized gas content that is not seen in Hα.

Extragalactic scattering

Scatter broadening manifests as a frequency-dependent temporal asymmetry in the burst dynamic spectrum, which is typically modelled as a one-sided exponential pulse broadening function convolved with a Gaussian pulse. Some bursts from FRB 20190520B have a burst structure that is suggestive of scattering: leading edges that are aligned across the radio frequency band coupled with temporal asymmetries that broaden at lower observing frequencies. However, a number of bursts are symmetric across the radio frequency band and vary significantly in burst width. When the burst S/N is low, it is also difficult to distinguish scattering from other burst substructure, such as the time–frequency drift of intensity islands9 (often referred to as the ‘sad trombone’ effect). Examples of bursts with and without evidence of pulse broadening are shown in Extended Data Fig. 7.

To measure the mean burst width and scattering timescale of FRB 20190520B, we first integrate the dynamic spectrum of each burst along the frequency axis and normalize the resulting burst profile. Each burst profile is fit with a model composed of a Gaussian component convolved with a one-sided exponential. Out of the 79 fitted scattering times, all scattering times with fractional uncertainties >50% are excised, leaving a subset of 26 bursts. The scattering times of these 26 bursts are then compared with their dynamic spectra and one-dimensional burst profiles to verify that the burst temporal asymmetry is broadly consistent with a v−4 frequency scaling. Supplementary Table 1 shows the scattering times re-scaled to 1.25 GHz assuming a v−4 frequency dependence. The remaining 53 bursts, whose scattering times have fractional uncertainties >50%, are then re-fit with a Gaussian-only model for the one-dimensional burst profile. The full-width at half-maximum burst widths are also shown in Supplementary Table 1. This provisional approach yields a mean scattering time and burst width of 10 ± 2.0 ms and 13.5 ± 1.2 ms at 1.25 GHz, respectively. A more detailed study investigating the burst widths in separate frequency subbands is in progress, and corroborates the scattering interpretation presented here.

The observed mean scattering time of 10 ± 2 ms at 1.25 GHz is probably too small for the FRB’s host galaxy to lie behind our proposed host-galaxy association. In this alternative scenario, the host galaxy would lie at a redshift zh > 0.241 that depends on the foreground galaxy’s contribution to the total DM budget. The FRB would pass through the putative intervening galaxy at a redshift zl > 0.241 at an impact parameter of about 4 kpc (based on the observed offset of the FRB localization in the optical images). The scattering contribution of the intervening galaxy lens is related to its DM contribution by52,53

$$tau ({rm{DM}},nu ,z)approx 48.03,{rm{mu }}{rm{s}}times frac{tilde{F}{G}_{{rm{scatt}}}{{rm{DM}}}_{{rm{l}}}^{2}}{{(1+{z}_{{rm{l}}})}^{3}{nu }^{4}},$$


where DMl is the contribution of the intervening galaxy in pc cm−3 in its rest frame, v is the observing frequency in GHz, zl is the intervening galaxy redshift, and (tilde{F}=zeta {{epsilon }}^{2}/f{({l}_{{rm{o}}}^{2}{l}_{{rm{i}}})}^{1/3}) quantifies the electron density fluctuations in the scattering layer for (tilde{F}) in units of pc−2/3 km−1/3, where (zeta ), ϵ2 and f describe the density fluctuation statistics and filling factor, l0 is the outer scale of turbulence and li is the inner scale. The fluctuation parameter (tilde{F}) in the Milky Way varies from ~10−3 pc−2/3 km−1/3 in the thick disk to (gtrsim )102 pc−2/3 km−1/3 near the inner Galaxy, and is typically about 0.1–1 pc−2/3 km−1/3 in the thin disk, but it is not generally known for other galaxies52. The geometric factor Gscatt is unity for scattering in a host galaxy but it can be very large for an intervening galaxy, in which case ({G}_{{rm{scatt}}}approx 2{d}_{{rm{sl}}}{d}_{{rm{lo}}}/{d}_{{rm{so}}},L), where dsl, dso and dso are the angular diameter distances between the source and lens, lens and observer, and source and observer, respectively, and L is the path length through the lens. The numerical pre-factor in equation (8) is for all angular diameter distances in Gpc and L in Mpc.

For the nominal DM contribution implied by the Hα emission, DMl ≈ 392 pc cm−3, and a path length through the intervening galaxy L ≈ 5 kpc, the implied scattering time is τ ≈ 20 s at 1.25 GHz, orders of magnitude larger than the observed mean scattering time. For a smaller gas temperature T4 ≈ 0.5, yielding a DM contribution DMl ≈ 290 pc cm−3, the implied scattering time is still τ ≈ 18 s at 1.25 GHz, assuming (tilde{F}) 0.1 ≈ pc−2/3 km−1/3 and L ≈ 5 kpc. However, there is significant latitude in this estimate, depending on the assumed gas temperature, path length sampled by the FRB and the value of (tilde{F}). Although we have assumed that the FRB traces the entire DMl implied by the Hα emission, the FRB line of sight probably only traces a fraction of the Hα-emitting gas. If DMl is as small as 50 pc cm−3, then τ could be as small as 0.6 s for (tilde{F}) 0.1 ≈ pc−2/3 km−1/3 and L ≈ 5 kpc. Although the implied scattering time could be two orders of magnitude smaller or larger depending on the combination of (tilde{F}), the path length L and the parameters used to estimate the Hα DM, our fiducial estimates suggest that the observed scattering time is significantly smaller than the scattering expected from an intervening galaxy.

Rotation measure search from FAST observations

A search for RM was performed with the FAST data. The polarization was calibrated by correcting for differential gains and phases between the receptors through separate measurements of a noise diode injected at an angle of 45° from the linear receptors. We searched for the RM from −30 × 105 rad m−2 to 3.0 × 105 rad m−2. No significant peak was found in the Faraday spectrum. The observed lack of polarization could be due to the intrinsic low linear polarization, a depolarization process within the source or from intra-channel Faraday rotation.

Source link