The Curious Case of Twin Fast Radio Bursts: Evidence for Neutron Star Origin?

Fast radio bursts (FRBs) are brilliant short-duration flashes of radio emission originating at cosmological distances. The vast diversity in the properties of currently known FRBs, and the fleeting nature of these events make it difficult to understand their progenitors and emission mechanism(s). Here we report high time resolution polarization properties of FRB 20210912A, a highly energetic event detected by the Australian Square Kilometre Array Pathfinder (ASKAP) in the Commensal Real-time ASKAP Fast Transients (CRAFT) survey, which show intra-burst PA variation similar to Galactic pulsars and unusual variation of Faraday Rotation Measure (RM) across its two sub-bursts. The observed intra-burst PA variation and apparent RM variation pattern in FRB 20210912A may be explained by a rapidly-spinning neutron star origin, with rest-frame spin periods of ~1.1 ms. This rotation timescale is comparable to the shortest known rotation period of a pulsar, and close to the shortest possible rotation period of a neutron star. Curiously, FRB 20210912A exhibits a remarkable resemblance with the previously reported FRB 20181112A, including similar rest-frame emission timescales and polarization profiles. These observations suggest that these two FRBs may have similar origins.


Introduction
Fast radio bursts (FRBs; e.g.Lorimer et al. 2007;Thornton et al. 2013) are intense short-lived radio signals of cosmological origin, the progenitors of which remain unknown to date (e.g.Petroff et al. 2022).There have been a plethora of FRB observations since their discovery (e.g.Shannon et al. 2018;CHIME/FRB Collaboration et al. 2021;Law et al. 2023), which have revealed a vast diversity of burst profiles (e.g.Pleunis et al. 2021), polarization properties (e.g.Day et al. 2020a), host galaxies (e.g.Bhandari et al. 2022), and local magnetoionic environments (e.g.Mannings et al. 2021;Mckinven et al. 2023).This diversity makes it difficult to infer progenitor properties, especially when allowing for selection biases (Macquart & Ekers 2018), effects of propagation through ionized media on the observed burst-properties (e.g.Petroff et al. 2022), and the possibility of multiple progenitor populations (e.g.Caleb et al. 2018).Time resolved analysis of the bursts, with full polarization information, provides key insights to the nature of the FRB progenitors, since changes on sub-millisecond timescales can only be attributed to the progenitor itself, or the magneto-ionic environment in the immediate vicinity of the progenitor (e.g.Luo et al. 2020).Such studies require very high signal-to-noise ratio (S/N) polarization profiles of FRBs at microsecond time resolution, which are relatively rare for non-repeating FRBs (see also Pandhi et al. 2024).
In this work, we present high time resolution polarization properties of FRB 20210912A, which are remarkably similar to those of the previously reported FRB 20181112A (Cho et al. 2020;Prochaska et al. 2019).These observations suggest that these two apparently non-repeating FRBs may have nearidentical progenitors, possibly rapidly rotating neutron stars with similar spin periods.We briefly describe the observations and data analysis methods in Section 2. High time resolution properties of FRB 20210912A are presented in Section 3 and their similarities with those of FRB 20181112A are described in Section 4. We present a possible interpretation of the observations in Section 5 and further discussion on the proposed interpretation in Section 6, and conclude with a summary of the results in Section 7.

Observation and Data Processing
Both FRB 20210912A and FRB 20181112A were detected by the Commensal Real-time ASKAP Fast Transients (CRAFT; Bannister et al. 2017) survey on the Australian Square Kilometre Array Pathfinder (ASKAP; Hotan et al. 2021), by passing an incoherent sum of total power from all antennas to the Fast Real-time Engine for Dedispersing Amplitudes (FREDDA; Qiu et al. 2023).Details of detection and localization are listed in Table 1.The real-time search pipeline, upon detection of FRBs, triggers recording of the raw voltage streams from each ASKAP antenna which are used for detailed offline analysis.Post-processing of the FRB data was carried out using the CRAFT Effortless Localization and Enhanced Burst Inspection pipeline (CELEBI; Scott et al. 2023).Offline correlation of voltage data and interferometric imaging of FRBs, as a part of post-processing, enabled phase-coherent beam-forming at the FRB sky location.The beam-formed data were used to estimate the optimum DM through structure maximization (Sutinjo et al. 2023).Polarization calibration was applied as part of post-processing, using the Vela pulsar (PSR J0835−4510) as the calibrator (Scott et al. 2023;Dial et al. in preparation).
Coherently de-dispersed and polarization-calibrated complex voltage data for two orthogonal linear receptors (X and Y, in the coordinate system defined by the antenna dipoles) were used to construct dynamic spectra of all Stokes' parameters (I, Q,U and V ) adopting the following convention (1) The choice of this convention was driven by the handedness of the ASKAP phased array feed (details are discussed in Dial et al. in preparation).The observed position angle of linear polarization is given by which is measured relative to the coordinate system defined by the receiver dipoles.We do not convert this to absolute position angles.
Full Stokes (I, Q,U,V ) dynamic spectra were constructed at different time and frequency resolutions; however, for most of the analysis in this work, we used 64 frequency channels  (Sutinjo et al. 2023) (channel width ≈ 5.25 MHz).The sensitivity of the system drops sharply at both edges of the observing band.Hence, 5% of the channels at either end of the band (i.e.10% of the channels in total) were excluded, and the effective bandwidth for all subsequent analysis is ≈ 300 MHz.

Measurement of Faraday Rotation
Linearly polarized electromagnetic waves propagating through magnetized plasma with a parallel (to the direction of propagation) component of the magnetic field undergo wavelength (λ) dependent rotation of the PA, which is known as Faraday Rotation.The rotation angle (δPA) is given by where λ 0 is the wavelength corresponding to a reference frequency and PA λ0 is the PA at the reference frequency.The proportionality constant, RM, is known as the Faraday rotation measure.We estimated the RM of the FRBs using two different methods -a linear fit to the variation of PA obs with λ 2 (Equation 6) and the technique of RM synthesis (Burn 1966;Brentjens & de Bruyn 2005;Heald 2009).Linear fits were carried out using 64-channel spectra, and it was verified that using lower resolution spectra yields consistent results.As the change in PA obs across the observing band is less than 90 • (see Appendix B), no corrections were needed to account for wrapping of angles.
Results obtained from linear fits are quoted as the estimated values of RM.Independent estimates of RM from the RMsynthesis method, obtained using the publicly available package RM Tools (Purcell et al. 2020), were used to validate the results from linear fits.In all cases, estimates of RM obtained from these two methods agree well within the uncertainties.It was verified that using finer frequency resolutions (up to 8192channel spectra) does not significantly change the results from RM synthesis.
For both FRBs, the average RM (i.e.RM avg ) was measured from the time-averaged spectra over the entire burst.Additionally, we also estimated RM over smaller time bins to probe the RM variation across the bursts.

Correction for Faraday Rotation
The observed Q,U dynamic spectra were 'corrected' (derotated) in order to remove the effect of Faraday rotation, using the average RM for each FRB, applying the wavelengthdependent transformation is the wavelength-dependent de-rotation angle.The reference wavelength (λ 0 ) for this de-rotation was chosen to be the wavelength corresponding to the central frequency of the observing band (not infinite frequency, as is sometimes chosen).We emphasize that for each FRB, de-rotation to Q,U dynamic spectra was applied for the average RM only.No correction or de-rotation was applied for short time-scale intra-burst RM variations.

Polarization Time Profiles and Spectra
Time profiles for all four Stokes parameters were constructed by averaging the dynamic spectra over all frequency channels.The de-rotated Q,U dynamic spectra were averaged over frequency to generate corrected Q,U time profiles.The PA profiles were then generated from the corrected Q,U profiles using the relation while the linearly polarized flux density profiles were generated using A further bias correction was applied to L, to account for unpolarized noise (see Everett & Weisberg 2001;Day et al. 2020b).The total polarized flux density profiles were generated using Estimates of the PA were discarded if L < 3 σ I (where σ I is the RMS noise in the total intensity profile), or if the uncertainty ∆PA > 5 • .Note that as the de-rotation was done with respect to the central frequency of the observing band, rather than the position angle of linear polarization at an infinite frequency.Spectra for all four Stokes parameters were generated by averaging the corresponding dynamic spectra (corrected dynamic spectra for Q,U) over a specified time range (or the entire burst duration).Spectra for PA, L and P were generated from the spectra of the Stokes parameters using the same equations mentioned above.

FRB 20210912A
Phase-coherent beam-forming of FRB 20210912A revealed two prominent sub-bursts: a strong primary sub-burst (A) followed by a weaker secondary one (B), separated by ∆T = 1.27 ± 0.11 ms in the observer frame, as shown in Figure 1 (see also Marnoch et al. 2023).Details of the measurements of sub-burst separation are given in Appendix A. Full Stokes (I, Q,U,V ) dynamic spectra and time profiles of FRB 20210912A are shown in Figure A1.The high detection signal to noise ratio (S/N ≈ 500) of this FRB facilitates time-resolved analysis across individual sub-bursts.Each sub-burst of FRB 20210912A is composed of multiple components with different spectral shape (see Figures A1 and A2) while low-intensity emission is present between the two prominent sub-bursts.
Despite a deep optical search with the Very Large Telescope (VLT), the host galaxy of FRB 20210912A remains hitherto undetected (Marnoch et al. 2023).Optical limits imply a distant host galaxy at a redshift of z ≳ 0.7, assuming a galaxy at least as luminous as the dwarf host galaxy of FRB 20121102A, the least luminous known FRB host galaxy (Tendulkar et al. 2017).Including these constraints with multiparameter fits to the cosmological redshift-DM relation and uncertainties therein (the 'Macquart relation';Macquart et al. 2020;James et al. 2022) yields a redshift estimate of z = 1.18 ± 0.24 (Marnoch et al. 2023).

Faraday Rotation Measure
The average RM of FRB 20210912A, measured from Q−U spectra time-averaged over the entire burst profile, is RM avg = 4.55 ± 0.49 rad m −2 .The Galactic RM in the direction of this FRB, 8 ± 4 rad m −2 (Hutschenreuter et al. 2022;Prochaska et al. 2023), is poorly constrained.Hence it was not possible to obtain a reliable estimate of the extragalactic component of the RM, however it is unlikely to be large.
However, the RMs of the two sub-bursts A and B are significantly different from each other, with RM A = −2.33 ± 0.37 rad m −2 and RM B = 11.32 ± 0.75 rad m −2 , as shown in Figure 1 and listed in Table 2. Polarization spectra of sub-bursts A and B (before correcting the Q − U dynamic spectra for RM avg ) show a clear difference between the slopes of the PA vs λ 2 curves for the two sub-bursts (see Appendix B).The absolute difference between the RMs of the two sub-bursts is |RM A − RM B | = 13.7 ± 0.8 rad m −2 .The high S/N of FRB 20210912A allows us to probe the temporal variation of RM within each sub-burst, at time-scales of tens of µs.Both subbursts of FRB 20210912A exhibit short time-scale (∼ 10 µs) variation of RM across them, as shown in Figure 2, with RM varying monotonically on either side of an extremum in each sub-burst.The extrema, which occur close to the sub-burst peaks, have opposite natures (minimum and maximum) in the two sub-bursts.
RM-synthesis yielded entirely consistent values of RM with those obtained from linear fits (described above), as listed in Table 2 and shown in Appendix B.Here we emphasize that the measured RM represents ∂PA/∂λ 2 , i.e. the local slope of the PA vs λ 2 curve, and may not be associated with the phenomenon of Faraday rotation.We note that the PA spectra, for some time bins, show hints of deviation from a linear variation of PA with λ 2 (see Appendix B).
The circular polarization fraction (V /I) shows weak (but measurable) dependence on wavelength.The (local) slope of the V /I vs λ 2 curve, κ[= ∂(V /I)/∂λ 2 ], for the spectrum integrated over the entire burst, is 10.5 ± 1.1 m −2 .The values of κ, for spectra integrated over each of the sub-bursts A and B, are consistent within the errors.However, κ shows temporal variation across each of the individual sub-bursts, with generally steeper values close to the centres of the sub-bursts and shallower at the edges, following broadly the same pattern as the apparent RM variation (though without any sign reversal).Correlated variation of apparent RM and κ may arise from Generalized Faraday Effects (see Kennett & Melrose 1998;Ilie et al. 2019;Noutsos et al. 2009;Kumar et al. 2022), in which case the λ-dependence of PA would not follow the Faraday law.

Polarization Time Profile
After correcting for RM avg = 4.55 rad m −2 , both sub-bursts of FRB 20210912A are found to be highly polarized, with total polarization fractions ≳ 70%.The fractional linear and circular polarization, as well as the position angle (PA) of linear polarization, vary across the sub-bursts, as shown in Figure 3.
As described in Sections 2.2 and 2.3, the PA of linear polarization was calculated at the central frequency of the observing band, after correcting for the average RM.Extrapolation of PA to infinite frequency (i.e.PA λ=0 ) has not been performed.For time independent Faraday rotation, as expected from the inter-stellar and the inter-galactic media (which are not likely to significantly change on timescales of ∼ms), PA at the central observing frequency has a constant offset from PA λ=0 .As discussed in the previous sub-section, the apparent short timescale RM variation may not be associated with Faraday rotation, in which case a linear (with respect to λ 2 ) extrapolation of PA to infinite frequency would not be meaningful.
The two sub-bursts show opposite signs of PA evolution near the peaks, with the PA rotating clockwise at the peak of the first sub-burst, and counter-clockwise at the peak of the second, as shown in Figure 3.For both sub-bursts, the fastest rate of PA rotation temporally coincides with the intensity peak within the estimated uncertainties (see Section 5.3 and Appecndix D).

Similarities with FRB 20181112A
FRB 20181112A, a FRB also detected and localized by ASKAP in the CRAFT survey (Cho et al. 2020;Prochaska et al. 2019), has a total intensity profile qualitatively similar to that of FRB 20210912A.Detailed analysis revealed further surprising similarities between the time scales and polarization properties of these two apparently unrelated events.Quantifying the similarity between these two FRBs is difficult, as discussed in Appendix E, due to the small number of high-S/N FRBs with time-resolved polarisation properties, and the lack of an appropriate null hypothesis of FRB behaviour against which to test.For the ease of comparison, high time resolution data for FRB 20181112A have been re-analysed using the same methods that were used for FRB 20210912A in this work and the results of the re-analysis are entirely consistent with the previously published (Cho et al. 2020).

Burst Profile and Emission Timescale
We used a time resolution (observer frame) of 3.8 µs to study the high time resolution properties FRB 20181112A, chosen such that the fine structures are resolved while keeping S/N sufficiently high.FRB 20181112A also exhibits a bright primary sub-burst (A) followed by a relatively faint secondary sub-burst (B), with ∆T = 0.81 ± 0.06 ms in the observer frame.Unlike FRB 20210912A, FRB 20181112A exhibits two more faint components (see Cho et al. 2020).
The total intensity profile of FRB 20181112A, when scaled to the same peak intensity and temporal separation between sub-bursts, has a remarkable correspondence to that of FRB 20210912A, as shown in Figure 4. We find (see Appendix A) that the ratio of the widths of sub-bursts A and B (i.e.FWHM A /FWHM B ) -which is 0.31 ± 0.02 for FRB 20181112A and 0.32 ± 0.01 for FRB 20210912A -is the same for these two FRBs within 5% (and 1σ).The width of the primary sub-burst relative to the separation between subburst peaks (FWHM A /T) -which is 0.046 ± 0.003 for FRB 20181112A and 0.052 ± 0.005 for FRB 20210912A -agrees within 1σ.The width of the secondary sub-burst relative to the separation between sub-burst peaks (FWHM B /T) -which is 0.148 ± 0.007 for FRB 20181112A and 0.16 ± 0.01 for FRB 20210912A -also agrees within 1σ.Table 3 summarises the temporal properties of the bursts.This means, although the absolute (observed) timescales of the two FRBs are different, their relative emission timescales are surprisingly similar.The above observed time scales have been modified by cosmic expansion, and hence redshift measurements of the FRB host galaxies are crucial to infer the intrinsic timescales.Optical follow-up observations have revealed that FRB 20181112A originates from a galaxy at z = 0.4755 (Prochaska et al. 2019), implying a rest-frame sub-burst separation of ∆T = 0.55 ± 0.04 ms.The rest-frame emission time-scales of these two FRBs would be identical if the host galaxy of FRB 20210912A is at a redshift of z = 1.35, which is entirely plausible given the redshift estimate in Section 3.

Intra-burst Variation of Rotation Measure
The average RM of FRB 20181112A, measured over the entire burst profile, is RM avg = 13.2 ± 1.0 rad m −2 .The RMs of the two sub-bursts A and B are significantly different from each other, by ∆RM AB = 15.2 ± 3.7 rad m −2 , as shown in Figure 1 and listed in Table 2.Although the average RM of FRB 20181112A -as well as the RMs of its two sub-bursts are different from those of FRB 20210912A -the absolute difference between the RMs of sub-bursts A and B are surprisingly similar (and formally consistent within the uncertainties) for these two FRBs.
The Galactic RM in the direction of FRB 20181112A is 16 ± 6 rad m −2 (Hutschenreuter et al. 2022;Prochaska et al. 2023).The sightline to FRB 20181112A through the Galactic interstellar medium (ISM) cannot change appreciably over timescales of ∼ ms.Hence the excess RM -i.e. the RM of an FRB after subtracting the Galactic contribution -follows the same variation pattern as that of the total RM (which are shown in Figures 2 and 5) with a constant offset equal to the Galactic RM in the direction of the FRB.
We note that the observed wavelength is longer than the emitted wavelength (in the rest-frame of the FRB host galaxy) by a factor of (1 + z).Assuming ∆RM AB has an origin within the FRB host galaxy (including regions close to the FRB source), the intrinsic value of ∆RM AB is larger than its observed value by a factor of (1 + z) 2 .This would imply different values of intrinsic difference between the RMs of the subbursts in these two FRBs, ≈ 65 rad m −2 for FRB 20210912A and ≈ 33 rad m −2 for FRB 20181112A.
Only sub-burst A of FRB 20181112A has sufficient S/N to probe the temporal variation of RM within the sub-burst.The RM profile is qualitatively similar to that of sub-burst A of FRB 20210912A, with a minimum close to the peak, as shown

Polarization Profiles
As mentioned earlier, the polarization time profiles were obtained by averaging the dynamic spectra over the frequency band, after correcting for the average RM.For both FRB 20210912A and FRB 20181112A, the fractional linear and circular polarization vary across the sub-bursts, as seen in Figures 3 and 6; fractional circular polarization shows weak frequency dependence (see Figures A5 and A9).
FRB 20181112A exhibits PA evolution across its primary sub-burst (A) similar to that of FRB 20210912A, as shown in Figure 6.The fastest rate of PA variation occurs near the peak of the sub-burst, as is observed for FRB 20210912A (see Appendix D for details).However, the lack of sufficient S/N does not allow probing the temporal variation of the PA across the secondary sub-burst (B) of FRB 20181112A.

Possible Interpretation
Several pieces of circumstantial evidence suggest that the progenitors of at least some FRBs are likely to be compact objects, possibly neutron stars (e.g.Petroff et al. 2022;Farah et al. 2018;Luo et al. 2020).The observed properties of FRB 20210912A and FRB 20181112A exhibit features qualitatively similar to those observed in Galactic pulsars -including high polarization fraction, intra-burst variation of fractional linear and circular polarization, variation of the position angle (PA) of linear polarization, short timescale appparent RM variation (e.g.Mitra et al. 2023;Smits et al. 2006;Yan et al. 2011;Dai et al. 2015;Noutsos et al. 2009) -supporting a neutron-star origin of these two events.Based on these qualitative similarities with the Galactic pulsars, here we propose a possible interpretation for the observed properties of FRB 20210912A and its striking similarities with FRB 20181112A.However, we acknowledge that alternate interpretations of the observations remain possible and may lead to completely different conclusions about the progenitor of these two FRBs.

Short Timescale RM Variation
As shown in Figuress 1, 2 and 5, sub-bursts of both FRB 20210912A and FRB 20181112A show significantly different RMs, while the observed RM is also found to vary across individual sub-bursts at timescales of ∼ 10 µs.The observed variation of RM is unlikely to be associated with changes in the degree of Faraday rotation in the inter-stellar or the intergalactic plasma, as magneto-ionic properties of these media are not expected to change on such short time scales.
Previous studies on RM variation in Galactic pulsars suggest that such short timescale 'apparent' RM variation may arise from scatter broadening of the pulse due to propagation through inhomogeneous and turbulent media, incoherent addition of quasi-orthogonally polarized emission modes with different spectral behaviour, or magnetospheric/generalized Faraday effects (e.g.Dai et al. 2015;Noutsos et al. 2009Noutsos et al. , 2015;;Ramachandran et al. 2004;Ilie et al. 2019).We reiterate that in all these cases, the wavelength dependence of PA is not governed by the Faraday law, and hence the apparent RM only represents the local slope of the variation of PA with λ 2 (i.e.∂PA/∂λ 2 ).The apparent RM hence cannot be used to estimate the value of PA at infinite frequency, which is the rationale behind our choice of normalization in Equation 6 and the reference frequency for de-rotation (see Section 2.2).
The hints of deviation from the Faraday law observed in the polarization spectra of FRB 20210912A (see Appendix B) cannot distinguish among the possible reasons behind the apparent RM variation.However, the correlated variation of apparent RM and κ (see Figure 2) suggests that the observed behaviour is likely to originate from 'generalized Faraday effects' in dense ionized media close to the source -possibly in the magnetosphere or near wind region of a neutron star (e.g.Cho et al. 2020;Kennett & Melrose 1998;Ilie et al. 2019;Lyutikov 2022).In this case, the RM variation pattern and its associated timescale is expected to be related to the magnetic field geometry near the emission source (e.g.Wang et al. 2011;Lyutikov 2022).
The difference between the measured RMs of the two subbursts of FRB 20210912A and FRB 20181112A, and the qualitative similarity in the variation of the apparent RM across sub-burst A of both FRBs at comparable timescales, indicate that the physical reasons behind the RM variation are likely to be same in both FRBs.The observed reversal in the nature of RM variation between sub-bursts A and B of FRB 20210912A may be associated with the reversal of the magnetic field geometry near opposite poles of a compact magnetized progenitor, as is expected from generalized Faraday effects in the magnetosphere or in the inner wind region (see Wang et al. 2011;Lyutikov 2022).

PA Evolution across Sub-bursts
Pulsar-like PA evolution has been observed for other FRBs (e.g.Pandhi et al. 2024;Mckinven et al. 2024).The PA 'swing' in pulsars (albeit typically for average profiles) is generally described by the 'rotating vector model' (RVM; Radhakrishnan & Cooke 1969;Johnston & Kramer 2019), where the PA traces the projection of the magnetic field at the emission site onto the sky plane as the neutron star rotates.The 'swing' of PA is attributed to the change in viewing geometry of the magnetic field around the neutron star.The fastest rate of PA rotation is expected to coincide with the 'centre' of the emission beam in this model, as is observed for FRB 20210912A and FRB 20181112A (see also e.g.Blaskiewicz et al. 1991).
The opposite signs of PA evolution in the two sub-bursts of FRB 20210912A can be qualitatively explained by a scenario where they are associated with emission from opposite magnetic poles of a neutron star and the line of sight intersects the emission beams from opposite poles at either side of the beam centre, e.g. from above and below.Such behaviour has been observed in some Galactic pulsars that show inter-pulse emission, albeit for average profiles (e.g.Johnston & Kramer 2019;Kramer & Johnston 2008).
In a simple RVM, the magnetic field structure of a neutron star is assumed to be di-polar and the radio emission is assumed to originate near the 'polar cap' region.In this simple geometry, the PA evolution across a pulse is given by where Θ is the magnetic obliquity (i.e. the angle between the rotation axis and the magnetic axis), α is the inclination (i.e.angle between the rotation axis and the line-of-sight), φ is the rotation phase and PA 0 is the observed PA at a reference rotation phase φ 0 .However, in reality, many pulsars exhibit complex temporal variation of the PA with pulse phase (e.g.Mitra et al. 2023;Smits et al. 2006).Significant deviations from a simple RVM may occur due to a number of factors, including relativistic effects, wobbling of the neutron star, complex magnetic field structures (deviations from a di-polar geometry), and presence of orthogonal polarization modes (e.g.Cordes et al. 1978;Blaskiewicz et al. 1991;Hibschman & Arons 2001).This makes quantitative fits of the RVM difficult for many sources.We also note that single pulses of pulsars often exhibit significant deviations from the PA trends of their average profiles (e.g.Singh et al. 2024).

A Rotating Vector Model for FRB 20210912A
Assuming a simple RVM with di-polar magnetic field (Equation 12), the fastest rate of PA evolution is given by in the observer frame, where T obs is the observed rotation period and β is the 'impact angle' (= α − Θ).Assuming that subbursts A and B are associated with opposite magnetic poles, we have from geometry.Neglecting the difference in emission heights at the two poles (e.g.Johnston & Kramer 2019), the time difference between locations of the fastest rate of PA evolution in sub-bursts A and B is approximately equal to half of the rotation period.
The rate of PA change (dPA/dt) was calculated by fitting local tangents to the PA curves, details of which are described in Appendix D. The fastest PA evolution rate in sub-burst A is 0.424 ± 0.016 deg µs −1 while sub-burst B has a fastest PA swing rate of −0.177 ± 0.006 deg µs −1 .The locations of the fastest rates of PA evolution in sub-bursts A and B are separated by 1.24 ± 0.03 ms, implying a rotation period (in the observer frame) of 2.48 ± 0.06.Using these estimates we infer an inclination of α = 76.2 • ± 1.7 • , magnetic obliquity of Θ A = 59.1 • ± 1.4 • (primary pole), Θ B = 120.9• ± 1.4 • (secondary pole) and impact angles of The half opening angle of the emission beam is given by (e.g.Johnston & Kramer 2019), where W is the width of the pulse.Using the FWHM of subburst A, we infer a half opening angle of ρ e = 44.8• ± 1.7 • for the emission beam.We note that this estimate relies on a number of simplifying assumptions which may not always hold, and hence the uncertainties are underestimated.Using the RVM-derived observed rotation speed, and incorporating redshift uncertainty, we infer an intrinsic spin period of 1.14 ± 0.13 ms for the progenitor of FRB 20210912A.
The PA evolution for an RVM with the inferred viewing geometry and rotation period is shown in Figure 3.While the observed PA evolution of FRB 20210912A near the peaks of the two sub-bursts is well-described by a di-polar RVM with the estimated parameters, significant deviations occur away from the peaks.In particular, these deviations are most apparent trailing the primary sub-burst and leading the secondary sub-burst, and coincide with contributions from faint extended emission components that exhibit different spectral properties (see Appendix A).This faint extended emission appears to have a flat PA profile.Nonetheless, as the steepest derivative of PA has been used to estimate the RVM parameters, deviations from the model far from this point of steepest rate of PA change do not contribute appreciably to our estimates of the uncertainties on the inferred parameters, and hence these uncertainties may be underestimated.We also cannot rule out the possibility that the PA evolution near the peaks are attributable to some physical mechanism other than the RVM, largely because the detailed physics of the FRB emission mechanism is not yet understood.Thus other interpretations of our measurements may yield different conclusions on the properties of the progenitor of FRB 20210912A.

Similar Progenitors for Two FRBs?
The striking resemblance between FRB 20210912A and FRB 20181112A -including profile shape, differential RM and short timescale RM variation pattern, polarization properties and evolution of PA across sub-bursts -suggests similar origin for these two apparently non-repeating FRBs.Their near-identical rest-frame emission time-scales -which would be exactly the same if the host galaxy of FRB 20210912A is at a redshift of z = 1.35 -may be attributable to (near-)identical physical conditions of their progenitors.This opens up the intriguing possibility of the existence of a class of transients with the same characteristic rest-frame emission time scales -cosmological "standard clocks".
The observed properties of both FRB 20210912A and FRB 20181112A appear broadly consistent with emission from rotating compact magnetized objects with rotation periods of ≈ 1.1 ms.This inferred rotation speed is higher than that of the fastest known millisecond pulsar (period = 1.4 ms), and close to the maximum allowed rotation speed for neutron stars (Hessels et al. 2006;Haskell et al. 2018).These two FRBs could hence be associated with impulsive radio emission from nearmaximally-rotating neutron stars.The hypothesis of a millisecond neutron star progenitor would naturally explain the intrinsic similarities between these two FRBs, due to the physical limit of maximum rotation speed.
The lack of significant time lag between the peak of the total intensity profile and the point of steepest PA variationassuming that the total intensity peak coincides with the centre of the emission beam and the time lag could be caused by aberration and retardation effects (e.g.Blaskiewicz et al. 1991;Johnston & Kramer 2019) -suggests that the observed radio emission originates close to the neutron star surface, at emission heights of ≲ 10% of the radius of the light cylinder (see Appendix D).However, this estimate critically relies on several assumptions which may not be valid in these cases.

A Possible Sub-class of FRBs?
The existence of two near-identical FRBs does not imply that all FRBs have similar origins.Some FRBs have been observed to exhibit quasi-periodicity which does not appear related to a spin period (Chime/Frb Collaboration et al. 2022;Pastor-Marazuela et al. 2023).Observations of pulsars and magnetars have shown that quasi-periodic temporal structures can originate with frequencies orders of magnitude higher than the spin period (e.g.Kramer et al. 2023), but such bursts present very differently in the polarization domain, showing flat PA curves in stark contrast to FRB 20210912A and FRB 20181112A.However, the existence of two remarkably similar FRBs suggests that at least a sub-class of FRBs may originate in near-maximally rotating neutron stars, although identification of such events may not always be possible due to various possible reasons discussed in Appendix F.
We do not find any other event in the current CRAFT FRB sample with high time resolution data available (Shannon et al. in preparation) that has observed properties and rest-frame timescales similar to FRB 20210912A and FRB 20181112A.Based on this fact, we estimate a 68% confidence limit on the fraction of such FRBs detected by ASKAP/CRAFT of 0.06-0.43(see Appendix F).This compares to the small fraction (≈ 2%) of known pulsars that show inter-pulse emissionevidence for emission from both poles -but with faster spinning pulsars having a greater prevalence for inter-pulses (e.g.Weltevrede & Johnston 2008;Keith et al. 2010;Kramer et al. 1998).A thorough search for similar events in other FRB surveys (e.g.CHIME/FRB Collaboration et al. 2021;Law et al. 2023) is beyond the scope of this work.

Further Implications
Our interpretation of FRB 20210912A and FRB 20181112A does not explain the physics of the FRB emission mechanism -in particular, why the emission is observable for only a short duration.This is the case for the vast majority of FRBs detected to date, as well as rotating neutron stars in the Galaxy that exhibit bright yet isolated radio pulses (Rotating Radio Transients;e.g. McLaughlin et al. 2006)).However, the neutron star magnetosphere-based interpretation presented here does not preclude the later detection of a repeat burst from one of these sources.If a repeat burst was detected, the rapid spin-down of these objects should be detectable: assuming spin-down is governed by magnetic dipole radiation, a spin-down of 0.1 ms would be expected within 90 days if FRB 20210912A behaves like the Crab pulsar (period P and its derivative Ṗ obey P Ṗ = 1.4 × 10 −14 s (Lyne et al. 2015)), and within 36 minutes if it behaves like young magnetars such as SGR J1935+2154 (P Ṗ = 4.6 × 10 −11 s (Israel et al. 2016)).While a relationship between period and redshift (analogous to the Macquart relation between DM and redshift) would be apparent regardless of the lifetime of the progenitors, the relationship would be tightest in the instance where progenitor lifetimes are short and detectable bursts are most commonly observed while the progenitor is still nearmaximally rotating.This scenario is consistent with emission near the light cylinder (e.g.Cognard et al. 1996).Confirmation of a periodicity-redshift relation for FRBs showing similar polarization properties as FRB 20210912A and FRB 20181112A would thus enable a redshiftless tool for FRB cosmology.

Summary
In this work, we present high-time-resolution polarization properties of FRB 20210912A, which shows remarkable resemblance with the previously reported FRB 20181112A.These two apparently non-repeating FRBs have similar burst structures, near-identical rest-frame emission timescales, and similar PA evolution and similar variation of (apparent) RM across the bursts.The observed PA swing and apparent RM variation pattern in these two FRBs may be explained by a rapidly-spinning-neutron-star origin, with rest-frame spin periods of ∼ 1.1 ms -comparable to the shortest known period of a pulsar and close to the shortest possible rotation period of a neutron star.We emphasize that other interpretations of these observations remain possible, which may lead to completely different conclusions.Nevertheless, the observed properties of these two FRBs provide a unique opportunity to probe the progenitors of such energetic events and hint at the existence of a class of cosmological transients with the same characteristic rest-frame emission time scales.gion) Noongar people and pay our respects to elders past, present and emerging.CWJ and MG acknowledge support by the Australian Government through the Australian Research Council's Discovery Projects funding scheme (project DP210102103).RMS acknowledges support through ARC Future Fellowship FT190100155.RMS and ATD acknowledge support by the Australian Government through the Aus-

A.2. FRB 20181112A
Each of the two sub-bursts, A and B, of FRB 20181112A is found to be optimally described by 4 components, as shown in Figure A4.The details of the measurements are listed in Table 3.We note that the optimum model does not capture the faint emission component at t ≈ 1.2 ms (see also Cho et al. 2020).However, this faint component does not have any significant overlap with the two prominent sub-bursts, and hence does not affect the measurement of the sub-burst widths or the separation between them.
The scattering timescale, which was kept independent for each sub-burst, was found to be consistent in the two sub-bursts with best-fit values of τ A = 19 ± 4 µs and τ B = 24 ± 9 µs.These estimates are also consistent with the scattering time-scales measured by Cho et al. (2020) and Prochaska et al. (2019).
The FWHM of sub-bursts A and B are FWHM A = 37 ± 2 µs and FWHM B = 120 ± 6 µs, respectively.The peaks of the two sub-bursts are separated by T AB = 0.809 ± 0.063 ms.

B. Frequency dependence of polarization
The polarization spectra of FRB 20210912A and FRB 20181112A, before any correction for RM, are shown in Figures A5, A6, A7, A9 and A10, integrating over the entire bursts, each sub-burst as well as over shorter time ranges.The wavelength dependence of PA was fitted with a linear relation between PA and λ 2 , as described in Section 2.1, and the measured slope is quoted as the estimate of RM.The non-Gaussian statistics of the PA errors (e.g.Ilie et al. 2019) have not been taken into account.We normalized the relation between PA and λ 2 at the centre of the observing band (ν 0 ) using the functional form where λ 0 is the wavelength corresponding to ν 0 and PA λ0 represents the value of PA at this wavelength.Normalization at infinite frequency (i.e.λ = 0), using a form PA = PA λ=0 + RMλ 2 , was not chosen because in case of a non-linear relation between PA and λ 2 , PA λ=0 does not carry any physical significance.Note that the choice of normalization point does not affect the estimate of the slope of the relation (RM) and its uncertainties.
The RM estimates obtained from this method are entirely consistent with the RM estimates obtained from RM-synthesis.Note that both these methods estimate the local slope of the PA with respect to λ 2 , i.e.RM ≡ ∂PA/∂λ 2 .In case the PA has a non-linear dependence on λ 2 , the measured RM is an estimate of the co-efficient of the linear term in the Taylor series expansion at the centre of the observing band.
The observed frequency dependence of the fractional circular polarization was quantified by the (local) slope of V /I with respect to λ 2 , i.e. κ ≡ ∂(V /I)/∂λ 2 .Note that this does not imply assumption of a linear relation between V /I and λ 2 .For a non-linear relation, κ is an estimate of the co-efficient of the linear term in the Taylor series expansion at the centre of the observing band.We estimated the value of κ from polarization spectra integrated over the same time ranges as RM measurements.

B.1. FRB 20210912A
The two sub-bursts, A and B, of FRB 20210912A have different RMs as evident from the slopes of PA with respect to λ 2 in Figure A5.Sub-burst A (middle panel) shows deviation from faraday law (a linear relation between PA and λ 2 ), which leads to a relatively poor fit.Reduced chi-squared (χ 2 r ) of the corresponding fits are mentioned in the lower panels of the Figures.Both sub-bursts, however, have the same value of κ within the errors.Within each sub-burst, the polarization spectra show significant temporal variation with RM and κ varying at timescales of ∼ 10 µs, as shown in Figures A6 and A7.Both RM and κ have extreme values close to the sub-burst peaks.In sub-burst A, the deviation from Faraday law is larger close to the peak.Such deviations are not apparent in sub-burst B.
The estimates of RM obtained from the RM-synthesis method are shown in Figure A8, which agree with the estimates from the fit (shown in Figure 2) within the errors.The values of PA λ=0 (at infinite frequency) obtained from RM-synthesis are also shown in Figure A8 (lower panels), which show correlated variation with the RM estimates.This also indicates to a deviation from the Faraday law in the wavelength dependence of PA.

B.2. FRB 20181112A
The two sub-bursts, A and B, of FRB 20210912A also show different RMs as evident from the slopes of PA with respect to λ 2 in Figure A9.Within sub-burst A RM varies at timescales of ∼ 10 µs, as shown in Figures A10 with extreme value close to the peak.The PA spectra do not show as large deviation from Faraday law as seen in FRB 20210912A.The value of κ is consistent in the two sub-bursts and shows no measurable temporal variation within sub-burst A.  The estimates of RM obtained from the RM-synthesis method are shown in Figure A11, which agree with the estimates from the fit (shown in Figure 5) within the errors.The values of PA λ=0 (at infinite frequency) obtained from RM-synthesis show correlated variation with the RM estimates, similar to FRB 20210912A.

C. Rotation Measure Variation
Short timescale (∼ 10 µs) variation of RM is observed across the sub-bursts of FRB 20210912A and FRB 20181112A, with RM varying monotonically on either side of extrema close to the peaks.To characterize the RM variation pattern, we empirically fit the RM profile of each sub-burst with a Gaussian of the form where RM max , t 0 and w RM are the peak (positive or negative), centre and the characteristic width of the Gaussian, respectively.
C.1.FRB 20210912A RM 0 was set (not fit) to the arithmetic mean of the RMs of the two sub-bursts, i.e.RM 0 = (RM A + RM B )/2 = 4.5 rad m −2 .The peaks of the best-fit Gaussians are RM max A = −17.75± 0.6 rad m −2 and RM max B = 15.5 ± 1.4 rad m −2 .The centres are located at t A 0 = 0.043 ± 0.002 ms and t B 0 = 1.289 ± 0.005 ms.The FWHMs of the best-fit Gaussians are 80 ± 3 µs and 99 ± 13 µs for sub-bursts A and B, respectively.We note that the ratio of the FWHMs is significantly different from the ratio of the burst widths.

C.2. FRB 20181112A
Sub-burst A shows an RM variation pattern qualitatively similar to those of the sub-bursts of FRB 20210912A, while sub-burst B of FRB 20181112A does not have sufficient S/N to probe RM variation across it.RM 0 was set (not fit) to the arithmetic mean of the RMs of the two sub-bursts, i.e.RM 0 = (RM A + RM B )/2 = 18.1 rad m −2 .The peak of the best-fit Gaussian is RM max A = −9.5 ± 0.6 rad m −2 and it is located at t A 0 = 0.011 ± 0.002 ms.The FWHM of the best-fit Gaussian is 42 ± 7 µs.The ratio of the FWHM of this best-fit Gaussian (associated with sub-burst A) to T AB is consistent within the uncertainties for FRB 20181112A (0.052 ± 0.009) and FRB 20210912A (0.061 ± 0.002).not have enough signal to ratio to probe any PA variation across it.The relation in Equation D6 implies an emission height of h em = (0.04 ± 0.02)R lc for sub-burst A.

E. Quantifying Similarities between Two FRBs
We would like to quantify the degree of similarity between FRB 20181112A and FRB 20210912A, to answer the question: what is the likelihood that two FRBs appear so similar due purely to random behaviour, rather than a common underlying physical mechanism?A statistical analysis of similarities between pulse profiles of two FRBs -taking into account different (and possibly unknown) redshifts, different scattering timescales, different S/N etc. -is difficult.An intuitive way to quantify similarity between any two FRBs would be to maximize Pearson's correlation coefficient (e.g.Freedman et al. 2007) between the total intensity profiles I(t) by varying three parameters: the relative amplitude ∆A, the relative start time ∆t 0 , and a time compression factor ∆τ .Doing this for a large sample of FRBs would produce a distribution of maximized correlation coefficients against which the value for FRB 20181112A and FRB 20210912A could be evaluated.
The first problem encountered with such an approach -or any approach trying to quantify similarity -is that no suitable model of FRB time-frequency profiles exist to form a null hypothesis against which to test, e.g. by generating synthetic FRBs.Relying on data however encounters several issues: many FRBs are sufficiently scattered that their intrinsic structure is unresolved, either due to instrumental time resolution, or because their shape is dominated by an exponential scattering tail.Clearly, any two such FRBs, when varying ∆A, ∆t 0 , and ∆τ , will correlate extremely well.Similarly, even FRBs with complicated structure are often dominated by a strong primary peak, which will dominate any tests of correlation regardless of the details of fine structure.Therefore, in any such analysis, it seems reasonable to exclude each FRB's primary peak, and to examine secondary structures for the degree of correlation.
It is at this point that the small number of high S/N, time-resolved FRBs with secondary structures becomes a problem.Of the FRB-detecting facilities with significant numbers of unbiased detections (as opposed to follow-up observations of known objects), only UTMOST, DSA-110, and ASKAP have data for which useful comparisons can be made.Excluding those FRBs with simple, single-pulse structures, and requiring a S/N of 50 (so that a secondary peak at the 10% power level would be detected at 5 sigma), leaves a total of 13 (2 UTMOST (Farah et al. 2018(Farah et al. , 2019)), 2 DSA (Sherman et al. 2023), and 9 ASKAP (Scott et al. 2023)) events.FRB 20181112A and FRB 20210912A are 'obviously' the most similar of these, but we do not consider this sample size sufficient to determine if the similarity is extraordinary.Furthermore, these two FRBs not only have similar rest-frame intensity profiles but also exhibit similar RM variation and PA swing.Capturing all this information in a single statistic is even more challenging and will be attempted as part of a separate work.

F. Number of such FRBs
We consider what fraction, f class , of FRBs detected by CRAFT could plausibly belong to the same class as FRBs 20181112A and 20210912A.However, because this class is currently defined only by these two members, rather than a large analysis of a population of bursts (e.g.Pleunis et al. 2021), the precise class definitions are ambiguous.
Using a strict definition of this class as having a broadband, narrow initial pulse, followed by a dimmer secondary pulse of amplitude less than 50% of the primary pulse, only FRBs 20181112A and 20210912A of the 22 FRBs detected by CRAFT in incoherent sum mode are class members (Cho et al. 2020;Day et al. 2020a;Bhandari et al. 2020Bhandari et al. , 2023;;Marnoch et al. 2023;Scott et al. 2023).Ten have large scattering tails however (τ scat > 0.5 ms), which would make the detection of a small secondary peak very difficult.Such scattering is most likely to arise either from the FRB host galaxy's interstellar medium (ISM) or circum-galactic medium (CGM), and is thus not intrinsic to the FRB emission mechanism (Sammons et al. 2023).Therefore, these FRBs could also be members of the same fundamental class, although a measurement of a rotation period for them would be unlikely.
FRBs 20190102C and 20190611B both exhibit two sub-pulses, but with different relative powers to FRBs 20181112A and 20210912A.FRB 20190102C has a small precursor burst offset from the main pulse by ∼ 0.4 ms and ∼10% of its amplitude, while FRB 20190611B has two bright peaks of almost equal magnitude, separated by ∼ 1.0 ms.If these peaks are associated with emission from opposite poles of a neutron star, the implied rotation periods would be 0.62 ms and 1.5 ms at their respective redshifts of 0.29 and 0.378 (Macquart et al. 2020).The former is excluded on causality considerations (Rhoades & Ruffini 1974;Haskell et al. 2018); the latter remains a plausible candidate.
The remaining eight FRBs do not appear to exhibit structures consistent with that of FRB 20181112A and FRB 20210912A.Our observation of two of 12 weakly scattered FRBs sets a 68% confidence limit on f class of 0.06-0.34;allowing FRB 20190611B to be a potential class member yields the range 0.12-0.43.

Figure 1 .
Figure 1.Burst profile and Faraday rotation measure (RM = ∂PA/∂λ 2 ) in FRB 20210912A (left) and FRB 20181112A (right).The frequencyaveraged Stokes-I profiles of the FRBs are shown in red at a time resolution of 3.8 µs (in normalized flux density units not shown in the plots).The x-errorbars represent the time-range for the corresponding RM measurements.The absolute difference between the RMs of the sub-bursts is ≈ 15 rad m −2 for both FRBs (in the observer frame).See Sections 3 and 4 for details.

Figure 2 .
Figure 2. Short timescale variation of RM [= ∂PA/∂λ 2 ; upper panel] and κ [= ∂(V /I)/∂λ 2 ; lower panel] in FRB 20210912A for sub-bursts A (left) and B (right).The frequency-averaged Stokes-I profile is shown in red at a time resolution of 3.8 µs (in normalized flux density units not shown in the plots).The uncertainties in the abscissa are the time range for the corresponding measurements.The cyan dashed lines show best fit Gaussian profiles to the RM variation (see Appendix C for details).The red lines show the total intensity for each component.

Figure 3 .
Figure 3.Time resolved polarization of FRB 20210912A.[A, upper panel] The frequency-averaged normalized total intensity (I), linearly (L) and circularly (V) polarized intensity at a time resolution of ≈ 3.8 µs.[B, lower panel] Position angle (PA) of linear polarization.Corrections for the average rotation measure (RMavg = 4.55 rad m −2 ) have been applied.The dashed curve shows the PA profile corresponding to a rotating vector model with inclination of α = 76.2 • and magnetic obliquity of Θ = 59.1 • .See text for details.

Figure 4 .
Figure 4. Scaled burst profiles of FRB 20181112A and FRB 20210912A.The frequency-averaged Stokes-I (total intensity) profiles of FRB 20210912A (blue) and FRB 20181112A (red) are plotted against time normalized by the separation between the two sub-bursts (∆T) for each FRB, at a time resolution of ≈ 9.5 µs.Flux densities are normalized by the peak of the profile.

Figure 5 .
Figure 5. Short timescale variation of RM [= ∂PA/∂λ 2 ; upper panel] and κ [= ∂(V /I)/∂λ 2 ; lower panel] in FRB 20181112A for sub-bursts A. The frequency-averaged Stokes-I profile is shown in red at a time resolution of 3.8 µs (in normalized flux density units not shown in the plots).The x-errorbars represent the time range for the corresponding measurements.The cyan dashed line shows the best fit Gaussian profile to the RM variation.

Figure 6 .
Figure 6.Time resolved polarization of FRB 20181112A.[A, upper panel] The frequency-averaged normalized total intensity (I), linearly (L) and circularly (V) polarized intensity at a time resolution of ≈ 3.8 µs.[B, lower panel] Position angle (PA) of linear polarization.Corrections for the average rotation measure (RMavg = 13.15 rad m −2 ) have been applied.

VFigure A4 .
Figure A1.Full Stokes time profile and dynamic spectra (I, Q,U,V ) of FRB 20210912A at a time resolution of 3.8 µs.The flux densities have been normalized by the peak of the total intensity profile.Q and U have been corrected for RMavg = 4.55 rad m −2 .

Figure A5 .
Figure A5.Polarization spectra of FRB 20210912A integrated over the entire burst (left), sub-burst A (middle) and sub-burst B (right) before correcting for the average RM.The time-ranges for averaging are mentioned in each panel.Polarization fractions weakly vary with frequency.Slope of the V /I vs λ 2 curve, κ, is mentioned in the top right corner of each panel.Frequency dependence of PA, especially in sub-burst A, shows hints of deviation from the Faraday law.Reduced chi-squared (χ 2 r ) of the fits are mentioned in the lower panels.Data points are plotted after averaging two adjacent channels of the 64-channel spectra on which the fits were performed.

Figure
Figure A6.Time-resolved polarization spectra of FRB 20210912A in sub-burst A before correcting for the average RM.The time-ranges for averaging and the slope of the V /I vs λ 2 curve, κ, are mentioned in each panel.Frequency dependence of PA shows hints of deviation from the Faraday law in time bins close the the intensity peak.Reduced chi-squared (χ 2 r ) of the fits are mentioned in the lower panels.Data points are plotted after averaging two adjacent channels of the 64-channel spectra on which the fits were performed.

Figure
Figure A7.Time-resolved polarization spectra of FRB 20210912A in sub-burst B before correcting for the average RM.The time-ranges for averaging and the slope of the V /I vs λ 2 curve, κ, are mentioned in each panel.Reduced chi-squared (χ 2 r ) of the corresponding fits are mentioned in the lower panels.See Appendix B for details.Data points are plotted after averaging two adjacent channels of the 64-channel spectra on which the fits were performed.

Figure A8 .Figure A9 .
Figure A8.Rotation measure (RM) of FRB 20210912A estimated using the RM-synthesis method, in sub-bursts A (left) and B (right).The frequency-averaged Stokes-I profile of the FRB is shown in red at a time resolution of 3.8 µs (in normalized flux density units not shown in the plots).The x-errorbars represent the time-range for the corresponding measurements.The lower panels show the corresponding PA at infinite frequency, assuming that PA has a linear dependence on λ 2 .See Appendix B for details.

Figure A11 .
Figure A11.Rotation measure (RM) of FRB 20181112A estimated using the RM-synthesis method, in sub-burst A. The frequency-averaged Stokes-I profile of the FRB is shown in red at a time resolution of 3.8 µs (in normalized flux density units not shown in the plots).The x-errorbars represent the time-range for the corresponding measurements.The lower panel shows the corresponding PA at infinite frequency, assuming that PA has a linear dependence on λ 2 .See Appendix B for details.

FigureFigure
Figure A12.PA swing in sub-bursts A (left) and B (right) of FRB 20210912A.The dashed black lines show the best-fit Gaussians to the rate of change of PA.See text for details.

Table 2 .
RMs of FRB 20181112A and FRB 20210912A * Absolute difference between RM of sub-bursts A and B

Table 3 .
Time scales of FRB 20181112A and FRB 20210912A 2 ] is observed for FRB 20181112A, although we can not rule out such a variation as the measurements are of low (≲ 2σ) significance.The relatively lower S/N of Stokes V (compared to FRB 20210912A), due to a combination of relative faintness and a lower degree of circular polarization, do not allow more accurate measurement of the temporal variation of κ in FRB 20181112A.
tralian Research Council's Discovery Projects funding scheme (project DP220102305).LM acknowledges the receipt of an MQ-RES scholarship from Macquarie University.KG acknowledges support through Australian Research Council Discovery Project DP200102243.This scientific work uses data obtained from the Australian Square Kilometre Array Pathfinder (ASKAP), located at Inyarrimanha Ilgari Bundara / the Murchison Radio-astronomy Observatory.We acknowledge the Wajarri Yamaji People as the Traditional Owners and native title holders of the Observatory site.ASKAP uses the resources of the Pawsey Supercomputing Research Centre.CSIRO's ASKAP radio telescope is part of the Australia Telescope National Facility.Establishment of ASKAP, Inyarrimanha Ilgari Bundara, the CSIRO Murchison Radioastronomy Observatory and the Pawsey Supercomputing Research Centre are initiatives of the Australian Government, with support from the Government of Western Australia and the Science and Industry Endowment Fund.Operation of ASKAP is funded by the Australian Government with support from the National Collaborative Research Infrastructure Strategy.This work was performed on the OzSTAR national facility at Swinburne University of Technology.The OzSTAR program receives funding in part from the Astronomy National Collaborative Research Infrastructure Strategy (NCRIS) allocation provided by the Australian Government, and from the Victorian Higher Education State Investment Fund (VHESIF) provided by the Victorian Government.This research has made use of NASA's Astrophysics Data System Bibliographic Services.