• Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.
  • Timing results for the black-widow pulsar J2051-0827 are presented, using a 21-year dataset from four European Pulsar Timing Array telescopes and the Parkes radio telescope. This dataset, which is the longest published to date for a black-widow system, allows for an improved analysis that addresses previously unknown biases. While secular variations, as identified in previous analyses, are recovered, short-term variations are detected for the first time. Concurrently, a significant decrease of approx. 2.5x10-3 cm-3 pc in the dispersion measure associated with PSR J2051-0827 is measured for the first time and improvements are also made to estimates of the proper motion. Finally, PSR J2051-0827 is shown to have entered a relatively stable state suggesting the possibility of its eventual inclusion in pulsar timing arrays.
  • We present evidence for a small glitch in the spin evolution of the millisecond pulsar J0613$-$0200, using the EPTA Data Release 1.0, combined with Jodrell Bank analogue filterbank TOAs recorded with the Lovell telescope and Effelsberg Pulsar Observing System TOAs. A spin frequency step of 0.82(3) nHz and frequency derivative step of ${-1.6(39) \times 10^{-19}\,\text{Hz} \ \text{s}^{-1}}$ are measured at the epoch of MJD 50888(30). After PSR B1821$-$24A, this is only the second glitch ever observed in a millisecond pulsar, with a fractional size in frequency of ${\Delta \nu/\nu=2.5(1) \times 10^{-12}}$, which is several times smaller than the previous smallest glitch. PSR J0613$-$0200 is used in gravitational wave searches with pulsar timing arrays, and is to date only the second such pulsar to have experienced a glitch in a combined 886 pulsar-years of observations. We find that accurately modelling the glitch does not impact the timing precision for pulsar timing array applications. We estimate that for the current set of millisecond pulsars included in the International Pulsar Timing Array, there is a probability of $\sim 50$% that another glitch will be observed in a timing array pulsar within 10 years.
  • We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.
  • We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.
  • In this paper we present results from flux density measurements for 21 pulsars over a wide frequency range, using the Giant Metrewave Radio Telescope (GMRT) and the Effelsberg telescope. Our sample was a set of mostly newly discovered pulsars from the selection of candidates for gigahertz-peaked spectra (GPS) pulsars. Using the results of our observations along with previously published data, we identify two new GPS pulsars. One of them, PSR J1740+1000, with dispersion measure of 24 pc cm$^{-3}$, is the first GPS pulsar with such a low DM value.We also selected several strong candidates for objects with high frequency turnover in their spectra which require further investigation.We also revisit our source selection criteria for future searches for GPS pulsars.
  • We present new results of fitting 108 spectra of radio pulsars with the flicker noise model proposed by Loehmer et al. (2008) and compare them with the spectral indices of power-law fits published by Maron et al. (2000). The fits to the model were carried out using the Markov chain Monte Carlo (MCMC) method appropriate for the non-linear fits. Our main conclusion is that pulsar radio spectra can be statistically very well described by the flicker noise model over wide frequency range from a few tens of MHz up to tens of GHz. Moreover, our dataset allows us to conduct statistical analysis of the model parameters. As our results show, there is a strong negative correlation between the flicker noise spectrum model parameters log $S_0$ and $n$ and a strong positive relationship between n and the power-law spectral index $\alpha$. The latter implies that their physical meaning is similar, however the flicker noise model has an advantage over broken power-law model. Not only it describes the spectra in higher frequency range with only two parameters, not counting scaling factor $S_0$, but also it shows smooth transition from flat to steep behaviour at lower and higher frequencies, respectively. On the other hand, there are no correlations of the flicker noise model parameters $S_0$, $\tau$ and $n$ with any of pulsar physical properties.
  • We show the results of microsecond resolution radio data analysis focused on flux measurements of single pulses of PSR B1133+16. The data were recorded at 4.85 GHz and 8.35 GHz with 0.5 GHz and 1.1 GHz bandwidth, respectively, using Radio Telescope Effelsberg (MPIfR). The most important conclusion of the analysis is, that the strongest single pulse emission at 4.85 GHz and 8.35 GHz contributes almost exclusively to the trailing part of the leading component of the pulsar mean profile, whereas studies at lower frequencies report that the contribution is spread almost uniformly covering all phases of the pulsar mean profile. We also estimate the radio emission heights to be around 1%--2% of the light cylinder radius which is in agreement with previous studies. Additionally these observations allowed us to add two more measurements of the flux density to the PSR B1133+16 broadband radio spectrum covering frequencies from 16.7 MHz up to 32 GHz. We fit two different models to the spectrum: the broken power law and the spectrum based on flicker noise model, which represents the spectrum in a simpler but similarly accurate way.
  • Aims. To investigate the flux density modulation from pulsars and the existence of specific behaviour of modulation index versus frequency. Methods. Several pulsars have been observed with the Effelsberg radio telescope at 8.35 GHz. Their flux density time series have been corrected for interstellar scintillation effects. Results. We present the measurement of modulation indices for 8 pulsars. We confirm the presence of a critical frequency at ~1 GHz for these pulsars (including 3 new ones from this study). We derived intrinsic modulation indices for the resulting flux density time series. Our data analysis revealed strong single pulses detected from 5 pulsars.
  • Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.
  • The mode switching phenomenon of PSR B0329+54 is investigated based on the long-term monitoring from September 2003 to April 2009 made with the Urumqi 25m radio telescope at 1540 MHz. At that frequency, the change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratios between the leading and trailing components are measured for the individual profiles averaged over a few minutes. It is found that the ratios follow normal distributions, where the abnormal mode has a wider typical width than the normal mode, indicating that the abnormal mode is less stable than the normal mode. Our data show that 84.9% of the time for PSR B0329+54 was in the normal mode and 15.1% was in the abnormal mode. From the two passages of eight-day quasi-continuous observations in 2004, and supplemented by the daily data observed with 15 m telescope at 610 MHz at Jodrell Bank Observatory, the intrinsic distributions of mode timescales are constrained with the Bayesian inference method. It is found that the gamma distribution with the shape parameter slightly smaller than 1 is favored over the normal, lognormal and Pareto distributions. The optimal scale parameters of the gamma distribution is 31.5 minutes for the abnormal mode and 154 minutes for the normal mode. The shape parameters have very similar values, i.e. 0.75^{+0.22}_{-0.17} for the normal mode and 0.84^{+0.28}_{-0.22} for the abnormal mode, indicating the physical mechanisms in both modes may be the same. No long-term modulation of the relative intensity ratios was found for both the modes, suggesting that the mode switching was stable. The intrinsic timescale distributions, for the first time constrained for this pulsar, provide valuable information to understand the physics of mode switching.
  • We investigate a high frequency turn-over effect in radio spectra for pulsars with positive or flat spectral index. Using GMRT and Effelsberg observatory, we estimated the flux density to reconstruct pulsar spectra. We find objects that have a maximum flux in their spectrum above 1 GHz and whose spectral indices are positive at lower frequencies. Some pulsars with a turn-over in their spectrum at high frequencies are found to exist in very interesting environments. We call these objects gigahertz-peaked spectra pulsars.
  • Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^{\alpha}$, where $f$ is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude $A$ as a function of the unknown spectral slope $\alpha$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $\alpha=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.
  • We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, {\mu}_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.
  • The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.
  • We present results from the high precision timing analysis of the pulsar-white dwarf (WD) binary PSR J1012+5307 using 15 years of multi-telescope data. Observations were performed regularly by the European Pulsar Timing Array (EPTA) network, consisting of Effelsberg, Jodrell Bank, Westerbork and Nan\c{c}ay. All the timing parameters have been improved from the previously published values, most by an order of magnitude. In addition, a parallax measurement of $\pi = 1.2(3)$ mas is obtained for the first time for PSR J1012+5307, being consistent with the optical estimation from the WD companion. Combining improved 3D velocity information and models for the Galactic potential the complete evolutionary Galactic path of the system is obtained. A new intrinsic eccentricity upper limit of $e<8.4\times 10^{-7}$ is acquired, one of the smallest calculated for a binary system and a measurement of the variation of the projected semi-major axis also constrains the system's orbital orientation for the first time. It is shown that PSR J1012+5307 is an ideal laboratory for testing alternative theories of gravity. The measurement of the change of the orbital period of the system of $\dot{P}_{b} = 5(1)\times 10^{-14}$ is used to set an upper limit on the dipole gravitational wave emission that is valid for a wide class of alternative theories of gravity. Moreover, it is shown that in combination with other binary pulsars PSR J1012+5307 is an ideal system to provide self-consistent, generic limits, based only on millisecond pulsar data, for the dipole radiation and the variation of the gravitational constant $\dot{G}$.
  • We have used the 76-m Lovell, 94-m equivalent WSRT and 100-m Effelsberg radio telescopes to investigate the simultaneous single-pulse properties of the radio emitting magnetar AXP XTE J1810-197 at frequencies of 1.4, 4.8 and 8.35 GHz during May and July 2006. We study the magnetar's pulse-energy distributions which are found to be very peculiar as they are changing on time-scales of days and cannot be fit by a single statistical model. The magnetar exhibits strong spiky single giant-pulse-like subpulses, but they do not fit the definition of the giant pulse or giant micropulse phenomena. Measurements of the longitude-resolved modulation index reveal a high degree of intensity fluctuations on day-to-day time-scales and dramatic changes across pulse phase. We find the frequency evolution of the modulation index values differs significantly from what is observed in normal radio pulsars. We find that no regular drifting subpulse phenomenon is present at any of the observed frequencies at any observing epoch. However, we find a quasi-periodicity of the subpulses present in the majority of the observing sessions. A correlation analysis indicates a relationship between components from different frequencies. We discuss the results of our analysis in light of the emission properties of normal radio pulsars and a recently proposed model which takes radio emission from magnetars into consideration.
  • As part of a European Pulsar Network (EPN) multi-telescope observing campaign, we performed simultaneous multi-frequency observations at 1.4, 4.9 and 8.4 GHz during July 2006 and quasi-simultaneous multi-frequency observations from Decem- ber 2006 until July 2007 at 2.7, 4.9, 8.4, 14.6 and 32 GHz, in order to obtain flux density measurements and spectral features of the 5.5-sec radio-emitting magnetar AXP J1810-197. We monitored the spectral evolution of its pulse shape which consists of a main pulse (MP) and an interpulse (IP). We present the flux density spectrum of the average profile and of the separate pulse components of this first-known radio-emitting transient anomalous X-ray pulsar. We observe a decrease of the flux density by a factor of 10 within 8 months and follow the disappearance of one of the two main components. Although the spectrum is generally flat, we observe large fluctuations of the spectral index with time. For that reason we have made some measurements of modulation indices for individual pulses in order to also investigate the origin of these fluctuations.
  • PSR B0031-07 is well known to exhibit three different modes of drifting sub-pulses (mode A, B and C). It has recently been shown that in a multifrequency observation, consisting of 2700 pulses, all driftmodes were visible at low frequencies, while at 4.85 GHz only mode-A drift or non-drifting emission was detected. This suggests that modes A and B are emitted in sub-beams, rotating at a fixed distance from the magnetic axis, with the mode-B sub-beams being closer to the magnetic axis than the mode-A sub-beams. Diffuse emission between the sub-beams can account for the non-drifting emission. Using the results of an analysis of simultaneous multifrequency observations of PSR B0031-07, we set out to construct a geometrical model that includes emission from both sub-beams and diffuse emission and describes the regions of the radio emission of PSR B0031-07 at each emission frequency for driftmodes A and B. Based on the vertical spacing between driftbands, we have determined the driftmode of each sequence of drift. To restrict the model, we calculated average polarisation and intensity characteristics for each driftmode and at each frequency. The model reproduces the observed polarisation and intensity characteristics, suggesting that diffuse emission plays an important role in the emission properties of PSR B0031-07. The model further suggests that the emission heights of this pulsar range from a few kilometers to a little over 10 kilometers above the pulsar surface. We also find that the relationships between height and frequency of emission that follow from curvature radiation and from plasma-frequency emission could not be used to reproduce the observed frequency dependence of the width of the average intensity profiles.
  • Pulse-to-pulse intensity variations are a common property of pulsar radio emission. For some of the objects single pulses are often 10-times stronger than their average pulse. The most dramatic events are so-called giant radio pulses (GRPs). They can be thousand times stronger than the regular single pulses from the pulsar. Giant pulses are a rare phenomenon, occurring in very few pulsars which split into two groups. The first group contains very young and energetic pulsars like the Crab pulsar, and its twin (PSR B0540-69) in the Large Magellanic Cloud (LMC), while the second group is represented by old, recycled millisecond pulsars like PSR B1937+21, PSR B1821-24, PSR B1957+20 and PSR J0218+4232 (the only millisecond pulsar detected in gamma-rays). We compare the characteristics of GRPs for these two pulsar groups. Moreover, our latest findings of new features in the Crab GRPs are presented. Analysis of our Effelsberg data at 8.35 GHz shows that GRPs do occur in all phases of its ordinary radio emission, including the phases of the two high frequency components (HFCs) visible only between 5 and 9 GHz.
  • Here we present the results from an analysis of a multifrequency simultaneous observation of PSR B0031$-$07. We have constructed a geometrical model, based on an empirical relationship between height and frequency of emission, that reproduces many of the observed characteristics. The model suggests very low emission altitudes for this pulsar of only a few kilometers above the star's surface.
  • Topics of the Colloquium: a) Interplanetary scintillation b) Interstellar scintillation c) Modeling and physical origin of the interplanetary and the interstellar plasma turbulence d) Scintillation as a tool for investigation of radio sources e) Seeing through interplanetary and interstellar turbulent media Ppt-presentations are available on the Web-site: http://www.prao.ru/conf/Colloquium/main.html
  • Individual giant radio pulses (GRPs) from the Crab pulsar last only a few microseconds. However, during that time they rank among the brightest objects in the radio sky reaching peak flux densities of up to 1500 Jy even at high radio frequencies. Our observations show that GRPs can be found in all phases of ordinary radio emission including the two high frequency components (HFCs) visible only between 5 and 9 GHz (Moffett & Hankins, 1996). This leads us to believe that there is no difference in the emission mechanism of the main pulse (MP), inter pulse (IP) and HFCs. High resolution dynamic spectra from our recent observations of giant pulses with the Effelsberg telescope at a center frequency of 8.35 GHz show distinct spectral maxima within our observational bandwidth of 500 MHz for individual pulses. Their narrow band components appear to be brighter at higher frequencies (8.6 GHz) than at lower ones (8.1 GHz). Moreover, there is an evidence for spectral evolution within and between those structures. High frequency features occur earlier than low frequency ones. Strong plasma turbulence might be a feasible mechanism for the creation of the high energy densities of ~6.7 x 10^4 erg cm^-3 and brightness temperatures of 10^31 K.
  • Recent work has made it clear that the ``standard model'' of pulsar radio emission cannot be the full answer. Some fundamental assumptions about the magnetic field and plasma flow in the radio-loud region have been called into question by recent observational and theoretical work, but the solutions to the problems posed are far from clear. It is time to formulate and carry out new observational campaigns designed to address these problems; sensitive low-frequency observations will an important part of such a campaign. Because pulsars are strong at low frequencies, we believe there will be a good number of candidates even for high-time-resolution single pulse work, as well as mean profile and integrated spectrum measurements. Such data can push the envelope of current models, test competing theories of the radio loud region, and possibly provide direct measures of the state of the emitting plasma.
  • Rough estimates and upper limits on current and particle densities form the basis of most of the canonical pulsar models. Whereas the surface of the rotating neutron star is capable of supplying sufficient charges to provide a current that, given the polar cap potential, could easily fuel the observed energy loss processes, observational and theoretical constraints provide strict upper limits to the charge densities. The space charge of a current consisting solely of particles having only one sign creates a compensating potential that will make the maximum current dependent on potential and distance. In the non-relativistic case this fact is expressed in the familiar Child-Langmuir law. Its relativistic generalization and subsequent application to the inner pulsar magnetosphere provides clear limits on the strength and radial extension of charged currents originating on the polar cap. Violent Pierce-type oscillations set in, if one attempts to inject more current than the space charge limit into a given volume. These considerations apply wherever there is a significant amount of charged current flow, in particular in the gap regions. There they can be used to derive limits on the size of such gaps and their stability.