Search

The Parkes 20 cm Multibeam pulsar surveys have discovered nearly half of the known pulsars and revealed many distant pulsars with high dispersion measures. Using a sample of 1,301 pulsars from these surveys, we have explored the spatial distribution and birth rate of normal pulsars. The pulsar distances used to calculate the pulsar surface density are estimated from the YMW16 electron-density model. When estimating the impact of the Galactic background radiation on our survey, we projected pulsars in the Galaxy onto the Galactic plane, assuming that the flux density distribution of pulsars is uniform in all directions, and utilized the most up-to-date background temperature map. We also used an up-to-date version of the ATNF Pulsar Catalogue to model the distribution of pulsar flux densities at 1400 MHz. We derive an improved radial distribution for the pulsar surface density projected on to the Galactic plane, which has a maximum value at $\sim$4 kpc from the Galactic Centre. We also derive the local surface density and birthrate of pulsars, obtaining 47 $\pm$ 5 $\mathrm{kpc^{-2}}$ and $\sim$ 4.7 $\pm$ 0.5 $\mathrm{kpc^{-2}\ Myr^{-1}}$, respectively. For the total number of potentially detectable pulsars in the Galaxy, we obtain (1.1 $\pm$ 0.2) $\times$ $10^{4}$ and (1.1 $\pm$ 0.2) $\times$ $10^{5}$ before and after applying the TM98 beaming correction model. The radial distribution function is used to estimate the proportion of pulsars in each spiral arm and the Galactic centre.

Understanding the jitter noise resulting from single-pulse phase and shape variations is important for the detection of gravitational waves using pulsar timing array. We presented measurements of jitter noise and single-pulse variability of 12 millisecond pulsars that are part of the International Pulsar Timing Array sample using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We found that the levels of jitter noise can vary dramatically among pulsars. A moderate correlation with a correlation coefficient of 0.57 between jitter noise and pulse width is detected. To mitigate jitter noise, we performed matrix template matching using all four Stokes parameters. Our results revealed a reduction in jitter noise ranging from 6.7\% to 39.6\%. By performing longitude-resolved fluctuation spectrum analysis, we identified periodic intensity modulations in 10 pulsars. In PSR J0030+0451, we detected single-pulses with energies more than 10 times the average pulse energy, suggesting the presence of giant pulses. We also observed a periodic mode-changing phenomenon in PSR J0030+0451. We examined the achievable timing precision by selecting a sub-set of pulses with a specific range of peak intensity, but no significant improvement in timing precision is achievable., Comment: 14 pages, 9 figures, Accepted for publication in ApJ

We present the discovery and timing results of four pulsars discovered in a pilot survey at intermediate Galactic latitudes with the Five-hundred Aperture Spherical Telescope (FAST). Among these pulsars, two belong to the category of millisecond pulsars (MSPs) with spin periods of less than 20 ms. The other two fall under the classification of "mildly recycled" pulsars, with massive white dwarfs as companions. Remarkably, this small survey, covering an area of 4.7 $deg^2$ , led to the discovery of four recycled pulsars. Such success underscores the immense potential of future surveys at intermediate Galactic latitudes. In order to assess the potential yield of MSPs, we conducted population simulations and found that both FAST and Parkes new phased array feed surveys, focusing on intermediate Galactic latitudes, have the capacity to uncover several hundred new MSPs., Comment: 7 pages, 4 figures, 2 tables, accepted to ApJ

Black widows are millisecond pulsars ablating their companions. The material blown from the companion blocks the radio emission, resulting in radio eclipses. The properties of the eclipse medium are poorly understood. Here, we present direct evidence of the existence of magnetic fields in the eclipse medium of the black widow PSR J2051$-$0827 using observations made with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We detect a regular decrease in rotation measure (RM) in the egress of eclipse, changing from $60\,\rm rad\,m^{-2}$ to $-28.7\,\rm rad\,m^{-2}$. The RM gradually changes back to normal when the line-of-sight moves away from the eclipse. The estimated line-of-sight magnetic field strength in the eclipse medium is $\sim 0.1$ G. The RM reversal could be caused by a change of the magnetic field strength along the line of sight due to binary orbital motion. The RM reversal phenomenon has also been observed in some repeating fast radio bursts (FRBs), and the study of spider pulsars may provide additional information about the origin of FRBs., Comment: 7 pages, 3 figures, accept for publication in ApJ

Flux densities are basic observation parameters to describe pulsars. In the most updated pulsar catalog, 24% of the listed radio pulsars have no flux density measurement at any frequency. Here, we report the first flux density measurements, spectral indices, pulse profiles, and correlations of the spectral index with pulsar parameters for 19 pulsars employing the Ultra-Wideband Low (UWL) receiver system installed on the Parkes radio telescope. The results for spectral indices of 17 pulsars are in the range between -0.6 and -3.10. The polarization profiles of thirteen pulsars are shown. There is a moderate correlation between the spectral index and spin frequency. For most pulsars detected, the S/N ratio of pulse profile is not high, so DM, Faraday rotation measure (RM), and polarization can not be determined precisely. Twenty-nine pulsars were not detected in our observations. We discuss the possible explanations for why these pulsars were not detected., Comment: 19 pages, 6 figures, 5 tables, accepted for publication in RAA

We report the phase-connected timing ephemeris, polarization pulse profiles, Faraday rotation measurements, and Rotating-Vector-Model (RVM) fitting results of twelve millisecond pulsars (MSPs) discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal radio Astronomy FAST survey (CRAFTS). The timing campaigns were carried out with FAST and Arecibo over three years. Eleven of the twelve pulsars are in neutron star - white dwarf binary systems, with orbital periods between 2.4 and 100 d. Ten of them have spin periods, companion masses, and orbital eccentricities that are consistent with the theoretical expectations for MSP - Helium white dwarf (He WD) systems. The last binary pulsar (PSR J1912$-$0952) has a significantly smaller spin frequency and a smaller companion mass, the latter could be caused by a low orbital inclination for the system. Its orbital period of 29 days is well within the range of orbital periods where some MSP - He WD systems have shown anomalous eccentricities, however, the eccentricity of PSR J1912$-$0952 is typical of what one finds for the remaining MSP - He WD systems., Comment: 11 pages, 5 figures, MNRAS accepted

We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $\tau\propto\Delta{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering due to multi-path propagation effects in the eclipsing material. We show that Kolmogorov turbulence in the eclipsing material with an inner scale of $\sim100$ m and an outer scale of the size of the eclipse region can naturally explain the observation. Our results show that the eclipsing material in such systems can be highly turbulent and suggest that scattering is one of the main eclipsing mechanisms at around 1.4 GHz., Comment: 8 pages, 4 figures, MNRAS accepted

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^{\alpha}$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $\alpha = -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $\alpha = -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability., Comment: 17 pages, 12 figures, accepted in MNRAS

We report on an {unusually} bright observation of PSR J1720$-$0533 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulsar is in a black widow system that {was discovered by the Commensal Radio Astronomy FAST Survey (CRAFTS). By coincidence, a bright scintillation maximum was simultaneous with the eclipse in our observation which allowed for precise measurements of flux density variations, as well as dispersion measure (DM) and polarization.} We found that there are quasi-periodic pulse emission variations with a modulation period of $\sim$ {22\,s} during the ingress of the eclipse, which could be caused by plasma lensing. {No such periodic modulation was found during the egress of the eclipse. } {The linear polarization of the pulsar disappears before the eclipse, even before there is a visually obvious change in DM. We also found that the pulse scattering maybe play an important role in the eclipse of PSR J1720$-$0533.}, Comment: 7 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

We present a novel methodological framework for quantum spatial search, generalising the Childs & Goldstone ($\mathcal{CG}$) algorithm via alternating applications of marked-vertex phase shifts and continuous-time quantum walks. We determine closed form expressions for the optimal walk time and phase shift parameters for periodic graphs. These parameters are chosen to rotate the system about subsets of the graph Laplacian eigenstates, amplifying the probability of measuring the marked vertex. The state evolution is asymptotically optimal for any class of periodic graphs having a fixed number of unique eigenvalues. We demonstrate the effectiveness of the algorithm by applying it to obtain $\mathcal{O}(\sqrt{N})$ search on a variety of graphs. One important class is the $n \times n^3$ rook graph, which has $N=n^4$ vertices. On this class of graphs the $\mathcal{CG}$ algorithm performs suboptimally, achieving only $\mathcal{O}(N^{-1/8})$ overlap after time $\mathcal{O}(N^{5/8})$. Using the new alternating phase-walk framework, we show that $\mathcal{O}(1)$ overlap is obtained in $\mathcal{O}(\sqrt{N})$ phase-walk iterations.

A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5-year data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models $\unicode{x2014}$ one with a common-spectrum process and one without $\unicode{x2014}$ we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings-Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be $A = 2.2^{+0.4}_{-0.3} \times 10^{-15}$ in units of gravitational-wave strain at a frequency of $1\, \text{yr}^{-1}$. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background., Comment: 10 pages, 5 figures

The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017$-$7156, J1125$-$6014, J1545$-$4550, and J1732$-$5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125$-$6014 implies a neutron star mass of $M_p = 1.5 \pm 0.2 M_\odot$ (68% credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600$-$3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M$_\odot$ (68% credibility interval). We measure the distance to PSR J1909$-$3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments., Comment: 44 pages, 32 figures, 7 tables. Accepted for publication in MNRAS; v2 includes corrected typos in parameter tables and reference added

We present radio observation of a millisecond pulsar PSR J0621+1002 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulsar shows periodic pulse intensity modulations for both the first and the third pulse components. The fluctuation spectrum of the first pulse component has one peak of 3.0$\pm$0.1 pulse periods, while that of the third pulse component has two diffused peaks of 3.0$\pm$0.1 and 200$\pm$1 pulse periods. The single pulse timing analysis is carried out for this pulsar and the single pulses can be divided into two classes based on the post-fit timing residuals. We examined the achievable timing precision using only the pulses in one class or bright pulses. However, the timing precision improvement is not achievable., Comment: accepted for publication in ApJ. 7 pages, 8 figures

This paper examines the performance of spatial search where the Grover diffusion operator is replaced by continuous-time quantum walks on a class of interdependent networks. We prove that for a set of optimal quantum walk times and marked vertex phase shifts, a deterministic algorithm for structured spatial search is established that finds the marked vertex with 100% probability. This improves on the original spatial search algorithm on the same class of graphs, which we show can only amplify to 50% probability. Our method uses $\left\lceil\frac{\pi}{2\sqrt{2}}\sqrt{N}\right\rceil$ marked vertex phase shifts for an $N$-vertex graph, making it comparable with Grover's algorithm for unstructured search. It is expected that this new framework can be readily extended to deterministic spatial search on other families of graph structures., Comment: 10 pages, 8 figures

Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no convincing detections we place upper limits on the GW amplitude for different event durations. From these bounds we place limits on the cosmic string tension of G mu ~ 10^{-5}, and highlight that this bound is independent from those obtained using other techniques. We discuss the physical implications of our results and the prospect of probing cosmic strings in the era of Square Kilometre Array (SKA)., Comment: Accepted for publication by MNRAS

This paper proposes a highly efficient quantum algorithm for portfolio optimisation targeted at near-term noisy intermediate-scale quantum computers. Recent work by Hodson et al. (2019) explored potential application of hybrid quantum-classical algorithms to the problem of financial portfolio rebalancing. In particular, they deal with the portfolio optimisation problem using the Quantum Approximate Optimisation Algorithm and the Quantum Alternating Operator Ansatz. In this paper, we demonstrate substantially better performance using a newly developed Quantum Walk Optimisation Algorithm in finding high-quality solutions to the portfolio optimisation problem.

Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes because incorrect models diminish sensitivity and may cause a spurious gravitational wave detection. Here we characterise noise processes for the 26 pulsars in the second data release of the Parkes Pulsar Timing Array using Bayesian inference. In addition to well-studied noise sources found previously in pulsar timing array data sets such as achromatic timing noise and dispersion measure variations, we identify new noise sources including time-correlated chromatic noise that we attribute to variations in pulse scattering. We also identify "exponential dip" events in four pulsars, which we attribute to magnetospheric effects as evidenced by pulse profile shape changes observed for three of the pulsars. This includes an event in PSR J1713$+$0747, which had previously been attributed to interstellar propagation. We present noise models to be used in searches for gravitational waves. We outline a robust methodology to evaluate the performance of noise models and identify unknown signals in the data. The detection of variations in pulse profiles highlights the need to develop efficient profile domain timing methods., Comment: 18 pages, 7 figures, 5 tables

We have observed PSR~B1534+12 (J1537+1155), a pulsar with a neutron star companion, using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We found that this pulsar shows two distinct emission states: a weak state with a wide pulse profile and a burst state with a narrow pulse profile. The weak state is always present. We cannot, with our current data, determine whether the pulse energy of the weak state follows a normal or a log-normal distribution. The burst state energy distribution follows a power-law. The amplitude of the single pulse emission in the burst state varies significantly; the peak flux intensity of the brightest pulse is 334 times stronger than that of the average pulse. We also examined the timing precision achievable using only bright pulses, which showed no demonstrable improvement because of pulse jitter and therefore quantified the jitter noise level for this pulsar., Comment: 6 pages, 6 figures, accepted for publication in The Astrophysical Journal Letters

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS)., Comment: 18 pages, 8 figures, 5 tables. Accepted for publication by MNRAS. Correcting error in the caption of Figure 8

We present long-term observations of two intermittent pulsars, PSRs~J1832+0029 and J1841$-$0500 using the Parkes 64\,m radio telescope. The radio emission for these pulsars switches "off" for year-long durations. Our new observations have enabled us to improve the determination of the on-off timescales and the spin down rates during those emission states. In general our results agree with previous studies of these pulsars, but we now have significantly longer data spans. We have identified two unexpected signatures in the data. Weak emission was detected in a single observation of PSR~J1832$+$0029 during an "off" emission state. For PSR~J1841$-$0500, we identified a quasi-periodic fluctuation in the intensities of the detectable single pulses, with a modulation period between 21 and 36 pulse periods., Comment: 7 pages, 7 figures, accepted for publication in ApJ

We have used the Nanshan 25-m Radio Telescope at Xinjiang Astronomical Observatory to obtain timing observations of 87 pulsars from 2002 July to 2014 March. Using the "Cholesky" timing analysis method we have determined positions and proper motions for 48 pulsars, 24 of which are improved positions compared to previously published values. We also present the first published proper motions for nine pulsars and improved proper motions for 21 pulsars using pulsar timing and position comparison method. The pulsar rotation parameters are derived and are more accurate than previously published values for 36 pulsars. Glitches are detected in three pulsars: PSRs J1722$-$3632, J1852$-$0635 and J1957+2831. For the first two, the glitches are large, with $\Delta\nu_g/\nu > 10^{-6}$, and they are the first detected glitches in these pulsars. PSR J1722$-$3632 is the second oldest pulsar with large glitch. For the middle-age pulsars ($\tau_c > 10^5$~yr), the calculated braking indices, $|n|$, are strongly correlated with $\tau_c$ and the numbers of positive and negative values of $n$ are almost equal. For young pulsars ($\tau_c < 10^5$~yr), there is no correlation between $|n|$ and $\tau_c$ and most have $n>0$., Comment: 21 pages, 8 figures, accepted for publication by ApJ

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches., Comment: 25 pages, 8 figures

Quantum walks (QW) are of crucial importance in the development of quantum information processing algorithms. Recently, several quantum algorithms have been proposed to implement network analysis, in particular to rank the centrality of nodes in networks represented by graphs. Employing QW in centrality ranking is advantageous comparing to certain widely used classical algorithms (e.g. PageRank) because QW approach can lift the vertex rank degeneracy in certain graphs. However, it is challenging to implement a directed graph via QW, since it corresponds to a non-Hermitian Hamiltonian and thus cannot be accomplished by conventional QW. Here we report the realizations of centrality rankings of both a three-vertex and four-vertex directed graphs with parity-time (PT) symmetric quantum walks. To achieve this, we use high-dimensional photonic quantum states, optical circuitries consisting of multiple concatenated interferometers and dimension dependent loss. Importantly, we demonstrate the advantage of QW approach experimentally by breaking the vertex rank degeneracy in a four-vertex graph. Our work shows that PT-symmetric quantum walks may be useful for realizing advanced algorithm in a quantum network.

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves., Comment: accepted by MNRAS

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science., Comment: Submitted to MNRAS and in review, 23 pages, 5 figures

The Parkes Pulsar Timing Array (PPTA) project monitors two dozen millisecond pulsars (MSPs) in order to undertake a variety of fundamental physics experiments using the Parkes 64m radio telescope. Since June 2017 we have been undertaking commensal searches for fast radio bursts (FRBs) during the MSP observations. Here, we report the discovery of four FRBs (171209, 180309, 180311 and 180714). The detected events include an FRB with the highest signal-to-noise ratio ever detected at the Parkes observatory, which exhibits unusual spectral properties. All four FRBs are highly polarized. We discuss the future of commensal searches for FRBs at Parkes., Comment: 8 pages, accepted for publication in MNRAS

In this paper we show how using complex valued edge weights in a graph can completely suppress the flow of probability amplitude in a continuous time quantum walk to specific vertices of the graph when the edge weights, graph topology and initial state of the quantum walk satisfy certain conditions. The conditions presented in this paper are derived from the so-called chiral quantum walk, a variant of the continuous time quantum walk which incorporates directional bias with respect to site transfer probabilities between vertices of a graph by using complex edge weights. We examine the necessity to break the time reversal symmetry in order to achieve zero transfer in continuous time quantum walks. We also consider the effect of decoherence on zero transfer and suggest that this phenomena may be used to detect decoherence in the system.

A quantum walk whose continuous limit coincides with Dirac equation is usually called a Dirac Quantum Walk (DQW). A new systematic method to build DQWs coupled to electromagnetic (EM) fields is introduced and put to test on several examples of increasing difficulty. It is first used to derive the EM coupling of a well-known $3D$ walk on the cubic lattice. Recently introduced DQWs on the triangular and honeycomb lattice are then re-derived, showing for the first time that these are the only DQWs that can be defined with spinors living on the vertices of these lattices. As a third example of the method's effectiveness, a new $3D$ walk on a parallelepiped lattice is derived. As a fourth, negative example, it is shown that certain lattices like the rhombohedral lattice cannot be used to build DQWs. The effect of changing representation in the Dirac equation is also discussed.

Timing analysis of PSR J1705$-$1906 using data from Nanshan 25-m and Parkes 64-m radio telescopes, which span over fourteen years, shows that the pulsar exhibits significant proper motion, and rotation instability. We updated the astrometry parameters and the spin parameters of the pulsar. In order to minimize the effect of timing irregularities on measuring its position, we employ the Cholesky method to analyse the timing noise. We obtain the proper motion of $-$77(3) \,mas\,yr$^{-1}$ in right ascension and $-$38(29) \,mas\,yr$^{-1}$ in declination. The power spectrum of timing noise is analyzed for the first time, which gives the spectral exponent $\alpha=-5.2$ for the power-law model indicating that the fluctuations in spin frequency and spin-down rate dominate the red noise. We detect two small glitches from this pulsar with fractional jump in spin frequency of $\Delta \nu/\nu\sim2.9\times10^{-10}$ around MJD~55199 and $\Delta \nu/\nu\sim2.7\times10^{-10}$ around MJD~55953. Investigations of pulse profile at different time segments suggest no significant changes in the pulse profiles around the two glitches., Comment: 10 pages, 8 figures, 7 tables, Astrophysics and Space Science, Volume 363, Issue 5, article id. 96

In this paper we present an efficiently scaling quantum algorithm which finds the size of the maximum common edge subgraph for a pair of arbitrary graphs and thus provides a meaningful measure of graph similarity. The algorithm makes use of a two-part quantum dynamic process: in the first part we obtain information crucial for the comparison of two graphs through linear quantum computation. However, this information is hidden in the quantum system with vanishingly small amplitude that even quantum algorithms such as Grover's search are not fast enough to distill the information efficiently. In order to extract the information we call upon techniques in nonlinear quantum computing to provide the speed-up necessary for an efficient algorithm. The linear quantum circuit requires $\mathcal{O}(n^3 \log^3 (n) \log \log (n))$ elementary quantum gates and the nonlinear evolution under the Gross-Pitaevskii equation has a time scaling of $\mathcal{O}(\frac{1}{g} n^2 \log^3 (n) \log \log (n))$, where $n$ is the number of vertices in each graph and $g$ is the strength of the Gross-Pitaveskii non-linearity. Through this example, we demonstrate the power of nonlinear quantum search techniques to solve a subset of NP-hard problems.

Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves., Comment: Accepted for publication by the Monthly Notices of the Royal Astronomical Society

We demonstrate coupling between the atomic spin and orbital-angular-momentum (OAM) of the atom's center-of-mass motion in a Bose-Einstein condensate (BEC). The coupling is induced by Raman-dressing lasers with a Laguerre-Gaussian beam, and creates coreless vortices in a $F=1$ $^{87}$Rb spinor BEC. We observe correlations between spin and OAM in the dressed state and characterize the spin texture; the result is in good agreement with the theory. In the presence of the Raman field our dressed state is stable for 0.1~s or longer, and it decays due to collision-induced relaxation. As we turn off the Raman beams, the vortex cores in the bare spin $|m_F=1\rangle$ and $|-1\rangle$ split. These spin-OAM coupled systems with the Raman-dressing approach have great potential for exploring new topological textures and quantum states.

Pulsar positions can be measured with high precision using both pulsar timing methods and very-long-baseline interferometry (VLBI). Pulsar timing positions are referenced to a solar-system ephemeris, whereas VLBI positions are referenced to distant quasars. Here we compare pulsar positions from published VLBI measurements with those obtained from pulsar timing data from the Nanshan and Parkes radio telescopes in order to relate the two reference frames. We find that the timing positions differ significantly from the VLBI positions (and also differ between different ephemerides). A statistically significant change in the obliquity of the ecliptic of $2.16\pm0.33$\,mas is found for the JPL ephemeris DE405, but no significant rotation is found in subsequent JPL ephemerides. The accuracy with which we can relate the two frames is limited by the current uncertainties in the VLBI reference source positions and in matching the pulsars to their reference source. Not only do the timing positions depend on the ephemeris used in computing them, but also different segments of the timing data lead to varying position estimates. These variations are mostly common to all ephemerides, but slight changes are seen at the 10$\mu$as level between ephemerides., Comment: accepted by MNRAS,22 pages, 5 figures

Pulse arrival time measurements over fourteen years with the Nanshan 25-m and Parkes 64-m radio telescopes have been used to determine the average profile and timing properties for PSR J1757$-$2421 (B1754$-$24). Analysis of the radio profile data shows a large variation of spectral index across the profile and an unusual increase in component separation with increasing frequency. The timing observations show that PSR B1754$-$24 underwent a large glitch with a fractional increase in spin frequency of $\Delta\nu_{\rm g}/\nu\sim 7.8\times10^{-6}$ in May 2011. At the time of the glitch there was a large permanent jump in the rotational frequency accompanied by two exponential recovery terms with timescales of 15(6) and 97(15) days, respectively., Comment: 9 pages, 8 figures, 5 tables, Accepted for publication in MNRAS

In this paper, we investigate the simulation of continuous-time quantum walks on specific classes of graphs, for which it is possible to fast-forward the time-evolution operator to achieve constant-time simulation complexity and to perform the simulation exactly, while maintaining $\mbox{poly}(\mbox{log}(n))$ efficiency. In particular, we discuss two classes of composite graphs, commuting graphs and Cartesian product of graphs, that contain classes of graphs which can be simulated in this fashion. This allows us to identify new families of graphs that we can efficiently simulate in a quantum circuit framework, providing practical and explicit means to explore quantum-walk based algorithms in laboratories.

A major advantage in using Szegedy's formalism over discrete-time and continuous-time quantum walks lies in its ability to define a unitary quantum walk on directed and weighted graphs. In this paper, we present a general scheme to construct efficient quantum circuits for Szegedy quantum walks that correspond to classical Markov chains possessing transformational symmetry in the columns of the transition matrix. In particular, the transformational symmetry criteria do not necessarily depend on the sparsity of the transition matrix, so this scheme can be applied to non-sparse Markov chains. Two classes of Markov chains that are amenable to this construction are cyclic permutations and complete bipartite graphs, for which we provide explicit efficient quantum circuit implementations. We also prove that our scheme can be applied to Markov chains formed by a tensor product. We also briefly discuss the implementation of Markov chains based on weighted interdependent networks. In addition, we apply this scheme to construct efficient quantum circuits simulating the Szegedy walks used in the quantum Pagerank algorithm for some classes of non-trivial graphs, providing a necessary tool for experimental demonstration of the quantum Pagerank algorithm.

Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 months, and results in a permanently modified pulse shape. Profile variations such as these may be the origin of timing noise observed in other millisecond pulsars. The sensitivity of pulsar-timing observations to gravitational radiation can be increased by accounting for this variability., Comment: Accepted to ApJ Letters

Circulant matrices are an important family of operators, which have a wide range of applications in science and engineering related fields. They are in general non-sparse and non-unitary. In this paper, we present efficient quantum circuits to implement circulant operators using fewer resources and with lower complexity than existing methods. Moreover, our quantum circuits can be readily extended to the implementation of Toeplitz, Hankel, and block circulant matrices. Efficient quantum algorithms to implement the inverses and products of circulant operators are also provided., Comment: An example application in solving the equation of motion for a vibrating system with cyclic symmetry is added in Section 7

Various quantum-walk based algorithms have been proposed to analyse and rank the centrality of graph vertices. However, issues arise when working with directed graphs --- the resulting non-Hermitian Hamiltonian leads to non-unitary dynamics, and the total probability of the quantum walker is no longer conserved. In this paper, we discuss a method for simulating directed graphs using PT-symmetric quantum walks, allowing probability conserving non-unitary evolution. This method is equivalent to mapping the directed graph to an undirected, yet weighted, complete graph over the same vertex set, and can be extended to cover interdependent networks of directed graphs. Previous work has shown centrality measures based on the CTQW provide an eigenvector-like quantum centrality; using the PT-symmetric framework, we extend these centrality algorithms to directed graphs with a significantly reduced Hilbert space compared to previous proposals. In certain cases, this centrality measure provides an advantage over classical algorithms used in network analysis, for example by breaking vertex rank degeneracy. Finally, we perform a statistical analysis over ensembles of random graphs, and show strong agreement with the classical PageRank measure on directed acyclic graphs.

We present proper motions for 15 pulsars which are observed regularly by the Nanshan 25-m radio telescope. Two methods, the frequentist method (Coles et al.2011) and the Bayesian (Lentati et al. 2014) method, are used and the results are compared. We demonstrate that the two methods can be applied to young pulsar data sets that exhibit large amounts of timing noise with steep spectral exponents and give consistent results. The measured positions also agree with very-long-baseline interferometric positions. Proper motions for four pulsars are obtained for the first time, and improved values are obtained for five pulsars., Comment: 8 pages, 2 figures and 7 tables. Accepted for publication in MNRAS

Toeplitz and Hankel matrices have been a subject of intense interest in a wide range of science and engineering related applications. In this paper, we show that quantum circuits can efficiently implement sparse or Fourier-sparse Toeplitz and Hankel matrices. This provides an essential ingredient for solving many physical problems with Toeplitz or Hankel symmetry in the quantum setting with deterministic queries.

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only., Comment: 29 pages. 16 figures. Accepted for publication in MNRAS

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used., Comment: 25 pages, 6 tables, 5 figures. Accepted for publication in MNRAS

Detailed structures of the Vela pulsar (PSR B0833-45, with a period of $89.33$ milliseconds) are predicted by adopting a recently-constructed unified treatment of all parts of neutron stars: the outer crust, the inner crust and the core based on modern microscopic Brueckner-Hartree-Fock calculations. To take the pulsar mass ranging from $1.0M_{\odot}$ to $2.0M_{\odot}$, we calculate the central density, the core/crust radii, the core/crust mass, the core/crustal thickness, the moment of inertia, and the crustal moment of inertia. Among them, the crustal moment of inertia could be effectively constrained from the accumulated glitch observations, which has been a great debate recently, known as "glitch crisis". Namely, superfluid neutrons contained in the inner crust, which are regarded as the origin of the glitch in the standard two-component model, could be largely entrained in the nuclei lattices, then there may not be enough superfluid neutrons ($\sim 4/5$ less than the previous value) to trigger the large glitches ($\Delta \nu/\nu_0 \sim 10^{-6}$) in the Vela pulsar. We then provide the first analysis of the crisis based on the microscopic basis, by confronting the glitch observations with the theoretical calculations for the crustal moment of inertia. We find that despite some recent opposition to the crisis argument, the glitch crisis is still present, which means that besides the crust superfluid neutrons, core neutrons might be necessary for explaining the large glitches of the Vela pulsar., Comment: 28 pages, 6 figures, 4 tables