The Relativistic Spin Precession in the Compact Double Neutron Star System PSR J1946+2052

We observe systematic profile changes in the visible pulsar of the compact double neutron star system PSR J1946+2052 using observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The interpulse of PSR J1946+2052 changed from a single-peak to a double-peak shape from 2018–2021. We attribute this evolution as the result of the relativistic spin precession of the pulsar. With the high sensitivity of FAST, we also measure significant polarization for the first time, allowing us to model this with the precessional rotating vector model. Assuming, to the first order, a circular hollow-cone-like emission beam pattern and taking the validity of general relativity (GR), we derive the binary’s orbital inclination angle ( 63∘−3∘+5∘ ) and pulsar’s spin geometry. The pulsar’s spin vector and the orbital angular momentum vector are found to be only slightly misaligned ( 0.∘21−0.∘10+0.∘28 ). The quoted uncertainties do not reflect the systematic uncertainties introduced by our model assumptions. By simulating future observations of profile and polarization evolution, we estimate that we could constrain the precession rate within a 43% uncertainty in 9 yr. Hence, we suggest that the system’s profile evolution could be combined with precise pulsar timing to test GR in the future.

Relativistic spin precession is another relativistic effect caused by the curvature of space-time in the presence of a large mass object.This precession causes the pulsar's spin vector to precess around the total angular momentum vector and moves the pulsar's emission beam with respect to our line of sight.The relativistic spin precession is usually not directly detectable through pulsar timing, but sometimes can be observed through the temporal evolution of the pulse profile and its polarization.
The first indication for profile changes caused by relativistic spin precession was given by Weisberg et al. (1989), observing the relative amplitude of the two prominent components of PSR B1913+16's profile changing over time.A systematic decrease in the pulse width change was later detected, and the profile changes were analyzed and modelled by Kramer (1998).A quantitative test of the precession rate was not possible, but the changes were found to be consistent with the expectation from GR.After PSR B1913+16, the relativistic spin precession has been detected through profile change in more pulsar binary systems, such as PSR J1537+1155 (Stairs et al. 2004), PSR J1141−6545 (Hotan et al. 2005), PSR J0737−3039B (Burgay et al. 2005), PSR J1906+0746 (Lorimer et al. 2006;Desvignes et al. 2019) and, more recently, in PSR J1757−1854 (Cameron et al. 2023).Some pulsars in DNSs, such as PSR J0737−3039A, show no sign of temporal profile evolution.This could be due to a small misalignment angle between the pulsar's spin vector and the orbital angular momentum, as inferred for PSR J0737−3039A (Ferdman et al. 2013).As the polarization position angle (PPA) swing is a marker of pulsar geometry, the temporal evolution of the PPA is also a potential (and powerful) indicator of relativistic spin precession (Stairs et al. 2004;Desvignes et al. 2019).
Assuming GR, the rate of spin precession of the pulsar is determined from the system's Keplerian parameters and two objects' masses (Barker & O'Connell 1975): (1) where G is the gravitational constant, c is the velocity of light, P b is the orbital period, e is the orbital eccentricity, m 1 and m 2 are masses of the pulsar and companion measured in units of solar mass.The Keplerian parameters such as P b and e can be derived through pulsar timing, and masses can be measured if more than two post-Keplerian parameters are detected.If the spin precession rate is derived independently through the profile and PPA evolution, then it can be used to test GR (Stairs et al. 2004;Fonseca et al. 2014;Desvignes et al. 2019).In the case of PSR J1906+0746, the measured spin precession rate agrees with GR to within 5% (Desvignes et al. 2019).In the particular case of PSR J0737−3039A/B, the spin precession rate is determined through modelling pulsar B's eclipse morphology evolution and that leads to a rate in agreement with GR to within 13% in Breton et al. (2008) and recently 6% in Lower et al. (2024).At last, if the spin precession rate could not be directly measured, one still can use the relativistic spin precession to model the binary system's geometry by assuming GR is correct (Kramer 1998;Manchester et al. 2010;Ferdman et al. 2013).Stovall et al. (2018) discovered the DNS system PSR J1946+2052, which has the shortest orbital period among known DNS systems at 1.88 h.The eccentricity of the orbit is 0.064.They detected only one post-Keplerian parameter, the periastron advance, to be 25.6 • yr −1 , indicating the system's total mass is 2.50(4) M ⊙ under the validity of GR.The pulsar was considered to be the primary neutron star of the system due to the spin period of 16.9 ms and the magnetic field strength at the surface of 4 × 10 9 G.They suggested that this system resembles the double pulsar system PSR J0737-3039A/B in configuration and binary evolution.Therefore, they proposed that the pulsar's spin vector could be nearly aligned with the orbital angular momentum vector and thus spin precession should be hard to observe.In their work, they did not detect any polarization.
The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is a single-dish radio telescope (Nan 2008;Nan et al. 2011) and is operating since 2019 (Jiang et al. 2020).Owing to the high sensitivity of FAST, we report here the detection of polarization and significant profile changes for PSR J1946+2052 for the first time.In section 2, we describe the data reduction and specifically demonstrate the removal of the ionospheric Faraday rotation effect.In section 3, we show the results of our analysis of pulse profile and polarization and model the geometry of this pulsar system.In section 4 we discuss the evolution of this DNS and predict how the measurement of spin precession rate would advance in future observations.

OBSERVATION AND DATA PROCESSING
We observed PSR J1946+2052 with the centre beam of FAST's 19-beam receiver at a central frequency of 1250MHz using a bandwidth of 500MHz for 11 epochs.The first session, conducted on March 30, 2019, was observed in 2 polarization (AA, BB) mode, while the remaining sessions were carried out in full Stokes mode in 2021 and 2022.Details of the length and date of observations are shown in Table 1.We recorded all data in PSRFITS (Hanisch et al. 2001) format with a sampling time of 49.152 µs and 8-bit digitization.To calibrate each observation's polarization, we pointed the telescope off the pulsar and turned on the noise diode in a 0.2-s cycle for 1 minute before observing the pulsar.We used the DSPSR (van Straten & Bailes 2011) analysis program and the PSRCHIVE (Hotan et al. 2004) software package to fold and de-disperse the pulsar data, to eliminate radio frequency interference (RFI), and to calibrate the polarization.We obtained the ephemeris from PSRCAT 1 (Manchester et al. 2005) and updated it with our own time-of-arrival (TOA) measurements using TEMPO2 (Nice et al. 2015).With the local solution of the ephemeris derived from TEMPO, for each epoch, we folded the data and scrunched time.To avoid effects due to uncertainties in the long-term ephemeris, we aligned pulses from different observations by using the maximum of the normalized pulse profile.

Spin phase
Figure 1.This figure presents the pulse profile and polarization information.In the bottom panels, we present the profile of each observation's relative intensity, linear polarization, and circular polarization with black, red, and blue lines respectively.PPAs and precessional RVM curves are shown in the top panels.Data from 2019-03-30 do not have polarization information and position angle because they are formed by 2 polarization.The grey points indicate the part of PPA that we consider to have orthogonal polarization mode.
was measured using rmfit in PSRCHIVE.However, due to the low Signal-to-Noise Ratio (SNR) and low degree of linear polarization, we could not derive a precise RM for the whole set of observations.Instead, we used the observation with the highest SNR to obtain RM=−27(11) rad m −2 and removed the Faraday rotation effect in polarization with a fixed RM −27 rad m −2 to analyze the relative change of the geometrical parameters (see Sec. 3).Then we computed the Stokes parameters I, Q, U, and V and exported the archive data into a numpy (Perez & Granger 2007) readable format using the Python package of PSRCHIVE.
RM is a crucial parameter in studying pulsar polarization.For example, 1 rad m −2 shifts the PPA by almost 14 • at 1250 MHz.As explained in Sotomayor-Beltran et al. ( 2013): where RM int , RM ion , RM ism , RM igm denote the intrinsic RM and RM of the ionosphere, interstellar medium, intergalactic medium respectively, where the latter is not applicable in this case.We can also expect a change in RM ism to be small, based on typical line-of-sights with a relatively uniform magnetic field along those and conservative estimates of a δDM of approximately 10 −4 pc cm −3 .Hence, we consider mainly the influence of ionospheric rotation measure (RM ion ) and its variation on the measurement of the PPA.We calculated RM ion of each epoch by using the ionFR package (Sotomayor-Beltran et al. 2013) and global ionospheric map (GIM) products.The GIM products were generated by different analysis centres using total electron content (TEC) gathered from global positioning system (GPS) stations around the world, combined with a specific mathematical model (Li et al. 2017).For each observation, dozens of GIM analysis centres can provide measurements of the ionospheric RM, and their estimates may differ slightly due to model dependence.Figure 2 shows an example of the RM ion of PSR J1946+2052 derived from GIM products on March 9, 2021, which are available on NASA's website. 3To remove the bias in the GIM products, we interpolated values of RM ion and their errors during the time span of our observations and obtained the weighted average and error for each analysis centre.Then we weighted these RM ion from different analysis centres with their errors and derived the average and the uncertainty of each day's ionospheric RM.The RM ion values for each observation epoch are presented in Table 1.
We have detected polarization in observations with FAST lasting longer than one hour.The profiles at 1.25 GHz are presented in Fig. 1 with 1024 phase bins.We consider the emission feature near pulse phases ∼ 0.36 − 0.40 to be the interpulse, complementing the stronger main pulse at pulse phases ∼ 0.77 − 0.93.The observation in 2019 has no polarization information so we show only total intensity (black line) for this epoch in Fig. 1.For the other epochs, we show linear polarization (red line), circular polarization (blue line), and the measured PPAs in the upper part of the plots.3. RESULTS As shown in Fig. 1, we have observed significant changes in the pulse profiles of PSR J1946+2052 over a 3-year time span, which can be attributed to relativistic spin precession.The changes manifest themselves most prominently in a change in the relative amplitudes of both interpulse and main pulse over time.To obtain a smooth pulse profile with less noise, we downsample the spin phase to 256 bins to extract the double-peak structure of both the interpulse and main pulse.We normalize each integrated profile to the peak of the main pulse's leading component, and the relative amplitudes of other components are shown in panels (b) and (c) of Fig. 3.The uncertainties in the relative amplitudes are calculated as the root mean squares of each profile's off-pulse region.
Relativistic spin precession can cause changes in the pulse profile and the PPAs.To reveal these effects, we analyze the full polarization profile of PSR J1946+2052.We use a simple circular hollow-cone-shape emission beam to model the profile change.For polarization, the PPA swings of PSR J1946+2052 are not perfect S-shape, and the points shown in grey are not easily consistent with a typical RVM swing.They may be caused by an unresolved orthogonal jump, as their location coincides with a minimum in linear polarisation and a change in the sign of circular polarisation.In our analysis, we then ignore this small range of PA values (displayed as grey points in Fig. 1).We present details of our analysis in the following sections.The most important profile evolution is that the interpulse has only one peak on 2019 March 30 but then becomes doublepeaked in observations taken after 2021 January 13.Assuming a circular hollow-cone emission beam, this profile variation can be modelled by the movement of our line of sight from the edge of the emission beam to the interior of the hollow cone.This suggests that initially, the line of sight was near the edge of the beam, but our assumption of the beam shape may be wrong (cf.Desvignes et al. 2019).
To estimate the separations between the two components of the interpulse and the main pulse, we use two Gaussian functions with identical widths to fit these two pulse components separately in each observation.The component separation is determined as the distance between the centres of the two Gaussians.The separation between the interpulse and the main pulse does not enter the analysis at this stage.For the interpulse seen on 2019-03-30 with only one peak, this method cannot be applied.Here, we fix the widths of the two Gaussian functions to the average value obtained from later observations.This fixed width is consistent with the distribution of the widths we derive from the 2021-2022 epochs.We also make the two Gaussian functions have the same amplitude.With these additional constraints, we can obtain an estimated peak component separation for the 2019 observation.We display the Gaussian functions that are used to fit the interpulse and the main pulse in panel (a) of Fig. 3.The peak component separations obtained using these methods are shown in panels (d) and (e) of Fig. 3. Figure 4 displays the stacked PPA measurements obtained from PSR J1946+2052, which shows no clear evidence of evolution.The relatively large uncertainties and the absence of significant changes in the swing shape suggest that the expected spin precession has not yet reached a level of detectability through PPA measurements.This may be caused by the geometrical configuration, a special precession phase, or oversized uncertainties in our measurements.

Convention and definitions
In this study, we adopt the PSR/IEEE convention (van Straten et al. 2010) to define the geometrical angles and their relationships, which is commonly referred to as "the observer's convention".These angles are illustrated in Fig. 5.It should be noted that this convention differs from the convention used in previous studies such as Damour & Taylor (1992) and Kramer & Wex (2009) (hereafter referred to as the DT convention) in which the line of sight is defined as pointing out from the sky plane.The angle between the line of sight and the spin vector (λ) and the inclination angle of the orbit (i) are complementary in these two conventions.Apart from those changes in convention, the formalism below follows Kramer & Wex (2009).
As the spin angular momentum of the pulsar is much smaller than its total angular momentum, we can assume that the direction of the orbital angular momentum is constant with time and is the same as the total angular momentum.This means that the inclination angle i remains constant with time.As the pulsar's spin vector precesses around the orbital angular momentum, the line of sight will cut through different latitudes of the emission beam, causing the angle λ to change with time according to the following function: where δ is the misalignment angle of the spin vector from the orbital angular momentum vector.Φ SO (t) is the precession phase and is described as: where Φ 0 is the reference precession phase, T 0 is the time when the precession phase is Φ 0 and Ω SO is the rate of the precession.In our work, we define T 0 as the date of the first observation (2019-03-30).
The projected position angle of the spin vector (Ψ 0 ) on the celestial sphere changes due to spin precession.In our convention, the position angle equals 0 in the North of the sky plane (I 0 ) and increases counter-clockwise as turning to the East (J 0 ).Ψ 0 is described as where Ω asc is the ascending node of the orbital plane, ∆ A is the angle shift caused by aberration, ∆ F is the change of the Faraday rotation effect which is illustrated in Eq. ( 2), and η(t) is the longitude of the pulsar's spin precession.Because it is difficult to distinguish the specific values of Ω asc , ∆ψ A and ∆ F , we simply define Ψ 0 (t) as and η(t) can be derived from

Modeling pulse profile separations and PPAs
We can model the separation evolution by assuming a circular hollow-cone-shape emission beam.For the dipole field emission, the separation between the interpulse and the main pulse could deviate from 180 • due to the different emission heights, while the shape of these two pulses could be the hollow-cone shape.We assume each hollow-cone beam to be symmetrical so that the separation between two peak components in one pulse can be expressed as (Kramer 1998): where W(t) is the separation, ρ is the opening angle of the emission beam, α is the magnetic inclination angle and β(t) is the angle between the magnetic axis and the line of sight.β(t) is given by: Once we determine α for the main pulse, the interpulse's magnetic inclination angle is given by π − α.We have observed the interpulse changing from a one-peak to a two-peak shape, but the short time span covered by our observations alone is insufficient to constrain the six geometrical parameters that we need to model the precession.
For modelling the PPAs, we continue to adopt the precessional RVM (Kramer & Wex 2009) that considers relativistic spin precession.The conventional RVM is used to model the PPA that increases clockwise.However, the PPA we use increases counter-clockwise as we derive the Stokes parameters from PSRCHIVE, which follows the PSR/IEEE convention (van Straten et al. 2010).To modify the RVM equation accordingly, we refer to Everett & Weisberg (2001).As Ψ 0 (t) and β(t) vary with time due to precession, the precessional RVM can be described by Eq. ( 11), where ϕ is the spin phase, ϕ 0 is the reference phase of the spin axis and is relative to the main pulse.Due to the limited SNR of our polarization data, we cannot consider the difference in emission height between the interpulse and the main pulse through PPAs, as done in Johnston & Kramer (2019).
Although the reference phase of the spin axis position angle is not expected to change due to relativistic spin precession, we fit ϕ 0 for every epoch because the profiles are aligned based on the maximum of the main pulse.We exclude the part of the PPA where apparent orthogonal polarization modes are present.The selected PPA points are shown in Fig. 1 as black points.
In order to constrain the geometry using profile separations and PPAs, we combine the profile separation model and precessional RVM and reconstruct the likelihood as given by Eq. ( 11). .We illustrate the geometrical angles in our convention.I0 represents the North of the sky plane, and J0 represents the East.The fundamental reference frame consists of I0, J0, and K0. S denotes the direction of the pulsar's spin vector, Sp is the direction of the projection of the spin vector, and k is the direction of the orbital angular momentum.The orbital reference frame consists of i, j, and k.The magnetic axis m, the main pulse's opening angle ρmain, and the interpulse's opening angle ρinter describe the emission beams.The definition of other angles can be found in the text.
where i and j indicate each day and each PPA point in day i's observation.As shown in Eq. 11, we weigh each component of the likelihood with the number of data points as follows: where C inter , C main and C PPA,i are the weighted coefficients we used for the component separations of the interpulse, the main pulse and each day's PPA swing.N inter is the number of data points of the interpulse's separation.Note that N inter = N main because we obtain one data point for the component separation of the main pulse and interpulse in each observation.N PPA,i indicates the number of data points of each day's PPA swing.This likelihood enables us to fit the profile separations with RVM information of β and ensures that β < ρ.We utilize the Markov Chain Monte Carlo (MCMC) method based on the Python packages emcee (Foreman-Mackey et al. 2013) and corner (Foreman-Mackey 2016) to obtain the distribution of geometrical parameters.
Since we have no prior information on the geometry, we search for the optimal geometry parameters with uniform priors over the allowed range while assuming that GR is correct.The individual masses of the binary system are not accurately determined yet.Still, we can use the mass function (0.268184(12) M ⊙ ) and total mass (2.50(4) M ⊙ ) of this sys-tem from Stovall et al. (2018), to express the precession rate as a function of the orbital inclination angle i.
We display the MCMC result of geometrical angles in Fig. 6 and ϕ 0 of each epoch in Fig. 7.We present in Fig. 4 the best fit of the PPA swings, while the best fit of separations and 1-σ samples are displayed in Fig. 3. Our results show that δ is very small, and is consistent with the expectation by Stovall et al. (2018).The value of α is related to the main pulse and indicates that the pulsar is a nearly orthogonal rotator.The observed significant difference in profile separations between the interpulse and the main pulse is because the spin axis and the magnetic axis are not entirely orthogonal.Given the small value of δ, we cannot derive a precise precession phase Φ 0 .Assuming negligible effects from additional Faraday rotation and aberration, the longitude of the ascending node Ω can be determined to a small range.The inclination angle i has been constrained by the fit to i = 63 •+5 • −3 • , which is consistent with the updated timing result (Freire, private communication) where i is estimated to be between 65 • and 80 • .Based on our estimate of i, we can predict the companion mass and the precession rate to be 1.33(4) M ⊙ and 7.96(23) • yr −1 respectively.Thus the mass of PSR J1946+2052 is predicted to be 1.17(6) M ⊙ .The distribution of these two parameters is shown in Fig. 8.By calculating β with the best fit geometrical parameters (the dark red curve in Fig. 9), we predict that the separation of the interpulse will disappear in 2032 because our = 86.68+1.96 2.30 line of sight will move out of the hollow-cone shape emission beam.While due to the small amplitude of β, the separation of the main pulse will stay visible all the time.We display the evolution of β in Fig. 9.With improved data and longer time baselines, in the future, we may be able to combine the timing analysis with profile and polarization studies to derive the geometry of this system precisely and even test GR.PSR J0737-3039A, thus the variation of the peak separation is significant.The trend of the separation change is different for the interpulse and the main pulse.The separation of the interpulse is increasing, while that of the main pulse appears to be decreasing.The variation in the interpulse from one peak to two peaks suggests that our line of sight is moving towards the magnetic axis of the interpulse if the beam is a hollow cone.In this work, we make the assumption that the two radio beams of the pulsar are symmetrical.Based on our fitting result, the line of sight is near the edge of two emission beams and moving toward the centre of the beams.Therefore we would expect the separation between the main pulse peaks to increase as well, which is, however, the opposite.This could be due to the complexity of the emission pattern in the interior of the main pulse's emission beam, as shown in PSR J1906+0746 (Desvignes et al. 2019).
In this work, we have assumed the validity of GR and used a simple beam model to constrain binary parameters.To determine when we can utilize PSR J1946+2052 for testing GR with spin precession, we conduct simulations of future observations of profile separations and PPA swings, assuming that our current best-estimated geometry is accurate.We produce a set of simulated profile peak separations and PPA swings every two months, based on the projections from our bestfit parameters, and incorporate the precession rate we derive from these parameters.To account for measurement errors on the actual profile peak separations and position angle swings, we introduce noise to our projected data.The error of each data point is determined by using the bootstrap method among corresponding existing errors.As the number of new simulated PPA swings grows faster than the number of simulated profile component separations, PPA swings will dominate the global fit.To counterbalance this discrepancy, we resample the simulated PPA swings and average the values for each year.This should not impact our fit outcome, as the position angle swing varies very slowly with time.
We perform simulations for additional 3-year, 6-year, and 9-year data sets, and present the results in Table 2, except ϕ 0 of each epoch.We fit for the value of Ω SO without assuming GR and find that it can be constrained to a relative precision of only 43% in 9 years (total of 11 years).However, these simulations do not incorporate pulsar timing.With 11-year FAST timing, we could possibly derive the DNS's masses and inclination angle i as well as potentially the system's longitude angle of the ascending node Ω asc .The timing measurement could be used as the prior in our analysis of the profile and PPA, significantly improving the constraint on the other geometrical parameters and the spin precession rate under the assumption of GR.This improvement is important for the pulsar emission model and binary evolution.
In summary, we observed PSR J1946+2052 for 11 times and find significant profile evolution and attribute this evolution to relativistic spin precession.We use the circular hollow-cone emission beam to derive the geometry of this system.We derive a global fit model by combining the separation fit and PPA fit to estimate the geometrical angles.In the fit, we follow the results of PSR J1946+2052's mass function and total mass in Stovall et al. (2018) and assume GR to be correct.The result shows that the spin vector of this pulsar is nearly aligned with the orbital angular momentum, indicating that the second SN should be fast and symmetric, as suggested by Stovall et al. (2018).We also constrain the companion mass and the precession rate under the scenario of GR.We emphasize that the derived uncertainties do not take into account possible systematic shortcomings of the assumed beam model.The future timing of PSR J1946+2052 will improve our understanding of the pulsar emission model as well as binary evolution and constrain the spin precession rate under the assumption of GR.

Figure 3 .
Figure 3.We present the down-sampled pulse profile and the evolution in relative amplitude and separation of each pulse, providing insight into the temporal changes of PSR J1946+2052's pulse profile.(a) The intensity profile of each epoch.The date of the observation goes from bottom to top.The dashed light coral and sea blue lines represent the Gaussian components that we derive from the fits of the pulse profiles.(b) The evolution of the interpulse's relative amplitude.(c) The evolution of the main pulse's relative amplitude.Only trailing components are shown because the leading components are used for normalization.(d) The evolution of the interpulse's separation between two components.(e) The evolution of the main pulse's separation between two components.The black solid and light coral lines in panels (d) and (e) indicate the best fit of our global fit model and samples in the 1 − σ confidence interval.(see details in Sec.3.2.2).

Figure 4 .
Figure 4.In the upper panel, we display the PPA measurements at different epochs.The solid lines are the best-fit curves of the global fit model.In the lower panel, we display the residuals of the best-fit PPA swings based on the first curve on 2021-03-09 and also the residual between the PPA measurements and the best-fit PPA curves.We complete the points where the residual would jump for 180 • .
Figure 5.We illustrate the geometrical angles in our convention.I0 represents the North of the sky plane, and J0 represents the East.The fundamental reference frame consists of I0, J0, and K0. S denotes the direction of the pulsar's spin vector, Sp is the direction of the projection of the spin vector, and k is the direction of the orbital angular momentum.The orbital reference frame consists of i, j, and k.The magnetic axis m, the main pulse's opening angle ρmain, and the interpulse's opening angle ρinter describe the emission beams.The definition of other angles can be found in the text.

Figure 6 .
Figure 6.The posterior distribution of the geometrical parameters in the global fit.The angles are given in units of degree.The value of each parameter indicates the result in a 68% confidence level.

Figure 9 .
Figure9.The evolution of β derived from the best fit geometrical parameters of PSR J1946+2052.For the upper panel, the blue solid and dashed lines indicate the best fit ρ of the main pulse and its 1-σ interval.The dark blue and translucent blue curves indicate the best fit β of the main pulse and samples in the 1-σ interval.For the lower panel, the red solid and dashed lines indicate the best fit ρ of the interpulse and its 1-σ interval.The dark red and translucent red curves indicate the best fit β of the interpulse and samples in the 1-σ interval.The grey area is the time span of our observations.The red dotted line in the lower panel indicates the day when the separation of the interpulse disappears, which is predicted to be around 2032.

Table 1 .
Date, length, and ionosphere RM of each observation RMion values on 2021-03-09 derived from different GIM analysis centres are indicated by various colours.The blank parts without points indicate the time when PSR J1946+2052 is not visible to FAST.The grey area represents the time span of our observation.
Everett & Weisberg (2001)As, we consider the unbiased linear polarization L true as detailed inEverett & Weisberg (2001).As scintillation modifies the SNR from epoch to epoch, the number of well-defined PPA points varies accordingly.The shown PPA values are corrected for Faraday rotation due to the ionospheric RM ion contribution by subtracting λ 2 × RM ion from the PPA values (measured in radians), where λ is the observing wavelength.

Table 2 .
Geometrical parameters in simulations