![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-88441208-56ff-4dbc-8185-d74cf7b51516.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-cb489ed6-7846-4233-bcc2-b93d99ef7170.jpg)

## Results from the long-term Fermi-LAT observations of the white dwarf binary pulsar AR Scorpii

## Q. Kaplan , P. J. Meintjes , K. K. Singh 1 1 1 2 , , H. J. van Heerden , F. A. Ramamonjisoa , 1 1 1 ∗

- I. P. van der Westhuizen
- 1.) Department of Physics, University of the Free State, PO Box 339, Bloemfontein, 9300
- 2.) Astrophysical Sciences Division, Bhaba Atomic Research Centre, Mumbai, India, 400085

E-mail: KaplanQ@ufs.ac.za

The discovery of the white dwarf binary system AR Scorpii (AR Sco) with its fascinating nonthermal dominated multi-frequency emission has sparked renewed interest in potential high energy gamma-ray emission from white dwarf pulsars. The Spectral Energy Distribution (SED) below and above optical shows evidence of non-thermal synchrotron emission, with pulsed emission in optical and X-ray bands at the white dwarf spin period ( P ∗ = 117s) as well as a beat period ( P b = 118 1s) with the binary period. . From an energy perspective, the highly magnetic rotating white dwarf can accelerate particles to TeV energies. In this study, a search for high energy gamma-ray emission was conducted between 100 MeV - 500 GeV by analysing the newly available Fermi-LAT Pass 8 data with the new Fermi 1.0.1 Science Tools. Binned likelihood analysis was done using power law, broken power law and log parabola models. From the selected Region of Interest (ROI) centred on AR Sco's position, we calculated a significance of √ TS ≤ 3.87 σ for the integrated gamma-ray activity between 100 MeV - 500 GeV at a photon flux level of 0.486 ± 0.261 x 10 -8 photons cm -2 s -1 using the broken power law model. This resulted in a 3 σ upper-limit detection from the position of AR Sco. The location of AR Sco inside the Rho Ophiuchi (Rho Oph) molecular cloud complex combined with the poor spatial resolution of Fermi-LAT, complicates any positive identification of low-level gamma-ray activity at the location that coincides with the position of AR Sco.

36th International Cosmic Ray Conference - ICRC2019 July 24th - August 1st, 2019 Madison, Wisconsin, USA

∗ Speaker.

http://pos.sissa.it/

## 1. Introduction

Since the discovery of the peculiar close binary system AR Scorpii (AR Sco) in 2016, ([1]) which displays non-thermal dominated multi-frequency emission from radio to X-rays (see Fig. 1), it has been the focus of intensive observational and theoretical investigation. This system consists of a 0.8 - 1.29 M /circledot highly magnetic white dwarf with a surface field B ≤ 500MG (e.g. [2, 3]) orbiting an 0.28-0.45 M /circledot M5 dwarf star with an orbital period of P orb = 3 56hours. Assuming the . secondary star fills its Roche lobe, the distance to the system is estimated to be approximately d ∼ 116 pc [1].

Figure 1: The multi-frequency Spectral Energy Distribution (SED) of AR Scorpii from radio to X-ray frequencies. The red and blue curves represent model atmospheres of the red M5 secondary dwarf star (T = 3100 K) and white dwarf (T = 9750 K) respectively. Adopted from [1].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-2a174709-d719-4a4f-8617-a84270e85050.jpg)

In the optical part of the spectrum (see [1]) the system reveals strong brightness variations every 117 s (i.e. pulsations), which is the spin period of the white dwarf, superimposed on the orbital period as well as strong pulsations at 118.2 s, which is a beat period between the spin period and the binary orbital period (see Fig. 2). Both the spin and beat periods are clearly visible on a power spectrum (see Fig. 3), especially at the first harmonic where a clear differentiation can be distinguished between the spin and beat periods [1].

Figure 2: The orbital modulation from AR Scorpii in R-band, showing ellipsoidal variations that may be due to a Rochle lobe filling (or nearly lobe filling) secondary star. Superimposed on this orbital light curve, optical pulsations can be seen at the spin period (117 s) and beat period (118.2 s) with the binary system. Adopted from [1].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-9eac8ab3-0f29-4b24-817d-c23be25413e1.jpg)

Figure 3: Periodograms from radio to optical (HST), all revealing the strong double pulse beat period between the spin and orbital periods of the system. The beat pulse and the spin frequency are most prominent at the first harmonic, close to 17 mHz. The peak on the left represents the beat while the peak on the right represents the spin frequency. Adopted from [1].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-341ce114-4dcd-426f-9c6b-568778d6e5ce.jpg)

High-speed all-Stokes optical polarimetry utilizing the South African Astronomical Observatory (SAAO) HIPPO Polarimeter [4] on the 1.9 m telescope reveals strong linear polarization at levels up to 40 % with significantly lower levels of circular polarization ( &lt; 10%), see Fig. 4, [2, 3]. The nature of the linear polarization profile is consistent with synchrotron emission produced in the magnetic field of a rotating magnetic dipole, strengthening the notion that the white dwarf in AR Sco is a white dwarf pulsar.

Also, the notion that nearly the entire Spectral Energy Distribution (SED) may be dominated by non-thermal emission sparked enormous interest among theorists (e.g. [5, 6]). A favorite approach was to model the emission as being the result of the interaction between a highly magnetized white dwarf pulsar and an M-type secondary dwarf star (see Fig. 5). It is proposed that the WD is a nearly perpendicular rotator where both open field lines sweep past the M-type in each rotation period. The interaction between the particle beam streaming out from the open field line regions and the M-type wind would lead to the formation of a bow shock, where the relativistic electrons can be accelerated by the large magnetic field to produce possible synchrotron radiation in the magnetic field of the magnetosphere. However, it has also been shown that very high energy gamma-ray production through a hadronic channel like π ◦ is also possible [7], as well as inverse Compton scattering.

Recent studies of the pulsed nature of the emission in X-rays using XMM-Newton data [8] reveal pulsed emission consistent with the white dwarf spin period with a pulsed fraction of ∼ 14 %. Here it is suggested that the X-ray emission is produced mainly by a thermal plasma on the secondary star heated by magnetic pumping from the white dwarf's magnetic field lines that sweeps across it during every rotation cycle. The X-ray luminosity of Lx = 4 × 10 30 erg s -1 is

Figure 4: Time series polarimetry (10s bins) that shows the degree of linear polarization, position angle of linear polarization and the degree of circular polarization. Adapted from [2].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-d2263def-adb8-402a-9c79-8449bcd4204e.jpg)

Figure 5: Proposed schematic illustration of the WD/MD binary in AR Sco, where (a) shows the increase of emission from the white dwarf as the electron/positron beam approaches the M-dwarf wind (point A), and (b) shows the episode when the open field lines of the WD interacts with the M-dwarf to form a bow shock, where electrons are accelerated to relativistic energies. Adopted from [5].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-c4c25ea0-0d1e-467d-b573-4dae0bc4f920.jpg)

then produced mainly by the thermalized electrons on the face of the secondary star, with the nonthermal component, peaking in IR/Optical, produced by relativistic electrons accelerated in the magnetosphere of the white dwarf [8]. It is also suggested that the relativistic electrons are trapped within the white dwarf's magnetic field lines due to the magnetic mirror effect. Both the models presented by Geng et. al. [5] and Takata et. al. [8] suggest non-thermal emission inside this region.

For AR Sco the ratio of the X-ray luminosity to spin-down power, i.e. α = Lx / L ˙ ν s wd , ∼ 10 -3 , is similar to another cataclysmic variable AE Aquarii and place both AR Sco and AE Aquarii in the same class as spin-powered pulsars [9]. The high ratio of spin-down power to thermal X-ray emission provides a reservoir for particle acceleration that makes AR Sco a possible candidate for high energy gamma-ray emission like most other spin-powered pulsars [9].

## 2. Using 10 years of Fermi-LAT data to search for γ -ray emission

A Fermi-LAT dataset from the past decade (4Aug2008-18Mar2019) was extracted from the Fermi Science Support Center (FSSC) in the energy range of 100 MeV to 500 GeV. By using the Pass 8 data analysis pipeline, which provides a better determination of the diffuse galactic gammaray emission and a significant improvement in terms of energy resolution from previous Fermi-LAT pipelines, it was possible to do a standard Binned Likelihood Analysis on the Fermi-LAT dataset.

The analysis was performed using the recently released Fermi Tools software, where a number of energy ranges were experimented with to fit and produce likelihood statistics for AR Sco. The event files, i.e. the photon and spacecraft files that were extracted from the FSSC site, were chosen to have a Region of Interest (ROI) of 10 ◦ centred at ARSco (RA:16 21 h m 47 28 , . s Dec: -22 53 10 39 , J2000) in order to extract the high energy photons for analysis. ◦ ′ . ′′ During reduction, only SOURCE class events ( evclass=128 ) were considered, along with the relational filter ( DATA\_QUAL &gt; 0 and LAT\_CONFIG==1 ) to select good time intervals. All the point sources in the fourth Fermi-LAT catalogue (4FGL) located within the ROI were modelled in the spectral fits, including the latest isotropic background and galactic diffuse emission. All the point sources included in the background model file are the same as reported in the 4FGL catalogue.

As AR Sco is not listed in the 4FGL catalogue we have followed the same methodology as ([10]), which reported the detection of a gamma-ray flare from the high redshift blazar DA 193, by manually adding the source to the model centred at its coordinates. The associated parameters of spectral models of all the sources within a 2 ◦ radius were chosen to be free to vary while performing the likelihood fitting. Since AR Sco is not part of the 4FGL catalogue, it was first added by parameterizing a power law spectral shape where the prefactor and photon index of the model were allowed to vary during optimization. A power law model was chosen since it is the easiest to model and is analogous with non-thermal emission. After the power law model was added, both broken power law and log parabola were also used to do a likelihood fit. The final results are summarized in Table 1.

Table 1: Spectral model parameters obtained for AR Sco using Binned Likelihood Analysis. Optimization of AR Sco was done using a power law ( Γ = - (2.7 ± 0.8)), broken power law and log parabola to determine the likelihood of detection of gamma-rays from the ROI centred on AR Sco.

| Spectral model   | IntFlux(photons/cm 2 /s)   |   N pred |    TS | Significance   |
|------------------|----------------------------|----------|-------|----------------|
| Power Law        | 0.565 ± 0.311 x 10 - 8     |  1092.78 |  8.53 | ≤ 2.92 σ       |
| Broken Power Law | 0.486 ± 0.261 x 10 - 8     |  1126.57 | 15.04 | ≤ 3.87 σ       |
| Log Parabola     | 0.251 0.179 x 10 - 8       |   737.39 | 13.15 | 3.62 σ         |

±

≤

Results from a previous Fermi-LAT study, shown in Table 2 ([11]), using the methodology of removing sources with a TS value less than 25 showed a possible gamma-ray emission from the region. However, during the analysis process one of the sources closest to AR Sco,the unidentified source 4FGL J1623.7-2315 nearly 0 6 . ◦ away, was not taken into account during optimization and may have contributed to emission attributed to AR Sco. These results revealed an excess at the √ TS ≤ 4.89 σ level from the position of AR Sco between the energy range of 0.1 - 500GeV. Since

all the sources closest to the ROI centre has to be taken into account, the results obtained from this method (see Table 2) needs to be adjusted by incorporating the associated parameters of spectral models of all the sources within a 2 ◦ radius mentioned earlier.

Table 2: Summary of model parameters from the power law fit for AR Sco reported in [11] excluding the nearby unidentified source 4FGL J1623.7-2315 . The significance levels reported by these authors ([11]) are therefore perhaps too optimistic due to contamination of this near-by ( 0.6 ◦ ) source !!

| Energy Bin     | Spectral Index   |   N pred | IntFlux(photons/cm 2 /s)   |    TS | Significance   |
|----------------|------------------|----------|----------------------------|-------|----------------|
| 100MeV-500 GeV | 2.70 ± 0.82      |  2041.91 | 1.076 ± 0.318 x 10 - 8     | 23.93 | ≤ 4.89 σ       |
| 100MeV-10 GeV  | 2.81 ± 0.22      |  1866.78 | 0.997 ± 0.317 x 10 - 8     | 19.35 | ≤ 4.41 σ       |
| 10GeV-500 GeV  | 6.03 1.04        |     0    | 5.681 1.160 x 10 - 16      |  0    | 0              |

±

±

These new updated results (see Table 1) of our present analysis in the energy range of 0.1 -500 GeV reveals low level emission at the ≤ 3.87 σ significance level from the position of AR Sco. After the likelihood analysis has been performed, SEDs were plotted using e.g. the power law model over 9 energy bins (see Fig. 6). The SED that was obtained using the Fermi build-in function bdlikeSED . The energy flux values can be considered as as 2 σ upper limits based on a limiting TS &lt; 4 value per bin. The distribution of TS values per energy bin (see Figure 7) also illustrates that the most significant emission is concentrated in the energy bins below 3 GeV.

Figure 6: Energy flux SED obtained using bdlikeSED. The flux values are 2 σ upper limits.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-e56197f4-6fc1-46fa-b3fc-db8d0a1989d3.jpg)

Figure 7: TS values vs energy of γ -ray photons.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00283v1-figure-a773ddbd-0f02-4cf4-9623-a5bade0563ab.jpg)

From the distribution of TS vs energy results, see Fig. 7, we noticed a higher TS value from the second energy bin compared to the total TS obtained, see Table 1. Likelihood analysis has been done between the 100MeV-3GeV energy range to investigate this higher TS value. The results within this energy range shows a total TS of 11.46 at a photon flux of 0.611 ± 0.482 x 10 -8 photons cm -2 s -1 , which corresponds to the maximum TS value obtained in Fig. 7. Due to Fermi's low spatial resolution at lower energies, it may be possible for other nearby sources to still contribute to the detected emission.

## 3. Conclusions

A search for possible gamma-ray emission from AR Sco was carried out by analysing the Fermi-LAT Pass 8 dataset using the Fermi Science Tools. The results from a binned-likelihood analysis on a ROI with a radius of 10 ◦ centred on the position of AR Sco, optimized with powerlaw, broken power law and log parabola models, revealed low-level gamma-ray flux between 100 MeV - 500 GeV at the 3.87 σ level by considering the broken power law as the best fit. This is below the 5 σ level heralding a possible detection. Due to the complexity of the surrounding region and the Rho Oph molecular cloud, as well as the poor spatial resolution of the Fermi-LAT telescope at lower energies, it is difficult to pinpoint the exact origin of this low level gamma-ray emission. Unbinned likelihood analysis is underway to try and maximise the probability of detection from the emitting region.

## References

- [1] T. R. Marsh et al. A radio-pulsing white dwarf binary star , Nature 537 (2016) 374
- [2] D.A.H. Buckley, P.J. Meintjes, S.B. Potter, T.R. Marsh &amp; B.T. Gänsicke Polarimetric evidence of a white dwarf pulsar in the binary system AR Sco , Nature Astronomy 1 (2017) 0029
- [3] D.A.H. Buckley, S.B Potter, P.J. Meintjes, T.R. Marsh and b B.T. Bänsicke Polarimetric evidence of the first white dwarf pulsar: The Binary System AR Scorpii , Galaxies 6 (2018) 14
- [4] S.B. Potter et al. Polarized QPOs from the INTEGRAL polar IGRJ14536-5522 (= Swift J 1453.4 - 5524) , MNRAS 402 (2010) 1161
- [5] J.-J. Geng, B. Zhang, Y.-F. Huang A model of white dwarf pulsar AR Sco , Astrophys. J. Lett. 831 (2016) L1
- [6] J.I. Katz AR Sco: Presenting a white dwarf synchronar , Astrophys. J. 835 (2017) 150
- [7] W. Bednarek Hadronic modelfor the non-thermal radiation from the binary system AR Scorpii , MNRAS 476 (2018) L10-L14
- [8] J. Takata, C. -P. Hu, L.C.C. Lin, P.H.T. Tam, P.S. Pal, C.Y. Hui, A.K.H Kong and K.S. Cheng A non-thermal pulsed X-ray emission of AR Scorpii , The Astrophysical Journal 853(2) (2018), 106
- [9] W. Becker &amp; J. Trümper The x-ray luminosity of rotation powered neutron stars , Astronomy &amp;Astrophysics 326 (1997) 682
- [10] V.S. Paliya et al. Detection of a Gamma-Ray Flare from the High-redshift Blazar DA 193 , The Astrophysical Journal 871 (2019) 211
- [11] Q. Kaplan, P.J. Meintjes, K.K. Singh, H.J. van Heerden and F.A. Ramamonjisoa, White Dwarf Pulsars as Possible Gamma-Ray Sources , PoS(HEASA 2018)037 , in Proceedings of High Energy Astrophysics in Southern Africa-HEASA-2018 , 1-3 August 2018, Parys (South Africa)
Q. Kaplan P. J. Meintjes K. K. Singh H. J. van Heerden F. A. Ramamonjisoa I. P. van der Westhuizen
Results from the long-term Fermi-LAT observations of the white dwarf binary pulsar AR Scorpii
The discovery of the white dwarf binary system AR Scorpii (AR Sco) with its
fascinating non-thermal dominated multi-frequency emission has sparked renewed
interest in potential high energy gamma-ray emission from white dwarf pulsars.
The Spectral Energy Distribution (SED) below and above optical shows evidence
of non-thermal synchrotron emission, with pulsed emission in optical and X-ray
bands at the white dwarf spin period ($P_{*} = 117 \, \mbox{s}$) as well as a
beat period ($P_{\rm b} = 118.1 \, \mbox{s}$) with the binary period. From an
energy perspective, the highly magnetic rotating white dwarf can accelerate
particles to TeV energies. In this study, a search for high energy gamma-ray
emission was conducted between 100 MeV - 500 GeV by analysing the newly
available Fermi-LAT Pass 8 data with the new Fermi 1.0.1 Science Tools. Binned
likelihood analysis was done using power law, broken power law and log parabola
models. From the selected Region of Interest (ROI) centred on AR Sco's
position, we calculated a significance of $\sqrt{TS}\leq$ 3.87$\sigma$ for the
integrated gamma-ray activity between 100 MeV - 500 GeV at a photon flux level
of 0.486 $\pm$ 0.261 x 10$^{-8}$ photons cm$^{-2}$ s$^{-1}$ using the broken
power law model. This resulted in a 3$\sigma$ upper-limit detection from the
position of AR Sco. The location of AR Sco inside the Rho Ophiuchi (Rho Oph)
molecular cloud complex combined with the poor spatial resolution of Fermi-LAT,
complicates any positive identification of low-level gamma-ray activity at the
location that coincides with the position of AR Sco.