Introduction

The variations of Faraday rotation measure (RM)1, an integral of the line-of-sight component of magnetic field strength and electron density, and dispersion measure (DM)1, defined as the line-of-sight integral over the free electron density, of repeating FRBs contain essential information about their physical nature and central engine, which are as yet unknown1,2,3,4. Recently, the rotation measure of FRB 20201124A shows dramatically oscillations5. FRB 20201124A was detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME)6 on 24 November 20207. It has a dispersion measure of 413.52 pc cm−3, larger than the Galactic DM contribution. This FRB is extremely active, with burst rate up to 50 per hour5. Thanks to its repetition property, it was localized to a galaxy named SDSS J050803.48 + 260338.0 at z = 0.098 by Australian Square Kilometre Array Pathfinder (ASKAP)8, the Five-hundred-meter Aperture Spherical radio Telescope (FAST)5, Very Large Array9, and European Very Long Baseline Interferometry Network10.

FAST detected 1863 independent bursts in a total of 88 hours from 1 April to 11 June 2021, covering 1.0 GHz to 1.5 GHz5. FAST observations discovered three previously unseen characteristics of repeating FRBs. First, this FRB shows dramatic Faraday rotation measure variations on a day time scale. The maximum variation of RM is about 500 rad m−2. A small variation of RM (about 1.8%) was also found from 20 bursts observed by the Effelsberg Radio Telescope11. Second, it is the first repeating FRB shows circular polarization. Third, it is the repeating FRB that has the widest mean burst width12. These characteristics cannot be explained by current theoretical models. FRB 20201124A can be analogy to the pulsed emissions from the Galactic binary system PSR B1259-63/LS 288313, containing PSR B1259-63 and the Be star LS 2883. The pulsed emission of PSR B1259-63 also shows variations in its RM, DM, scattering timescale and polarized intensity on short time-scales when it approaches periastron. It has been argued that a highly magnetized neutron star in an interacting high-mass X-ray binary system can account for all of the observational phenomena of FRB 20180916B14. A Be/X-ray binary system, which is composed of a neutron star and a Be star with a decretion disk can well explain the frequency-dependent active window of FRB 20180916B15.

In this work, we propose that FRB 20201124A is produced by a magnetar orbiting around a Be star companion with a decretion disk. The interaction between radio bursts and the disk of the Be star can naturally explain the varying RM, depolarization, large scattering timescale, and Faraday conversion of FRB 20201124A.

Results

PSR B1259-63/LS 2883

PSR B1259-63/LS 2883 is a binary system hosting a radio pulsar and a high-mass Be star LS2883. The spin period of the pulsar is 47.76 ms and the spin-down luminosity is 8 × 1035 erg/s16. Its distance is 2.6 kpc17. The orbital period is 1237 days with eccentricity e = 0.916,17,18. LS2883 is an early Be type star that possesses a powerful stellar wind and an equatorial decretion disk19. There is some evidences that the Be-star disk is highly inclined to the orbital plane both from the pulsed emission20, the unpulsed emission21, and timing measurements22.

This system shows unpulsed emissions in multiple bands, such as in radio16, X-rays23, GeV24,25, and TeV γ-rays26. The light curves show double-peak profiles, which can be attributed to the interaction between the pulsar radiation and the disk27.

Four pulsed radio observations near periastron have been made in the year of 199413, 199720, 200028 and 200429. When the pulsar approaches the disk, the pulsed emissions of PSR B1259-63 show large variations in flux, degree of linear polarization, DM and RM13,28,29. These variations can be explained as the pulsed emissions interacting with the disk13,30. The properties of pulsed emissions, such as degree of polarization, DM and RM, are similar for the four periastron passages29.

Rotation measure variation

We attempt to explain the observed properties of FRB 20201124A using a physical model for the decretion disk of a Be star (see methods, subsection The disk model). Our model is shown in Fig. 1. The Be star locates at the disk center. For the coordinate system, we set the Be star as the origin and fix the x-axis and y-axis on the disk. The orbit has an inclination ϕ with respect to the disk. Generally, a magnetar is formed by a supernova explosion. A supernova explosion is usually asymmetric and gives a large kick to the newly formed magnetar31. Due to the natal kick of the magnetar, the orbital plane and the disk plane are misaligned32. From observations, the inclination angle has a wide range in Be star/X-ray binary systems, from 25 to 7033. In our model, after the companion dies as a supernova, a magnetar is formed, leading to FRBs generated in the magnetar magnetosphere3,34, which is supported by the FAST observations5. As the orbital motion of the binary, radio bursts from the magnetar pass through different components of the disk to reach the observer, which causes the variation of RM. When the pathway of radio bursts does not intersect the disk, however, a non-evolving RM is expected.

Fig. 1: A schematic diagram of the magnetar/Be star binary model.
figure 1

a The Be star locates at the center of the disk. The magnetar is shown as the red point. The stellar disk plane (shown as the black dash line) is inclined to the orbital plane (shown as purple dash lines) with the angle ϕ = π/12. The purple dash line is the major axis of the magnetar's orbit. When the pathways of radio bursts pass through the disk, the interaction between bursts and disk can reproduce the observed variable rotation measure, depolarization, large scattering timescale, and Faraday conversion. b The face-on view of the Be star's disk. The magnetic field shown as arrows is assumed as azimuthal (or toroidal) in the disk. This model predicts that the RM contribution from the disk changes sign when the magnetar passes in front of the Be star. The RM variation discovered by FAST supports this scenario.

For the binary system PSR B1259-63/LS 2883, significant RM changes were observed both in magnitude and in sign before and after the 2004 periastron passage29. Compared to the Galactic PSR B1259-63, the RM variation of FRB 20201124A is small. A thin disk is considered. We assume the density of disk at the stellar surface ρ0 = 3 × 10−14 g cm−3 and β = 4, which are the typical value for disks of Be stars35. The orbital period is assumed to be 80 days with orbital eccentricity e = 0.75. Using the disk model in equation (1) and the magnetic field profile in equation (2), we can reproduce the overall RM variation with time. For the spherical coordinate origin at the Be star, the observer’s direction is θ = 3.08 rad and ϕ = 2.18 rad. The derived RM variation is shown as the red line with magnetic field B0 = 10 G at the stellar surface in Fig. 2. This magnetic field strength is consistent with those derived from RM variation and Faraday conversion5. The blue points is the mean observed RM variation ΔRM = RM-RMc in a day. RMc is taken to be constant, as observed after MJD 59350. The blue region shows the range of ΔRM in each day.

Fig. 2: The fit of RM evolution of FRB 20201124A using the binary model.
figure 2

The blue scatter is the mean of observed ΔRM in day, and the red line is the predicted evolution of our model. The blue shadow region is the range of observed RM. In the calculation, the disk density is ρ0 = 3 × 10−14 g cm−3 close to the star and has a steep-decay index β = 4. The orbital period is taken to be 80 days with orbital eccentricity e = 0.75. Our model can reproduce the main structure of RM evolution. However, there are some minor structures can not be reproduced. These minor structures may be caused by the clumps in the disk. Source data are provided as a Source Data file.

Some dips and peaks of ΔRM are not well reproduced. The possible reason is that the disk is very clumpy both in density and/or magnetic field. From observations, a clumpy disk was required to explain the multi-wavelength observations of the Galactic binary system PSR B1259-63/LS 288336,37. The clumpy disk is also supported by numerical simulations. For example, a two-armed spiral density enhancement caused by the tidal interaction was found using three-dimensional smoothed particle hydrodynamics simulation38. From observations, a stellar wind from a massive star is clumpy, with the filling factor of fc = 0.139. At present, there are no constraints on disk clumps. For simplicity, we assume that they are similar to the clumps in a stellar wind, without magnetic field clumps. The density inside the clumps is 1/fc larger that of the smooth disk. Because the interspace between clumps is not empty, the clumping factor can be expressed as \({f}_{c}=\langle {n}_{e}^{2}\rangle /{\langle {n}_{e}\rangle }^{2}\). The spatial distribution of clumps remains unknown. Using ad hoc clumps to fit the RM variation is less productive. So we discuss the effect of clumps on RM quantitatively. For the RM peak near MJD 59340, the observed value is about ten times larger than the predication of the smooth disk. This implies the value of fc is about 0.1 if only density clumps are considered, being similar to those in a stellar wind.

The observations show that some bursts have low polarization before MJD 59350. The fluctuations of the electron density in the disk can cause differential Faraday rotation along different paths, which can account for the observed burst depolarization (see methods, subsection Depolarization in the disk). The mean burst width of FRB 20201124A is larger than the widths of other known repeating FRBs. The fluctuations of the electron density also can cause scatter broadening of the pulses (see methods, subsection Scattering timescale).

The disk also affects the DM and detection of radio bursts. In Fig. 3, we show the DM and the optical depth due to the free-free absorption contributed by the disk. The DM change shown as the blue line is small. The maximal change is 0.12 pc cm−3. FAST found no significant DM variation in a 95%-confidence-level5, with an upper limit ΔDM ≤ 4.9 pc cm−3. Therefore, the model prediction is consistent with FAST observations. Next, we consider the free-free absorption as a function of orbital phase. The optical depth τ due to free-free absorption shown as the red line is less than unity. The DM and free-free absorption caused by the stellar wind can be safely ignored (see methods, subsection The effect of a stellar wind). The presence of disk clumps also affects the interaction with the radio bursts. The DM and free-free absorption are enhanced by a factor of 1/fc = 10 with respect to the case of a smooth disk40. From Fig. 3, the maximum DM contribution is 4 pc cm−3 in the clumps, which is less than the observed DM variation5. After considering the disk clump, the optical depth to free-free absorption is less than 1. So the clumpy disk does not affect the detection of radio bursts. For FRBs, the induced Compton scattering should be considered41. The induced Compton scattering due to the disk is also estimated to be small, and does not affect FRB detection. Therefore, the disk is transparent to radio emission, and has no effect on the detection of FRBs.

Fig. 3: The DM and free-free absorption optical depth τ contributed by the disk.
figure 3

The same parameters as Fig. 2 are used. The DM variation shown as the blue line is small, which is consistent with the FAST observation. The optical depth τ due to free-free absorption by the disk denoted as the red line is much less than 1. So the existence of the disk has no effect on the detection of radio bursts.

The circular polarization from Faraday conversion

For non-repeating FRBs, variations in both linear and circular polarization have been observed42. However, previous observations show that repeating FRBs only have high degrees of linear polarization. FRB 20201124A is the first one showing high circular polarization. For the Galactic binary system PSR B1259-63/LS 2883, the linear and circular polarizations varied rapidly in an hour timescale near the periastron29. Physically, besides Faraday rotation, the propagation of FRBs through a magneto-ionic environment can also lead to Faraday conversion. In addition, a fraction of circular polarization of FRB 20201124A is produced by radiation mechanism5. There are two mechanisms to generate Faraday conversion. The first case is that radio signals propagate into a relativistic plasma43. The second case is that the magnetic field quasi-perpendicular to the wave vector is large, which had been used to interpret the circular polarization observed in solar radio bursts44. For the first case, termination shocks can accelerate electron and positron pairs from the magnetar wind to relativistic speeds. Because electrons and positrons contribute with the opposite sign to the circularly polarized component, the different distributions of electrons and positrons are required to generate circular polarization. In our model, a weak magnetar wind is considered. Therefore, the termination shock may be absent. For the binary system PSR B1259-63/LS 2883, the minimum Lorentz factor of the shocked electrons is estimated as45 105. In Fig. 4, the Faraday conversion coefficient as a function of minimum Lorentz factor \({\gamma }_{\min }\) is shown. From this figure, we can see that the Faraday conversion coefficient is found to approach zero when \({\gamma }_{\min }\) is about 200. From the accurate calculations, the exact Faraday conversion is found to approach zero with the increasing electron energy43, because ultra-relativistic electrons hardly interact with a radiation field. When the magnetar is in front of the Be star, FRB signals cannot interact with the plasma shocked by the termination shock, leading to no conversion. So, the Faraday conversion due to relativistic plasma is hard to account for the observed circular polarization (see methods, subsection Origin of the circular polarization).

Fig. 4: Faraday conversion coefficient ρQ distribution as a function of minimum Lorentz factor \({\gamma }_{\min }\).
figure 4

An electron column density of Ne = 1 pc cm−3 is assumed. The plot uses ν = 1.25 GHz, θ = π/4 and p = 2.3. Faraday conversion only takes effect for \({\gamma }_{\min } \; < \; 200\).

Next we use the magnetic field quasi-perpendicular to the wave vector to explain the origin of circular polarization. When the cyclotron frequency is approaching the observation frequency, the Faraday conversion will take effect. From the FAST observations, the required magnetic field B is about 7 Gauss for Faraday conversion5, which is the same order of the parallel magnetic field B. This mechanism naturally explains the circular polarization of FRB 20201124A. FAST observations indicate that the Faraday conversion occurs in the RM variation stage (e.g., bursts 779 and 926). When the RM variation stops, the circularly polarized bursts lack qusai-periodic structures (e.g. Burst 1472). There is no Faraday conversion. This is consistent with the prediction of our model. For the above orbital parameters, the observed Faraday conversion appears when bursts’ paths penetrate the disk. When the pathway of radio bursts does not intersect the disk, the RM variation stops. Meanwhile, there is no perpendicular magnetic field required for Faraday conversion.

Candidates of a magnetar/Be star binary

The total energy emitted by radio bursts during FAST observation is about 3 × 1041 erg in 88 h5. For dipolar magnetic field Bd = 1013 G, the total magnetic energy of a neutron star (NS) is Emag = 3 × 1043 erg. In order to explain the energy budget of FRB 20201124A, about 3% NS magnetic energy is transferred to radio bursts in 88 h, which is impossible. Moreover, the Galactic FRB 20200428 suggests FRBs are accompanied by X-ray bursts (XRBs). The energy ratio between a radio burst and its associated XRB is about 10−546. If this value is valid for FRB 20201124A, the energy of XRBs is EX = 3 × 1041/10−5 ~ 1046 erg. For Emag > EX, the magnetic field Bd is larger than 1015 G. So it should be a magnetar. The bursting activity may be related to the crustal magnetic energy of young magnetars47. If the converted magnetic energy comes from the toroidal component, it is at least two orders of magnitude larger than the dipolar magnetic field of the neutron star. From numerical simulations48 and observations49, the toroidal component could be ten times larger than the dipolar component. So, if the required magnetic field is toroidal (>1015 G), the dipolar magnetic field is larger than 1014 G, which also supports a magnetar. Therefore, if the central engine is a compact stellar object, it should be a magentar from energy budget.

There are two candidates of a magnetar/Be star binary. The first one is LS I + 61 30350. The second questionably one is SGR 0755-2933, which is also referred to as 2SXPS J075542.5-29335351. LS I + 61 303 is a Be gamma-ray binary, with period about 26.496 d and a high eccentric orbit (e = 0.537)52. The mass of the Be star in this binary is about 13 M. It has been found that some magnetar-like bursts are possibly from this binary50,53. More recently, transient radio pulsed emissions from its direction were detected by FAST54, with a period of 269 ms. This supports the presence of a rapidly rotating neutron star in this binary54. So the magnetar/Be binary used in this paper is close to LS I + 61 303. They have similar orbital period and eccentricity. Searching for FRBs from LS I + 61 303 is attractive. From the FAST observation, it has been found a rapid change of the plasma properties along the line of sight, which may be caused by the variation of the stellar wind properties54. Because the brightness temperatures of FRBs are much higher than that of pulsations of neutron star55, the induced Compton scattering may be significant41,56.

However, no FRB-like signal from the two systems has been reported. The possible reasons are as follows. First, FRBs are rare phenomena. They may be highly beamed57,58, and have narrow spectra, making them hard to discover. For example, more than 20 magnetars have been discovered in our Galaxy. Only FRB 200428 emitted from SGR 1935+2154 has been discovered59,60. A bright X-ray burst is associated with FRB 200428, which has a strikingly different spectrum from other X-ray bursts61,62. FAST observed SGR 1935+2154 for eight hours in its active phase57. 29 X-ray bursts are detected, but no FRBs coincident with the bursts. Second, the high brightness temperatures of FRBs require that the emission process mechanism must be coherent1,3,4. The physical conditions required to produce coherent radiation in magnetars may be hard to satisfy57. Third, the burst rate of FRBs deviates from Poissonian and varies significantly63,64,65. For FRB 20201124A, FAST found a sudden quenching of burst activity5. If similar active properties of FRBs are produced in magnetar/Be star binaries, enough and proper-time observations should be performed to search FRB signals from them.

Discussion

Considering the interaction between a Be star disk and FRBs, we can naturally explain the mysterious features of FRB 20201124A. This model can well produce the variation of RM with a magnetic field of 10 Gauss at the stellar surface. This magnetic field and the strong electron density fluctuations can generate the observed depolarization of the FRBs by differential Faraday rotation along different paths. The magnetic field perpendicular to the wave vector in the disk, having the same order of the magnetic field along the line of sight, can account for the observed circular polarization.

Recent observations show that the RM of FRB 20190520B is rapidly varying with two sign changes66,67, with similar magnitude as that of PSR B1259-63/LS 2883. The sinusoidal-like form of RM variation emerges. A dense magnetized screen near the FRB source can account for the RM variation and depolarization66. In our model, the distance of the screen can be of order the separation of the two stars. For orbital period P = 600 d, our model can well explain the RM variation of FRB 20190520B, which is shown in Fig. 5. This figure gives the total RM in the source frame. The RM contributed besides the disk is assumed to be RM\({}_{{{{{{{{\rm{ob}}}}}}}}}{(1+{z}_{{{{{{{{\rm{ob}}}}}}}}})}^{2}=-10301.92\) rad m−2, where RMob = − 6700 rad m−2 is the first observed RM by Green Bank Telescope (GBT) at MJD 5910666 and zob = 0.241 is the redshift of FRB 20190520B. The semimajor axis can be derived as \(a={(G{M}_{T}{P}^{2}/2\pi )}^{1/3}=8.5\) AU, where MT = 30M is the total mass of the binary. Our model predicts the RM evolution would be quasi-periodic. If a large number of RM detections spanning a long timescale are accumulated, the orbital period could be derived from these RM data.

Fig. 5: The fit of RM evolution for FRB 20190520B using the binary model.
figure 5

The blue points are source RM calculated from RM\({}_{{{{{{{{\rm{ob}}}}}}}}}\times {(1+{z}_{{{{{{{{\rm{ob}}}}}}}}})}^{2}\), where RMob is the observed value and zob = 0.241 is the redshift. The red line is the fit from the binary model. The disk density is ρ0 = 8 × 10−13 g/cm3 with a steep-decay index of − 3. The mass of the Be star is assumed to be 30M with radius R = 10R. The magnetic field at the stellar surface is 100 G. The orbital period is taken as 600 d with orbital eccentricity e = 0.75. Source data are provided as a Source Data file.

Methods

The disk model

There are some evidences for a disk around a Be star68 and especially for LS 288319. The disk opening angle is found to be less than 2069. The disk size ranges from a few stellar radius to a few tens stellar radius from interferometry observations35. The density of disk can be well described by a radial power-law form, with surface density between about 10−13 to a few times 10−10 g cm−335. From the DM and RM evolution of PSR B1259-63/LS 2883, ρ0 = 10−16 g cm−3 was used30.

The classical formation channel of a Be/X-ray binary (or a magnetar/Be star binary) is as follows. In an interacting binary, the B type star accretes matter from the companion star. Therefore, it is spun up to nearly critical velocity, leading to the formation of a decretion disk, which is supported by observations70. At this stage, the orbit and disk are almost in the same plane71. In a magnetar/Be star binary, the companion is a magnetar, which is usually formed by a supernova explosion. The process of supernova explosions is asymmetric and gives a very large kick to the newly formed magnetar31,72. Due to the kick of the magnetar, the orbital plane and the disk is misaligned. Besides, the kick also makes the orbit be eccentric33,73.

From observations, some Be/X-ray binaries have been found to be high misaligned33. For the Galactic binary PSR B1259-63/LS 2883, the inclination angle between the orbital plane and the disk plane are measured using several methods from observations. The inclination angle is found to be 10 to 40 by fitting to the observed variations in DM and RM of PSR B1259-6330. Using the X-ray data of PSR B1259-63, the derived inclination angel is about 7074. The inclination angle of PSR B1259-63/LS 2883 system is about 25 from the optical spectra of LS 288375. Using the 23 yr of pulsar timing data, the inclination angel is found to be 3518.

We consider the rotationally supported, geometrically thin disk model with isothermal hydrogen gas. The density profile of disk can be described as35

$$\rho (r,z)={\rho }_{0}{\left(\frac{r}{{R}_{\star }}\right)}^{-\beta }\exp [-{z}^{2}/{H}^{2}(r)],$$
(1)

where r is the distance from the star, ρ0 is the disk density at the stellar radius R, β is the density slope, z is the height, and H(r) is the scale height of the disk. Observational properties of Be stars are well described by this disk model and the density slope β is from 2 to 435.

The magnetic field, B, in the disk is poorly known. We assume that it is azimuthal and has a form

$${{{{{{{\bf{B(r)}}}}}}}}={{{{{{{{\bf{B}}}}}}}}}_{{{{{{{{\bf{0}}}}}}}}}({R}_{\star }/r),$$
(2)

with B0 is the magnetic field strength at the stellar surface. This model can explain the RM variation of PSR B1259-63/LS 288313. Given the above model, we can calculate the RM as a function of orbital phase. The RM contributed by the free electrons in the disk is

$${{{{{{{\rm{RM}}}}}}}}=8.1\times 1{0}^{5}{\int}_{{{{{{{{\rm{LOS}}}}}}}}}{n}_{e}{{{{{{{{\bf{B}}}}}}}}}_{\parallel }dl\,{{{{{{{\rm{rad\,{m}}}}}}}^{-2}}},$$
(3)

where ne is the free electron density in unit of cm−3, B is the magnetic field along the line of sight in unit of Gauss, andlis in unit of parsec. For the azimuthal magnetic field described above, we expect the RM contributed by the disk to change sign when the magnetar passes in front of the Be star, which is supported by FAST data.

The effect of tidal forces on the disk

The tidal force of the magnetar could affect the disk76. The process of disk truncation in Be/X-ray binaries was first discussed in 1997 from the orbital period-Hα equivalent width diagram77. From observations, 26 Be stars show disk truncation signature78. Among them, six Be stars are known to have close binary companions in circular orbits78. In the classical viscous disk model79, the criterion for disk truncation can be written as

$${T}_{{{{{{{{\rm{vis}}}}}}}}}+{T}_{{{{{{{{\rm{res}}}}}}}}}\le 0,$$
(4)

where Tvis and \({T}_{{{{{{{{\rm{res}}}}}}}}}\) are the viscous torque and the resonant torque, respectively. Using the standard formulae of Tvis and \({T}_{{{{{{{{\rm{res}}}}}}}}}\)80, \({T}_{{{{{{{{\rm{vis}}}}}}}}}+{T}_{{{{{{{{\rm{res}}}}}}}}}\le 0\) is met for a small value of α (viscosity parameter). The drift time-scale is τdrift = Δr/vr, where vr is the radial velocity, and Δr is the gap size between the truncation radius and the radius where the gravity by the magnetar begins to dominate. The spread of the disk can make the truncation inefficient for a system with τdrift/Porb = Δr/0.1csPorb < 1. For a typical value α = 0.1, τdrift/Porb = 0.1 in our model, which supports the truncation effect is inefficient. This conclusion was also confirmed by numerical simulations. For misaligned Be/X-ray binary systems, smoothed particle hydrodynamics simulations show that Be star disks are not truncated in highly eccentric binaries38. In addition, the disk may be warped by the tidal force38, leading to non-smooth distributions of density and magnetic field in the disk. This is the possible reason to cause the clumps in the disk. The spread of disk could cause the smaller density of the disk’s outer part, compared to last periastron. Therefore, the variation of RM will become smaller in magnitude over time. In our model, the Be/magnetar binary has a highly misaligned disk, eccentric orbit and ignorable disk truncation. These three requirements are self-consistent.

Besides the tidal force, the magnetar wind may also affect the disk. In our model, the radio bursts can be powered by the magnetic-energy release of the magnetar. The magnetar wind may be weak, similar to that of Galactic magnetars (i.e., SGR 1935+2154). Strong magnetar winds can produce a bright, compact persistent radio source associated with FRBs81,82, which has been observed in FRB 20121102A83 and FRB 20190520B84. However, there is no compact persistent radio emission associated FRB 20201124A10. The extended radio emission around FRB 20201124A is attributed to star-formation activity8,9,85.

Fitting RM using the binary model

Compared to the Galactic PSR B1259-63, the RM variation of FRB 20201124A is small. We consider a thin disk by assuming ρ0 = 3 × 10−14 g cm−3, and β = 4. The radius and mass of Be star is taken to be R = 5R and M = 8M. Typical velocities of the plasma in the disk are86 vd = 150−300 km s−1. The vertical scale height is derived from

$$H(r)={c}_{s}{\left(\frac{r}{G{M}_{\star }}\right)}^{1/2}r,$$
(5)

where cs is the sound speed. For an isothermal gas, the sound speed is \({c}_{s}={(kT/\mu {m}_{H})}^{1/2}\), where T is the isothermal temperature, μ is the mean molecular weight, and mH is the mass of hydrogen. We assume μ = 1 and T = 104 K. In order to prevent that the disk is significantly truncated by tidal forces, we consider a large eccentric orbit with e = 0.75. In our model, the RM variation is mainly affected by the orbital motion. We assume the orbit period is P = 80 days. Then the semi-major axis can be derived from

$${a}_{0}={\left[\frac{G({M}_{\star }+{M}_{NS}){P}^{2}}{4{\pi }^{2}}\right]}^{1/3},$$
(6)

where MNS = 1.4M is the magnetar mass. To match the observed RM evolution, we assume that the inclination angle is ϕ = π/12 and the observer’s direction is θ = 3.08 rad and ϕ = 2.18 rad. We assume the magnetic field B0 = 10 G. From observations, the surface magnetic field of a Be star has been measured. The magnetic field in the longitudinal direction is between 2 and 300 Gauss for 15 Be stars87. From a large sample of 60 Be stars, the maximal surface magnetic field 4000 Gauss was derived88. Therefore, the value of B0 = 10 G is well in the allowed range.

The DM contributed by the disk is

$${{{{{{{\rm{DM}}}}}}}}={\int}_{{{{{{{{\rm{LOS}}}}}}}}}{n}_{e}dl,$$
(7)

where ne is the free electron density along the line of sight. In the observations, FAST can detect FRB signals. Therefore, the optical depth

$$\tau =8.2\times 1{0}^{-2}{T}^{-1.35}{\nu }^{-2.1}\int\nolimits_{0}^{d}{n}_{e}^{2}dl$$
(8)

of free-free absorption must be less than unity. In the above equation, T is the temperature in Kelvin, ν is the observing frequency in GHz, and d is the distance to the magnetar (in parsec). The same parameters for the disk is used to calculate the DM and optical depth.

The induced Compton scattering is crucial for high brightness temperature sources. In our model, the scattering occurs in the nearest vicinity of the source. For a single narrow beam of uniform brightness and solid angle ΔΩ, the effective optical depth is41,56

$${\tau }_{{{{{{{{\rm{ind}}}}}}}}}=\frac{k{T}_{B}}{{m}_{e}{c}^{2}}({{\Delta }}{{{\Omega }}}^{2}){\tau }_{T},$$
(9)

where TB is the brightness temperature, and τT = ne(r)σTr is the Thompson scattering optical depth. For radio bursts with duration w and flux S occurring at redshift z = 0.098, the brightness temperature is

$${T}_{B}=S{d}^{2}/2k{(\nu w)}^{2} \sim 6.6\times 1{0}^{30}\,{{{{{{{\rm{K}}}}}}}}\,\left(\frac{S}{{{{{{{{\rm{Jy}}}}}}}}}\right){\left(\frac{\nu }{{{{{{{{\rm{GHz}}}}}}}}}\right)}^{-2}{\left(\frac{w}{{{{{{{{\rm{ms}}}}}}}}}\right)}^{-2}{\left(\frac{d}{450{{{{{{{\rm{Mpc}}}}}}}}}\right)}^{2}.$$
(10)

The radiation subtends a solid angle \({{\Delta }}{{\Omega }}={(cw)}^{2}/{r}_{d}^{2}=1{0}^{-10}\) with the distance rd = 50R between the magnetar and the disk. If we assume radio bursts travel at distance 2R away from the Be star, using the base density 5 × 10−14 g cm−3 of the disk and equation (9), we can estimate τind is 0.025, which is much less than one. Therefore, radio bursts can escape from the disk.

The effect of a stellar wind

A massive star usually has a large wind mass loss, which strongly depends on the stellar effective temperature89. The effect of a stellar wind on DM and free-free absorption should be considered. For an isotropic wind, the number density of the wind is given by

$${\rho }_{{{{{{{{\rm{wind}}}}}}}}}={\rho }_{0,{{{{{{{\rm{w}}}}}}}}ind}{(r/{R}_{\star })}^{-2},$$
(11)

where r is the radial distance from the star. The stellar surface density is \({\rho }_{0,{{{{{{{\rm{w}}}}}}}}ind}=\dot{M}/4\pi {R}_{\star }^{2}{v}_{w}\) with mass loss rate \(\dot{M}\) and wind velocity vw = 3 × 108 cm s−1. From observations, the mass loss rate of Be stars has been estimated to be90,91 10−11−10−8M yr−1. A typical value ρ0,wind = 4.18 × 10−17 g cm−3 is adopted, which corresponds to the mass loss rate 3 × 10−10M yr−1. The ratio of the mass flux density in the disk over the stellar wind is

$$\frac{{\rho }_{0}{v}_{d}}{{\rho }_{0,wind}{v}_{w}} \sim 80,$$
(12)

which is well in range of 30–10092.

Because of adiabatic cooling, the wind temperature decreases slowly with radial distance from the star. The temperature of the stellar wind follows a power-law evolution due to adiabatic cooling93

$${T}_{{{{{{{{\rm{wind}}}}}}}}}={T}_{0,{{{{{{{\rm{w}}}}}}}}ind}{(r/{R}_{\star })}^{-{\beta }_{1}}.$$
(13)

In our calculation, we use the effective temperature of star T0,wind = 3 × 104 K and β1 = 2/3. In this case, the maximum DMw and τw contributed by the stellar wind are 0.54 pc cm−3 and 0.005, respectively. Based on these results, the contribution of the stellar wind can be safely ignored. We are not aware of magnetic field measurements in the wind of a massive star. For the Galactic PSR B1259-63, the RM variation cannot be explained by the contribution from the stellar wind30. The geometry of the magnetic field in the stellar wind may be toroidal30,94. Therefore, the radio burst propagates orthogonally to the toroidal magnetic field in a significant fraction of the orbit, which contradicts the observed RM evolution.

Below, we compare the magnetic field in the pulsar wind and in stellar disk. Since the magnetic field is dipolar inside the magnetar light cylinder, rLC = 5 × 109Pm,s cm, and toroidal ( r−1) in the wind, the magnetic field strength at a distance r is

$${B}_{{{{{{{{\rm{wind}}}}}}}}}=230{B}_{14}{P}_{s}^{-3}\left(\frac{5\times 1{0}^{9}\,{{{{{{{\rm{cm}}}}}}}}}{r}\right)\,{{{{{{{\rm{G}}}}}}}},$$
(14)

for the period of magnetar Pm = 2 s. We make the comparison at the same radius r = 2R = 10R. Bwind = 0.9 G is less than the magnetic field in the disk. Therefore, the magnetic field in the disk dominates.

Depolarization in the disk

The observed pulse depolarization is attributed to random fluctuations in the magnetic field and/or electron density in the disk, which can cause differential Faraday rotation along different paths of the light rays95. If the Faraday dispersion function is well described by a Gaussian distribution with variance σ2 = k2δ(RM)2〉, the measured polarized intensity can be parameterized as

$$P({\lambda }^{2})={P}_{i}\exp [-2{k}^{2}\langle \delta {({{{{{{{\rm{R}}}}}}}}M)}^{2}\rangle {\lambda }^{4}],$$
(15)

where Pi is the intrinsic polarized intensity, λ is the observing wavelength, k = 0.81 is a constant, and \({\langle \delta {({{{{{{{\rm{R}}}}}}}}M)}^{2}\rangle }^{1/2}\) is the variance in ΔRM. For simplicity, we assume that \({\langle \delta {({{{{{{{\rm{R}}}}}}}}M)}^{2}\rangle }^{1/2}\) arises from fluctuations in free electron density ne. Furthermore, we assume that the variance of electron density ne satisfies \({\langle \delta {{{{{{{{{\rm{n}}}}}}}}}_{{{{{{{{\rm{e}}}}}}}}}}^{2}\rangle }^{1/2}\propto {n}_{e}\) in the disk30, which is similar to the small-scale turbulence in the interstellar medium. Letting \({\langle \delta {({{{{{{{\rm{R}}}}}}}}M)}^{2}\rangle }^{1/2}\) = fΔRM, the stochastic Faraday rotation leads to the depolarization when ΔRM > f−1λ−2. From the FAST observation, the value of ΔRM is about 200 rad m−2 at frequency 1.25 GHz. Therefore, f is about 0.1. We conclude that stochastic Faraday rotation may be responsible for the observed depolarization.

Scattering timescale

From the CHIME observation, the largest scattering timescale of FRB 20201124A is 14 ms7. The fluctuations of the electron density in the disk can cause scatter broadening of the pulses. The electron plasma with density fluctuations of δne and spatial scale a produces a typical scattering angle at wavelength λ of

$${\theta }_{s}=\frac{1}{2\pi }{\left(\frac{L}{a}\right)}^{1/2}{r}_{0}{\lambda }^{2}\delta {{{{{{{{\rm{n}}}}}}}}}_{{{{{{{{\rm{e}}}}}}}}},$$
(16)

where r0 is the classical electron radius and thickness of the screen L is assumed to be roughly equal to its distance d from the magnetar96. If the scattering is due to a thin screen at a distance rs from the magnetar, the scattering angle θs is \({(2{\tau }_{s}c/{r}_{s})}^{1/2}\), with the scattering time τs. We use the burst 20210405B with scattering time τs = 12.17 ms detected by CHIME7 to estimate the spatial scale a. The distance between the burst and the screen is about 42 solar radius. Therefore, the scattering angle θs is about 0.9. From equation (16), the density fluctuation of electron plasma in the disk is δne = 80a1/2 cm−3. If δne is set to be equal to the total electron density ne in the scattering region, from the disk model, then δne = 3.9 × 106 cm−3 is found, implying that a = 2.4 × 109 cm.

Origin of circular polarization

FRB 20201124A is the first repeating FRB with significant circular polarization as discovered by FAST5 and Parkes97. Two possible mechanisms can cause Faraday conversion: linearly polarized light is converted to circularly polarized one and vice versa. The polarization measurement of FRB 121102 has been used to constrain the local magneto-ionic environment using Faraday conversion98,99.

First, we consider that the circular polarization arises from the propagation through mildly relativistic plasma. For the PSR B1259-63/LS 2883, X-ray, GeV gamma-ray and Tev gamma-ray radiations have been observed. The natural physical mechanism for these high-energy emissions are as follows. An interaction between the relativistic wind of the pulsar and the low-velocity wind from the Be star can form a termination shock. The electron and positron pairs from the pulsar wind are accelerated to relativistic velocities by the terminal shock and emit broadband non-thermal emission through synchrotron radiation and inverse Compton scattering. In our model, a termination shock may not exist, due to the weak magnetar wind.

If only the relativistic plasma is considered, the change of polarization angle is

$${{\Delta }}\phi =\frac{1}{2}\,{{{{{{{\rm{RRM}}}}}}}}\,{\lambda }^{3},$$
(17)

where RRM is the relativistic rotation measure (RRM), and λ is the observed wavelength. We assume a power-law distribution for accelerated particles

$$N(\gamma )={N}_{0}{\gamma }^{-p},\,\,{\gamma }_{\min }\le \gamma \le {\gamma }_{\max },$$
(18)

where \({\gamma }_{\min }\) and \({\gamma }_{\max }\) are the minimum and maximum Lorentz factors, respectively. For the power-law index p ≠ 1, the normalization factor is \({N}_{0}={n}_{r}(p-1){({\gamma }_{1}^{1-p}-{\gamma }_{2}^{1-p})}^{-1}\) with the number density nr of the relativistic particles. The definition of RRM is100

$${{{{{{{\rm{RRM}}}}}}}}=\frac{{e}^{4}}{4{\pi }^{3}{\epsilon }_{0}{m}_{e}^{3}{c}^{4}}\frac{p-1}{p-2}\int {n}_{r}{({{{{{{{\bf{B}}}}}}}}\sin \theta )}^{2}{\gamma }_{\min }ds.$$
(19)

The fraction of circular polarization is Πc = RRMλ3 for the relativistic plasma dominated case. For relativistic plasma, the definition of DMr is DM\(\times {{{{{{{\mathcal{D}}}}}}}}\)101, with Doppler factor \({{{{{{{\mathcal{D}}}}}}}}={[\gamma (1-v/c\cos \theta )]}^{-1}\), where θ is the angle between the line of sight and the direction of plasma motion. The RRM effect should require \({n}_{e}\ll {n}_{r}{\gamma }_{1}\ln (\nu /{f}_{c}{\gamma }_{1})\).

For the 1994 periastron passage of PSR B1259-63, the changes of DM and RM cannot be explained in the termination shock model30. The observed changes of DM in four periastron passages of PSR B1259-63 are roughly consistent with the contributions from the disk and wind of Be star, with the disk component dominated102. Therefore, the cold plasma dominates the wave properties. The RRM effect is very weak. When the magnetar is in front of the Be star, radio bursts cannot interact with the plasma shocked by the termination shock. Therefore, there is no Faraday conversion, which is disfavored by the FAST observations.

Numerically, the approximated Faraday conversion coefficient is43,103

$${\rho }_{Q}=8.5\times 1{0}^{-3}\frac{2}{p-2}\left[{\left(\frac{\omega }{{f}_{c}\sin \theta {\gamma }_{\min }^{2}}\right)}^{(p-2)/2}-1\right]\frac{p-1}{{\gamma }_{\min }^{1-p}}{\left(\frac{{f}_{c}\sin \theta }{\omega }\right)}^{\frac{p+2}{2}}\frac{2\pi }{\omega },$$
(20)

where ω = 2πν, fc = eB/mec is the cyclotron frequency, and θ is the angle between wave vector and the line of sight. This approximation approaches zero if \({\gamma }_{\min }\ge 200\). In order to explain the high-energy emission from the Galactic binary system PSR B1259-63/LS 2883, an extreme large \({\gamma }_{\min }=1{0}^{5}\) is required45,102. Therefore, the relativistic plasma hardly generates the observed circular polarization, unless the minimum Lorentz factor is less than 100.

Next we consider the Faraday conversion occurring when the magnetic field is perpendicular to the wave propagation direction. The required magnetic field for Faraday conversion is about 7 G5, which is the same order of the magnetic field strength along the line of sight. The Be star disk can provide the required vertical magnetic field. This is the most natural way to explain the production of circular polarization.