Free Access
Volume 633, January 2020
Article Number A84
Number of page(s) 12
Section Astronomical instrumentation
Published online 15 January 2020

© ESO 2020

1. Introduction

The Argentine Institute of Radio astronomy (IAR; Instituto Argentino de Radioastronomía1) was founded in 1962 as a pioneer radio observatory in South America with two 30 m parabolic single-dish radio antennas (Fig. 1). Antenna 1 (A1) saw its first light in 1966, whereas Antenna 2 (A2) was built later in 19772. The IAR’s initial purpose was to perform a high sensitivity survey of neutral hydrogen (λ = 21 cm) in the southern hemisphere; this survey ended satisfactorily in the year 2000 with high-impact publications in collaboration with German and Dutch institutions (Testori et al. 2001; Bajaja et al. 2005; Kalberla et al. 2005).

thumbnail Fig. 1.

View of IAR antennas, A2 (left) and A1 (right).

Although the IAR has been a center of intense scientific and technological activity since it was founded, the radio antennas have not been employed in any scientific project since 2001. For the first time in over fifteen years, the IAR antennas are being upgraded to conduct high-quality radio astronomy. The Pulsar Monitoring in Argentina3 (PuMA) team is a collaboration of scientists and technicians from the IAR and the Rochester Institute of Technology (RIT). The collaboration has been working for two years with both antennas, including the implementation of a dedicated backend, the construction of a brand new frontend for A2, in addition to the formation of human resources for observations, data analysis, and pulsar astrophysics. This project represents the first systematic pulsar timing observations in South America and the beginning of pulsar science in Argentina.

In this work, we present the new hardware developments and observations of IAR’s radio antennas from the last two years. In Sect. 2, we give an overview of the current state of the radio observatory, its radio interference environment (RFI), the atomic clock availability, and the new developments for the acquisition software. In Sect. 3, we describe the observational techniques and capabilities, both hardware and observational cadence, as well as their automation. The calibration of total flux densities and polarization are part of our current developments.

In Sect. 4, we give a description of the various scientific projects that are being carried out or will be in the near future with further hardware improvements. One of the major goals of the PuMA collaboration involves high-cadence monitoring of millisecond pulsars (MSPs). In particular, monitoring of PSR J0437−4715 (one of the closest MSPs) is of great importance for gravitational wave detection using pulsar timing array techniques. Other close-by MSPs are also a target of the LIGO-Virgo collaboration for the search of continuous gravitational waves. In addition, interstellar scintillation can be studied with the same kind of data used at different time scales. Also, transient phenomena, such as fast radio bursts (FRBs), magnetars, and pulsar glitches are recent additions to the goals of our observations. Lastly, we discuss in Sect. 5 the impact that our contributions to radio observations in the southern hemisphere may have on those areas of current astrophysical research, and the potential for near-term improvement in both hardware and software.

2. The IAR observatory

Located in the provincial park Pereyra Iraola near the city of La Plata, Buenos Aires, the IAR itself is located at (latitude) and (longitude), with local time UTC-3. The IAR observatory has two 30 m diameter single-dish antennas, A1 and A2, aligned on a north-south direction (Fig. 1), separated by 120 m. These radio telescopes cover a declination range of −90° < δ <  −10° and an hour angle range of two hours east-west, −2h <  t <  2h. The angular resolution at 1420 MHz is ∼30′. The block diagram of Fig. 2 represents the connections of the antennas with their different modules.

thumbnail Fig. 2.

Current setup of IAR antennas.

2.1. Antenna 1

Since 2004, several updates and repairs were made on A1, including a complete front-end repair in 2009 and a new set of positional encoders installed in 2014 to keep the tracking system up to date. In 2015, we installed a software defined radio (SDR) module to perform pulsar observations. The characteristics of the current frontend of A1 are listed in Table 1.

Table 1.

Parameters of the antennas and their receivers (frontend and software configuration).

The frontend has as a feeder system a orthomode transducer (OMT) with a quadridge type wave guide to coax transducer, using a 90° hybrid coupler it gives both circular polarization products (LHC and RHC) from the linear products. Currently only one circular polarization is used. We have inserted new band pass radio frequency (RF) filters into the receivers of both antennas; the rest of the analog chain corresponds to a standard heterodyne receiver. A1 benefits from a lower insertion loss filter than A2 in the range 1150–1450 MHz, model ZX75BP-1280+ from Mini-circuits4.

The backend or acquisition module is based on two SDRs model B205 from Ettus5 using a Xilinx Spartan-6 XC6SLX75 field-programmable gate array (FPGA). This allows us to acquire raw samples from the front-end intermediate frequency as voltages in time series with a maximum rate of 56 MHz per board, and a universal serial bus (USB) 3.0 for connectivity. The sample rate is the same as the analog bandwidth due to the internal frontend in the SDR module. As each receiver has two digitizer boards, each with 56 MHz of bandwidth; we can use them in two modes: (i) as consecutive bands, giving a total of 112 MHz of bandwidth for a single polarization; (ii) by adding the two polarizations of 56 MHz bandwidth in order to obtain total power. A1 currently uses the first configuration.

The surface of A1 consists of a solid area at the center of the parabolic surface, while the rest is made of perforated aluminum sheets. This configuration gives an aperture efficiency of 32.8% (Testori et al. 2001).

2.2. Antenna 2

A2 has a newly developed receiver that is fully operational since November 2018. The receiver, different from the one in A1, has the digitization stage directly in the RF band, having as a result less RF components and more RF bandwidth available at lower cost. The RF feeder uses a turnstile feed with two orthogonal circular polarization outputs, each output is connected to the SDRs through their low noise amplifiers and filters.

The backend of A2 has the same scheme as A1 but with a different configuration. For A2, both SDR boards take data of each circular polarization at the same time, frequency and bandwidth. The processing software adds both polarizations to obtain total power. At the present time, we are not processing the polarization products as Stokes parameters, though this will be implemented in the processing software shortly.

A2 has a wider range of sensitivity, including higher frequencies up to 1600 MHz, with the use of a filter model VBFZ-1400-S+, also from Mini-circuits6 (Fig. 3). The surface of A2 presents a grid mesh for all its surface, giving a worse figure for its aperture efficiency, therefore resulting in a different Gain; the estimated aperture efficiency for A2 is 30.0%. The characteristics of the current frontend in A2 are listed in Table 1.

thumbnail Fig. 3.

RF filters installed in A1 (ZX75BP-1280-S+) and A2 (VBFZ-1400-S+).

2.3. Clock synchronization

Clock synchronization of the digital boards is performed with the one pulse per second (PPS) signal of a global positioning system (GPS) disciplined oscillator with an accuracy of 1.16 × 10−12 (one-day average). In order to get a precise data time stamp, the PPS signal is used in the SDR boards to synchronize the first time sample with the exact second of the GPS time. Then, the acquisition software reads the computer clock, which is synchronized with a Network Time Protocol (NTP) and a PPS signal from the GPS directly connected to the kernel OS through a serial port.

At less than one kilometer away from the IAR the Argentine-German Geodetic Observatory is located, AGGO7, dedicated, among other research, to Very Large Baseline Interferometry observations of quasars for geodetic purposes. This kind of measurements requires a precise clock for the synchronization with other observatories. This is achieved using an hydrogen maser clock with a short time stability of 10−15 (Allan Variance). With a locally developed RF-Over-Fiber device (Mena et al. 2013), we receive the 10 MHz signal from the AGGO’s hydrogen maser using a fiber optic cable that connects both institutions. This signal is used to synchronize a test unit backend with the aim to compare results from the different time bases used at the moment. In the future, this shall be used to synchronize the data acquisition.

2.4. Radio-Frequency-Intereference environment

The IAR is located in a rural area outside La Plata, Buenos Aires. Although this is not a radio frequency interferences (RFI) quiet zone, the radio band from 1 GHz to 2 GHz has a low level of RFI activity that ensures the capability to do radio astronomy8, as confirmed by the latest RFI measurement campaign from December 2017 (Gancio et al. 2014) for over a month in which the 90 % of RFIs are detected below the −160 power spectral density (PSD) [dBW m−2 Hz−1] (Fig. 4). Moreover, the IAR has implemented a new protocol to provide a clean RFI local environment that is compatible with the current research and technical activity of the nearby antennas.

thumbnail Fig. 4.

Wide bandwidth Power Spectrum Density obtained from a one-month average of measurements to characterize the RFI environment at IAR.

Observations with A1 have shown that RFIs affect ∼10% of the daily observing time, although during night-time the RFI activity reduces significantly and compromises less than 1% of the observing time. In average, observations with A2 have a significant less amount of RFI due to its distance to the administrative offices and laboratories.

Narrow-band mitigation is performed with the rfifind package from PRESTO9. First, the software parses the data in pieces of a certain time width (1 s in our case) per frequency interval. Then, it identifies in each one of these pieces weather the total power is too high, the data have an abnormal standard deviation, or the average of the data is above some given threshold. In that case, a mask is applied to the data before it is processed. In Fig. 5 we present an example mask for A1 and A2 that shows the flagged data as a function of time and frequency channel. We note that A1 is more affected by random RFIs due to its local environment. These RFIs are usually mitigated during night time, out of office hours. In the case of A2, the RFIs are predominantly monochromatic and their impact can be mitigated by using a larger number of frequency channels. We also note that these RFIs were proven to be polarized, so it is not straightforward to compare the masks from A1 and A2. Moreover, the mask criteria applies differently in each antenna given their different sensibility. A more detailed analysis of this RFI environment is ongoing.

thumbnail Fig. 5.

Example of a 3.5 h RFI mask for A1 (top) and A2 (bottom) for a simultaneous observation using rfifind with -time 1 and default parameters; colored sections are masked out. The mask criteria acts differently in each antenna given their different sensibility, so the plots are not in the same scale. RFIs affect only 1.3% of this observation for A1 and less than 6.3% for A2. We note that the persistent RFIs in A2 are monochromatic and affect individual channels only (2 of 64 in this case). The frequency channels on the borders of the bandwidth in A2 are removed due to the design of the receiver.

2.5. Future upgrades

We plan to increase the bandwidth of the receivers. This implies: (i) an upgrade of the frontend to be able to operate in a frequency band from 1 GHz up to 2 GHz with insertion losses below 1 dB; (ii) a new backend based on the CASPER boards like the SNAP board10 with bandwidths up to 500 MHz for each polarization. Moreover, the new receiver will benefit from state-of-the-art low-noise-amplifiers and electronics that will allow to reduce the system temperature to Tsys <  50 K (using its cryogenic capability).

3. Observations

3.1. Software infrastructure

The acquisition software was developed entirely at the IAR in C language. It processes the raw voltage samples at the desired rate without losses, while transforming them into a time series of RF channels. The software uses a scheme of synchronized threads in order to read the time samples from the different boards, process the fast Fourier transform (FFT) products, do the time average, and separate the final data product into channels before writing to disk. This is performed while keeping the synchronization with the PPS signal from the GPS Disciplined Oscillator (GPSDO). The final product is a file in Filterbank (SIGPROC) format11, compatible with standard pulsar reduction software like PRESTO (Ransom 2001; Ransom et al. 2002, 2003). Figure 6 represents the software diagram used in the C code, and Table 1 summarizes the main parameters of the digital receiver and the configuration used on each antenna.

thumbnail Fig. 6.

Pulsar software data acquisition block.

Tracking and pointing systems of A1 and A2 are controlled remotely through the IAR server. A weather station and a video-camera help to control the IAR environment. Data acquisition is performed by a different computer, connected to the backends. After acquisition, Filterbank files are saved and processed in the IAR storage system. The raw data files are also transferred to a data center at the Rochester Institute for Technology (RIT-PuMA-DEN Lab) for backup.

3.2. Automation

We are developing a distributed software architecture to control both IAR radio telescopes. Our goal is to generate a modern, dynamic and heterogeneous system in which modularity is an essential part, both in the development and in the expansion of the tools available for the observatory.

We are working in upgrading the current client-server architecture by means of building a scalable and dynamic control software that consists on a series of simple modules that perform specific tasks. These modules can be orchestrated by states, events and messages passed to a controller software that has enough privileges to make decisions upon the running modules. This common communication interface/application programming interface (API) allows the use of graphical user interfaces (GUIs), command line interfaces (CLIs), and simple viewers. A scheme of the software architecture is shown in Fig. 7.

thumbnail Fig. 7.

Software architecture for the control of the radio telescopes.

Finally, we plan to develop a scheduler to fully automate observations in order to offer the whole observation pipeline to the scientific community. In addition, the IAR is preparing a public proposer’s interface together with online tools to assess the technical aspects of a requested observation and a remote monitoring during its performance.

3.3. Data processing

As we mentioned before, we apply PRESTO to process the Filterbank files acquired. First, we employ the rfifind routine to generate a mask, which allows us to remove the RFIs. Then, using that mask, we fold the data with prepfold. For this step, we use the ATNF catalog12 (Manchester et al. 2005), data as input. The outputs of the last routine are a set of prepfold (.PFD) files.

For the moment we are working on three different projects: detection of glitches, pulse time of arrival (ToA) extraction, and flux density measurement. In the case of a glitch search, we analyze if the pulsar observed period matches with the expected topocentric period for the observing date derived from the latest reported ephemeris in the ATNF catalog.

If it is not the case, we reprocess the Filterbank file doing a fitting for the new period (see Sect. 4.4). The extraction of ToAs is handled with the PSRCHIVE (Hotan et al. 2004) package pat. These ToAs are processed with TEMPO2, using a suitable template for the pulse profile, to compute residuals needed for our scientific goals (see Sect. 4.1). The third project is the calibration of pulsar flux densities. For that, we employ the diode tube in each of our radio telescopes and calibration sources, such as Hydra A. Figure 8 shows a summary of the whole reduction process of the data.

thumbnail Fig. 8.

Data processing summary.

Further post-processing is carried out with the freely available PyPulse package13. For instance, we make use of this package to compute the signal-to-noise ratio of our observations defined as the ratio between the mean pulse peak and the rms of the noise. The latter is calculated from an off-pulse window of size 1/8 of the total phase bins in which the integrated flux density is minimal. We refer to Lam et al. (2016) for further details.

3.4. Observational capabilities and testing

We have been observing with both antennas since December 2018 to test and calibrate them. Both radio telescopes are capable of observing sources for almost four hours on a daily basis. We note that bright pulsars with sharp-peaked profiles are the easiest to detect; this is straightforward from the standard formula for the expected signal-to-noise ratio (S/N; Lorimer & Kramer 2012):


where Smean, P, and W are the mean flux density, period and equivalent width of the pulses, respectively, G, B, np, and Tsys are the antenna gain, bandwidth, number of polarizations, and system temperature, respectively, and tobs is the effective observing time. Using Eq. (1) we can estimate which pulsars can be observed with IAR antennas. We introduce the parameters given in Table 1 in the equation and we set the maximum observing time to t = 220 min, from where we obtain a value for the expected S/N. Considering that a reliable detection can be achieved for those pulsars with a S/N >  10, we make use of the python tool psrqpy (Pitkin 2018) to select the ones that IAR’s antennas are capable to detect14.

Our main target for testing purposes is the bright MSP J0437−4715 (Fig. 9). This pulsar allows us to test the timing quality of both antennas thanks to its high timing stability, high brightness, and short spin period. In Fig. 9, we compare our observations of J0437−4715 in the pre-upgrade and post-upgrade configuration. The upgrade consisted in an increase in bandwidth from 20 to 112 MHz, the incorporation of a better band-pass filter, and the use of both digitizer boards.

thumbnail Fig. 9.

Top panel: signal-to-noise distribution as a function of the integration time for different observations of J0437−4715 before and after the upgrade in A1. Bottom panel: typical pulse profiles in the pre-upgrade and post-upgrade configuration.

We note that this pulsar shows an important variation of the flux due to scintillation (Osłowski et al. 2014), most notably in observations after the upgrade which have a larger bandwidth. We show the preliminary timing results in the next section.

3.5. Collaborations with other radio observatories in Argentina

One spin-off from the development of the digital back-end receiver is its use as a stand alone unit for radio astronomical observations. This unit has been tested successfully in a 35 m Deep Space Antenna (DSA) named CLTC-CONAE-Neuquen15, installed in Neuquen, Argentina. Test observations were carried out in S and X bands targeting both pulsars (Vela and J0437−4715) and continuum sources. This will enable the possibility to do simultaneous observations between the IAR radio telescopes and the DSA antenna at different frequencies for radio astronomical research.

4. Enabling science projects

In this section we describe several scientific projects that we are –or will be– able to perform with the refurbished IAR antennas. We present the current state of each project, some preliminary results, and projections with future hardware improvements.

4.1. Pulsar timing and gravitational waves

Millisecond pulsars show a remarkable rotational stability. This allows one to predict the TOA of their pulses with high precision over long periods of time. Given a physical model, we can compare the predicted TOAs, tP, with the actual observed TOAs in a certain reference frame, tO, and compute the residuals δt = tP − tO. The residuals contain information of the astrophysical system and small effects due to different processes that can be incorporated in the timing model. For millisecond pulsars in particular, the residuals are dominated by white-noise when a large-enough bandwidth is used to resolve accurately the dispersion measure and can reach the order δt <  1 μs.

With this idea, pulsar timing arrays (PTAs), consisting usually of tens of precisely-timed MSPs, can be used as a Galactic scale detector of low-frequency gravitational waves. The main goal is to detect a “stochastic background” of such low-frequency gravitational waves, originated from an ensemble of unresolved supermassive black hole binaries (SMBHBs). Specifically, the effect of this background would appear in the PTA as a particular spatial angular correlation of the ToAs from different pulsars given by the Hellings-Downs curve (Hellings & Downs 1983). Several physical effects need to be modeled and corrected in order to find the effect of gravitational waves on the ToAs; chiefly, the pulsar dynamics and intrinsic instabilities, timing delays due to the interstellar medium and solar wind (see Hobbs & Dai 2017, for a recent review). There are three main PTA collaborations, NANOGrav from North America (Arzoumanian et al. 2018), EPTA from Europe (Desvignes et al. 2016), and PPTA from Australia (Reardon et al. 2016), together with an international PTA consortium (IPTA; Verbiest et al. 2016) that coordinates common efforts. Observatories in China, South Africa, India, and Argentina plan to join the IPTA soon.

Bright MSPs, such as PSR J0437−4715, are excellent targets for IAR’s antennas. We are currently performing an almost four hours per day monitoring of PSR J0437−4715 at 1400 MHz. These observations are projected to increase the sensitivity of pulsar timing arrays by increasing the observing cadence by a factor 20−30 and hence be sensitive to closer to merger (or less massive) SMBHB systems, and to reach potentially detectable SMBHB (Zhu et al. 2015) in their host galaxies (e.g., Fig. 2 in Burt et al. 2011). These observations also provide an overall sky coverage together with Parkes and MeerKAT observatories on the southern hemisphere. Moreover, the almost four hours per day of data that IAR can provide of J0437−4715 allows to minimize statistical uncertainties due to jitter of the pulses, improving the timing quality (see Shannon et al. 2014; Lam 2018). In our current set-up, the most limitating factor is the bandwidth. Full details of the IAR’s contribution to J0437−4715 timing and future projections will be presented in an upcoming work.

Current observations of J0437−4715 at IAR lead to residuals root mean squares (rms) of 0.55 μs for A1, 0.81 μs for A2, and 0.78 μs for A1+A2 (introducing a jump for matching both data sets), as displayed in Fig. 10. These values of rms ≲ 1 μs are compatible with Table 1 of Burt et al. (2011) expectations, but still far from PTA reported 0.1 μs in Perera (2019) or the 0.04 μs optimal reachable (see Osłowski et al. 2011). The timing precision (i.e., rms statistics or the single TOA precision) will improve with a larger bandwidth (up to 1 GHz). In addition, the precision of the timing parameters will improve with the continuation of daily observations to accumulate long term data (several years). Other MSPs within the reach of IAR observatory are given in Table 2 (extracted from the ATNF pulsar catalog16, Manchester et al. 2005), detailing their barycentric periods, mean flux densities at 1400 MHz, binary models, dispersion measures, and numbers of glitches observed.

thumbnail Fig. 10.

Residuals of the ToAs of the pulsar J0437−4715 measured at IAR with A1 and A2. The residuals and rms have been computed with the TEMPO2 package (Hobbs et al. 2006) and error bars correspond only to template-fitting errors (i.e., no systematics considered).

Table 2.

Observable MSPs from IAR with the current setup (top) and with future upgrades (bottom).

The detection of gravitational waves from compact binaries, along with their electromagnetic counterparts, notably enhances our comprehension of astrophysical processes. In particular, the case of the binary neutron star merger observed by LIGO-Virgo (Abbott et al. 2017a) gave birth to full fledged multi-messenger astronomy. In connection with PTA detections of single sources, the modeling of accreting matter around merging SMBHBs and the characteristic features of their electromagnetic spectra is currently an extremely active research area (Bowen et al. 2018, 2019; d’Ascoli et al. 2018). Likewise, the gravitational waves from merging supermassive black holes can present distinctive features in each polarization, which can inform us about the strong precession of the binary systems (Lousto & Healy 2019). In addition to our contributions to PTAs collaboration, we plan to search (or at least place constraints) for continuous gravitational waves from individual SMBHBs (Zhu et al. 2015; Kelley et al. 2018). The cadence of daily rather than monthly observations leads to sensitivity to sources of gravitational waves being produced closer to the merger of the supermassive black holes by a factor 302/3 (see e.g., Blanchet et al. 1996), hence an order of magnitude increase in its amplitude h. In particular, for sources of a few billion solar masses at z = 1, we expect to reach Earth with gravitational strains oscillations of up to h ∼ 10−14 (see the case of QSO 3C 186 in Lousto et al. 2017). Detection of gravitational waves from individual SMBHB (Detweiler 1979; Hellings & Downs 1983; Zhu et al. 2015) (as opposed to the stochastic background) may lead to important clues about the formation and evolution of such sources and can be performed by studying a few very well timed pulsars, like J0437−4715.

Another interesting millisecond pulsar to study with the next-generation backend is PSR J2241−5236 (Keith et al. 2011). This pulsar is currently observed by Parkes and MeerKAT and will be essential for the IPTA because of its excellent timing quality (Reardon, priv. comm.). However, since it is a black widow pulsar with an orbital period of ∼3.5 h, it shows orbital noise that reduces sensitivity to gravitational waves. Such noise can be characterized and modeled with high-cadence observations. This makes this pulsar an excellent target for IAR as it can monitor the full orbit on a daily basis. A phase-resolved analysis is expected to be viable with the future improvements in bandwidth and sensitivity.

4.2. Targeted pulsar studies for continuous gravitational waves detection from laser interferometry

In addition to the most remarkable detection of merging binary black holes and binary neutron stars (The LIGO Scientific Collaboration et al. 2019), the LIGO-Virgo collaboration monitors over 200 pulsars, looking for continuous gravitational waves coming from any time-varying (quadrupolar or higher) deformation of the spinning neutron stars (Abbott et al. 2017b, 2019a). The criteria to choose those pulsars is that their periods be less than 0.1 s, so the frequency of the emitted waves is at the start frequency of the LIGO-Virgo sensitivity curve, that is, currently with a frequency above 20 Hz. An important benchmark is given by the spin-down limit, obtained by equating the (radio) observed slowdown spinning rate to the expected rate due to the loss of energy by gravitational waves (Palomba 2005). For the Crab (PSR J0534+2200) and Vela (PSR J0835−4510) pulsars, current upper gravitational waves bounds show that this limit is now surpassed by nearly an order of magnitude, while for other six studied pulsars this spin-down limit was also recently reached (Abbott et al. 2019a).

The timing of those targeted pulsars is very important as their ToAs are used to construct ephemeris to search at specific frequencies in the LIGO-Virgo data. It is also important to know if and when those pulsars have glitches. In fact, overlooking glitches has a negative impact on standard CW analyses (Ashton et al. 2017, 2018). While identifying glitches can help to better model (Prix et al. 2011) potential emitters of continuous gravitational waves (see Abadie et al. 2011, for the 2006 Vela glitch). During the O1/O2 LIGO observing runs, several of the following pulsars have displayed glitches, like PSR J0205+6449 during O1 and five others during O2, including Crab and Vela pulsars (Keitel et al. 2019); the others are PSR J1028−5819, PSR J1718−382, and PSR J0205+6449 (see Abbott et al. 2017b, 2019a).

The third generation of laser interferometer detectors, with increased sensitivity in the low frequency band (starting at a few Hertz), are potentially sensitive enough to hold a chance to observe continuous gravitational waves from selected MSPs or from younger glitching ones (Glampedakis & Gualtieri 2018).

In particular, the pulsar J0711−6830 (at a distance of 0.11 kpc) is within a factor of 1.3 of the spin-down limit (assuming a 1038 kg m2 standard moment of inertia), and is one of our targeted pulsars for future observation at IAR. Another close-to-Earth recycled MSP, and close to its spin-down limit, is PSR J0437−4715, which is already being daily followed up at IAR. Other pulsars of interest for LIGO-Virgo (Table 2 of Abbott et al. 2019a) that are on the reach of IAR’s (future) observation capabilities (see Sect. 3.4) are pulsars J1744−1134, J1643−1224, J2241−5236, J2124−3358, J1603−7202, J0900−3144, and J1730−2304 (as displayed in Table 2).

The criteria for pulsar selection for a direct detection of continuous gravitational waves is similar to that for pulsar timing arrays, since both require (preferable non-glitching) MSPs in order to extract signals from observations over years (Woan et al. 2018). But while pulsar timing prefers more stable, recycled pulsars, young pulsars with larger asymmetries would be stronger gravitational waves emitters. Pulsars located close to Earth are preferred both for direct gravitational waves, due to the larger amplitude of the waveform strain (inversely proportional to the Earth-pulsar distance), and for the pulsar timing arrays, in order to obtain a better signal (pulse profile) to noise (instrumental and interstellar media) ratio. These conditions suit well with IAR capabilities, with an added value of daily observations that allow for shorter time scales studies than most observatories.

4.3. Magnetars

Magnetars are isolated young neutron stars with very large magnetic fields (on the order of 1015 G). About 30 magnetars have been reported, but a much larger population is expected, given their transient nature. Fermi and Swift satellites observe them in soft gamma-rays and in X-rays associated with soft-gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs). These are characterized by energetic winds, intense radiation, and a decaying magnetic field on scales from days to months. For a recent overview of the observational (including radio) properties of magnetars we refer to Esposito et al. (2018) and, for its physical modeling, to Turolla et al. (2015).

A possible connection with superluminous supernovae (SL-SNe) has been speculated. The explosions of these particular SNe are an order of magnitude more luminous than standard SNe and may lead to the formation of a highly-magnetized and fast-spinning magnetar, which in turn may energize the supernova remnants (Inserra et al. 2013). Magnetars may also be related to FRBs (Eftekhari et al. 2019). LIGO has searched (Abbott et al. 2016) for coincident FRB and gravitational waves signals in its first generation runs (2007–2013), and for magnetar bursts during the advanced LIGO-Virgo observations (Abbott et al. 2019b).

A few magnetars can be detected pulsating in radio wavelengths. They tend to display numerous glitches, and sometimes anti-glitches (i.e., spin-down; Archibald et al. 2013). Therefore, to study their early behavior, they need to be monitored with high cadence. Magnetars appear in the upper-right corner of the P diagram, namely, they have large periods and period derivatives (Lorimer & Kramer 2012). Under the assumption that those pulsars only brake due to dipole radiation emission (), they show the highest magnetic fields known, which gives them their name. Their spectrum is roughly Sν ∝ ν−0.5, thus harder than that of regular pulsars, which have Sν ∝ ν−1.8 (see for instance Dai et al. 2015).

Another characteristic feature of magnetars is that their pulses do not stabilize in shape, as opposed to regular pulsars that reach stability after average of a few hundred pulses (although they might switch between a few of them on minutes to hours scales; Esposito et al. 2018). Also, magnetars can be very bright (reaching 10 Jy in the L-band) in radio and also have an on-off-on behavior (Esposito et al. 2018).

The magnetar XTE J1810−197 has experienced periods of activity in X-rays (Ibrahim et al. 2004) and in radio frequencies, being the first magnetar in which radio pulsations were detected (Camilo et al. 2006). After being in a radio-quiet state for several years (Camilo et al. 2016), this magnetar has recently experienced another outburst (Lyne et al. 2018). As an exploratory study we dedicated observing time to this object from Dec. 14, 2018 to Mar. 1, 2019 (del Palacio et al. 2018). Single-polarization observations with a bandwidth of 56 MHz centered at 1420 MHz revealed significant pulsating radio emission from XTE J1810−197 with a barycentric spin period of P = 5.54137(3) s on MJD 58466.615, consistent with the values reported in Lyne et al. (2018). Unfortunately, we could not derive polarization angles and calibrated flux densities with these measurements. The pulse profiles from Dec. 14 showed a complex structure of a short, strong peak preceded by a less intense and longer in duration precursor, as reported at other frequencies in Levin et al. (2019). In turn, the precursor peak is not visible on subsequent observations as shown in Fig. 11.

thumbnail Fig. 11.

Pulse profiles of the magnetar J1810−197 at 1400 MHz measured by IAR antennas at different epochs.

Magnetars can display strong linear polarization (Camilo et al. 2007; Torne et al. 2017) and their study will benefit from those kind of measurements at IAR. Daily observations allow to follow glitches and the search for a binary magnetar. The increase in sensitivity of the antennas will allow us to investigate possible patterns of stabilization of individual pulses. This search would benefit of an appropriate machine learning clustering algorithm running on our data bank.

4.4. Glitches and young pulsars

Although pulsars have extremely stable periods over time, some young pulsars are prone to have glitches: sudden changes in their period due to interior changes in the star. Discovered 50 years ago, nowadays almost 200 pulsars are known to glitch (Manchester 2018). Southern (Yu et al. 2013) and northern (Espinoza et al. 2011; Fuentes et al. 2017) based surveys provide comprehensive catalogs17.

Magnetars present the largest relative glitches in frequency (ν = P−1) with Δν/ν ∼ 10−6 while, for young pulsars Δν/ν ∼ 10−7 − 10−8, and for MSPs Δν/ν ∼ 10−11. The increase in frequency is generally followed by an exponential decrease that, lasting 10–100 days, tries to recover the pre-glitch period, though a permanent change remains.

While the steady slow down of the pulsar spin is most likely produced by magnetic braking taking place outside the neutron star, glitches are thought to be produced by the sudden coupling of a fast rotating superfluid core with the crust, transferring to it some of the core’s angular momentum and hence producing the decrease of the pulsar period. Details of the modeling of such coupling have been challenged (Andersson et al. 2012; Chamel 2013; Piekarewicz et al. 2014) and are a matter of current research (see Haskell & Melatos 2015, for a review on models of pulsar glitches).

The Vela Pulsar (PSR B0833−45/J0835−4510) is one of the most active pulsars in terms of glitching, counting 20 in the last 50 years. The latest glitch occurred recently, around MJD 58515, and was reported by Sarkissian et al. (2019). We briefly summarize the radio timing observations performed at IAR of this event, which we first reported in Lopez Armengol et al. (2019). As part of the commissioning and developing stage, regular observations of Vela with both antennas were restarted by the end of January 2019 after one month of inactivity. We observed Vela on Jan. 29 (MJD 58512.14) and obtained Pbary = 89.402260(7) ms, consistent with the available ephemeris before the glitch. After the new glitch was reported, we started a daily follow-up of the event starting on Feb. 04 (MJD 58518.15). The monitoring initially consisted of a combination of short (10–15 min) and long (60–220 min) observations. The reconstruction of the post-glitch ephemeris, shown in Fig. 12, yields a period jump of ΔP ∼ −0.241 μs, equivalent to a frequency jump of 3.0 × 10−5 Hz, that is consistent with the value estimated by Kerr (2019), within 7% error.

thumbnail Fig. 12.

Vela’s glitch on February 1st. 2019 as measured at IAR in the 1400 MHz band. Ephemeris for the pre-glitch epoch were taken from Sarkissian et al. (2017) and the glitch date from Kerr (2019). The errorbars at one σ level are taken directly from the output of PRESTO and are likely to represent an underestimation. The dispersion is larger for short observations (10–15 min) than for long ones (> 60 min). We highlight that at MJD 58487 the data points from A1 and A2 overlap and are hardly distinguishable.

IAR’s program of pulsar observations considers their follow up for up to four hours per day. Hence, there is a chance that during this collected data a glitch could be observed “live”. In the case of the very bright Vela pulsar it will be possible to observe single pulses. In order to do so, we need to achieve a higher sensitivity. Ongoing tests suggest that IAR antennas are currently capable of detecting the Vela pulsar with an integration time as small as 0.4 s (i.e., five pulses added) with a significance greater than 5σ. With the future improvements in the antennas receivers (Sect. 2.5), which include a combination of broader bandwidth and reduction of system temperature, it will be possible to study the dynamical spectra of single pulses. Other glitching pulsars currently within the reach of IAR observatory and their properties are given in Table 3 (adapted from ATNF18 Catalog, Manchester et al. 2005).

Table 3.

Potentially observable glitching pulsars from IAR.

4.5. Fast-radio-burst observations

Fast radio bursts are intense bursts of radio emission that can reach flux densities of 10’s of Jy, have a duration of milliseconds and exhibit the same characteristic dispersion sweep in frequency as radio pulsars. The dispersion measure of these radio bursts together with models of the intergalactic medium suggest that the FRBs may come from as distant as 1 Gpc. The first FRB 010724 was reported by Lorimer et al. (2007) looking at archival survey data from the Parkes observatory. The single dispersed pulse, followed a ∼ν−2-law, had a width of less than 5 ms and a dispersion measure consistent with a distance of D ∼ 400 Mpc, which clearly placed the source at extragalactic distances. Since then over 72 FRBs have been reported (Ravi 2019), but most notably there are cases in which they have been observed to repeat. The first one, FRB 121102, was observed by the Arecibo and then Green bank and other radio telescopes, and also in X-rays (Spitler et al. 2016; Scholz et al. 2016), in 2015. The second repeating FRB 180814.J0422+73, has recently been reported by the CHIME/FRB team (CHIME/FRB Collaboration 2019) (while most of the FRBs have been found in the 1400 MHz band, CHIME found its 13 FRBs in the 400–800 MHz band). FRB repeaters observations are very important in order to better understand the physics underneath, and possibly to rule out models leading to single catastrophic events. Statistical analysis presented by Ravi (2019) suggests that FRBs are more likely produced by repeating sources; thus, radio monitoring for repeating FRBs seems feasible and it is crucial to probe the nature of these sources. Notably, during the writing of this manuscript, eight new FRB repeaters have been reported by CHIME/FRB (The CHIME/FRB Collaboration 2019).

Crucial to understand the sources of such FRBs is the identification of an optical or X-ray counterpart. For FRB 121102 a supernova type I in a dwarf galaxy at z = 0.19 has been found in a coincidence position (Chatterjee et al. 2017; Bassa et al. 2017; Eftekhari et al. 2019) using VLA observations. There are other models for FRBs ranging from SGRs, merging of white dwarfs or neutron stars, collapsing supra-massive neutron stars to evaporating primordial black holes and cusps of superconducting cosmic strings. For a thorough current review of the FRB field see Cordes & Chatterjee (2019).

There are essentially two strategies to follow to observe FRBs with IAR antennas. One is to observe areas of the sky we expect to produce FRBs, like nearby galaxy clusters. The second strategy is to follow known FRBs to look for a repeater. A southern hemisphere search for repeating FRBs has been started by the Australian SKA pathfinder (Bhandari et al. 2019); a new FRB 180924 event was reported already in Bannister et al. (2019) to occur at ∼4 kpc from the center of a luminous Galaxy at z = 0.32.

With rates up to 10 000 FRB/sky/day as reported in Table 3 of Petroff et al. (2019) and two observing antennas with large available observing time, we can do a first estimation of the rate of potentially detectable FRBs by IAR as IARFRB ×ΔΩ/(4π), with . This first guess suggests a likely detection after ∼500 h of accumulated observations. However, considering the lower sensitivity of IAR antennas with respect to other observatories such as Parkes, a more realistic perspective is to have a serendipitous detection after ∼5000 h of observations, although a large uncertainty could be added due the unknown luminosity function of the FRBs.

We note that observations of transient events at IAR will have a great synergy with other local observational facilities, such as the 47 cm optical telescope Telescopio Rafael Montemayor (TRM), to be installed in Argentina in early 2020. This telescope will be fully automated and it will have a 1° ×1° field of view, which makes it highly suited to monitor target of opportunity sources, such as counterparts of FRBs. Also note that the MASTER-OAFA robotic telescope (OAFA observatory of San Juan National University, Argentina)19 is already working on fast response FRB follow ups.

The biggest problem faced when looking for FRBs is that we do not know when and where one will flash next. The radio sky has to be observed until one of those events is captured; it is thought that thousands of FRBs occur per day, so serendipitous detections are feasible. This implies that large amounts of raw radio-astronomical observations have to be accumulated and analyzed. Observations of 400 MHz of bandwidth, scanned every 100 μs and in frequency bins of 0.2 MHz lead to a typical amount of 300 GB per hour of observation, that is potentially 7.2 TB per day and 2.6 PB per year. A minimal requirement of bandwidth of 200 MHz and frequency bins of 0.4 MHz would add up to 0.65 PB per year. Considering that with the IAR observatory we can track objects in the sky for up to four hours per day, any single observation can contain FRB data and one can think of applying a machine learning algorithm to look for those kind of patterns in the accumulated raw data. As an example of the potential of such archival searches, recently Zhang et al. (2019) found a new FRB (FRB 010312) in the original archival data set of the first (FRB 010724) detection, almost two decades later.

4.6. Interstellar medium scintillation

Scintillation of radio signals from pulsars is significantly affected by scattering in the turbulent interstellar medium. Understanding and accurately modeling this source of noise is crucial to improve the sensitivity of the pulsar timing arrays to detect gravitational waves from SMBHBs.

A classic study based on a daily analysis of the flux density at 610 MHz of 21 pulsars for five years by Stinebring et al. (2000) yielded valuable information of the flux variations caused by interstellar effects through the mechanism of “refractive scintillation”. Recent detailed studies (Reardon et al. 2019; Reardon 2018) of the long-term changes in the diffractive scintillation pattern of the binary pulsar PSR J1141−6545 during six years allowed to improve the determination of its distance, D = 156.79 ± 0.25 pc, the binary parameters, and to give an estimate of its proper motion.

We can perform similar studies from IAR in the southern hemisphere, either by systematic studies of a selected set of pulsars via daily monitoring or by use of the information already collected for targeted pulsars for other projects, like for J0437−4715. Our observations would supplement those in other southern observatories, in particular those from Australia, providing a double cadence, every nearly 12 h.

4.7. Tests of gravity with pulsar timing

Pulsars in a binary system can be used as tests of alternative theories of gravity. They have been extensively used in the post-Newtonian parametrization approach to weak gravitational fields (see Will 2018, for a review). In particular one can test General Relativity versus modified theories (see Will 2018; Renevey 2019, for a current review).

Notably, in a recent paper Yang et al. (2017) proposed to use scintillation measurements of PSR B0834+06 to test predictions of alternative theories of gravity. This work could be made extensive to J0437−4715 and other accurately measured MSPs. Pulsar timing have also been proposed as a mean to detect and to shed some light on the constituents of dark matter (e.g., Khmelnitsky & Rubakov 2014).

Pulsar timing arrays can constrain alternative theories of gravity via its strong sensitivity to gravitational waves polarizations, particularly the longitudinal polarization of vector and scalar modes (Lee et al. 2008; Alves & Tinto 2011; Chamberlin & Siemens 2012; Cornish et al. 2018). SKA, with its extraordinary sensitivity will be able to be at the forefront of those tests of gravity (e.g., Weltman et al. 2018).

5. Conclusions

We have shown that with IAR’s upgrades we can now perform four hours of continuous observations of bright pulsars of the southern hemisphere with a daily cadence for very long periods of time. This gives an estimate of 1000 h year−1 pulsar−1 antenna−1. These capabilities allow us to contribute to the international pulsar timing array efforts, and to follow up targeted MSPs of the LIGO-Virgo collaboration. We have also been able to already monitor magnetar activity and pulsar glitches. Notably IAR’s location provides a 12 h complementary window with respect to Australian Parkes’ and over five hours of MeerKAT observations, and thus can uniquely cover transient phenomena such as FRBs or live-glitches.

In order to increase the number of pulsars and the accuracy of measurements for those projects we are developing a tighter calibration of the intensity and polarization observations (particularly important for magnetars measurements), automation of observations (for applications to interstellar scintillation measurements and glitches surveys, for instance), and increase of the RF bandwidth by an order of magnitude (with new/additional USRP plates matching IAR site RFI frequency windows) and data downloading bandwidth (RF through fiber optics) for applications to FRB searches, for instance.

We plan to store all our raw data and make it available to the astronomical community for archival post-processing and further analysis, both for continuous gravitational waves searches and interstellar scintillation studies as well as for transient phenomena such as FRBs, (mini-)glitches, magnetars, and other unexpected astrophysical phenomena.

Future perspectives, in addition to the upgrades to IAR’s radio antennas, include observation time in the two 35 m antennas located in Mendoza20 and in Neuquen21, Argentina. Regarding optical counterparts, we plan to perform simultaneous observations with the facilities already present at San Juan22, Argentina.

The detection of gravitational waves from supermassive black hole mergers can yield fundamental insight into their astrophysical formation and growth scenarios as well as the structure and evolution of the universe, and, ultimately, provide crucial tests for classical theories of gravity. IAR is now ready to make its modest contributions toward this goal.


In 2019, A1 and A2 were baptized “Varsavsky” and “Bajaja”, respectively, in honor to their contributions to the IAR.


Argentina is a member of the International Telecommunication Union that protects radio bands for astronomical observations.


We note that this estimated way to calculate the signal-to-noise corresponds to the S/N as computed with a cross-correlation function instead of the peak to off-peak rms S/B, see Lam et al. (2016) to the mathematical relation between the two quantities.


The authors thank various members of the IAR’s technical staff for their work, as well as numerous members of the NANOGrav and LIGO-Virgo collaborations for very valuable discussions, in particular D. Keitel and J. Romano for reading our original manuscript. Noteworthy have also been the additions and corrections suggested by an anonymous referee. Part of this work was supported by the National Science Foundation (NSF) from Grants No. PHY-1912632, No. PHY-1607520, and No. PHY-1726215. FG and JAC acknowledge support by PIP 0102 (CONICET). This work received financial support from PICT-2017-2865 (ANPCyT). JAC was also supported by the Agencia Estatal de Investigación grant AYA2016-76012-C3-3-P from the Spanish Ministerio de Economía y Competitividad (MINECO) and by the Consejería de Economía, Innovación, Ciencia y Empleo of Junta de Andalucía under research group FQM-322, as well as FEDER funds.


  1. Abadie, J., Abbott, B. P., Abbott, R., et al. 2011, Phys. Rev. D, 83, 042001 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, Phys. Rev. D, 93, 122008 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017a, ApJ, 848, L12 [Google Scholar]
  4. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017b, ApJ, 839, 12 [NASA ADS] [CrossRef] [Google Scholar]
  5. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2019a, ApJ, 879, 10 [NASA ADS] [CrossRef] [Google Scholar]
  6. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2019b, ApJ, 874, 163 [NASA ADS] [CrossRef] [Google Scholar]
  7. Alves, M. E. D. S., & Tinto, M. 2011, Phys. Rev. D, 83, 123529 [NASA ADS] [CrossRef] [Google Scholar]
  8. Andersson, N., Glampedakis, K., Ho, W. C. G., & Espinoza, C. M. 2012, Phys. Rev. Lett., 109, 241103 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  9. Archibald, R. F., Kaspi, V. M., Ng, C.-Y., et al. 2013, Nature, 497, 591 [NASA ADS] [CrossRef] [Google Scholar]
  10. Arzoumanian, Z., Brazier, A., Burke-Spolaor, S., et al. 2018, ApJS, 235, 37 [Google Scholar]
  11. Ashton, G., Prix, R., & Jones, D. I. 2017, Phys. Rev. D, 96, 063004 [NASA ADS] [CrossRef] [Google Scholar]
  12. Ashton, G., Prix, R., & Jones, D. I. 2018, Phys. Rev. D, 98, 063011 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bajaja, E., Arnal, E. M., Larrarte, J. J., et al. 2005, A&A, 440, 767 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Bannister, K. W., Deller, A. T., Phillips, C., et al. 2019, Science, 365, 565 [NASA ADS] [CrossRef] [Google Scholar]
  15. Bassa, C. G., Tendulkar, S. P., Adams, E. A. K., et al. 2017, ApJ, 843, L8 [NASA ADS] [CrossRef] [Google Scholar]
  16. Bhandari, S., Bannister, K. W., James, C. W., et al. 2019, MNRAS, 486, 70 [NASA ADS] [CrossRef] [Google Scholar]
  17. Blanchet, L., Iyer, B. R., Will, C. M., & Wiseman, A. G. 1996, Class. Quant. Grav., 13, 575 [NASA ADS] [CrossRef] [Google Scholar]
  18. Bowen, D. B., Mewes, V., Campanelli, M., et al. 2018, ApJ, 853, L17 [NASA ADS] [CrossRef] [Google Scholar]
  19. Bowen, D. B., Mewes, V., Noble, S. C., et al. 2019, ApJ, 879, 76 [NASA ADS] [CrossRef] [Google Scholar]
  20. Burt, B. J., Lommen, A. N., & Finn, L. S. 2011, ApJ, 730, 17 [NASA ADS] [CrossRef] [Google Scholar]
  21. Camilo, F., Ransom, S. M., Halpern, J. P., et al. 2006, Nature, 442, 892 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  22. Camilo, F., Reynolds, J., Johnston, S., et al. 2007, ApJ, 659, L37 [NASA ADS] [CrossRef] [Google Scholar]
  23. Camilo, F., Ransom, S. M., Halpern, J. P., et al. 2016, ApJ, 820, 110 [NASA ADS] [CrossRef] [Google Scholar]
  24. Chamberlin, S. J., & Siemens, X. 2012, Phys. Rev. D, 85, 082001 [NASA ADS] [CrossRef] [Google Scholar]
  25. Chamel, N. 2013, Phys. Rev. Lett., 110, 011101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  26. Chatterjee, S., Law, C. J., Wharton, R. S., et al. 2017, Nature, 541, 58 [NASA ADS] [CrossRef] [Google Scholar]
  27. CHIME/FRB Collaboration (Amiri, M., et al.) 2019, Nature, 566, 235 [NASA ADS] [CrossRef] [Google Scholar]
  28. Cordes, J. M., & Chatterjee, S. 2019, ARA&A, 57, 417 [NASA ADS] [CrossRef] [Google Scholar]
  29. Cornish, N. J., O’Beirne, L., Taylor, S. R., & Yunes, N. 2018, Phys. Rev. Lett., 120, 181101 [NASA ADS] [CrossRef] [Google Scholar]
  30. Dai, S., Hobbs, G., Manchester, R. N., et al. 2015, MNRAS, 449, 3223 [NASA ADS] [CrossRef] [Google Scholar]
  31. d’Ascoli, S., Noble, S. C., Bowen, D. B., et al. 2018, ApJ, 865, 140 [NASA ADS] [CrossRef] [Google Scholar]
  32. del Palacio, S., Garcia, F., Combi, L., et al. 2018, ATel, 12323, 1 [NASA ADS] [Google Scholar]
  33. Desvignes, G., Caballero, R. N., Lentati, L., et al. 2016, MNRAS, 458, 3341 [NASA ADS] [CrossRef] [Google Scholar]
  34. Detweiler, S. 1979, ApJ, 234, 1100 [NASA ADS] [CrossRef] [Google Scholar]
  35. Eftekhari, T., Berger, E., Margalit, B., et al. 2019, ApJ, 876, L10 [NASA ADS] [CrossRef] [Google Scholar]
  36. Espinoza, C. M., Lyne, A. G., Stappers, B. W., & Kramer, M. 2011, MNRAS, 414, 1679 [NASA ADS] [CrossRef] [Google Scholar]
  37. Esposito, P., Rea, N., & Israel, G. L. 2018, ArXiv e-prints [arXiv:1803.05716] [Google Scholar]
  38. Fuentes, J. R., Espinoza, C. M., Reisenegger, A., et al. 2017, A&A, 608, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Gancio, G., Larrarte, J. J., Diaz, E., Aquino, F., & Spagnolo, S. 2014, International VLBI Service for Geodesy and Astrometry 2014 General Meeting Proceedings: VGOS: The New VLBI Network, 210 [Google Scholar]
  40. Glampedakis, K., & Gualtieri, L. 2018, in Astrophysics and Space Science Library, eds. L. Rezzolla, P. Pizzochero, D. I. Jones, N. Rea, & I. Vidaña, 457, 673 [NASA ADS] [CrossRef] [Google Scholar]
  41. Haskell, B., & Melatos, A. 2015, Int. J. Mod. Phys. D, 24, 1530008 [NASA ADS] [CrossRef] [Google Scholar]
  42. Hellings, R. W., & Downs, G. S. 1983, ApJ, 265, L39 [NASA ADS] [CrossRef] [Google Scholar]
  43. Hobbs, G., & Dai, S. 2017, Natl. Sci. Rev., 4, 707 [NASA ADS] [CrossRef] [Google Scholar]
  44. Hobbs, G. B., Edwards, R. T., & Manchester, R. N. 2006, MNRAS, 369, 655 [NASA ADS] [CrossRef] [Google Scholar]
  45. Hotan, A. W., van Straten, W., & Manchester, R. N. 2004, PASA, 21, 302 [NASA ADS] [CrossRef] [Google Scholar]
  46. Ibrahim, A. I., Markwardt, C. B., Swank, J. H., et al. 2004, ApJ, 609, L21 [NASA ADS] [CrossRef] [Google Scholar]
  47. Inserra, C., Smartt, S. J., Jerkstrand, A., et al. 2013, ApJ, 770, 128 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Keitel, D., Woan, G., Pitkin, M., et al. 2019, Phys. Rev. D, 100, 064058 [NASA ADS] [CrossRef] [Google Scholar]
  50. Keith, M. J., Johnston, S., Ray, P. S., et al. 2011, MNRAS, 414, 1292 [NASA ADS] [CrossRef] [Google Scholar]
  51. Kelley, L. Z., Blecha, L., Hernquist, L., Sesana, A., & Taylor, S. R. 2018, MNRAS, 477, 964 [NASA ADS] [CrossRef] [Google Scholar]
  52. Kerr, M. 2019, ATel, 12481, 1 [NASA ADS] [Google Scholar]
  53. Khmelnitsky, A., & Rubakov, V. 2014, JCAP, 2, 019 [NASA ADS] [CrossRef] [Google Scholar]
  54. Lam, M. T. 2018, ApJ, 868, 33 [NASA ADS] [CrossRef] [Google Scholar]
  55. Lam, M. T., Cordes, J. M., Chatterjee, S., et al. 2016, ApJ, 819, 155 [NASA ADS] [CrossRef] [Google Scholar]
  56. Lee, K. J., Jenet, F. A., & Price, R. H. 2008, ApJ, 685, 1304 [NASA ADS] [CrossRef] [Google Scholar]
  57. Levin, L., Lyne, A. G., Desvignes, G., et al. 2019, MNRAS, 488, 5251 [NASA ADS] [CrossRef] [Google Scholar]
  58. Lopez Armengol, F. G., Lousto, C. O., del Palacio, S., et al. 2019, ATel, 12482, 1 [NASA ADS] [Google Scholar]
  59. Lorimer, D. R., & Kramer, M. 2012, Handbook of Pulsar Astronomy (Cambridge: Cambridge University Press) [Google Scholar]
  60. Lorimer, D. R., Bailes, M., McLaughlin, M. A., Narkevic, D. J., & Crawford, F. 2007, Science, 318, 777 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  61. Lousto, C. O., & Healy, J. 2019, Phys. Rev. D, 99, 064023 [NASA ADS] [CrossRef] [Google Scholar]
  62. Lousto, C. O., Zlochower, Y., & Campanelli, M. 2017, ApJ, 841, L28 [NASA ADS] [CrossRef] [Google Scholar]
  63. Lyne, A., Levin, L., Stappers, B., et al. 2018, ATel, 12284, 1 [NASA ADS] [Google Scholar]
  64. Manchester, R. N. 2018, in Pulsar Astrophysics the Next Fifty Years, eds. P. Weltevrede, B. B. P. Perera, L. L. Preston, & S. Sanidas, IAU Symp., 337, 197 [NASA ADS] [Google Scholar]
  65. Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993 [NASA ADS] [CrossRef] [Google Scholar]
  66. Mena, J., Bandura, K., Cliche, J. F., et al. 2013, J. Instrum., 8, T10003 [CrossRef] [Google Scholar]
  67. Osłowski, S., van Straten, W., Hobbs, G. B., Bailes, M., & Demorest, P. 2011, MNRAS, 418, 1258 [NASA ADS] [CrossRef] [Google Scholar]
  68. Osłowski, S., van Straten, W., Bailes, M., Jameson, A., & Hobbs, G. 2014, MNRAS, 441, 3148 [NASA ADS] [CrossRef] [Google Scholar]
  69. Palomba, C. 2005, MNRAS, 359, 1150 [NASA ADS] [CrossRef] [Google Scholar]
  70. Perera, B. B. P., DeCesar, M. E., Demorest, P. B., et al. 2019, MNRAS, 490, 4666 [NASA ADS] [CrossRef] [Google Scholar]
  71. Petroff, E., Hessels, J. W. T., & Lorimer, D. R. 2019, A&ARv, 27, 4 [NASA ADS] [CrossRef] [Google Scholar]
  72. Piekarewicz, J., Fattoyev, F. J., & Horowitz, C. J. 2014, Phys. Rev. C, 90, 015803 [NASA ADS] [CrossRef] [Google Scholar]
  73. Pitkin, M. 2018, J. Open Source Softw., 3, 538 [NASA ADS] [CrossRef] [Google Scholar]
  74. Prix, R., Giampanis, S., & Messenger, C. 2011, Phys. Rev. D, 84, 023007 [NASA ADS] [CrossRef] [Google Scholar]
  75. Ransom, S. M. 2001, in American Astronomical Society Meeting Abstracts, BASS, 33, 1484 [NASA ADS] [Google Scholar]
  76. Ransom, S. M., Eikenberry, S. S., & Middleditch, J. 2002, AJ, 124, 1788 [Google Scholar]
  77. Ransom, S. M., Cordes, J. M., & Eikenberry, S. S. 2003, ApJ, 589, 911 [NASA ADS] [CrossRef] [Google Scholar]
  78. Ravi, V. 2019, Nat. Astron., 3, 928 [NASA ADS] [CrossRef] [Google Scholar]
  79. Reardon, D. J. 2018, PhD Thesis, Monash University, Victoria 3800, Australia [Google Scholar]
  80. Reardon, D. J., Hobbs, G., Coles, W., et al. 2016, MNRAS, 455, 1751 [NASA ADS] [CrossRef] [Google Scholar]
  81. Reardon, D. J., Coles, W. A., Hobbs, G., et al. 2019, MNRAS, 485, 4389 [NASA ADS] [CrossRef] [Google Scholar]
  82. Renevey, C. 2019, ArXiv e-prints [arXiv:1905.13720] [Google Scholar]
  83. Sarkissian, J. M., Reynolds, J. E., Hobbs, G., & Harvey-Smith, L. 2017, PASA, 34, e027 [NASA ADS] [CrossRef] [Google Scholar]
  84. Sarkissian, J., Hobbs, G., Reynolds, J., Palfreyman, J., & Olney, S. 2019, ATel, 12466, 1 [NASA ADS] [Google Scholar]
  85. Scholz, P., Spitler, L. G., Hessels, J. W. T., et al. 2016, ApJ, 833, 177 [NASA ADS] [CrossRef] [Google Scholar]
  86. Shannon, R. M., Osłowski, S., Dai, S., et al. 2014, MNRAS, 443, 1463 [NASA ADS] [CrossRef] [Google Scholar]
  87. Spitler, L. G., Scholz, P., Hessels, J. W. T., et al. 2016, Nature, 531, 202 [NASA ADS] [CrossRef] [Google Scholar]
  88. Stinebring, D. R., Smirnova, T. V., Hankins, T. H., et al. 2000, ApJ, 539, 300 [NASA ADS] [CrossRef] [Google Scholar]
  89. Testori, J. C., Reich, P., Bava, J. A., et al. 2001, A&A, 368, 1123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. The CHIME/FRB Collaboration (Andersen, B. C., et al.) 2019, ApJ, 885, L24 [NASA ADS] [CrossRef] [Google Scholar]
  91. The LIGO Scientific Collaboration, the Virgo Collaboration, Abbott, B. P., et al. 2019, Phys. Rev. X, 9, 031040 [Google Scholar]
  92. Torne, P., Desvignes, G., Eatough, R. P., et al. 2017, MNRAS, 465, 242 [NASA ADS] [CrossRef] [Google Scholar]
  93. Turolla, R., Zane, S., & Watts, A. L. 2015, Rep. Prog. Phys., 78, 116901 [Google Scholar]
  94. Verbiest, J. P. W., Lentati, L., Hobbs, G., et al. 2016, MNRAS, 458, 1267 [NASA ADS] [CrossRef] [Google Scholar]
  95. Weltman, A., Bull, P., Camera, S., et al. 2018, PASA, submitted [arXiv:1810.02680] [Google Scholar]
  96. Will, C. M. 2018, Theory and Experiment in Gravitational Physics (Cambridge: Cambridge University Press) [CrossRef] [Google Scholar]
  97. Woan, G., Pitkin, M. D., Haskell, B., Jones, D. I., & Lasky, P. D. 2018, ApJ, 863, L40 [NASA ADS] [CrossRef] [Google Scholar]
  98. Yang, H., Nishizawa, A., & Pen, U.-L. 2017, Phys. Rev. D, 95, 084049 [NASA ADS] [CrossRef] [Google Scholar]
  99. Yao, J. M., Manchester, R. N., & Wang, N. 2017, ApJ, 835, 29 [NASA ADS] [CrossRef] [Google Scholar]
  100. Yu, M., Manchester, R. N., Hobbs, G., et al. 2013, MNRAS, 429, 688 [NASA ADS] [CrossRef] [Google Scholar]
  101. Zhang, S.-B., Hobbs, G., Dai, S., et al. 2019, MNRAS, 484, L147 [NASA ADS] [CrossRef] [Google Scholar]
  102. Zhu, X.-J., Wen, L., Hobbs, G., et al. 2015, MNRAS, 449, 1650 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Parameters of the antennas and their receivers (frontend and software configuration).

Table 2.

Observable MSPs from IAR with the current setup (top) and with future upgrades (bottom).

Table 3.

Potentially observable glitching pulsars from IAR.

All Figures

thumbnail Fig. 1.

View of IAR antennas, A2 (left) and A1 (right).

In the text
thumbnail Fig. 2.

Current setup of IAR antennas.

In the text
thumbnail Fig. 3.

RF filters installed in A1 (ZX75BP-1280-S+) and A2 (VBFZ-1400-S+).

In the text
thumbnail Fig. 4.

Wide bandwidth Power Spectrum Density obtained from a one-month average of measurements to characterize the RFI environment at IAR.

In the text
thumbnail Fig. 5.

Example of a 3.5 h RFI mask for A1 (top) and A2 (bottom) for a simultaneous observation using rfifind with -time 1 and default parameters; colored sections are masked out. The mask criteria acts differently in each antenna given their different sensibility, so the plots are not in the same scale. RFIs affect only 1.3% of this observation for A1 and less than 6.3% for A2. We note that the persistent RFIs in A2 are monochromatic and affect individual channels only (2 of 64 in this case). The frequency channels on the borders of the bandwidth in A2 are removed due to the design of the receiver.

In the text
thumbnail Fig. 6.

Pulsar software data acquisition block.

In the text
thumbnail Fig. 7.

Software architecture for the control of the radio telescopes.

In the text
thumbnail Fig. 8.

Data processing summary.

In the text
thumbnail Fig. 9.

Top panel: signal-to-noise distribution as a function of the integration time for different observations of J0437−4715 before and after the upgrade in A1. Bottom panel: typical pulse profiles in the pre-upgrade and post-upgrade configuration.

In the text
thumbnail Fig. 10.

Residuals of the ToAs of the pulsar J0437−4715 measured at IAR with A1 and A2. The residuals and rms have been computed with the TEMPO2 package (Hobbs et al. 2006) and error bars correspond only to template-fitting errors (i.e., no systematics considered).

In the text
thumbnail Fig. 11.

Pulse profiles of the magnetar J1810−197 at 1400 MHz measured by IAR antennas at different epochs.

In the text
thumbnail Fig. 12.

Vela’s glitch on February 1st. 2019 as measured at IAR in the 1400 MHz band. Ephemeris for the pre-glitch epoch were taken from Sarkissian et al. (2017) and the glitch date from Kerr (2019). The errorbars at one σ level are taken directly from the output of PRESTO and are likely to represent an underestimation. The dispersion is larger for short observations (10–15 min) than for long ones (> 60 min). We highlight that at MJD 58487 the data points from A1 and A2 overlap and are hardly distinguishable.

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.