Issue |
A&A
Volume 672, April 2023
|
|
---|---|---|
Article Number | A151 | |
Number of page(s) | 7 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/202245214 | |
Published online | 14 April 2023 |
New upper limits on low-frequency radio emission from isolated neutron stars with LOFAR
1
Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, PO Box 94249 1090 GE Amsterdam, The Netherlands
e-mail: ines.pastormarazuela@uva.nl
2
ASTRON, the Netherlands Institute for Radio Astronomy, PO Box 2 7790 AA Dwingeloo, The Netherlands
3
NYU Abu Dhabi, PO Box 129188 Abu Dhabi, UAE
Received:
14
October
2022
Accepted:
2
January
2023
Neutron stars that show X-ray and γ-ray pulsed emission must generate electron-positron pairs somewhere in the magnetosphere. Pairs like this are also required for radio emission, which poses the question why a number of these sources appear to be radio quiet. We carried out a deep radio search toward four such neutron stars that are isolated X-ray or γ-ray pulsars, but for which no radio pulsations have been detected so far. These sources are 1RXS J141256.0+792204 (Calvera), PSR J1958+2846, PSR J1932+1916, and SGR J1907+0919. A search at lower radio frequencies, where the radio beam is thought to be wider, increases the chances of detecting these sources compared to the earlier higher-frequency searches. We thus carried out a search for periodic and single-pulse radio emission with the LOFAR radio telescope at 150 MHz. We used the known periods and searched a wide range of dispersion measures because the distances are only poorly constrained. We did not detect pulsed emission from any of the four sources. However, we place highly constraining upper limits on the radio flux density at 150 MHz, of ≲1.4 mJy.
Key words: stars: neutron / stars: magnetars
© The Authors 2023
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Through their spin and magnetic field, neutron stars act as powerful cosmic dynamos that can generate a wide variety of electromagnetic emission. There thus exist many subclasses of neutron stars, with different observed behaviors. The evolutionary links between some of the classes are established, while for others, these connections are currently unknown. The largest group in this varied population is formed by the regular rotation-powered radio pulsars. The fast-spinning, high magnetic field population is generally constituted by young pulsars. They show a high spin-down energy-loss rate Ė and a number of energetic phenomena such as radio giant pulse (GP) emission. The most extreme of these fast-spinning and/or high-field sources might also power fast radio bursts (FRBs; e.g., Pastor-Marazuela et al. 2022). On the long-period outskirts of the P–Ṗ diagram, slowly rotating pulsars (e.g., Young et al. 1999; Tan et al. 2018) and magnetars (e.g., Caleb et al. 2022; Hurley-Walker et al. 2022) sometimes continue to radiate.
Some neutron stars, however, only emit intermittently at radio frequencies. Rotating radio transients (RRATs) burst very irregularly, and most are found near the death line in the P–Ṗ diagram (Keane et al. 2011), between the canonical radio pulsars and magnetars. The exact evolutionary connection between RRATs and the steadily radiating normal pulsars is unclear, but studies suggest the presence of an evolutionary link between these different classes (e.g., Burke-Spolaor 2012).
Finally, populations of neutron stars exist that do not appear to emit in radio at all: radio-quiet magnetars such as most anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs), X-ray dim isolated neutron stars (XDINSs; Haberl 2007), and γ-ray pulsars (e.g., Abdo et al. 2013). These are able to produce high-energy emission, but are often radio quiet. Gençali & Ertan (2018) proposed that RRATs can evolve into XDINSs through a fallback accretion disk, thus becoming radio quiet. However, the magnetar SGR 1935+2154 was recently seen to emit a bright radio burst that bridged the gap in radio luminosities between regular pulsars and FRBs (CHIME/FRB Collaboration 2020; Bochenek et al. 2020; Maan et al. 2022b). This suggests that magnetars might explain the origin of some, if not all, extragalactic FRBs.
Some of these magnetars might produce radio emission that is only visible at low radio frequencies. Detections of radio pulsations of the γ- and X-ray pulsar Geminga, PSR J0633+1746, have been claimed at and below the 100 MHz observing frequency range (Malofeev & Malov 1997; Malov et al. 2015; Maan 2015), although a very deep search using the low-frequency array (LOFAR; van Haarlem et al. 2013) failed to detect any such pulsations (Coenen 2013). These low-frequency detections offer an intriguing possibility to better understand the radio emission mechanism of these enigmatic objects. Radio detections of a magnetar with LOFAR, complementary to higher-frequency studies such as Camilo et al. (2006) and Maan et al. (2022a) for XTE J1810−197, might offer insight into emission mechanisms and propagation in ultra-strong magnetic fields.
The XDINSs feature periods that are as long as those in magnetars, but their magnetic field strength is less extreme. The XDINSs form a small group of seven isolated neutron stars that show thermal emission in the soft X-ray band. Since their discovery with ROSAT in the 1990s, several attempts failed to detect these sources at radio frequencies (e.g., Kondratiev et al. 2009). As these campaigns operated above 800 MHz, a sensitive lower-frequency search may be opportune. It has been proposed (e.g., Komesaroff 1970; Cordes 1978) and observed (e.g., Chen & Wang 2014) that pulsar profiles are usually narrower at higher frequencies and become broader at lower radio frequencies. This suggests that the radio emission cone is broader at low frequencies, and sweeps across a larger fraction of the sky as seen from the pulsar. Additionally, radio pulsars often present negative spectral indices, and are thus brighter at lower frequencies (Bilous et al. 2016). If all neutron star radio beams are broader and brighter at lower frequencies, the probability of detecting radio emission from γ- and X-ray isolated neutron stars (INSs) increases at the lower radio frequencies offered through LOFAR. The earlier observations that resulted in nondetections might then just have missed the narrower high-frequency beam, where the wider lower-frequency beam may, in contrast, actually enclose Earth. In this situation, LOFAR might detect the source.
Recently, a number of radio pulsars were discovered that shared properties with XDINSs and RRATs, such as soft X-ray thermal emission, a similar position in the P–Ṗ diagram, and a short distance to the Solar System. These sources, PSR J0726−2612 (Rigoselli et al. 2019) and PSR J2251−3711 (Morello et al. 2020), support the hypothesis that XDINS are indeed not intrinsically radio quiet, but that their radio beam points away from us. These shared properties might reflect a potential link between the radio and X-ray emitting pulsars with XDINSs and RRATs. A firm low-frequency radio detection of INSs would thus tie these observationally distinct populations of neutron stars together.
In this work we present LOFAR observations of four INSs that brightly pulsate at X-ray or γ-ray energies, but have not been detected in radio. These sources are listed in Sect. 2, and their parameters are presented in Table 1.
Parameters of the observed pulsars and observational setup of the observations in the LC3_036 proposal.
2. Targeted sources
2.1. J1412+7922
The INS 1RXS J141256.0+792204, called “Calvera” and hereafter J1412+7922, was first detected with ROSAT (Voges et al. 1999) as an X-ray point source. It was subsequently detected with Swift and Chandra (Rutledge et al. 2008; Shevchuk et al. 2009). X-ray observations confirmed its neutron star nature through the detection of P ≃ 59 ms pulsations by Zane et al. (2011), and allowed for the determination of its spin-down luminosity Ė ∼ 6 × 1035 erg s−1, characteristic age τc ≡ P/2Ṗ ∼ 3 × 105 years, and surface dipole magnetic field strength Bs = 4.4 × 1011 G by Halpern et al. (2013). Although these values are not unusual for a rotation-powered pulsar, the source is not detected in radio (Hessels et al. 2007; Zane et al. 2011) or γ-rays (Mereghetti et al. 2021). The X-ray emission can be modeled with a two-temperature blackbody spectrum (Zane et al. 2011), similar to other XDINS (Pires et al. 2014). However, the spin period of J1412+7922 is much faster than typically observed in XDINS. Since the source is located at high Galactic latitudes and its inferred distance is relatively small (∼3.3 kpc; Mereghetti et al. 2021), the path through the interstellar medium is not long enough to explain the radio nondetections by high dispersion measure (DM) or scattering values.
2.2. J1958+2846
Discovered by Abdo et al. (2009) through a blind frequency search of Fermi-LAT γ-ray data, INS PSR J1958+2846, hereafter J1958+2846, has shown no X-ray or radio continuum emission counterpart so far (Ray et al. 2011; Frail et al. 2016). Arecibo observations have placed very constraining upper limits of 0.005 mJy at 1510 MHz (Ray et al. 2011). Searches for pulsations from the source using the single international LOFAR station FR606 by Grießmeier et al. (2021) also found no periodic signal.
The double-peaked pulse profile of J1958+2846 can be interpreted as a broad γ-ray beam. The earlier higher-frequency radio nondetections might be due to a narrower radio beam and to an unfavorable rotation geometry with respect to the line of sight. If the radio beam is indeed wider at lower frequencies, LOFAR would have higher chances of detecting it. In this case, a setup more sensitive than the Grießmeier et al. (2021) single-station search is required.
Modeling by Pierbattista et al. (2015) indicated that the γ-ray pulse profile of J1958+2846 can be well fit by one-pole caustic emission (OPC; Romani & Watters 2010; Watters et al. 2009) or an outer-gap model (OG; Cheng et al. 2000). In both cases, the γ-rays are generated at high altitudes above the NS surface. Each model constrains the geometry of the pulsar. For the OPC model, the angle between the rotation and magnetic axes α = 49°, while the angle between the observer line of sight and the rotational axis ζ = 85°. The OG model reports similarly large angles, with the NS equator rotating in the plane that also contains Earth, and an oblique dipole: α = 64° , ζ = 90°. If this model is correct, the low-frequency radio beam would thus need to be wider than ∼30° to encompass the telescope. This is uncommonly wide; only 8 out of the 600 pulsars in the ATNF catalog that are not recycled and have a published 400 MHz flux have a duty cycle suggestive of a beam wider than 30% (Manchester et al. 2005). Because a width like this is unlikely, a total-intensity detection would thus to first order suggest a geometry where α and ζ are closer than follows from Pierbattista et al. (2015), even if this suggestion would only be qualitative. Subsequent follow-up measurements of polarization properties throughout the pulse, and a fit of these properties to the rotating vector model (RVM; Radhakrishnan & Cooke 1969), can quantify the allowed geometries to within a relatively precise combinations of α and ζ. In a similar study on radio-loud γ-ray pulsars, Rookyard et al. (2015) already reported that RVM fits suggest that the magnetic inclination angles α are much smaller than predicted by the γ-ray light curve models. This in turn confirms that deep radio searches can lead to detections even when the γ-ray light curves suggest that the geometry is unfavorable.
2.3. J1932+1916
The INS PSR J1932+1916, hereafter J1932+1916, was discovered in Fermi-LAT data through blind searches with the Einstein at Home volunteer computing system (Clark et al. 2017). J1932+1916 is the youngest and γ-ray brightest of the four γ-ray pulsars presented from that effort in Pletsch et al. (2013). The period is 0.21 s, and the characteristic age is 35 kyr. Frail et al. (2016) reported no continuum 150 MHz source at this position with GMRT at a flux density upper limit of 27 mJy beam−1, with 1σ errors. If the flux density they found at the position of the pulsar is in fact the pulsed emission from J1932+1916, then a LOFAR periodicity search as described here should detect the source at an S/N of 15 if the duty cycle is 10%. Karpova et al. (2017) reported a potential pulsar wind nebula (PWN) association from Swift and Suzaku observations. However, no X-ray periodicity searches have been carried out before.
2.4. J1907+0919
The soft gamma repeater J1907+0919, also known as SGR 1900+14, was detected through its bursting nature by Mazets et al. (1979). Later outbursts were detected in 1992 (Kouveliotou et al. 1993), 1998 (Hurley et al. 1999), and 2006 (Mereghetti et al. 2006). The August 1998 outburst allowed the detection of an X-ray period of ∼5.16 s, and thus confirmed the nature of the source as a magnetar (Hurley et al. 1999; Kouveliotou et al. 1999). Frail et al. (1999) detected a transient radio counterpart that appeared simultaneously to the 1998 outburst, and they identified the radio source as a synchrotron-emitting nebula. Shitov et al. (2000) claimed to have found radio pulsations at 111 MHz from four to nine months after the 1998 burst, but the number of trials involved in the search, the small bandwidth of the system, and the low S/N of the presented plots lead us to conclude that the confidence level for these detections is low. No other periodic emission has been found at higher radio frequencies (Lorimer & Xilouris 2000; Fox et al. 2001; Lazarus et al. 2012).
This paper is organised as follows: in Sect. 3 we explain how we used LOFAR (van Haarlem et al. 2013) to observe the sources mentioned above. In Sect. 4 we detail the data reduction procedure, including the periodicity and the single-pulse searches that we carried out. In Sect. 5 we present our results, including the upper limit that we set on the pulsed emission. In Sect. 6 we discuss the consequences of these nondetections for the radio-quiet pulsar population, and in Sect. 7 we conclude.
3. Observations
We observed the four sources with the largest possible set of high band antennas (HBAs) that LOFAR can coherently beam form. Each observation thus added 22 HBA core stations, covering 78.125 MHz bandwidth in the 110 MHz–190 MHz frequency range (centered on 148.92 MHz), with 400 channels that were 195 kHz wide. The LOFAR beam-forming abilities allow us to simultaneously observe different regions of the sky (van Leeuwen & Stappers 2010; Stappers et al. 2011; Coenen et al. 2014). For our point-source searches of INSs, we used three beams per observation; one beam pointed to the source of interest, the second beam on a nearby known pulsar, and the third beam as a calibrator blank-sky beam to cross-check potential candidates as possibly arising from radio frequency interference (RFI). We carried out observations between 16 January 2015 and 15 February 2015 under project ID LC3_0361. We integrated for 3 h on each of our sources. The data were taken in Stokes I mode. Since the periods of the γ-ray pulsars are known, the time resolution of each observation was chosen such as to provide good coverage of the pulse period at a sampling time between 0.16−1.3 ms. The observation setup is detailed in Table 1.
4. Data reduction
The data were preprocessed by the LOFAR pulsar pipeline after each observation (Alexov et al. 2010; Stappers et al. 2011) and stored on the LOFAR Long Term Archive2 in PSRFITS format (Hotan et al. 2004). The 1.5 TB of data was then transferred to one of the nodes of the Apertif real-time FRB search cluster ARTS (van Leeuwen 2014; van Leeuwen et al. 2022).
We performed a periodicity search as well as a single-pulse search using PRESTO3 (Ransom 2001). The data were cleaned of RFI using first rfifind, and then impulsive and periodic signals at DM = 0 pc cm−3 were removed. Next we searched the clean data for periodic signals and single pulses. We searched for counterparts around the known P and Ṗ of each pulsar. Additionally, we performed a full blind search in order to search for potential pulsars in the same field of view because many new pulsars are found at low frequencies (Sanidas et al. 2019) and serendipitous discoveries happen regularly (e.g., Oostrum et al. 2020). Because the DM of our sources is unknown, we searched over a range of DMs going from 4 pc cm−3 to 400 pc cm−3. The DM-distance relation is not precise enough to warrant a much smaller DM range, even for sources for which a distance estimate exists; and a wider DM range allows for discovery of other pulsars contained in our field of view. The highest DM pulsar detected with LOFAR has a DM = 217 pc cm−3 (Sanidas et al. 2019). We thus searched up to roughly twice this value to ensure that any detectable sources were covered. We determined the optimal dedispersion parameters with DDplan from PRESTO. The sampling time variation between some of the four observations had a slight impact on the exact transitions of the step size, but the data were generally dedispersed in steps of 0.01 pc cm−3 up to DM = 100 pc cm−3, then by 0.03 pc cm−3 steps up to 300 pc cm−3, and finally, using 0.05 pc cm−3 steps.
We manually inspected all candidates down to σ = 4, resulting in ∼1400 candidates per beam. To verify our observational setup, we used the same blind-search technique to our test pulsars B1322+83 and B1933+16, which we detected. The test pulsar B1953+29 was not detected because the sampling time of the observation of J1958+2846 was not adapted to its ∼6 ms period. However, we were able to detect B1952+29 (Hewish et al. 1968) in this same pointing. Even though it is located at > 1° from the targeted coordinates, it is bright enough to be visible as a side-lobe detection.
The candidates from the PRESTO single-pulse search were further classified using the deep-learning classification algorithm developed by Connor & van Leeuwen (2018), which has been verified and was successful in the Apertif surveys (e.g., Connor et al. 2020; Pastor-Marazuela et al. 2021). This reduced the number of candidates significantly by sifting out the remaining RFI. The remaining candidates were inspected visually.
5. Results
In our targeted observations, we were unable to detect any plausible astronomical radio pulsations or single pulses. We determined new 150 MHz flux upper limits by computing the sensitivity limits of our observations. To establish these sensitivity limits, we applied the radiometer equation adapted to pulsars, as detailed below. We determined the telescope parameters that are input to this equation by following the procedure4 described in Kondratiev et al. (2016) and Mikhailov & van Leeuwen (2016). This approach takes the system temperature (including the sky temperature), the projection effects governing the effective area of the fixed tiles, and the amount of time and bandwidth removed due to RFI into account to produce the overall observation system-equivalent flux density (SEFD).
For the sensitivity limit on the periodic emission, we used the following equation (see, e.g., Dewey et al. 1985):
where β ≲ 1 is a digitization factor, Tsys (K) is the system temperature, G (K Jy−1) is the telescope gain, Δν (Hz) is the observing bandwidth, and tobs (s) is the observation time. P (s) represents the spin period, and W (s) gives the pulsed width assuming a pulsar duty cycle of 10%. To facilitate direct comparison of the periodic emission limits to values reported in Ray et al. (2011) and Grießmeier et al. (2021), for example, we used a minimum signal-to-noise ration S/Nmin = 5. A more conservative option, given the high number of candidates per beam, would arguably be to use a limit of S/N = 8. We did, however, review all candidates with S/N > 4 by eye and the reported sensitivity limits can be scaled to a different S/N value.
The sensitivity limit on the single-pulse emission, Slim,sp, is computed as follows:
where all variables are the same as in Eq. (2). We searched for single pulses down to a signal-to-noise ratio S/Nmin = 7.
We report these periodic and single-pulse sensitivity limits, computed at the coordinates of the central beam of each observation, in Table 1. Even though all observations are equally long, the estimated Slim,p values are different. This is mostly due to the strong dependence of the LOFAR effective area, and hence the sensitivity, on the elevation.
In Fig. 1 we compare our upper limits to those established in previous searches, mostly using the same techniques. Our upper limit on the flux of J1907+0919 is ∼50× deeper than the claimed 1998–1999 detections at the same 3m wavelength with BSA (Shitov et al. 2000). Other searches were generally undertaken at higher frequencies (Hessels et al. 2007; Zane et al. 2011; Ray et al. 2011; Pletsch et al. 2013; Grießmeier et al. 2021). When we assume that the radio spectra of these four pulsars are described by a single power-law Sν ∝ να with a spectral index of α = −1.4 (Bates et al. 2013; Bilous et al. 2016), the upper limits we present here for J1412+7922 and J1932+1916 are the most stringent for any search so far. The upper limits on J1958+2846 (Arecibo; Ray et al. 2011) and J1907+0919 (GBT; Lazarus et al. 2012) are more sensitive than ours by a factor of 2–3. However, pulsars present a broad range of spectral indices. When we take the mean ±2σ measured by Jankowski et al. (2018), spectral indices can vary from −2.7 to −0.5. The flux upper limits we measure would be the deepest assuming a −2.7 spectral index, but the shallowest at −0.5.
Fig. 1. Flux density upper limits of this work at 150 MHz (filled symbols) with S/N = 5 for comparison to earlier searches of the same sources (empty symbols). Solid lines going through our upper limit estimates with spectral index α = −1.4 are overlaid to show the scaling of our sensitivity limits. Our limits are plotted slightly offset from the 150 MHz observing frequency (dashed line) for better visibility. The faded green marker for SGR J1907+0919 represents the detection claimed by Shitov et al. (2000). |
6. Discussion
6.1. Comparison to previous limits
For J1958+2846 and J1932+1916, we can make a straightforward relative comparison between the results we presented here and the existing limit at 150 MHz from the single-station LOFAR campaign by Grießmeier et al. (2021). Our 22 core stations are each one-quarter of the area of the FR606 station and are coherently combined, leading to a difference of a factor in area A for the radiometer equation and Slim. The integration time t of 3 h is shorter than the FR606 total of 8.3 h (J1958+2846) and 4.1 h (J1932+1916), leading to a factor in the radiometer equation. Other factors such as the sky background and the influence of zenith angle on the sensitivity are probably mostly the same for both campaigns. Our Slim is thus times deeper than the Grießmeier et al. (2021) upper limit for J1958+2846, and 4.7 times for J1932+1916. These factors agree well with the actual limits listed in Table 1.
Bilous et al. (2016) measured the mean flux density Smean of 158 pulsars detected with LOFAR, where . Compared to these LOFAR detections, our upper limit on J1412+7922 is deeper than all 158 sources (100%), J1958+2846 is deeper than 156 sources (99%), J1932+1916 is deeper than 144 sources (93%), and J1907+0919 is deeper than 109 sources (69%). The flux upper limits we have set on each of the sources in our sample are some of the deepest compared to other LOFAR radio pulsar detections. Longer observing times are thus unlikely to result in a detection or improve our flux upper limits. Additional follow-up would only be constraining with more sensitive radio telescopes.
6.2. Emission angles and intensity
Different pulsar emission mechanism models exist that predict that radio and γ-ray emission is formed simultaneously in the pulsar magnetosphere. The emission sites are not necessarily co-located, however. The periodic radio emission is generally thought to be formed just above the polar cap. The high-energy polar cap (PC) model next assumes that the γ-ray emission is also produced near the surface of the NS and near the magnetic polar caps. In the models of the outer magnetosphere emission, such as the OG or the OPC models, on the other hand, the γ-ray emission is produced high up in the magnetosphere of the NS, within the extent of the light cylinder.
Of the sources in our sample, specific high-energy geometry models have only been proposed for J1958+2846 (Pierbattista et al. 2015). A detection could have confirmed one of these (Sect. 2.2). Conclusions can be drawn from the nondetections for our sample in general, however. The two general high-energy model classes mentioned above predict different testable beam widths. Our radio nondetections, when attributed to radio beams that are not wide enough to encompass Earth, favor models of the outer magnetosphere (see, e.g., Romani & Watters 2010). The reason is that in the OG and OPC models, the γ-ray beam (which is detected for our sources) is much broader than the radio beam. The radio beam, being much narrower, is unlikely to cut through our line of sight. This model class is thus more applicable than a class in which the radio and high-energy beam are of similar angular size, such as the PC model (or, to a lesser extent, the slot-gap model; Muslimov & Harding 2003; Pierbattista et al. 2015). In this case, detections in both radio and high energy would more often be expected. Our results thus favor OG and OPC models over PC models for high-energy emission.
While it is instructive to discuss the coverage of the radio pulsar beam in binary terms – it either hits or misses Earth – this visibility is not very unambiguous in practice. The beam edge is not sharp. In a beam-mapping experiment enabled by the geometric precession in PSR J1906+0745 (van Leeuwen et al. 2015), the flux at the edge of the beam is over 100× dimmer than the peak, but it is still present and detectable (Desvignes et al. 2019). Deeper searches thus continue to have value, even if nondetections at the same frequency already exist. The detection of PSR J1732−3131 only at 327 and potentially even 34 MHz (Maan & Aswathappa 2014) shows that emission beam widening (or, possibly equivalently, a steep spectral index) at low frequencies is a real effect for γ-ray pulsars as well.
6.3. Emission mechanism and evolution
Most models explain the radio quietness of an NS through a chance beam misalignment, as above. It might also be a more intrinsic property, however. In at least two regions in the P–Ṗ diagram, radio emission may be increasingly hard to generate.
The first parameter space of interest is for sources close to the radio death line (Chen & Ruderman 1993). XDINSs are preferably found there, which suggests that these sources are approaching a state in their evolution in which radio emission ceases generally. In normal pulsars, the death line represents the transition into a state in which electron–positron pair formation over the polar cap ceases completely. When the pulsar rotates too slowly to generate a large enough potential drop over the polar cap, as required for this formation, the radio emission stops (Ruderman & Sutherland 1975). The high-energy emission also requires pair formation, but this could occur farther out. We note that polar-cap pair formation can continue at longer periods if the NS surface magnetic field is not a pure dipole. With this decreased curvature radius, the NS may continue to radiate. Evidence for these higher-order fields is present in a number of pulsars, for instance, PSR J0815+0939 (Szary & van Leeuwen 2017) and PSR B1839−04 (Szary et al. 2020). This would also influence the interpretation of any polarization information, as the RVM generally assumes a dipole field.
None of the sources in our sample are close to this death line (see Fig. 2), but SGR J1907+0919 is beyond a different purported boundary: the photon splitting line (Baring & Harding 2001). In pulsars in this second parameter space of interest, where magnetic fields are stronger than the quantum critical field, 4.4 × 1013 G (Fig. 2), pair formation cannot compete with magnetic photon splitting. These high-field sources might then be radio quiet, but X-ray or γ-ray bright. We mark the critical field line for a dipole in Fig. 2, but we note, as Baring & Harding (2001) did, that higher multipoles and general relativistic effects can subtly change the quiescence limit on a per-source basis. However, based on its spindown dipole magnetic field strength of 7 × 1014 G, our nondetection of SGR J1907+0919 supports the existence of this limit.
Fig. 2. P − Ṗ diagram showing the location of the sources presented in this work. All pulsars from the ATNF Pulsar Catalogue (Manchester et al. 2005) are shown as gray dots. Different pulsar classifications are encircled by different symbols. The sources discussed in this work are shown as black stars, from left to right: J1412+7922, J1932+1916, J1932+1916, and J1907+0919. The orange shaded region is delimited by the death line, while the green shaded region is delimited by the photon splitting line. Plot generated with psrqpy (Pitkin 2018). |
6.4. Propagation effects
While the emission-beam widening and the negative spectral index provide potential advantages in a search for pulsars at low frequencies, some propagation effects such as dispersion and scattering intensify there. This hinders the detection of certain sources. The largest pulsar DM detected with LOFAR is 217 pc cm−3, while many Galactic pulsars are known to have DM > 1000 pc cm−3. Although the sources studied in this work do not have radio detections and thus lack a known DM, we can estimate this DM if a hydrogen column density NH were measured from soft X-ray detections. He et al. (2013) reported the correlation between NH and DM as follows: .
While J1958+2846 and J1932+1916 have only been detected in γ-rays, J1412+7922 and J1907+0919 have soft X-ray detections for which NH has been measured. For J1907+0919, Kouveliotou et al. (1999) measured a high NH value of 3.4 − 5.5 × 1022 cm−2. The correlation suggests a DM of 1100−1800 pc cm−3. At a DM this large, the detection limit of LOFAR is severely impacted. Because J1907+0919 is a very slow rotator, the intrachannel dispersion delay still only is about 10% of the period, which means that peridiocity searches might still detect it in principle. The flux density per bin is much decreased when the pulse is smeared out over 100 s of time bins, however.
In contrast, Shevchuk et al. (2009) reported a measured NH = 3.1 ± 0.9 × 1020 cm−2 for J1412+7922. We thus estimate its DM to be in the range 5–15 pc cm−3. This low DM would easily have been detected with LOFAR.
7. Conclusion
We have conducted deep LOFAR searches of periodic and single-pulse radio emission from four isolated neutron stars. Although we validated the observational setup with the detection of the test pulsars, we did not detect any of the four targeted pulsars. This can be explained with an intrinsic radio-quietness of these sources, as was proposed previously. It might also be caused by a chance misalignment between the radio beam and the line of sight.
With the new upper limits, we can rule out the hypothesis that INSs were not previously detected at radio frequencies around 1 GHz because their spectrum were steeper than that of regular radio pulsars. Since radio emission from magnetars has been detected after high-energy outbursts (e.g., Maan et al. 2022b), additional radio observations of J1907+0919 if the source reactivates might be successful at detecting single-pulse or periodic emission in the future.
After we completed the current manuscript as Pastor-Marazuela (2022), Arias et al. (2022) posted a preprint presenting partly the same data.
Acknowledgments
This research was supported by the Netherlands Research School for Astronomy (‘NOVA5-NW3-10.3.5.14’), the European Research Council under the European Union’s Seventh Framework Programme (FP/2007–2013)/ERC Grant Agreement No. 617199 (‘ALERT’), and by Vici research programme ‘ARGO’ with project number 639.043.815, financed by the Dutch Research Council (NWO). We further acknowledge funding from National Aeronautics and Space Administration (NASA) grant number NNX17AL74G issued through the NNH16ZDA001N Astrophysics Data Analysis Program (ADAP) to SMS. This paper is based (in part) on data obtained with the International LOFAR Telescope (ILT) under project code LC3_036 (PI: van Leeuwen). LOFAR (van Haarlem et al. 2013) is the low frequency array designed and constructed by ASTRON. It has observing, data processing, and data storage facilities in several countries, that are owned by various parties (each with their own funding sources), and that are collectively operated by the ILT foundation under a joint scientific policy. The ILT resources have benefitted from the following recent major funding sources: CNRS-INSU, Observatoire de Paris and Université d’Orléans, France; BMBF, MIWF-NRW, MPG, Germany; Science Foundation Ireland (SFI), Department of Business, Enterprise and Innovation (DBEI), Ireland; NWO, The Netherlands; The Science and Technology Facilities Council, UK; Ministry of Science and Higher Education, Poland.
References
- Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, Science, 325, 840 [Google Scholar]
- Abdo, A. A., Ajello, M., Allafort, A., et al. 2013, ApJS, 208, 17 [Google Scholar]
- Alexov, A., Hessels, J., Mol, J. D., Stappers, B., & van Leeuwen, J. 2010, Astronomical Data Analysis Software and Systems XIX, 434 [Google Scholar]
- Arias, M., Botteon, A., Bassa, C. G., et al. 2022, A&A, 667, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Baring, M. G., & Harding, A. K. 2001, ApJ, 547, 929 [NASA ADS] [CrossRef] [Google Scholar]
- Bates, S. D., Lorimer, D. R., & Verbiest, J. P. W. 2013, MNRAS, 431, 1352 [Google Scholar]
- Bilous, A., Kondratiev, V., Kramer, M., et al. 2016, A&A, 591, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bochenek, C. D., Ravi, V., Belov, K. V., et al. 2020, Nature, 587, 59 [NASA ADS] [CrossRef] [Google Scholar]
- Bogdanov, S., Ho, W. C. G., Enoto, T., et al. 2019, ApJ, 877, 69 [Google Scholar]
- Burke-Spolaor, S. 2012, Proc. Int. Astron. Union, 8, 95 [NASA ADS] [CrossRef] [Google Scholar]
- Caleb, M., Heywood, I., Rajwade, K., et al. 2022, Nat. Astron., 6, 828 [NASA ADS] [CrossRef] [Google Scholar]
- Camilo, F., Ransom, S. M., Halpern, J. P., et al. 2006, Nature, 442, 892 [NASA ADS] [CrossRef] [Google Scholar]
- Chen, J. L., & Wang, H. G. 2014, ApJS, 215, 11 [NASA ADS] [CrossRef] [Google Scholar]
- Chen, K., & Ruderman, M. 1993, ApJ, 402, 264 [NASA ADS] [CrossRef] [Google Scholar]
- Cheng, K. S., Ruderman, M., & Zhang, L. 2000, ApJ, 537, 964 [NASA ADS] [CrossRef] [Google Scholar]
- CHIME/FRB Collaboration (Andersen, B. C., et al.) 2020, Nature, 587, 54 [Google Scholar]
- Clark, C. J., Wu, J., Pletsch, H. J., et al. 2017, ApJ, 834, 106 [Google Scholar]
- Coenen, T. 2013, PhD Thesis, University of Amsterdam, The Netherlands [Google Scholar]
- Coenen, T., van Leeuwen, J., Hessels, J. W. T., et al. 2014, A&A, 570, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Connor, L., & van Leeuwen, J. 2018, ApJ, 156, 256 [CrossRef] [Google Scholar]
- Connor, L., van Leeuwen, J., Oostrum, L. C., et al. 2020, MNRAS, 499, 4716 [CrossRef] [Google Scholar]
- Cordes, J. M. 1978, ApJ, 222, 1006 [NASA ADS] [CrossRef] [Google Scholar]
- Desvignes, G., Kramer, M., Lee, K., et al. 2019, Science, 365, 1013 [Google Scholar]
- Dewey, R. J., Taylor, J. H., Weisberg, J. M., & Stokes, G. H. 1985, ApJ, 294, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Fox, D. W., Kaplan, D. L., Kulkarni, S. R., & Frail, D. A. 2001, ApJ, submitted [arXiv:astro-ph/0107520] [Google Scholar]
- Frail, D. A., Kulkarni, S. R., & Bloom, J. S. 1999, Nature, 398, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Frail, D. A., Jagannathan, P., Mooley, K. P., & Intema, H. T. 2016, ApJ, 829, 119 [NASA ADS] [CrossRef] [Google Scholar]
- Gençali, A. A., & Ertan, Ü. 2018, MNRAS, 481, 244 [CrossRef] [Google Scholar]
- Grießmeier, J. M., Smith, D. A., Theureau, G., et al. 2021, A&A, 654, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Haberl, F. 2007, Ap&SS, 308, 181 [NASA ADS] [CrossRef] [Google Scholar]
- Halpern, J. P., Bogdanov, S., & Gotthelf, E. V. 2013, ApJ, 778, 120 [NASA ADS] [CrossRef] [Google Scholar]
- He, C., Ng, C.-Y., & Kaspi, V. M. 2013, ApJ, 768, 64 [NASA ADS] [CrossRef] [Google Scholar]
- Hessels, J. W. T., Stappers, B. W., Rutledge, R. E., Fox, D. B., & Shevchuk, A. H. 2007, A&A, 476, 331 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hewish, A., Bell, S. J., Pilkington, J. D. H., Scott, P. F., & Collins, R. A. 1968, Nature, 217, 709 [NASA ADS] [CrossRef] [Google Scholar]
- Hotan, A. W., van Straten, W., & Manchester, R. N. 2004, PASA, 21, 302 [Google Scholar]
- Hurley, K., Li, P., Kouveliotou, C., et al. 1999, ApJ, 510, 111 [Google Scholar]
- Hurley-Walker, N., Zhang, X., Bahramian, A., et al. 2022, Nature, 601, 526 [NASA ADS] [CrossRef] [Google Scholar]
- Jankowski, F., van Straten, W., Keane, E. F., et al. 2018, MNRAS, 473, 4436 [Google Scholar]
- Karpova, A., Shternin, P., Zyuzin, D., Danilenko, A., & Shibanov, Y. 2017, MNRAS, 466, 1757 [CrossRef] [Google Scholar]
- Keane, E. F., Kramer, M., Lyne, A. G., Stappers, B. W., & McLaughlin, M. A. 2011, MNRAS, 415, 3065 [Google Scholar]
- Komesaroff, M. M. 1970, Nature, 225, 612 [NASA ADS] [CrossRef] [Google Scholar]
- Kondratiev, V. I., McLaughlin, M. A., Lorimer, D. R., et al. 2009, ApJ, 702, 692 [NASA ADS] [CrossRef] [Google Scholar]
- Kondratiev, V. I., Verbiest, J. P. W., Hessels, J. W. T., et al. 2016, A&A, 585, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kouveliotou, C., Fishman, G. J., Meegan, C. A., et al. 1993, Nature, 362, 728 [NASA ADS] [CrossRef] [Google Scholar]
- Kouveliotou, C., Strohmayer, T., Hurley, K., et al. 1999, ApJ, 510, 115 [Google Scholar]
- Lazarus, P., Kaspi, V. M., Champion, D. J., Hessels, J. W. T., & Dib, R. 2012, ApJ, 744, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Lorimer, D. R., & Xilouris, K. M. 2000, ApJ, 545, 385 [NASA ADS] [CrossRef] [Google Scholar]
- Maan, Y. 2015, ApJ, 815, 126 [Google Scholar]
- Maan, Y., & Aswathappa, H. A. 2014, MNRAS, 445, 3221 [Google Scholar]
- Maan, Y., Surnis, M. P., Chandra Joshi, B., & Bagchi, M. 2022a, ApJ, 931, 67 [NASA ADS] [CrossRef] [Google Scholar]
- Maan, Y., van Leeuwen, J., Straal, S., & Pastor-Marazuela, I. 2022b, ATel, 15697, 1 [NASA ADS] [Google Scholar]
- Malofeev, V. M., & Malov, O. I. 1997, Nature, 389, 697 [NASA ADS] [CrossRef] [Google Scholar]
- Malov, O. I., Malofeev, V. M., Teplykh, D. A., & Logvinenko, S. V. 2015, Astron. Rep., 59, 183 [NASA ADS] [CrossRef] [Google Scholar]
- Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993 [Google Scholar]
- Mazets, E. P., Golenetskii, S. V., & Guryan, Y. A. 1979, Sov. Astron. Lett., 5, 343 [Google Scholar]
- Mereghetti, S., Esposito, P., Tiengo, A., et al. 2006, ApJ, 653, 1423 [NASA ADS] [CrossRef] [Google Scholar]
- Mereghetti, S., Rigoselli, M., Taverna, R., et al. 2021, ApJ, 922, 253 [NASA ADS] [CrossRef] [Google Scholar]
- Mikhailov, K., & van Leeuwen, J. 2016, A&A, 593, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Morello, V., Keane, E. F., Enoto, T., et al. 2020, MNRAS, 493, 1165 [NASA ADS] [CrossRef] [Google Scholar]
- Muslimov, A. G., & Harding, A. K. 2003, ApJ, 588, 430 [NASA ADS] [CrossRef] [Google Scholar]
- Oostrum, L. C., van Leeuwen, J., Maan, Y., Coenen, T., & Ishwara-Chandra, C. H. 2020, MNRAS, 492, 4825 [Google Scholar]
- Pastor-Marazuela, I. 2022, PhD Thesis, University of Amsterdam, The Netherlands [Google Scholar]
- Pastor-Marazuela, I., Connor, L., van Leeuwen, J., et al. 2021, Nature, 596, 505 [NASA ADS] [CrossRef] [Google Scholar]
- Pastor-Marazuela, I., van Leeuwen, J., Bilous, A., et al. 2022, A&A, submitted [arXiv:2202.08002] [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pierbattista, M., Harding, A. K., Grenier, I. A., et al. 2015, A&A, 575, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pires, A. M., Haberl, F., Zavlin, V. E., et al. 2014, A&A, 563, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pitkin, M. 2018, J. Open Source Softw., 3, 538 [NASA ADS] [CrossRef] [Google Scholar]
- Pletsch, H. J., Guillemot, L., Allen, B., et al. 2013, ApJ, 779, L11 [NASA ADS] [CrossRef] [Google Scholar]
- Radhakrishnan, V., & Cooke, D. J. 1969, ApJ, 3, 225 [Google Scholar]
- Ransom, S. M. 2001, PhD Thesis, Harvard University, USA [Google Scholar]
- Ray, P. S., Kerr, M., Parent, D., et al. 2011, ApJS, 194, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Rigoselli, M., Mereghetti, S., Suleimanov, V., et al. 2019, A&A, 627, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Romani, R. W., & Watters, K. P. 2010, ApJ, 714, 810 [NASA ADS] [CrossRef] [Google Scholar]
- Rookyard, S. C., Weltevrede, P., & Johnston, S. 2015, MNRAS, 446, 3367 [NASA ADS] [CrossRef] [Google Scholar]
- Ruderman, M. A., & Sutherland, P. G. 1975, ApJ, 196, 51 [Google Scholar]
- Rutledge, R. E., Fox, D. B., & Shevchuk, A. H. 2008, ApJ, 672, 1137 [NASA ADS] [CrossRef] [Google Scholar]
- Sanidas, S., Cooper, S., Bassa, C. G., et al. 2019, A&A, 626, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shevchuk, A. S. H., Fox, D. B., & Rutledge, R. E. 2009, ApJ, 705, 391 [NASA ADS] [CrossRef] [Google Scholar]
- Shitov, Y. P., Pugachev, V. D., & Kutuzov, S. M. 2000, Pulsar Astronomy - 2000 and Beyond, ASP Conf. Ser., 202, 685 [NASA ADS] [Google Scholar]
- Stappers, B. W., Hessels, J. W. T., Alexov, A., et al. 2011, A&A, 530, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Szary, A., & van Leeuwen, J. 2017, ApJ, 845, 95 [NASA ADS] [CrossRef] [Google Scholar]
- Szary, A., van Leeuwen, J., Weltevrede, P., & Maan, Y. 2020, ApJ, 896, 168 [NASA ADS] [CrossRef] [Google Scholar]
- Tan, C. M., Bassa, C. G., Cooper, S., et al. 2018, ApJ, 866, 54 [NASA ADS] [CrossRef] [Google Scholar]
- van Haarlem, M. P., Wise, M. W., Gunst, A. W., et al. 2013, A&A, 556, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Leeuwen, J. 2014, in The Third Hot-wiring the Transient UniverseWorkshop (HTU-III), Santa Fe, NM, 79 [Google Scholar]
- van Leeuwen, J., & Stappers, B. W. 2010, A&A, 509, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Leeuwen, J., Kasian, L., Stairs, I. H., et al. 2015, ApJ, 798, 118 [NASA ADS] [CrossRef] [Google Scholar]
- van Leeuwen, J., Kooistra, E., Oostrum, L., et al. 2022, ArXiv e-prints [arXiv:2205.12362] [Google Scholar]
- Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389 [NASA ADS] [Google Scholar]
- Watters, K. P., Romani, R. W., Weltevrede, P., & Johnston, S. 2009, ApJ, 695, 1289 [Google Scholar]
- Young, M. D., Manchester, R. N., & Johnston, S. 1999, Nature, 400, 848 [NASA ADS] [CrossRef] [Google Scholar]
- Zane, S., Haberl, F., Israel, G. L., et al. 2011, MNRAS, 410, 2428 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Parameters of the observed pulsars and observational setup of the observations in the LC3_036 proposal.
All Figures
Fig. 1. Flux density upper limits of this work at 150 MHz (filled symbols) with S/N = 5 for comparison to earlier searches of the same sources (empty symbols). Solid lines going through our upper limit estimates with spectral index α = −1.4 are overlaid to show the scaling of our sensitivity limits. Our limits are plotted slightly offset from the 150 MHz observing frequency (dashed line) for better visibility. The faded green marker for SGR J1907+0919 represents the detection claimed by Shitov et al. (2000). |
|
In the text |
Fig. 2. P − Ṗ diagram showing the location of the sources presented in this work. All pulsars from the ATNF Pulsar Catalogue (Manchester et al. 2005) are shown as gray dots. Different pulsar classifications are encircled by different symbols. The sources discussed in this work are shown as black stars, from left to right: J1412+7922, J1932+1916, J1932+1916, and J1907+0919. The orange shaded region is delimited by the death line, while the green shaded region is delimited by the photon splitting line. Plot generated with psrqpy (Pitkin 2018). |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.