Source author record

V. Ravi

V. Ravi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2022arXiv

New Tests of Millilensing in the Blazar PKS 1413+135

Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational millilensing, PKS 1413+135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the millilens itself. We discuss two possible candidates for the putative millilens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989-day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in August 2022, May 2025, and February 2028.

preprint2022arXiv

The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate

Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise'' of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6σ$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.

preprint2021arXiv

Discovering New Strong Gravitational Lenses in the DESI Legacy Imaging Surveys

We have conducted a search for new strong gravitational lensing systems in the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys' Data Release 8. We use deep residual neural networks, building on previous work presented in Huang et al. (2020). These surveys together cover approximately one third of the sky visible from the northern hemisphere, reaching a z band AB magnitude of ~22.5. We compile a training sample that consists of known lensing systems as well as non-lenses in the Legacy Surveys and the Dark Energy Survey. After applying our trained neural networks to the survey data, we visually inspect and rank images with probabilities above a threshold. Here we present 1210 new strong lens candidates.

preprint2021arXiv

Finding Strong Gravitational Lenses in the DESI DECam Legacy Survey

We perform a semi-automated search for strong gravitational lensing systems in the 9,000 deg$^2$ Dark Energy Camera Legacy Survey (DECaLS), part of the DESI Legacy Imaging Surveys (Dey et al.). The combination of the depth and breadth of these surveys are unparalleled at this time, making them particularly suitable for discovering new strong gravitational lensing systems. We adopt the deep residual neural network architecture (He et al.) developed by Lanusse et al. for the purpose of finding strong lenses in photometric surveys. We compile a training set that consists of known lensing systems in the Legacy Surveys and DES as well as non-lenses in the footprint of DECaLS. In this paper we show the results of applying our trained neural network to the cutout images centered on galaxies typed as ellipticals (Lang et al.) in DECaLS. The images that receive the highest scores (probabilities) are visually inspected and ranked. Here we present 335 candidate strong lensing systems, identified for the first time.

preprint2020arXiv

Radio observations of two intermittent pulsars: PSRs J1832+0029 and J1841-0500

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.

preprint2020arXiv

The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135

PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a "multizone" model is needed to account for the high-energy emission.

preprint2016arXiv

Are the distributions of Fast Radio Burst properties consistent with a cosmological population?

High time resolution radio surveys over the last few years have discovered a population of millisecond-duration transient bursts called Fast Radio Bursts (FRBs), which remain of unknown origin. FRBs exhibit dispersion consistent with propagation through a cold plasma and dispersion measures indicative of an origin at cosmological distances. In this paper we perform Monte Carlo simulations of a cosmological population of FRBs, based on assumptions consistent with observations of their energy distribution, their spatial density as a function of redshift and the properties of the interstellar and intergalactic media. We examine whether the dispersion measures, fluences, inferred redshifts, signal-to-noises and effective widths of known FRBs are consistent with a cosmological population. Statistical analyses indicate that at least 50 events at Parkes are required to distinguish between a constant co-moving FRB density, and a FRB density that evolves with redshift like the cosmological star formation rate density.

preprint2016arXiv

Fast Radio Transient searches with UTMOST at 843 MHz

We report the first radio interferometric search at 843 MHz for fast transients, particularly Fast Radio Bursts (FRBs). The recently recommissioned Swinburne University of Technology's digital backend for the Molonglo Observatory Synthesis Telescope array (the UTMOST) with its large collecting area (18,000 $\mathrm{m^2}$) and wide instantaneous field of view (7.80 $\mathrm{deg^2}$) is expected to be an efficient tool to detect FRBs. As an interferometer it will be capable of discerning whether the FRBs are truly a celestial population. We show that UTMOST at full design sensitivity can detect an event approximately every few days. We report on 2 preliminary FRB surveys at about 7% and 14% respectively of the array's final sensitivity. Several pulsars have been detected via single pulses and no FRBs were discovered with pulse widths ($W$), in the range 655.36 $μ$s $< W < 41.9$ ms and dispersion measures (DMs) in the range $100 < $DM$< 2000$ $\mathrm{pc\,cm^{-3}}$. This non-detection sets a 2$σ$ upper limit of the sky rate of not more than 1000 events $\mathrm{sky^{-1}}$ $\mathrm{day^{-1}}$ at 843 MHz down to a flux limit of 11 Jy for 1 ms FRBs. We show that this limit is consistent with previous survey limits at 1.4 GHz and 145 MHz and set a lower limit on the mean spectral index of FRBs of $α> -3.2$.

preprint2016arXiv

Radio detections of southern ultracool dwarfs

We report the results of a volume-limited survey using the Australia Telescope Compact Array to search for transient and quiescent radio emission from 15 southern hemisphere ultracool dwarfs. We detect radio emission from 2MASSW J0004348-404405 increasing the number of radio loud ultracool dwarfs to 22. We also observe radio emission from 2MASS J10481463-3956062 and 2MASSI J0339352-352544, two sources with previous radio detections. The radio emission from the three detected sources shows no variability or flare emission. Modelling this quiescent emission we find that it is consistent with optically thin gyrosynchrotron emission from a magnetosphere with an emitting region radius of (1 - 2)$R_*$, magnetic field inclination 20$^{\circ}$ - 80$^{\circ}$, field strength $\sim$10 - 200 G, and power-law electron density $\sim$10$^4$ - 10$^8$ cm$^{-3}$. Additionally, we place upper limits on four ultracool dwarfs with no previous radio observations. This increases the number of ultracool dwarfs studied at radio frequencies to 222. Analysing general trends of the radio emission for this sample of 15 sources, we find that the radio activity increases for later spectral types and more rapidly rotating objects. Furthermore, comparing the ratio of the radio to X-ray luminosities for these sources, we find 2MASS J10481463-3956062 and 2MASSI J0339352-352544 violate the Guedel-Benz relation by more than two orders of magnitude.

preprint2016arXiv

The disturbance of a millisecond pulsar magnetosphere

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.

preprint2016arXiv

The Fluence and Distance Distributions of Fast Radio Bursts

Fast radio bursts (FRB) are millisecond-duration radio pulses with apparent extragalactic origins. All but two of the FRBs have been discovered using the Parkes dish which employs multiple beams formed by an array of feed horns on its focal plane. In this paper, we show that (i) the preponderance of multiple-beam detections, and (ii) the detection rates for varying dish diameters, can be used to infer the index $α$ of the cumulative fluence distribution function (the log$N$-log$F$ function: $α=1.5$ for a non-evolving population in a Euclidean universe). If all detected FRBs arise from a single progenitor population, multiple-beam FRB detection rates from the Parkes telescope yield the constraint $0.52<α<1.0$ with $90$% confidence. Searches at other facilities with different dish sizes refine the constraint to $0.66<α<0.96$. Our results favor FRB searches with smaller dishes, because for $α<1$, the gain in field-of-view for a smaller dish is more important than the reduction in sensitivity. Further, our results suggest that (i) FRBs are not standard candles, and (ii) the distribution of distances to the detected FRBs is weighted towards larger distances. If FRBs are extragalactic, these results are consistent with a cosmological population, which would make FRBs excellent probes of the baryonic content and geometry of the Universe.

preprint2016arXiv

The magnetic field and turbulence of the cosmic web measured using a brilliant fast radio burst

Fast radio bursts (FRBs) are millisecond-duration events thought to originate beyond the Milky Way galaxy. Uncertainty surrounding the burst sources, and their propagation through intervening plasma, has limited their use as cosmological probes. We report on a mildly dispersed (dispersion measure 266.5+-0.1 pc cm^-3), exceptionally intense (120+-30 Jy), linearly polarized, scintillating burst (FRB 150807) that we directly localize to 9 arcmin^2. Based on a low Faraday rotation (12.0+-0.7 rad m^-2), we infer negligible magnetization in the circum-burst plasma and constrain the net magnetization of the cosmic web along this sightline to <21 nG, parallel to the line-of-sight. The burst scintillation suggests weak turbulence in the ionized intergalactic medium.

preprint2015arXiv

A fast radio burst in the direction of the Carina dwarf spheroidal galaxy

We report the real-time discovery of a fast radio burst (FRB 131104) with the Parkes radio telescope in a targeted observation of the Carina dwarf spheroidal galaxy. The dispersion measure of the burst is 779 cm$^{-3}$ pc, exceeding predictions for the maximum line-of-sight Galactic contribution by a factor of 11. The temporal structure of the burst is characterized by an exponential scattering tail with a timescale of 2.0$^{+0.8}_{-0.5}$ ms at 1582 MHz that scales as frequency to the power $-$4.4$^{+1.6}_{-1.8}$ (all uncertainties represent 95% confidence intervals). We bound the intrinsic pulse width to be $<0.64$ ms due to dispersion smearing across a single spectrometer channel. Searches in 78 hours of follow-up observations with the Parkes telescope reveal no additional sporadic emission and no evidence for associated periodic radio emission. We hypothesize that the burst is associated with the Carina dwarf galaxy. Follow-up observations at other wavelengths are necessary to test this hypothesis.

preprint2015arXiv

A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $λ^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

preprint2015arXiv

Gravitational waves from binary supermassive black holes missing in pulsar observations

Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.

preprint2015arXiv

Pulsar Observations of Extreme Scattering Events

Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular flux measurements of compact extragalactic sources. They are characterized by a flux variation over a period of weeks, suggesting the passage of a "diverging plasma lens" across the line of sight. Modeling the refraction of such a lens indicates that the structure size must be of order AU and the electron density of order 10s of cm^{-3}. Similar structures have been observed in measurements of pulsar intensity scintillation and group delay. Here we report observations of two ESEs showing increases in both intensity scintillation and dispersion made with the Parkes Pulsar Timing Array (PPTA). These allow us to make more complete models of the ESE, including an estimate of the "outer-scale" of the turbulence in the plasma lens. These observations show clearly that the ESE structure is fully turbulent on an AU scale. They provide some support for the idea that the structures are extended along the line of sight, such as would be the case for a scattering shell. The dispersion measurements also show a variety of AU scale structures which would not be called ESEs, yet involve electron density variations typical of ESEs and likely have the same origin.

preprint2015arXiv

Timing analysis for 20 millisecond pulsars in the Parkes Pulsar Timing Array

We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024-0719, J1045-4509, J1600-3053, J1603-7202, and J1730-2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437-4715 and J1909-3744 with $M_p=1.44\pm0.07$ $M_\odot$ and $M_p=1.47\pm0.03$ $M_\odot$ respectively. The improved orbital period-derivative measurement for PSR J0437-4715 results in a derived distance measurement at the 0.16% level of precision, $D=156.79\pm0.25$ pc, one of the most fractionally precise distance measurements of any star to date.

preprint2015arXiv

Versatile Directional Searches for Gravitational Waves with Pulsar Timing Arrays

By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be "phased-up" to form time series of the two polarisation modes of GWs coming from a particular direction of the sky. Our technique requires no assumptions regarding the time-domain behaviour of a GW signal. This method has already been used to place stringent bounds on GWs from individual SMBHBs in circular orbits. Here, we describe the methodology and demonstrate the versatility of the technique in searches for a wide variety of GW signals including bursts with unmodeled waveforms. Using the first six years of data from the Parkes Pulsar Timing Array, we conduct an all-sky search for a detectable excess of GW power from any direction. For the lines of sight to several nearby massive galaxy clusters, we carry out a more detailed search for GW bursts with memory, which are distinct signatures of SMBHB mergers. In all cases, we find that the data are consistent with noise.

preprint2014arXiv

An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set

We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from $5\times 10^{-9}$ to $2\times 10^{-7}$ Hz. Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude $h_0$ for a range of GW frequencies. For example, at $10^{-8}$ Hz our analysis has excluded with $95\%$ confidence the presence of signals with $h_0\geqslant 1.7\times 10^{-14}$. Our new limits are about a factor of four more stringent than those of Yardley et al. (2010) based on an earlier PPTA data set and a factor of two better than those reported in the recent Arzoumanian et al. (2014) paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of $10^{9} M_{\odot}$ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of $4 \times 10^{-3} {\rm{Mpc}}^{-3} {\rm{Gyr}}^{-1}$ at $95\%$ confidence on the coalescence rate of nearby ($z \lesssim 0.1$) supermassive binary black holes in circular orbits with chirp masses of $10^{10}M_{\odot}$.

preprint2014arXiv

Binary super-massive black hole environments diminish the gravitational-wave signal in the pulsar timing band

We assess the effects of super-massive black hole (SMBH) environments on the gravitational-wave (GW) signal from binary SMBHs. To date, searches with pulsar timing arrays for GWs from binary SMBHs, in the frequency band $\sim1-100$\,nHz, include the assumptions that all binaries are circular and evolve only through GW emission. However, dynamical studies have shown that the only way that binary SMBH orbits can decay to separations where GW emission dominates the evolution is through interactions with their environments. We augment an existing galaxy and SMBH formation and evolution model with calculations of binary SMBH evolution in stellar environments, accounting for non-zero binary eccentricities. We find that coupling between binaries and their environments causes the expected GW spectral energy distribution to be reduced with respect to the standard assumption of circular, GW-driven binaries, for frequencies up to $\sim20$\,nHz. Larger eccentricities at binary formation further reduce the signal in this regime. We also find that GW bursts from individual eccentric binary SMBHs are unlikely to be detectable with current pulsar timing arrays. The uncertainties in these predictions are large, owing to observational uncertainty in SMBH-galaxy scaling relations and the galaxy stellar mass function, uncertainty in the nature of binary-environment coupling, and uncertainty in the numbers of the most massive binary SMBHs. We conclude, however, that low-frequency GWs from binary SMBHs may be more difficult to detect with pulsar timing arrays than currently thought.

preprint2014arXiv

Limitations in timing precision due to single-pulse shape variability in millisecond pulsars

High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intra-observation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909-3744 shows the lowest levels of jitter noise, which we estimate to contribute $\sim$10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437-4715 decorrelates over a bandwidth of $\sim$2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.

preprint2014arXiv

Prospects for gravitational-wave detection and supermassive black hole astrophysics with pulsar timing arrays

Large-area sky surveys show that massive galaxies undergo at least one major merger in a Hubble time. Ongoing pulsar timing array (PTA) experiments are aimed at measuring the gravitational wave (GW) emission from binary supermassive black holes (SMBHs) at the centres of galaxy merger remnants. In this paper, using the latest observational estimates for a range of galaxy properties and scaling relations, we predict the amplitude of the GW background generated by the binary SMBH population. We also predict the numbers of individual binary SMBH GW sources. We predict the characteristic strain amplitude of the GW background to lie in the range $5.1\times10^{-16}<A_{\rm yr}<2.4\times10^{-15}$ at a frequency of $(1\,{\rm yr})^{-1}$, with 95% confidence. Higher values within this range, which correspond to the more commonly preferred choice of galaxy merger timescale, will fall within the expected sensitivity ranges of existing PTA projects in the next few years. In contrast, we find that a PTA consisting of at least 100 pulsars observed with next-generation radio telescopes will be required to detect continuous-wave GWs from binary SMBHs. We further suggest that GW memory bursts from coalescing SMBH pairs are not viable sources for PTAs. Both the GW background and individual GW source counts are dominated by binaries formed in mergers between early-type galaxies of masses $\gtrsim5\times10^{10}M_{\odot}$ at redshifts $\lesssim 1.5$. Uncertainties in the galaxy merger timescale and the SMBH mass - galaxy bulge mass relation dominate the uncertainty in our predictions.

preprint2014arXiv

Searching for gravitational wave memory bursts with the Parkes Pulsar Timing Array

Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected.

preprint2013arXiv

Gravitational-wave Limits from Pulsar Timing Constrain Supermassive Black Hole Evolution

The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be $Ω_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability.

preprint2012arXiv

Detection of 107 glitches in 36 southern pulsars

Timing observations from the Parkes 64-m radio telescope for 165 pulsars between 1990 and 2011 have been searched for period glitches. A total of 107 glitches were identified in 36 pulsars, where 61 have previously been reported and 46 are new discoveries. Glitch parameters were measured by fitting the timing residual data. Observed relative glitch sizes Δν_g/νrange between 10^-10 and 10^-5, where ν= 1/P is the pulse frequency. We confirm that the distribution of Δν_g/νis bimodal with peaks at approximately 10^-9 and 10^-6. Glitches are mostly observed in pulsars with characteristic ages between 10^3 and 10^5 years, with large glitches mostly occurring in the younger pulsars. Exponential post-glitch recoveries were observed for 27 large glitches in 18 pulsars. The fraction Q of the glitch that recovers exponentially also has a bimodal distribution. Large glitches generally have low Q, typically a few per cent, but large Q values are observed in both large and small glitches. Observed time constants for exponential recoveries ranged between 10 and 300 days with some tendency for longer timescales in older pulsars. Shorter timescale recoveries may exist but were not revealed by our data which typically have observation intervals of 2 - 4 weeks. For most of the 36 pulsars with observed glitches, there is a persistent linear increase in \dotνin the inter-glitch interval. Where an exponential recovery is also observed, the effects of this are superimposed on the linear increase in \dotν. In some cases, the slope of the linear recovery changes at the time of a glitch. The \ddotνvalues characterising the linear changes in \dotνare almost always positive and, after subtracting the magnetospheric component of the braking, are approximately proportional to the ratio of |\dotν| and the inter-glitch interval, as predicted by vortex-creep models.

preprint2012arXiv

Development of a pulsar-based timescale

Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.

preprint2012arXiv

Does a 'stochastic' background of gravitational waves exist in the pulsar timing band?

We investigate the effects of gravitational waves (GWs) from a simulated population of binary super-massive black holes (SMBHs) on pulsar timing array datasets. We construct a distribution describing the binary SMBH population from an existing semi-analytic galaxy formation model. Using realizations of the binary SMBH population generated from this distribution, we simulate pulsar timing datasets with GW-induced variations. We find that the statistics of these variations do not correspond to an isotropic, stochastic GW background. The "Hellings & Downs" correlations between simulated datasets for different pulsars are recovered on average, though the scatter of the correlation estimates is greater than expected for an isotropic, stochastic GW background. These results are attributable to the fact that just a few GW sources dominate the GW-induced variations in every Fourier frequency bin of a 5-year dataset. Current constraints on the amplitude of the GW signal from binary SMBHs will be biased. Individual binary systems are likely to be detectable in 5-year pulsar timing array datasets where the noise is dominated by GW-induced variations. Searches for GWs in pulsar timing array data therefore need to account for the effects of individual sources of GWs.

preprint2012arXiv

Measurement and correction of variations in interstellar dispersion in high-precision pulsar timing

Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparallelled measurement of the turbulent interstellar plasma at au scales. In this paper we demonstrate that without consideration of wavelength- independent red-noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimise this additional noise, and factors, such as scintillation, that can exacerbate the problem. Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603-7202 of ~2x10^{-3} cm^{-3}pc for about 250 days. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously.

preprint2012arXiv

Observations and modelling of pulsed radio emission from CU Virginis

We present 13 cm and 20 cm radio observations of the magnetic chemically peculiar star CU Virginis taken with the Australia Telescope Compact Array. We detect two circularly polarised radio pulses every rotation period which confirm previous detections. In the first pulse, the lower frequency emission arrives before the higher frequency emission and the ordering reverses in the second pulse. In order to explain the frequency dependence of the time between the two pulses, we construct a geometric model of the magnetosphere of CU Virginis, and consider various emission angles relative to the magnetic field lines. A simple electron cyclotron maser emission model, in which the emission is perpendicular to the magnetic field lines, is not consistent with our data. A model in which the emission is refracted through cold plasma in the magnetosphere is shown to have the correct pulse arrival time frequency dependence.

preprint2012arXiv

The Parkes Pulsar Timing Array Project

A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.

preprint2011arXiv

The magnetosphere of the ultracool dwarf DENIS 1048-3956

Ultracool dwarfs, the least-massive contributors to the stellar mass function, exhibit striking magnetic properties that are inconsistent with trends for more massive stars. Here, we present the widest-band radio observations to date of an ultracool dwarf, DENIS-P J104814.9-395604, in four 2 GHz bandwidths between wavelengths of 1 cm and 10 cm. These data were obtained with the Australia Telescope Compact Array using the new Compact Array Broadband Backend instrument. We detected a stable negatively-sloped power-law spectrum in total intensity, with spectral index alpha=1.71+-0.09. Circular polarization fractions between 0.25 and 0.4 were found at the low-frequency end of our detection band. We interpret these results as indicative of gyrosynchrotron emission. We suggest that the radio emission originates from beyond the co-rotation radius, R_C, of the star. Adopting this model, we find R_C between 1.2-2.9 R_*, and a non-thermal electron density and magnetic field strength between 10^(5)-10^(7.2) cm^(-3) and 70-260 G respectively at R_C. The model accounts for the violation of the Guedel-Benz relation between X-ray and radio luminosities of low-mass stars by DENIS-P J104814.9-395604.

preprint2011arXiv

The non-uniform, dynamic atmosphere of Betelgeuse observed at mid-infrared wavelengths

We present an interferometric study of the continuum surface of the red supergiant star Betelgeuse at 11.15 microns wavelength, using data obtained with the Berkeley Infrared Spatial Interferometer each year between 2006 and 2010. These data allow an investigation of an optically thick layer within 1.4 stellar radii of the photosphere. The layer has an optical depth of ~1 at 11.15 microns, and varies in temperature between 1900 K and 2800 K and in outer radius between 1.16 and 1.36 stellar radii. Electron-hydrogen atom collisions contribute significantly to the opacity of the layer. The layer has a non-uniform intensity distribution that changes between observing epochs. These results indicate that large-scale surface convective activity strongly influences the dynamics of the inner atmosphere of Betelgeuse, and mass-loss processes.

preprint2010arXiv

Observations of radio pulses from CU Virginis

The magnetic chemically peculiar star CU Virginis is a unique astrophysical laboratory for stellar magnetospheres and coherent emission processes. It is the only known main sequence star to emit a radio pulse every rotation period. Here we report on new observations of the CU Virginis pulse profile in the 13 and 20\,cm radio bands. The profile is known to be characterised by two peaks of 100$\%$ circularly polarised emission that are thought to arise in an electron-cyclotron maser mechanism. We find that the trailing peak is stable at both 13 and 20\,cm, whereas the leading peak is intermittent at 13\,cm. Our measured pulse arrival times confirm the discrepancy previously reported between the putative stellar rotation rates measured with optical data and with radio observations. We suggest that this period discrepancy might be caused by an unknown companion or by instabilities in the emission region. Regular long-term pulse timing and simultaneous multi-wavelength observations are essential to clarify the behaviour of this emerging class of transient radio source.

preprint2010arXiv

The many faces of Betelgeuse

The dynamics of the surface and inner atmosphere of the red supergiant star Betelgeuse are the subject of numerous high angular resolution and spectroscopic studies. Here, we present three-telescope interferometric data obtained at 11.15 microns wavelength with the Berkeley Infrared Spatial Interferometer (ISI), that probe the stellar surface continuum. We find striking variability in the size, effective temperature, and degree of asymmetry of the star over the years 2006-2009. These results may indicate an evolving shell of optically thick material close to the stellar photosphere.

preprint2010arXiv

Wide radio beams from gamma-ray pulsars

We investigate the radio and gamma-ray beaming properties of normal and millisecond pulsars by selecting two samples from the known populations. The first, Sample G, contains pulsars which are detectable in blind searches of gamma-ray data from the Fermi Large Area Telescope. The second, Sample R, contains pulsars detectable in blind radio searches which have spin-down luminosities Edot > 10^{34} erg/s. We analyse the fraction of the gamma-ray-selected Sample G which have detectable radio pulses and the fraction of the radio-selected Sample R which have detectable gamma-ray pulses. Twenty of our 35 Sample G pulsars have already observed radio pulses. This rules out low-altitude polar-cap beaming models if, as is currently believed, gamma-ray beams are generated in the outer magnetosphere and are very wide. We further find that, for the highest-Edot pulsars, the radio and gamma-ray beams have comparable beaming factors, i.e., the beams cover similar regions of the sky as the star rotates. For lower-Edot gamma-ray emitting pulsars, the radio beams have about half of the gamma-ray sky coverage. These results suggest that, for high-Edot young and millisecond pulsars, the radio emission originates in wide beams from regions high in the pulsar magnetosphere, probably close to the null-charge surface and to the gamma-ray emitting regions. Furthermore, it suggests that for these high-Edot pulsars, as in the gamma-ray case, features in the radio profile represent caustics in the emission beam pattern.