Draft version March 31, 2021 Preprint typeset using LATEX style emulateapj v. 12/16/11 REVEALING A NEW BLACK WIDOW BINARY 4FGL J0336.0+7502 Kwan-Lok Li1, Y. X. Jane Yap2, Chung Yue Hui3, and Albert K. H. Kong2 Draft version March 31, 2021 ABSTRACT We report on a discovery of a promising candidate as a black widow millisecond pulsar binary, 4FGL J0336.0+7502, which shows many pulsar-like properties in the 4FGL-DR2 catalog. Within the 95% error region of the LAT source, we identified an optical counterpart with a clear periodic- ity at Porb = 3:718178(9) hours using the Bohyunsan 1.8-m Telescope, Lulin One-meter Telescope, Canada{France{Hawaii Telescope, and Gemini-North. At the optical position, an X-ray source was marginally detected in the Swift/XRT archival data, and the detection was confirmed by our Chan- dra/ACIS DDT observation. The spectrum of the X-ray source can be described by a power-law +1:2 −14 −2 −1 model of Γx = 1:6 ± 0:7 and F0:3−7keV = 3:5−1:0 × 10 erg cm s . The X-ray photon index and the low X-ray-to-γ-ray flux ratio (i.e., < 1%) are both consistent with that of many known black widow pulsars. There is also a hint of an X-ray orbital modulation in the Chandra data, although the significance is very low (1:3σ). If the pulsar identity and the X-ray modulation are confirmed, it would be the fifth black widow millisecond pulsar binary that showed an orbitally-modulated emission in X-rays. Subject headings: binaries: close | gamma rays: stars | pulsars: general | X-rays: binaries 1. INTRODUCTION bands (see, e.g., Romani & Sanchez 2016; Yap et al. The Fermi Large Area Telescope (LAT) has been ob- 2019 for the details), although exceptions exist (e.g., serving the MeV/GeV γ-ray sky since June 2008. With 3FGL J0212.1+5320 that does not show any observable the ten years of data taken between 2008 and 2018, 5064 pulsar heating effect; Li et al. 2016). Some recently dis- sources are detected in the Fermi LAT 10-Year Point covered BW/RB candidates include 3FGL J0954.8−3948 Source Catalog (4FGL-DR2; Abdollahi et al. 2020; Bal- (Li et al. 2018), 4FGL J2333.1−5527 (Swihart let et al. 2020). While a major portion of the cataloged et al. 2020), 4FGL J0935.3+0901 (Wang et al. sources are known systems (e.g., active galaxies, pulsars, 2020), 4FGL J0407.7−5702 (Miller et al. 2020), and etc.), about one-fourth of them are unidentified at other 4FGL J0940.3−7610 (Swihart et al. 2021; see the Table wavelengths. Other than active galaxies, many of these 3 of Hui & Li 2019 and the references therein for more unidentified γ-ray sources are believed to be pulsar sys- candidates). tems. In this paper, we present a multi-wavelength study There have been multi-wavelength searching cam- for 4FGL J0336.0+7502, which is a new BW MSP paigns conducted for new candidates of γ-ray pulsars candidate identified by our unidentified Fermi-LAT from the list of the unidentified Fermi-LAT sources (e.g., sources observing campaign. The study includes (i) Hui et al. 2015; Braglia et al. 2020). Machine learning the Fermi-LAT γ-ray properties of 4FGL J0336.0+7502 techniques were also applied on the classification for pul- (§2); (ii) optical photometric observations taken by the sars based on only the γ-ray properties recently (e.g., Saz Bohyunsan 1.8-m Telescope, Lulin One-meter Telescope Parkinson et al. 2016; Luo et al. 2020; Hui et al. 2020). (LOT), Canada{France{Hawaii Telescope (CFHT), and These efforts have led to at least a dozen candidates for Gemini-North (§3); (iii) Swift/XRT and Chandra/ACIS- further radio/γ-ray pulsation searches. Many of these S X-ray observations (§4); and (iv) a discussion for candidates could be associated with two special pulsar 4FGL J0336.0+7502 based on its multi-wavelength prop- classes, black widow (BW) and redback (RB), which are erties (§5). millisecond pulsars in compact binaries (the orbital pe- 2. GAMMA-RAY PROPERTIES arXiv:2103.01557v2 [astro-ph.HE] 30 Mar 2021 riods are often less than a day). Besides the compact 4FGL J0336.0+7502 is a bright γ-ray source (detec- orbits, the two classes are characterised by the very low- tion significance of 31:1σ) located at a Galactic latitude mass companions (i.e., 0.1{0.4 M for RBs and < 0:1M of b = 15:5◦ (Ballet et al. 2020). It was first discovered for BWs; Roberts 2013; Chen et al. 2013) ablated by the in γ-rays by Fermi-LAT as 1FGL J0334.2+7501 in the strong radiations that originate from the primary pul- 1FGL catalog (Abdo et al. 2010), and was subsequently sars. The radiation would also heat up the tidally-locked cataloged in 2FGL, 3FGL, and 4FGL(-DR2) (Nolan et al. companion one-sided, and this so-called pulsar heating 2012; Acero et al. 2015; Abdollahi et al. 2020; Ballet et al. effect can result in orbital modulations in the optical 2020). In 4FGL-DR2 (using data taken from August 1 Department of Physics, National Cheng Kung University, 2008 to August 2018), 4FGL J0336.0+7502 was classified 70101 Tainan, Taiwan; [email protected] (KLL) as a steady source on a yearly time-scale with a variabil- 4 2 Institute of Astronomy, National Tsing Hua University, ity index of 9.8. The average γ-ray flux in 100 MeV{ Hsinchu 30013, Taiwan 3 Department of Astronomy and Space Science, Chungnam 4 A source with a variability index greater than 18.48 has less National University, Daejeon, Republic of Korea than 1% chance to be stable 2 Li et al. TABLE 1 The 4FGL-DR2 LogParabola parameters of 4FGL J0336.0+7502 and ten arbitrarily chosen black widow millisecond pulsars. Name Γγ βγ 4FGL J0336.0+7502 1:78 ± 0:10 0:39 ± 0:07 PSR J0023+0923 1:89 ± 0:11 0:46 ± 0:09 PSR J0610−2100 1:99 ± 0:09 0:51 ± 0:06 PSR J0952−0607 2:02 ± 0:18 0:68 ± 0:11 PSR J1124−3653 1:89 ± 0:06 0:34 ± 0:04 PSR J1301+0833 1:93 ± 0:08 0:41 ± 0:06 PSR J1544+4937 2:07 ± 0:10 0:20 ± 0:05 PSR J1731−1847 2:33 ± 0:24 0:71 ± 0:25 PSR J1805+0615 2:01 ± 0:14 0:50 ± 0:12 PSR J2047+1053 1:70 ± 0:15 0:39 ± 0:08 PSR J2214+3000 1:68 ± 0:04 0:49 ± 0:03 −12 −2 −1 100 GeV is Fγ = (7:4 ± 0:5) × 10 erg cm s that is Within the Fermi-LAT 95% error circle, we spotted a among the top 30% of all 4FGL-DR2 sources. According variable star that was undetected in the first 9 frames, to 4FGL-DR2, the γ-ray spectrum in the LAT energy and then became observable with the magnitude rais- domain is significantly curved with > 8σ, and can be ing from R ≈ 23 mag to ≈ 21 mag in just 40 min- described by a LogParabola model, utes. We stacked the first 9 images, and the vari- able was marginally detected with R ≈ 24 mag (Figure 3). The variable source is also cataloged in the Gaia dN E −(Γγ +βγ log(E=Eb)) / ; DR3 with α(J2000) = 03h36m10s:1811(1), δ(J2000) = dE Eb +75◦0301700: 268(1) (54.0424214(5), 75.0547967(3); Gaia where Γ is the photon index, β defines the degree of Collaboration et al. 2020), at which an X-ray source was γ γ marginally detected in archival Swift/XRT observations curvature, and Eb is a fixed scale parameter. The best-fit (see §4.1). spectral parameters are Γγ = 1:78±0:10 and βγ = 0:39± 0:07. In Table 1, we have listed the Γγ and βγ parameters 3.2. Lulin One-meter Telescope of ten arbitrarily chosen black widow millisecond pulsars in 4FGL-DR2 for comparison. We followed-up the variable source using LOT at Lulin The γ-ray properties of 4FGL J0336.0+7502, includ- Observatory. The observing dates are 2019 October ing the off-plane Galactic latitude, the low long-term 17/18, 2019 November 10, and 2020 January 8/9. As variability, and the significant spectral curvature, are the source is faint for a 1-m class telescope, the data exactly the characteristics that are commonly seen in were taken unfiltered with 600 seconds per frame. 174 γ-ray pulsars (see, e.g., Hui et al. 2015). In fact, frames were taken in total in the 5 days. Standard data 4FGL J0336.0+7502 has been suggested by Saz Parkin- reduction and analysis processes were carried out using son et al. (2016) as a good MSP candidate using machine IRAF as described in the previous section. Flux calibra- learning techniques. This motivated us to investigate tion was done using the Pan-STARRS catalog (PS1 DR2; the system further in multi-wavelength. With the high Flewelling et al. 2020), although the observations are un- quality of the LAT data, the 95% positional accuracy of filtered. 4FGL J0336.0+7502 is around 1:05, which is fine enough The target's brightness was varying in the LOT light to allow a feasible search for the counterpart in optical curve as we have seen in the BOAO data. About half of and X-rays (see Figure 1).