Free Access
Issue
A&A
Volume 647, March 2021
Article Number A191
Number of page(s) 6
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202037702
Published online 01 April 2021

© ESO 2021

1. Introduction

The discovery of pulsars in 1967 was associated with the detection of regular (periodic) pulses of emission in the meter range of wavelengths. These pulses were similar to “weak sporadic interference” (Hewish et al. 1968). In 2006, pulsars were discovered that emitted irregular pulses (McLaughlin et al. 2006). These pulsars were called rotating radio transients (RRATs). Time intervals between RRAT pulses can range from tens of seconds to hours. It is still not completely clear whether these pulsars are the same objects as canonical pulsars. The search for an evolutionary relationship between canonical pulsars and RRATs was considered in papers Keane & Kramer (2008) and Keane et al. (2011).

There are several hypotheses that attempt to explain the irregularity of RRAT emission. According to Weltevrede et al. (2006), RRATs can be ordinary pulsars with “extreme bursts of radio emission”. If the distance to a pulsar is great, the sensitivity of the radio telescope is not sufficient to detect periodic emission and to obtain an average profile. At the same time, the energy of individual giant pulses may be high enough to detect them. According to another hypothesis from Zhang et al. (2007), RRATs are pulsars with long (extreme) nullings. In the work of Burke-Spolaor & Bailes (2010), it is assumed that the duration of the “fading” window of pulsars increases in proportion to their lifetime. Thus, according to this hypothesis, RRATs represent a late evolutionary phase in the pulsar’s life. This work also shows that the hypothesis of giant pulses is not applicable for at least some of the observed fast radio transients. For example, the source J0941-39 is, for part of the time, observed as a RRAT, and part of the time as pulsar. In Wang et al. (2007), it is assumed that RRAT are an extreme type of mode-switching pulsars. In the papers Li (2006) and Luo & Melrose (2007), it is assumed that RRATs are extinct pulsars that are turned of the interaction of the pulsar’s magnetosphere with the incident matter or the interaction of the magnetosphere with the surrounding matter.

Thus, the only certain thing is that RRATs are a subclass of pulsars that are detected by the emission of their individual pulses, but the reason for this sporadic emission remains unclear. It is possible that different hypotheses are valid for different RRATs, and that RRAT is an intermediate class compared to different types of known pulsars. There are few comprehensive RRAT studies, and therefore it is not yet possible to discard any hypotheses.

Research into rotating radio transients is difficult. Because the appearance of pulses is unpredictable, and the time interval between detected pulses can reach up to several hours, a long-time observation is necessary for their statistic accumulation. It is known that to obtain an average pulsar profile, which is one of the basic characteristics of a pulsar, you need to accumulate several hundred, or better, several thousand pulses (Lorimer & Kramer 2004). If the RRAT emits an average of one pulse per minute, then it takes more than 16 hours of continuous observations to accumulate 1000 pulses and obtain an average profile. Until now, RRAT detection takes place on the world’s best telescopes, and the allocation of dozens of hours for a comprehensive study of even one RRAT is an expensive task.

Since August 2014 on the Large Phased Array of P.N.Lebedev Physical Institute (LPA LPI), round-the-clock observations have been carried out, including those used to search for pulsars (Tyul’bashev et al. 2016). The processing of semi-annual data revealed 33 new RRAT (Tyul’bashev et al. 2018a; Tyul’bashev et al. 2018b). For one of these rotating transients (RRAT J0139+33, also known as PSR J0139+3336 Michilli et al. 2020), a lot of pulses were detected in the blind search. J0139+33 has high peak flux densities (Speak) in individual pulses. Therefore, we chose this transient to test the hypotheses about the nature of RRATs.

2. Observations and data processing

The LPA LPI telescope is a Meridian-type radio telescope with a filled aperture and is a phased array. The telescope’s receiving elements consist of 16 384 half-wave dipoles. The size of the antenna is 187 m in the east-west direction and 384 m in the north-south direction. After a major upgrade of the antenna, its central frequency is 110.3 MHz, and its full band is 2.5 MHz. Several independent radio telescopes were created on the basis of a single antenna field. Details about the antenna upgrade and the main observation programs for its multibeam implementation are given in the works of Shishov et al. (2016) and Tyul’bashev et al. (2016).

Observations of the monitoring program are carried out simultaneously in 96 beams of the LPA LPI antenna. The data are recorded by a 32-channel digital recorder with a 78 kHz channel bandwidth and a readout time of τ = 12.5 ms. The size of the beam is approximately 0.5° ×1°. The instant observation area of LPA when operating in the monitoring mode is about 50 square degrees. An area of 17 000 square degrees is viewed per 24 hours, and each point in the sky is observed once in 24 hours for 3.5 min at half the power of the beam when passing through the meridian at a declination of zero degrees. The accumulated total observation time at each point in the sky for one year is over 21 hours. A calibration signal is recorded six times in 24 hours in an “OFF-ON-OFF” manner (calibration step) using a known temperature, allowing us to equalize the gain in the frequency channels (see details in Tyul’bashev et al. 2020).

While searching for new RRATs on LPA LPI, a candidate was recorded in the catalog for further verification if the signal-to-noise ratio (S/N) captured by the program was ≥7 after compensation of the dispersion delay (Tyul’bashev et al. 2018a). Initially, this S/N level was chosen based on the considerations that these signals are reliably detected visually on dynamic spectra. At the same time, for RRATs that have a defined period (P0) and a dispersion measure (DM), it is possible to reprocess the data and purposefully search for weaker pulses. These pulses may be hardly distinguishable or not at all visible on the dynamic spectra, but they are reliably detected when taking the known period into account.

There are two main problems incurred by the search for transient pulses. The first one is an interference removal. Most of interferences are of an impulse nature. We compared the S/N of detected pulses on DM = 0 pc ⋅ cm−3 and on the DM of the transient, which allowed us to remove interferences. The second problem is proper background subtraction. In the initial RRAT search, linear approximation was used to subtract the background signal. The time interval at which the subtraction was performed was approximately 15 time seconds. After removing the interference, the noise root dispersion was evaluated over the entire interval. The obtained estimate of the noise root-mean-square deviations was used to determine S/N ≥ 7 for observed peak signals. At the same time, if a fifteen-second recording still contains interference that lasts more than a few seconds or scintillations from background sources, the resulting noise dispersion estimate will be incorrect, as will the S/N estimate. In particular, the characteristic timescale of scintillations on interplanetary plasma at a frequency of 111 MHz is comparable in duration to the recording of the RRAT signal (if DM < 100 pc · cm−3) in the dynamic spectrum. If a single scintillation coincides with the time of arrival of the transient pulse, a situation may arise where the background signal must be subtracted at an interval of a few seconds.

The current implementation of the background subtraction algorithm is that the background signal is determined on an array of 32-channel data added with a zero dispersion measure after removing pulse interferences. A power polynomial was fit into the resulting array. It is believed that this polynomial fits the background signal. A check for polynomials of different powers showed that the most precise fit is for polynomials with a power greater than seven. An illustration of the background subtraction can be found in the Fig. 1.

thumbnail Fig. 1.

Vertical scale in the upper line of figures on the dynamic spectrum shows the frequencies used in 32-channel observations. The horizontal scale of the lower line of drawings shows the recording time. The upper part of the figure shows the dynamic spectrum, and the lower part of the figure shows the pulse profile with height in arbitrary units (a.u.). The blue color shows the pulse obtained by adding up the channels represented in the dynamic spectra after taking the dispersion measure into consideration, and the red color shows the signal added up with the zero dispersion measure. The detected pulse has S/N = 18.3.

At the next step, an array of 32-channel data was added, taking the dispersion measure of the studied transient into account, and a power polynomial was subtracted from the resulting array. The noise dispersion was determined at an interval of four seconds after excluding the signal from the transient itself.

There were two main criteria for selecting the pulses of the RRAT under consideration. First, the S/N of detected pulses should be as follows: S/N ≥ 7. Second, the peak height in an array added with a known dispersion measure of RRAT is 1.5 times higher than in an array added with a zero dispersion measure of variance. The second criterion allowed us to remove interference with a duration of more than a second. As an illustration of the search algorithms, we present three pulses from J0139+33 detected in the data of August 5, 2017 (see Fig. 2).

thumbnail Fig. 2.

From left to right: three pulses from RRAT J0139+33, detected in a single observation session, with S/N = 29; 10; 11. The horizontal and vertical axes are similar to the axes in Fig. 1.

We performed a pulse search at an S/N = 4. Extremely weak pulses that are still visible on dynamic spectra have S/N < 5 (see Fig. 3). However, with these S/N, false detections appear, which are practically absent for transients with S/N > 7. Therefore, pulses with S/N < 7 were excluded from further statistical analysis.

thumbnail Fig. 3.

One of the weak pulses from J0139+33. The horizontal and vertical axes are similar to the axes in Fig. 1. The detected pulse has an S/N = 4.9.

3. Results

The search for all pulses that fall into the beam of LPA is for RRAT J0139+33, with P0 = 1.2479 s, DM = 21.23 pc cm−3, D = 0.98 kpc (Sanidas et al. 2019), and for the canonical pulsar B0320+39, with P0 = 3.032 s, DM = 26.18 pc cm−3, D = 0.95 kpc (ATNF catalog)1. Both studied objects are second pulsars, have close dispersion measures, and are located at the same distances.

Because the LPA LPI is a phased array, the shape of its beam obeys the dependence (sin(x)/x)2. The consequence is that the transient pulse recorded at the top of the emission pattern (the moment the source passes through the meridian) and the same energy pulse recorded on the slope of the emission pattern will have different responses (different heights). The coordinates of the objects under study are known, therefore it is possible to compensate for the heights of the detected pulses, taking the shape of the LPA beam into account. However, the correction of the beam can be high in the area of the zeros of the beam. Additionally, the pulse coordinates for cases of particularly strong ionospheric scintillations can be shifted by tens of time seconds, and, therefore, the formal application of the correction to compensate for the height of the detected pulse can increase the pulse height many times and be incorrect. To draw the dependences, we used pulses found near the maximum of the emission pattern in the time interval ±1.5m from the top of the source. The correction of the height of the detected pulses does not exceed 40%, and therefore we can expect the resulting histograms not to be distorted. The ratio of peak flux densities of the strongest to the weakest (S/N = 7) detected pulses for RRAT J0139+33 is 22 (the strongest pulse S/N = 154), and for PSR B0320+39 it is 37.4 (the strongest pulse S/N = 262).

A total of 60830 (S/N ≥ 7) B0320+39 pulsar pulses were detected in 885 observation sessions and 1595 (S/N ≥ 7) J0139+33 transient pulses in 1669 sessions. On the interval ±1.5m from the time of passing through the Meridian of the studied objects, 39377 pulses for PSR B0320+39 were found, and for RRAT J0139+33, 1013 pulses were found. The average occurrence frequency of pulses from J0139+33 is pulse for 4.9 min of recording. The total loss of pulses for B0320+39 is 26%, and for J0139+33, it is 99.58%. A random check shows that in high-quality records, the number of registered B0320+39 pulsar pulses is almost the same as the number of periods. For illustration, Fig. 4 shows a typical records of one of the days with the highest number of registered pulses for B0320+39 and J0139+33.

thumbnail Fig. 4.

1.5 min of recording after subtracting the background signal and compensating for known dispersion measures for PSR B0320+39 and RRAT J0139+33. Only one pulse is not visible in the record for the pulsar, which is located between the points 109 000 and 110 000. For the transient, the recording for August 05, 2017 is shown. Three pulses are visible on the recording. The dynamic spectra of these pulses and the average profiles are shown in Fig. 1. The time interval between pulses labeled “1” and “3” is 67.44 s, which corresponds to 54 periods of J0139+33. The time interval between pulses labeled “2” and “3” is equal to one period. The vertical axis indicates the intensity in arbitrary units (a.u.) after removing the background signal. On the horizontal axis, the time is counted at points from the beginning of the observation session. One point corresponds to 12.5 ms. In both pictures, the intensity is normalized by the calibration step.

For all pulses detected in the ±1.5m interval of the maximum of the beam, distributions of peak flux densities expressed in signal-to-noise units were drawn (Fig. 5). It is clear that the distributions for B0320+39 and J0139+33 are different. The most precise fitting with a broken power-law dependence with exponents 0.4 and 1.8 was shown for J0139+33. Table 1 features some results of model testing using the standard deviation and the Akaike information criterion (AIC2) corrected for small sample size (AICc3) in the cases of lognormal and broken power dependence. The broken power dependence fits better for J0139+33 than lognormal dependence.

thumbnail Fig. 5.

Distribution of peak flux densities in S/N units on the logarithmic scale for the pulsar B0320+39 (left), for transient J0139+33 (right). The vertical axis shows the number of detected pulses, and the horizontal axis shows the S/N of these pulses. The histograms were drawn in increments of three on the abscissa axis. The figure shows that the distributions are different, despite the similar values of the dispersion measures, distances, and flux densities of the strongest observed pulses in the second pulsar B0320+39 and the second transient J0139+33. The lognormal distribution with probable power tail (see the fragment of the figure) was fit for B0320+39, and the broken power dependence was fit for J0139+33 (see the fragment of the figure). Standard deviations are shown in the upper corner of B0320+39.

Table 1.

Standard deviations and AICc of J0139+33.

The distribution of pulse energy for PSR B0320+39 was studied earlier at LPA LPI (Kazantsev & Potapov 2018). In Kazantsev & Potapov (2018), it was shown that giant pulses from this pulsar are unlikely. Among the data for RRAT J0139+33, the broken power distribution fits the most precisely, although a visible spread of points it seems to be associated with a low number of detected pulses.

4. Discussion of results

There are a number of signs for giant pulses. These include a high peak flux density in comparision with the average flux density of a pulsar, the power law of the pulse distribution, a short duration or the same time of an extremely short duration of a giant pulse, high-degree linear and circular polarization, localization at the longitude of an average pulse profile, and many more (see, e.g., Sutton et al. 1971; Kinkhabwala & Thorsett 2000; Soglasnov et al. 2004; Hankins & Eilek 2007; Kazantsev & Potapov 2018). Unfortunately, with our monitoring observations we can test only pulse distribution and estimate peak fluxes of pulses.

In early works, flux density distributions for RRAT pulses were already developed. The distributions for the four RRAT was obtained in a pioneering work by McLaughlin et al. (2006). For two of the four sources, the number of detected pulses was very low and amounted to 11 and 27 pulses, but the authors believe that the distribution of pulses by energy is described by a power function. For RRAT J1317-5759, the distribution is based on 108 pulses, and for RRAT J1819-1458, it is based on 229 pulses. The distribution of pulses by energy was also found to be a power function with an exponent of about one.

Other authors (Keane et al. 2010; Cui et al. 2017; Mickaliger et al. 2018; Meyers et al. 2019), however, when drawing pulse energy distributions for almost 30 RRATs, have found that in general the distributions are close to the lognormal or the sum of two lognormal distributions. Moreover, both transients defined in McLaughlin et al. (2006) as transients with a power distribution of pulses over energies are marked in Mickaliger et al. (2018) as transients of which the distribution consists of two lognormal distributions. The authors associate the lognormal distribution, which has a maximum at high flux densities, with possible flash activity in these RRATs. In Mickaliger et al. (2018), it is noted that there are known flashing pulsars with a similar (double-peak) distribution. In Cui et al. (2017), it is noted that for some of the transients the distribution is better described by the sum of the lognormal and power distributions.

The study of the energy distribution of pulses for canonical pulsars was carried out in the works Burke-Spolaor et al. (2012) and Mickaliger et al. (2018). In these papers, distributions were drawn for more than 300 pulsars, and it was shown that for most of them the dependence is well described by the lognormal distribution or the sum of the lognormal and power distributions. That is, in general, the distributions for RRATs and for pulsars are similar. At the same time, the distribution is lognormal on the side of low flux densities, and a power dependence is often added to the tail of the distribution.

The dependence mentioned in the previous section for J0139+33 is described by a power law. A power dependence in the pulse distribution is observed for a millisecond pulsar in the Crab nebula. The pulsar in the Crab nebula is the most famous pulsar with giant pulses. Mickaliger et al. (2012) presented 13 estimates of the exponent obtained by different authors for the distribution of pulses by energy at frequencies from 112 MHz to 4850 MHz. The estimation of the exponent varies from 2.1 to 4.2, and its median value falls to 2.8. For second pulsars, a similar behavior is demonstrated by the B0950+08 pulsar in which giant pulses are observed (Singal 2001; Kazantsev & Potapov 2018). The broken power-law dependence is also shown for three pulsars (B0531+21, B0950+08, B1237+25) with giant pulses (Smirnova 2012; Kazantsev & Potapov 2017).

The hypotheses mentioned in the introduction can be divided into two types. Either RRATs are pulsars that do not have periodical emission between the observed sporadic pulses, or RRATs are pulsars that have very weak periodical emission between the sporadic pulses.

The search for regular emission at a frequency of 111 MHz of 16 rotating transients discovered at 1400 MHz was done on the LPA LPI radio telescope in 2010–2012 (Losovsky & Dumsky 2014). Regular emission was not detected in individual sessions. However, a search using the fast folding algorithm (FFA) followed by the averaging of all observation sessions revealed weak regular emission of most of the studied RRATs. The estimated peak flux density of detected RRATs ranges from 20 mJy to 90 mJy. The average peak flux density was 45 mJy.

Is it possible to estimate the flux density of J0139+33? Indirect estimates of pertaining to the absence of periodic emission of J0139+33 can be derived from the work on the search for pulsars in the summed up power spectra at a frequency of 111 MHz (Tyul’bashev et al. 2017; Tyul’bashev et al. 2020). These works used six-channel data with a total observation bandwidth of 2.5 MHz and a sampling time of 0.1 s. The total increase in S/N over the four-year interval was approximately 30 times. No harmonics corresponding to the transient period were found in the summed power spectra. This means that there was no regular periodic emission from the direction of J0139+33. Because we know the exact coordinate of the transient and the pulse duration, we can give a rough estimate of the integral flux density, which is based on the assumption that the root-mean-square deviations of noise during observations are 0.1 Jy, and that the transient is a canonical pulsar with constant, albeit weak, emission. According to Table 1 from the work (Tyul’bashev et al. 2016) showing the expected sensitivity in a single session of observations on the LPA LPI as 15–20 mJy for the worst case at the optimal reading time and S/N > 7. The upper limit of the integral flux density from J0139+33 after summation of power spectra can be given as Sint < 2 mJy at the observation frequency of 111 MHz.

Rotating radio transient J0139+33 was also detected by studying its individual pulses in a survey on the search for pulsars on the LOw Frequency ARray (LOFAR) system (Sanidas et al. 2019). No periodic emission was detected for the transient when using standard search programs. According to Fig. 3 in this paper, the sensitivity of LOFAR survey observations for pulsars located in the extragalactic plane with a dispersion measure of less than 100 pc cm−3 is estimated as 1.2 mJy at the observation frequency of 135 MHz.

Thus, both observations on the LPA LPI and observations on LOFAR indicate that the integral flux density of J0139+33 does not exceed 2 mJy at the frequency of 111 MHz. We suppose that J0139+33 is a canonical pulsar, its flux density Sint = 2 mJy, and the half-width of the average profile coincides with the half-width of the observed transient pulses. According to Sanidas et al. (2019), the pulsar period P = 1.2479 ms, and the half-width of the profile w50 = 25 ms. Then the upper estimate of the peak flux density in the average profile Speak = (1.2479/0.025)×2 mJy ≤ 100 mJy. If we convert S/N for the strongest detected pulses into units of flux density using data on the fluctuation sensitivity of the LPA LPI antenna (Tyul’bashev et al. 2016), the peak flux density will be equal to 46 Jy (S/N = 154). Thus, the observed intensity of the strongest observed pulse is at least 460 times higher than the expected average pulse of J0139+33, if the rotating transient still has a weak regular emission. This excess of energy in strong pulses compared to the average pulsar ones is also characteristic of giant pulses. If weak regular emission is detected in the future (see work of Losovsky & Dumsky 2014), the conclusion about the detection of giant pulses J0139+33 will be strengthened.

Regular emission from J0139+33 has not yet been detected, and therefore it is impossible to be completely sure that J0139+33 is a weak pulsar with giant pulses. Detecting regular emission from J0139+33 will provide solid proof that all the observed transient pulses are gigantic.

5. Conclusion

Based on observations at a frequency of 111 MHz, we describe the distribution of peak flux densities for J0139+33 in units of S/N, which is in fact a broken power-law dependence. Regular emission of J0139+33 is not detected up to the integral flux densities of 1–2 mJy. Finally, we show that J0139+33 is a rotating transient with giant pulses.


Acknowledgments

The authors express their gratitude for the help in the course of this work to A. N. Kazantsev, V. A. Potapov, M. A. Kitaeva, and L. B. Potapova. We thank the anonymous referee for a number of helpful suggestions that improved the paper.

References

  1. Burke-Spolaor, S., & Bailes, M. 2010, MNRAS, 402, 855 [NASA ADS] [CrossRef] [Google Scholar]
  2. Burke-Spolaor, S., Johnston, S., Bailes, M., et al. 2012, MNRAS, 423, 1351 [NASA ADS] [CrossRef] [Google Scholar]
  3. Cui, B. Y., Boyles, J., McLaughlin, M. A., & Palliyaguru, N. 2017, ApJ, 840, 5 [Google Scholar]
  4. Hankins, T. H., & Eilek, J. A. 2007, ApJ, 670, 693 [NASA ADS] [CrossRef] [Google Scholar]
  5. Hewish, A., Bell, S., Pilkington, J., Scott, P., & Collins, R. 1968, Nature, 217, 709 [NASA ADS] [CrossRef] [Google Scholar]
  6. Kazantsev, A. N., & Potapov, V. A. 2017, Astron. Rep., 61, 747 [Google Scholar]
  7. Kazantsev, A. N., & Potapov, V. A. 2018, Res. Astron. Astrophys., 18, 097 [Google Scholar]
  8. Keane, E. F., & Kramer, M. 2008, MNRAS, 391, 2009 [NASA ADS] [CrossRef] [Google Scholar]
  9. Keane, E., Ludovici, D., Eatough, R., et al. 2010, MNRAS, 401, 1057 [NASA ADS] [CrossRef] [Google Scholar]
  10. Keane, E. F., Kramer, M., Lyne, A. G., Stappers, B. W., & McLaughlin, M. A. 2011, MNRAS, 415, 3065 [NASA ADS] [CrossRef] [Google Scholar]
  11. Kinkhabwala, A., & Thorsett, S. E. 2000, ApJ, 535, 365 [NASA ADS] [CrossRef] [Google Scholar]
  12. Li, X.-D. 2006, ApJ, 646, L139 [NASA ADS] [CrossRef] [Google Scholar]
  13. Lorimer, D. R., & Kramer, M. 2004, Handbook of Pulsar Astronomy, 4 [Google Scholar]
  14. Losovsky, B. Y., & Dumsky, D. V. 2014, Astron. Rep., 58, 537 [Google Scholar]
  15. Luo, Q., & Melrose, D. 2007, MNRAS, 378, 1481 [NASA ADS] [CrossRef] [Google Scholar]
  16. McLaughlin, M., Lyne, A., Lorimer, D., et al. 2006, Nature, 439, 817 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  17. Meyers, B. W., Tremblay, S. E., Bhat, N. D. R., et al. 2019, PASA, 36, e034 [Google Scholar]
  18. Michilli, D., Bassa, C., Cooper, S., et al. 2020, MNRAS, 491, 725 [Google Scholar]
  19. Mickaliger, M. B., McLaughlin, M. A., Lorimer, D. R., et al. 2012, ApJ, 760, 64 [NASA ADS] [CrossRef] [Google Scholar]
  20. Mickaliger, M. B., McEwen, A. E., McLaughlin, M. A., & Lorimer, D. R. 2018, MNRAS, 479, 5413 [Google Scholar]
  21. Sanidas, S., Cooper, S., Bassa, C. G., et al. 2019, A&A, 626, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Shishov, V., Chashei, I., Oreshko, V., et al. 2016, Astron. Rep., 60, 1067 [NASA ADS] [CrossRef] [Google Scholar]
  23. Singal, A. K. 2001, Ap&SS, 278, 61 [Google Scholar]
  24. Smirnova, T. V. 2012, Astron. Rep., 56, 430 [Google Scholar]
  25. Soglasnov, V. A., Popov, M. V., Bartel, N., et al. 2004, ApJ, 616, 439 [NASA ADS] [CrossRef] [Google Scholar]
  26. Sutton, J. M., Staelin, D. H., & Price, R. 1971, M., 46, 97 [Google Scholar]
  27. Tyul’bashev, S., Tyul’bashev, V., Oreshko, V., & Logvinenko, S. 2016, Astron. Rep., 60, 220 [NASA ADS] [CrossRef] [Google Scholar]
  28. Tyul’bashev, S., Tyul’bashev, V., Kitaeva, M., et al. 2017, Astron. Rep., 61, 848 [NASA ADS] [CrossRef] [Google Scholar]
  29. Tyul’bashev, S. A., Tyul’bashev, V. S., Malofeev, V. M., et al. 2018a, Astron. Rep., 62, 63 [NASA ADS] [CrossRef] [Google Scholar]
  30. Tyul’bashev, S. A., Tyul’bashev, V. S., & Malofeev, V. M. 2018b, A&A, 618, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Tyul’bashev, S. A., Kitaeva, M. A., Tyul’basheva, G. E., Malofeev, V. M., & Tyul’bashev, V. S. 2020, Astron. Rep., 64, 526 [Google Scholar]
  32. Wang, N., Manchester, R. N., & Johnston, S. 2007, MNRAS, 377, 1383 [NASA ADS] [CrossRef] [Google Scholar]
  33. Weltevrede, P., Stappers, B. W., Rankin, J. M., & Wright, G. A. E. 2006, ApJ, 645, L149 [NASA ADS] [CrossRef] [Google Scholar]
  34. Zhang, B., Gil, J., & Dyks, J. 2007, MNRAS, 374, 1103 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Standard deviations and AICc of J0139+33.

All Figures

thumbnail Fig. 1.

Vertical scale in the upper line of figures on the dynamic spectrum shows the frequencies used in 32-channel observations. The horizontal scale of the lower line of drawings shows the recording time. The upper part of the figure shows the dynamic spectrum, and the lower part of the figure shows the pulse profile with height in arbitrary units (a.u.). The blue color shows the pulse obtained by adding up the channels represented in the dynamic spectra after taking the dispersion measure into consideration, and the red color shows the signal added up with the zero dispersion measure. The detected pulse has S/N = 18.3.

In the text
thumbnail Fig. 2.

From left to right: three pulses from RRAT J0139+33, detected in a single observation session, with S/N = 29; 10; 11. The horizontal and vertical axes are similar to the axes in Fig. 1.

In the text
thumbnail Fig. 3.

One of the weak pulses from J0139+33. The horizontal and vertical axes are similar to the axes in Fig. 1. The detected pulse has an S/N = 4.9.

In the text
thumbnail Fig. 4.

1.5 min of recording after subtracting the background signal and compensating for known dispersion measures for PSR B0320+39 and RRAT J0139+33. Only one pulse is not visible in the record for the pulsar, which is located between the points 109 000 and 110 000. For the transient, the recording for August 05, 2017 is shown. Three pulses are visible on the recording. The dynamic spectra of these pulses and the average profiles are shown in Fig. 1. The time interval between pulses labeled “1” and “3” is 67.44 s, which corresponds to 54 periods of J0139+33. The time interval between pulses labeled “2” and “3” is equal to one period. The vertical axis indicates the intensity in arbitrary units (a.u.) after removing the background signal. On the horizontal axis, the time is counted at points from the beginning of the observation session. One point corresponds to 12.5 ms. In both pictures, the intensity is normalized by the calibration step.

In the text
thumbnail Fig. 5.

Distribution of peak flux densities in S/N units on the logarithmic scale for the pulsar B0320+39 (left), for transient J0139+33 (right). The vertical axis shows the number of detected pulses, and the horizontal axis shows the S/N of these pulses. The histograms were drawn in increments of three on the abscissa axis. The figure shows that the distributions are different, despite the similar values of the dispersion measures, distances, and flux densities of the strongest observed pulses in the second pulsar B0320+39 and the second transient J0139+33. The lognormal distribution with probable power tail (see the fragment of the figure) was fit for B0320+39, and the broken power dependence was fit for J0139+33 (see the fragment of the figure). Standard deviations are shown in the upper corner of B0320+39.

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.