Source author record

T. J. W. Lazio

T. J. W. Lazio 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

54works
9topics
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

54 published item(s)

preprint2022arXiv

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

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})^α$, 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 $α= -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 $α= -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.

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.

preprint2022arXiv

Toward Astrometric Constraints on a Supermassive Black Hole Binary in the Early-Type Galaxy NGC\,4472

The merger of two galaxies, each hosting a supermassive black hole (SMBH) of mass $10^6$\,M$_{\odot}$ or more, could yield a bound SMBH binary. For the early-type galaxy NGC\,4472, we study how astrometry with a next-generation Very Large Array (ngVLA) could be used to monitor the reflex motion of the primary SMBH of mass $M_{\rm pri}$, as it is tugged on by the secondary SMBH of mass $M_{\rm sec}$. Casting the orbit of the putative SMBH binary in terms of its period $P$, semimajor axis $a_{\rm bin}$, and mass ratio $q = M_{\rm sec} / M_{\rm pri} \le 1$, we find the following: (1) Orbits with fiducial periods of $P = 4$\,yr and 40\,yr could be spatially resolved and monitored. (2) For a 95\% accuracy of $2\,μ$as per monitoring epoch, sub-parsec values of $a_{\rm bin}$ could be accessed over a range of mass ratios notionally encompassing major ($q > \frac{1}{4}$) and minor ($q < \frac{1}{4}$) galaxy mergers. (3) If no reflex motion is detected for $M_{\rm pri}$ after 1(10)\,yr of monitoring, a SMBH binary with period $P = 4(40)$\,yr and mass ratio $q > 0.01(0.003)$ could be excluded. This would suggest no present-day evidence for a past major merger like that recently simulated, where scouring by a $q \sim 1$ SMBH binary formed a stellar core with kinematic traits like those of NGC\,4472. (4) Astrometric monitoring could independently check the upper limits on $q$ from searches for continuous gravitational waves from NGC\,4472.

preprint2021arXiv

A repeating fast radio burst source in a globular cluster

Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (Margalit & Metzger 2018, CHIME/FRB Collaboration et al. 2020). Recently, CHIME/FRB announced the discovery (Bhardwaj et al. 2021) of the repeating FRB 20200120E, coming from the direction of the nearby grand design spiral galaxy M81. Four potential counterparts at other observing wavelengths were identified (Bhardwaj et al. 2021) but no definitive association with these sources, or M81, could be made. Here we report an extremely precise localisation of FRB 20200120E, which allows us to associate it with a globular cluster (GC) in the M81 galactic system and to place it ~2pc offset from the optical center of light of the GC. This confirms (Bhardwaj et al. 2021) that FRB 20200120E is 40 times closer than any other known extragalactic FRB. Because such GCs host old stellar populations, this association strongly challenges FRB models that invoke young magnetars formed in a core-collapse supernova as powering FRB emission. We propose, instead, that FRB 20200120E is a highly magnetised neutron star formed via either accretion-induced collapse of a white dwarf or via merger of compact stars in a binary system (Margalit et al. 2019). Alternative scenarios involving compact binary systems, efficiently formed inside globular clusters, could also be responsible for the observed bursts.

preprint2021arXiv

Late-Time Evolution and Modeling of the Off-Axis Gamma-ray Burst Candidate FIRST J141918.9+394036

We present new radio and optical data, including very long baseline interferometry, as well as archival data analysis, for the luminous decades-long radio transient FIRST J141918.9+394036. The radio data reveal a synchrotron self-absorption peak around 0.3 GHz and a radius of around 1.3 mas (0.5 pc) 26 years post-discovery, indicating a blastwave energy $\sim5 \times 10^{50}$ erg. The optical spectrum shows a broad [OIII]$λ$4959,5007 emission-line that may indicate collisional-excitation in the host galaxy, but its association with the transient cannot be ruled out. The properties of the host galaxy are suggestive of a massive stellar progenitor that formed at low metallicity. Based on the radio light curve, blastwave velocity, energetics, nature of the host galaxy and transient rates we find that the properties of FIRST J1419+39 are most consistent with long gamma-ray burst (LGRB) afterglows. Other classes of (optically-discovered) stellar explosions as well as neutron star mergers are disfavored, and invoking any exotic scenario may not be necessary. It is therefore likely that FIRST J1419+39 is an off-axis LGRB afterglow (as suggested by Law et al. and Marcote et al.), and under this premise the inverse beaming fraction is found to be $f_b^{-1}\simeq280^{+700}_{-200}$, corresponding to an average jet half-opening angle $<θ_j>\simeq5^{+4}_{-2}$ degrees (68% confidence), consistent with previous estimates. From the volumetric rate we predict that surveys with the VLA, ASKAP and MeerKAT will find a handful of FIRST J1419+39-like events over the coming years.

preprint2020arXiv

A Demonstration of Extremely Low Latency $γ$-ray, X-Ray & UV Follow-Up of a Millisecond Radio Transient

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5σ$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3σ$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science.

preprint2020arXiv

A repeating fast radio burst source localised to a nearby spiral galaxy

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

preprint2020arXiv

Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

preprint2020arXiv

On Frequency-Dependent Dispersion Measures and Extreme Scattering Events

Radio emission propagating over an Earth-pulsar line of sight provides a unique probe of the intervening ionized interstellar medium (ISM). Variations in the integrated electron column density along this line of sight, or dispersion measure (DM), have been observed since shortly after the discovery of pulsars. As early as 2006, frequency-dependent dispersion measures have been observed and attributed to several possible causes. Ray-path averaging over different effective light-cone volumes through the turbulent ISM contributes to this effect as will DM misestimation due to radio propagation across compact lensing structures such as those caused by "extreme scattering events". We present methods to assess the variations in frequency-dependent dispersion measures due to the turbulent ISM versus these compact lensing structures along the line of sight. We analyze recent Low-Frequency Array (LOFAR) observations of PSR J2219+4754 to test the underlying physical mechanism of the observed frequency-dependent DM. Previous analyses have indicated the presence of strong lensing due to compact overdensities halfway between the Earth and pulsar. Instead we find the frequency dependence of the DM timeseries for PSR J2219+4754 is consistent with being due solely to ISM turbulence and there is no evidence for any extreme scattering event or small-scale lensing structure. The data show possible deviations from a uniform turbulent medium, suggesting that there may be an enhanced scattering screen near one of the two ends of the line of sight. We present this analysis as an example of the power of low-frequency observations to distinguish the underlying mechanisms in frequency-dependent propagation effects.

preprint2020arXiv

Radio Observations of HD80606 Near Planetary Periastron: II. LOFAR Low Band Antenna Observations at 30-78 MHz

All the giant planets in the solar system generate radio emission via the electron cyclotron maser instability, most notably giving rise to Jupiter's decametric emissions. An interaction with the solar wind is at least partially responsible for all of these solar system electron cyclotron masers. HD80606b is a giant planet with a highly eccentric orbit, leading to predictions that its radio emission may be enhanced substantially near periastron. This paper reports observations with the Low Frequency Array (LOFAR) of HD80606b near its periastron in an effort to detect radio emissions generated by an electron cyclotron maser instability in the planet's magnetosphere. The reported observations are at frequencies between 30 MHz and 78 MHz, and they are distinguished from most previous radio observations of extrasolar planets by two factors: (i) They are at frequencies near 50 MHz, much closer to the frequencies at which Jupiter emits (< 40 MHz) and lower than most previously reported observations of extrasolar planets; and (ii) Sensitivities of approximately a few millijanskys have been achieved, an order of magnitude or more below nearly all previous extrasolar planet observations below 100 MHz. We do not detect any radio emissions from HD80606b and use these observations to place new constraints on its radio luminosity. We also revisit whether the observations were conducted at a time when it was super-Alfvenic relative to the host star's stellar wind, which experience from the solar system illustrates is a state in which an electron cyclotron maser emission can be sustained in a planet's magnetic polar regions.

preprint2020arXiv

VLA/Realfast Detection of a Burst from FRB180916.J0158+65 and Tests for Periodic Activity

We report on the detection of a burst from FRB180916 by realfast/VLA and present software for interpreting fast radio bursts (FRB) periodicity. We demonstrate a range of periodicity analyses with bursts from FRB180916, FRB121102 and FRB180814. Our results for FRB180916 and FRB121102 are consistent with published results. For FRB180814, we did not detect any significant periodic episodes. The realfast-detected and other high-frequency bursts for FRB180916 tend to lie at the beginning of the activity window, indicating a possible phase-frequency relation. The python package $\texttt{frbpa}$ can be used to reproduce and expand on this analysis to test models for repeating FRBs.

preprint2019arXiv

A pulsar-based timescale from the International Pulsar Timing Array

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.

preprint2019arXiv

A Search for Late-Time Radio Emission and Fast Radio Bursts from Superluminous Supernovae

We present results of a search for late-time radio emission and Fast Radio Bursts (FRBs) from a sample of type-I superluminous supernovae (SLSNe-I). We used the Karl G. Jansky Very Large Array to observe ten SLSN-I more than 5 years old at a frequency of 3 GHz. We searched fast-sampled visibilities for FRBs and used the same data to perform a deep imaging search for late-time radio emission expected in models of magnetar-powered supernovae. No FRBs were found. One SLSN-I, PTF10hgi, is detected in deep imaging, corresponding to a luminosity of $1.2\times10^{28}$ erg s$^{-1}$. This luminosity, considered with the recent 6 GHz detection of PTF10hgi in Eftekhari et al (2019), supports the interpretation that it is powered by a young, fast-spinning ($\sim$ ms spin period) magnetar with $\sim$ 15 Msun of partially ionized ejecta. Broadly, our observations are most consistent with SLSNe-I being powered by neutron stars with fast spin periods, although most require more free-free absorption than is inferred for PTF10hgi. We predict that radio observations at higher frequencies or in the near future will detect these systems and begin constraining properties of the young pulsars and their birth environments.

preprint2019arXiv

Estimates of Fast Radio Burst Dispersion Measures from Cosmological Simulations

We calculate the dispersion measure (DM) contributed by the intergalactic medium (IGM) to the total measured DM for fast radio bursts (FRBs). We use the MareNostrum Instituto de Ciencias del Espacio (MICE) Onion Universe simulation (Fosalba et al. 2008) to track the evolution of the dark matter particle density over a large range of redshifts. We convert this dark matter particle number density to the corresponding free electron density and then integrate it to find the DM as a function of redshift. This approach yields an intergalactic DM of $\rm DM_{\rm IGM}(z = 1) = 800^{+7000}_{-170}$ pc cm$^{-3}$, with the large errors representative of the structure in the IGM. We place limits on the redshifts of the current population of observed FRBs. We also use our results to estimate the host galaxy contribution to the DM for the first repeater, FRB 121102, and show that the most probable host galaxy DM contribution, $\rm DM_{\rm host} \approx 310$ pc cm$^{-3}$, is consistent with the estimate made using the Balmer emission lines in the spectrum of the host galaxy, $\rm DM_{\rm Balmer} = 324$ pc cm$^{-3}$ (Tendulkar et al. 2017). We also compare our predictions for the host galaxy contribution to the DM for the observations of FRB 180924 (Bannister et al. 2019) and FRB 190523 (Ravi et al. 2019), both of which have been localized to a host galaxy.

preprint2019arXiv

The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science case and survey design

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$σ$ goal of 70 $μ$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

preprint2019arXiv

The NANOGrav 11-Year Data Set: Evolution of Gravitational Wave Background Statistics

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

preprint2019arXiv

The NANOGrav 11-Year Data Set: Limits on Gravitational Wave Memory

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $ηM \sim 2\times10^{10}M_\odot$ and inclination of $ι=π/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

preprint2016arXiv

Are we there yet? Time to detection of nanohertz gravitational waves based on pulsar-timing array limits

Decade-long timing observations of arrays of millisecond pulsars have placed highly constraining upper limits on the amplitude of the nanohertz gravitational-wave stochastic signal from the mergers of supermassive black-hole binaries ($\sim 10^{-15}$ strain at $f = 1/\mathrm{yr}$). These limits suggest that binary merger rates have been overestimated, or that environmental influences from nuclear gas or stars accelerate orbital decay, reducing the gravitational-wave signal at the lowest, most sensitive frequencies. This prompts the question whether nanohertz gravitational waves are likely to be detected in the near future. In this letter, we answer this question quantitatively using simple statistical estimates, deriving the range of true signal amplitudes that are compatible with current upper limits, and computing expected detection probabilities as a function of observation time. We conclude that small arrays consisting of the pulsars with the least timing noise, which yield the tightest upper limits, have discouraging prospects of making a detection in the next two decades. By contrast, we find large arrays are crucial to detection because the quadrupolar spatial correlations induced by gravitational waves can be well sampled by many pulsar pairs. Indeed, timing programs which monitor a large and expanding set of pulsars have an $\sim 80\%$ probability of detecting gravitational waves within the next ten years, under assumptions on merger rates and environmental influences ranging from optimistic to conservative. Even in the extreme case where $90\%$ of binaries stall before merger and environmental coupling effects diminish low-frequency gravitational-wave power, detection is delayed by at most a few years.

preprint2016arXiv

From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release

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.

preprint2016arXiv

Microarcsecond VLBI pulsar astrometry with PSRPI I. Two binary millisecond pulsars with white dwarf companions

Model-independent distance constraints to binary millisecond pulsars (MSPs) are of great value to both the timing observations of the radio pulsars, and multiwavelength observations of their companion stars. Very Long Baseline Interferometry (VLBI) astrometry can be employed to provide these model-independent distances with very high precision via the detection of annual geometric parallax. Using the Very Long Baseline Array, we have observed two binary millisecond pulsars, PSR J1022+1001 and J2145-0750, over a two-year period and measured their distances to be 700 +14 -10 pc and 613 +16 -14 pc respectively. We use the well-calibrated distance in conjunction with revised analysis of optical photometry to tightly constrain the nature of their massive (M ~ 0.85 Msun) white dwarf companions. Finally, we show that several measurements of their parallax and proper motion of PSR J1022+1001 and PSR J2145-0750 obtained by pulsar timing array projects are incorrect, differing from the more precise VLBI values by up to 5 sigma. We investigate possible causes for the discrepancy, and find that imperfect modeling of the solar wind is a likely candidate for the timing model errors given the low ecliptic latitude of these two pulsars.

preprint2016arXiv

The International Pulsar Timing Array: First Data Release

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.

preprint2016arXiv

The NANOGrav Nine-year Data Set: Observations, Arrival Time Measurements, and Analysis of 37 Millisecond Pulsars

We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.

preprint2016arXiv

Three-dimensional Tomography of the Galactic and Extragalactic Magnetoionic Medium with the SKA

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way.

preprint2015arXiv

NANOGrav Constraints on Gravitational Wave Bursts with Memory

Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.

preprint2015arXiv

Observing Radio Pulsars in the Galactic Centre with the Square Kilometre Array

The discovery and timing of radio pulsars within the Galactic centre is a fundamental aspect of the SKA Science Case, responding to the topic of "Strong Field Tests of Gravity with Pulsars and Black Holes" (Kramer et al. 2004; Cordes et al. 2004). Pulsars have in many ways proven to be excellent tools for testing the General theory of Relativity and alternative gravity theories (see Wex (2014) for a recent review). Timing a pulsar in orbit around a companion, provides a unique way of probing the relativistic dynamics and spacetime of such a system. The strictest tests of gravity, in strong field conditions, are expected to come from a pulsar orbiting a black hole. In this sense, a pulsar in a close orbit ($P_{\rm orb}$ < 1 yr) around our nearest supermassive black hole candidate, Sagittarius A* - at a distance of ~8.3 kpc in the Galactic centre (Gillessen et al. 2009a) - would be the ideal tool. Given the size of the orbit and the relativistic effects associated with it, even a slowly spinning pulsar would allow the black hole spacetime to be explored in great detail (Liu et al. 2012). For example, measurement of the frame dragging caused by the rotation of the supermassive black hole, would allow a test of the "cosmic censorship conjecture." The "no-hair theorem" can be tested by measuring the quadrupole moment of the black hole. These are two of the prime examples for the fundamental studies of gravity one could do with a pulsar around Sagittarius A*. As will be shown here, SKA1-MID and ultimately the SKA will provide the opportunity to begin to find and time the pulsars in this extreme environment.

preprint2015arXiv

Pulsar Observations Using the First Station of the Long Wavelength Array and the LWA Pulsar Data Archive

We present initial pulsar results from the first station of the Long Wavelength Array (LWA1) obtained during the commissioning period of LWA1 and early science results. We present detections of periodic emission from 44 previously known pulsars, including 3 millisecond pulsars (MSPs). The effects of the interstellar medium on pulsar emission are significantly enhanced at the low frequencies of the LWA1 band (10--88 MHz), making LWA1 a very sensitive instrument for characterizing changes in dispersion measures (DM) and other effects from the interstellar medium. Pulsars also often have significant evolution in their pulse profile at low frequency and a break in their spectral index. We report DM measurements for 44 pulsars, mean flux density measurements for 36 pulsars, and multi-frequency component spacing and widths for 15 pulsars with more than one profile component. For 27 pulsars, we report spectral index measurements within our frequency range. We also introduce the LWA1 Pulsar Data Archive, which stores reduced data products from LWA1 pulsar observations. Reduced data products for the observations presented here can be found on the archive. Reduced data products from future LWA1 pulsar observations will also be made available through the archive.

preprint2015arXiv

The structure function of Galactic \HI opacity fluctuations on AU scales based on MERLIN, VLA and VLBA data

We use MERLIN, VLA and VLBA observations of Galactic \HI absorption towards 3C~138 to estimate the structure function of the \HI opacity fluctuations at AU scales. Using Monte Carlo simulations, we show that there is likely to be a significant bias in the estimated structure function at signal-to-noise ratios characteristic of our observations, if the structure function is constructed in the manner most commonly used in the literature. We develop a new estimator that is free from this bias and use it to estimate the true underlying structure function slope on length scales ranging $5$ to $40$~AU. From a power law fit to the structure function, we derive a slope of $0.81^{+0.14}_{-0.13}$, i.e. similar to the value observed at parsec scales. The estimated upper limit for the amplitude of the structure function is also consistent with the measurements carried out at parsec scales. Our measurements are hence consistent with the \HI opacity fluctuation in the Galaxy being characterized by a power law structure function over length scales that span six orders of magnitude. This result implies that the dissipation scale has to be smaller than a few AU if the fluctuations are produced by turbulence. This inferred smaller dissipation scale implies that the dissipation occurs either in (i) regions with densities $\gtrsim 10^3 $cm$^-3$ (i.e. similar to that inferred for "tiny scale" atomic clouds or (ii) regions with a mix of ionized and atomic gas (i.e. the observed structure in the atomic gas has a magneto-hydrodynamic origin).

preprint2015arXiv

The Virtual Astronomical Observatory: Re-engineering Access to Astronomical Data

The U.S. Virtual Astronomical Observatory was a software infrastructure and development project designed both to begin the establishment of an operational Virtual Observatory (VO) and to provide the U.S. coordination with the international VO effort. The concept of the VO is to provide the means by which an astronomer is able to discover, access, and process data seamlessly, regardless of its physical location. This paper describes the origins of the VAO, including the predecessor efforts within the U.S. National Virtual Observatory, and summarizes its main accomplishments. These accomplishments include the development of both scripting toolkits that allow scientists to incorporate VO data directly into their reduction and analysis environments and high-level science applications for data discovery, integration, analysis, and catalog cross-comparison. Working with the international community, and based on the experience from the software development, the VAO was a major contributor to international standards within the International Virtual Observatory Alliance. The VAO also demonstrated how an operational virtual observatory could be deployed, providing a robust operational environment in which VO services worldwide were routinely checked for aliveness and compliance with international standards. Finally, the VAO engaged in community outreach, developing a comprehensive web site with on-line tutorials, announcements, links to both U.S. and internationally developed tools and services, and exhibits and hands-on training .... All digital products of the VAO Project, including software, documentation, and tutorials, are stored in a repository for community access. The enduring legacy of the VAO is an increasing expectation that new telescopes and facilities incorporate VO capabilities during the design of their data management systems.

preprint2014arXiv

A 24-Hour Global Campaign To Assess Precision Timing of the Millisecond Pulsar J1713+0747

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

preprint2014arXiv

Chasing the Identification of ASCA Galactic Objects (ChIcAGO) - An X-ray Survey of Unidentified Sources in the Galactic Plane I: Source Sample and Initial Results

We present the Chasing the Identification of ASCA Galactic Objects (ChIcAGO) survey, which is designed to identify the unknown X-ray sources discovered during the ASCA Galactic Plane Survey (AGPS). Little is known about most of the AGPS sources, especially those that emit primarily in hard X-rays (2-10 keV) within the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 X-ray flux range. In ChIcAGO, the subarcsecond localization capabilities of Chandra have been combined with a detailed multi-wavelength follow-up program, with the ultimate goal of classifying the >100 unidentified sources in the AGPS. Overall to date, 93 unidentified AGPS sources have been observed with Chandra as part of the ChIcAGO survey. A total of 253 X-ray point sources have been detected in these Chandra observations within 3' of the original ASCA positions. We have identified infrared and optical counterparts to the majority of these sources, using both new observations and catalogs from existing Galactic plane surveys. X-ray and infrared population statistics for the X-ray point sources detected in the Chandra observations reveal that the primary populations of Galactic plane X-ray sources that emit in the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 flux range are active stellar coronae, massive stars with strong stellar winds that are possibly in colliding-wind binaries, X-ray binaries, and magnetars. There is also a fifth population that is still unidentified but, based on its X-ray and infrared properties, likely comprise partly of Galactic sources and partly of active galactic nuclei.

preprint2014arXiv

NANOGrav Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries in Circular Orbits

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

preprint2014arXiv

The Very Large Array Low-frequency Sky Survey Redux (VLSSr)

We present the results of a recent re-reduction of the data from the Very Large Array (VLA) Low-frequency Sky Survey (VLSS). We used the VLSS catalog as a sky model to correct the ionospheric distortions in the data and create a new set of sky maps and corresponding catalog at 73.8 MHz. The VLSS Redux (VLSSr) has a resolution of 75 arcsec, and an average map RMS noise level of $σ\sim0.1$ Jy beam$^{-1}$. The clean bias is $0.66\timesσ$, and the theoretical largest angular size is 36 arcmin. Six previously un-imaged fields are included in the VLSSr, which has an unbroken sky coverage over 9.3 sr above an irregular southern boundary. The final catalog includes 92,964 sources. The VLSSr improves upon the original VLSS in a number of areas including imaging of large sources, image sensitivity, and clean bias; however the most critical improvement is the replacement of an inaccurate primary beam correction which caused source flux errors which vary as a function of radius to nearest pointing center in the VLSS.

preprint2013arXiv

An Asteroid Belt Interpretation for the Timing Variations of the Millisecond Pulsar B1937+21

Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.

preprint2013arXiv

Detection and Flux Density Measurements of the Millisecond Pulsar J2145-0750 below 100 MHz

We present flux density measurements and pulse profiles for the millisecond pulsar PSR J2145-0750 spanning 37 to 81 MHz using data obtained from the first station of the Long Wavelength Array. These measurements represent the lowest frequency detection of pulsed emission from a millisecond pulsar to date. We find that the pulse profile is similar to that observed at 102 MHz. We also find that the flux density spectrum between ~40 MHz to 5 GHz is suggestive of a break and may be better fit by a model that includes spectral curvature with a rollover around 730 MHz rather than a single power law.

preprint2013arXiv

Observations of Crab Giant Pulses in 20-84 MHz using LWA1

We report the detection and observed characteristics of giant pulses from the Crab Nebula pulsar (B0531+21) in four frequency bands covering 20-84 MHz using the recently-completed Long Wavelength Array Station 1 (LWA1) radio telescope. In 10 hours of observations distributed over a 72-day period in Fall of 2012, 33 giant pulses having peak flux densities between 400 Jy and 2000 Jy were detected. Twenty-two of these pulses were detected simultaneously in channels of 16 MHz bandwidth centered at 44 MHz, 60 MHz, and 76 MHz, including one pulse which was also detected in a channel centered at 28 MHz. We quantify statistics of pulse amplitude and pulse shape characteristics, including pulse broadening. Amplitude statistics are consistent with expectations based on extrapolations from previous work at higher and lower frequencies. Pulse broadening is found to be relatively high, but not significantly greater than expected. We present procedures that have been found to be effective for observing giant pulses in this frequency range.

preprint2012arXiv

A New Technique for Spectral Analysis of Ionospheric TEC Fluctuations Observed with the Very Large Array VHF System: From QP Echoes to MSTIDs

We have used a relatively long, contiguous VHF observation of a bright cosmic radio source (Cygnus A) with the Very Large Array (VLA) through the nighttime, midlatitude ionosphere to demonstrate the phenomena observable with this instrument. In a companion paper, we showed that the VLA can detect fluctuations in total electron content (TEC) with amplitudes of <0.001 TECU and can measure TEC gradients with a precision of about 0.0002 TECU/km. We detail two complementary techniques for producing spectral analysis of these TEC gradient measurements. The first is able to track individual waves with wavelengths of about half the size of the array (~20 km) or more. This technique was successful in detecting and characterizing many medium-scale traveling ionospheric disturbances (MSTIDs) seen intermittently throughout the night and has been partially validated using concurrent GPS measurements. Smaller waves are also seen with this technique at nearly all times, many of which move in similar directions as the detected MSTIDs. The second technique allows for the detection and statistical description of the properties of groups of waves moving in similar directions with wavelengths as small as 5 km. Combining the results of both spectral techniques, we found a class of intermediate and small scale waves which are likely the quasi-periodic (QP) echoes that have been observed to occur within sporadic-E (Es) layers. We find two distinct populations of these waves. The members of one population are coincident in time with MSTIDs and are consistent with being generated within Es layers by the E-F coupling instability. The other population seems more influenced by the neutral wind, similar to the predominant types of QP echoes found by the Sporadic-E Experiments over Kyushu (Fukao et al. 1998; Yamamoto et al. 2005).

preprint2012arXiv

Discovery of a Meter-Wavelength Radio Transient in the SWIRE Deep Field: 1046+59

We report the results of a low frequency radio variability and slow transient search using archival observations from the Very Long Array. We selected six 325 MHz radio observations from the spring of 2006, each centered on the Spitzer-Space-Telescope Wide-area Infrared Extragalactic Survey (SWIRE) Deep Field: 1046+59. Observations were spaced between 1 day to 3 months, with a typical single-epoch peak flux sensitivity below 0.2 \mjb near the field pointing center. We describe the observation parameters, data post-processing, and search methodology used to identify variable and transient emission. Our search revealed multiple variable sources and the presence of one, day-scale transient event with no apparent astronomical counterpart. This detection implies a transient rate of 1$\pm$1 event per 6.5 $°^2$ per 72 observing hours in the direction of 1046+59 and an isotropic transient surface density $Σ= 0.12 °^{-2}$ at 95% confidence for sources with average peak flux density higher than 2.1 mJy over 12 hr.

preprint2012arXiv

High-precision Measurements of Ionospheric TEC Gradients with the Very Large Array VHF System

We have used a relatively long, contiguous VHF observation of a bright cosmic radio source (Cygnus A) with the Very Large Array (VLA) to demonstrate the capability of this instrument to study the ionosphere. This interferometer, and others like it, can observe ionospheric total electron content (TEC) fluctuations on a much wider range of scales than is possible with many other instruments. We have shown that with a bright source, the VLA can measure differential TEC values between pairs of antennas (delta-TEC) with an precision of 0.0003 TECU. Here, we detail the data reduction and processing techniques used to achieve this level of precision. In addition, we demonstrate techniques for exploiting these high-precision delta-TEC measurements to compute the TEC gradient observed by the array as well as small-scale fluctuations within the TEC gradient surface. A companion paper details specialized spectral analysis techniques used to characterize the properties of wave-like fluctuations within this data.

preprint2012arXiv

Multiwavelength Constraints on Pulsar Populations in the Galactic Center

The detection of radio pulsars within the central few parsecs of the Galaxy would provide a unique probe of the gravitational and magneto-ionic environments in the Galactic Center (GC) and, if close enough to Sgr A*, precise tests of general relativity in the strong-field regime. While it is difficult to find pulsars at radio wavelengths because of interstellar scattering, the payoff from detailed timing of pulsars in the GC warrants a concerted effort. To motivate pulsar surveys and help define search parameters for them, we constrain the pulsar number and spatial distribution using a wide range of multiwavelength measurements. These include the five known radio pulsars within 15 arcmin of Sgr A*, radio and gamma-ray measurements of diffuse emission, non-detections in high frequency pulsar surveys of the central parsec, a catalog of radio point sources from an imaging survey, infrared observations of massive star populations in the central few parsecs, candidate pulsar wind nebulae in the inner 20 pc and estimates of the core-collapse supernova rate based on X-ray measurements. We find that under current observational constraints, the inner parsec of the Galaxy could harbor as many as ~10^3 active radio pulsars that are beamed towards Earth. Such a large population would distort the low-frequency measurements of both the intrinsic spectrum of Sgr A* and the free-free absorption along the line of sight of Sgr A*.

preprint2012arXiv

Prospects for probing strong gravity with a pulsar-black hole system

The discovery of a pulsar (PSR) in orbit around a black hole (BH) is expected to provide a superb new probe of relativistic gravity and BH properties. Apart from a precise mass measurement for the BH, one could expect a clean verification of the dragging of space-time caused by the BH spin. In order to measure the quadrupole moment of the BH for testing the no-hair theorem of general relativity (GR), one has to hope for a sufficiently massive BH. In this respect, a PSR orbiting the super-massive BH in the center of our Galaxy would be the ultimate laboratory for gravity tests with PSRs. But even for gravity theories that predict the same properties for BHs as GR, a PSR-BH system would constitute an excellent test system, due to the high grade of asymmetry in the strong field properties of these two components. Here we highlight some of the potential gravity tests that one could expect from different PSR-BH systems, utilizing present and future radio telescopes, like FAST and SKA.

preprint2012arXiv

Tiny scale opacity fluctuations from VLBA, MERLIN and VLA observations of HI absorption toward 3C 138

The structure function of opacity fluctuations is a useful statistical tool to study tiny scale structures of neutral hydrogen. Here we present high resolution observation of HI absorption towards 3C 138, and estimate the structure function of opacity fluctuations from the combined VLA, MERLIN and VLBA data. The angular scales probed in this work are ~ 10-200 milliarcsec (about 5-100 AU). The structure function in this range is found to be well represented by a power law S_tau(x) ~ x^{beta} with index beta ~ 0.33 +/- 0.07 corresponding to a power spectrum P_tau(U) ~ U^{-2.33}. This is slightly shallower than the earlier reported power law index of ~ 2.5-3.0 at ~ 1000 AU to few pc scales. The amplitude of the derived structure function is a factor of ~ 20-60 times higher than the extrapolated amplitude from observation of Cas A at larger scales. On the other hand, extrapolating the AU scale structure function for 3C 138 predicts the observed structure function for Cas A at the pc scale correctly. These results clearly establish that the atomic gas has significantly more structures in AU scales than expected from earlier pc scale observations. Some plausible reasons are identified and discussed here to explain these results. The observational evidence of a shallower slope and the presence of rich small scale structures may have implications for the current understanding of the interstellar turbulence.

preprint2011arXiv

An Arecibo Search for Pulsars and Transient Sources in M33

We report on a systematic and sensitive search for pulsars and transient sources in the nearby spiral galaxy M33, conducted at 1.4 GHz with the Arecibo telescope's seven-beam receiver system, ALFA. Data were searched for both periodic and aperiodic sources, up to 1000 pc cm^{-3} in dispersion measure and on timescales from $\sim$50 \us to several seconds. The galaxy was sampled with 12 ALFA pointings, or 84 pixels in total, each of which was searched for 2-3 hr. We describe the observations, search methodologies and analysis strategies applicable to multibeam systems, and comment on the data quality and statistics of spurious events that arise due to radio frequency interference. While these searches have not led to any conclusive signals of periodic or transient nature that originate in the galaxy, they illustrate some of the underlying challenges and difficulties in such searches and the efficacy of simultaneous multiple beams in the analysis of search output. The implied limits are $\la$5 uJy Mpc^2 in luminosity (at 1400 MHz) for periodic sources in M33 with duty cycles $\la$5%. For short-duration transient signals (with pulse widths $\la$100 \us), the limiting peak flux density is 100 mJy, which would correspond to a 5-$σ$ detection of bright giant pulses ($\sim$20 kJy) from Crab-like pulsars if located at the distance of M33. We discuss the implications of our null results for various source populations within the galaxy and comment on the future prospects to conduct even more sensitive searches with the upcoming next-generation instruments including the Square Kilometer Array and its pathfinders.

preprint2011arXiv

Probing the First Stars and Black Holes in the Early Universe with the Dark Ages Radio Explorer (DARE)

A concept for a new space-based cosmology mission called the Dark Ages Radio Explore (DARE) is presented in this paper. DARE's science objectives include (1) When did the first stars form? (2) When did the first accreting black holes form? (3) When did Reionization begin? (4) What surprises does the end of the Dark Ages hold (e.g., Dark Matter decay)? DARE will use the highly-redshifted hyperfine 21-cm transition from neutral hydrogen to track the formation of the first luminous objects by their impact on the intergalactic medium during the end of the Dark Ages and during Cosmic Dawn (redshifts z=11-35). It will measure the sky-averaged spin temperature of neutral hydrogen at the unexplored epoch 80-420 million years after the Big Bang, providing the first evidence of the earliest stars and galaxies to illuminate the cosmos and testing our models of galaxy formation. DARE's approach is to measure the expected spectral features in the sky-averaged, redshifted 21-cm signal over a radio bandpass of 40-120 MHz. DARE orbits the Moon for a mission lifetime of 3 years and takes data above the lunar farside, the only location in the inner solar system proven to be free of human-generated radio frequency interference and any significant ionosphere. The science instrument is composed of a three-element radiometer, including electrically-short, tapered, bi-conical dipole antennas, a receiver, and a digital spectrometer. The smooth frequency response of the antennas and the differential spectral calibration approach using a Markov Chain Monte Carlo technique will be applied to detect the weak cosmic 21-cm signal in the presence of the intense solar system and Galactic foreground emissions.

preprint2011arXiv

Prospects for Probing the Spacetime of Sgr A* with Pulsars

The discovery of radio pulsars in compact orbits around Sgr A* would allow an unprecedented and detailed investigation of the spacetime of the supermassive black hole. This paper shows that pulsar timing, including that of a single pulsar, has the potential to provide novel tests of general relativity, in particular its cosmic censorship conjecture and no-hair theorem for rotating black holes. These experiments can be performed by timing observations with 100 micro-second precision, achievable with the Square Kilometre Array for a normal pulsar at frequency above 15 GHz. Based on the standard pulsar timing technique, we develop a method that allows the determination of the mass, spin, and quadrupole moment of Sgr A*, and provides a consistent covariance analysis of the measurement errors. Furthermore, we test this method in detailed mock data simulations. It seems likely that only for orbital periods below ~0.3 yr is there the possibility of having negligible external perturbations. For such orbits we expect a ~10^-3 test of the frame dragging and a ~10^-2 test of the no-hair theorem within 5 years, if Sgr A* is spinning rapidly. Our method is also capable of identifying perturbations caused by distributed mass around Sgr A*, thus providing high confidence in these gravity tests. Our analysis is not affected by uncertainties in our knowledge of the distance to the Galactic center, R0. A combination of pulsar timing with the astrometric results of stellar orbits would greatly improve the measurement precision of R0.

preprint2011arXiv

PSRPI: A large VLBA pulsar astrometry program

Obtaining pulsar parallaxes via relative astrometry (also known as differential astrometry) yields distances and transverse velocities that can be used to probe properties of the pulsar population and the interstellar medium. Large programs are essential to obtain the sample sizes necessary for these population studies, but they must be efficiently conducted to avoid requiring an infeasible amount of observing time. This paper describes the PSRPI astrometric program, including the use of new features in the DiFX software correlator to efficiently locate calibrator sources, selection and observing strategies for a sample of 60 pulsars, initial results, and likely science outcomes. Potential applications of high-precision relative astrometry to measure source structure evolution in defining sources of the International Celestial Referent Frame are also discussed.

preprint2010arXiv

Identification of a Population of X-ray Emitting Massive Stars in the Galactic Plane

We present X-ray, infrared, optical and radio observations of four previously unidentified Galactic plane X-ray sources, AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. Detection of each source with the Chandra X-ray Observatory has provided sub-arcsecond localizations, which we use to identify bright infrared counterparts to all four objects. Infrared and optical spectroscopy of these counterparts demonstrate that all four X-ray sources are extremely massive stars, with spectral classifications Ofpe/WN9 (AX J163252-4746), WN7 (AX J184738-0156 = WR121a), WN7-8h (AX J144701-5919) and OIf+ (AX J144547-5931). AX J163252-4746 and AX J184738-0156 are both luminous, hard, X-ray emitters with strong Fe XXV emission lines in their X-ray spectra at ~6.7 keV. The multi-wavelength properties of AX J163252-4746 and AX J184738-0156 are not consistent with isolated massive stars or accretion onto a compact companion; we conclude that their X-ray emission is most likely generated in a colliding-wind binary system. For both AX J144701-5919 and AX J144547-5931, the X-ray emission is an order of magnitude less luminous and with a softer spectrum. These properties are consistent with a colliding-wind binary interpretation for these two sources also, but other mechanisms for the generation of X-rays cannot be excluded. There are many other as yet unidentified X-ray sources in the Galactic plane, with X-ray properties similar to those seen for AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. This may indicate a substantial population of X-ray-emitting massive stars and colliding-wind binaries in the Milky Way.

preprint2010arXiv

Pulsars with the Australian Square Kilometre Array Pathfinder

The Australian Square Kilometre Array Pathfinder (ASKAP) is a 36-element array with a 30-square-degree field of view being built at the proposed SKA site in Western Australia. We are conducting a Design Study for pulsar observations with ASKAP, planning both timing and search observations. We provide an overview of the ASKAP telescope and an update on pulsar-related progress.

preprint2010arXiv

Radio Observations of HD 80606 Near Planetary Periastron

This paper reports Very Large Array observations at 325 and 1425 MHz (90cm and 20cm) during and near the periastron passage of HD 80606b on 2007 November 20. We obtain flux density limits (3-sigma) of 1.7 mJy and 48 microJy at 325 and 1425 MHz, respectively, equivalent to planetary luminosity limits of 2.3 x 10^{24} erg/s and 2.7 x 10^{23} erg/s. These are well above the Jovian value (at 40 MHz) of 2 x 10^{18} erg/s. The motivation for these observations was that the planetary magnetospheric emission is driven by a stellar wind-planetary magnetosphere interaction so that the planetary luminosity would be elevated. Near periastron, HD 80606b might be as much as 3000 times more luminous than Jupiter. Recent transit observations of HD 80606b provide stringent constraints on the planetary mass and radius, and, because of the planet's highly eccentric orbit, its rotation period is likely to be "pseudo-synchronized" to its orbital period, allowing a robust estimate of the former. We are able to make robust estimates of the emission frequency of the planetary magnetospheric emission and find it to be around 60--90 MHz. We compare HD 80606b to other high-eccentricity systems and assess the detection possibilities for both near-term and more distant future systems. Of the known high eccentricity planets, only HD 80606b is likely to be detectable, as HD 20782B b and HD 4113b are both likely to have weaker magnetic field strengths. Both the forthcoming "EVLA low band" system and the Low Frequency Array may be able to improve upon our limits for HD 80606b, and do so at a more optimum frequency. If the low-frequency component of the Square Kilometre Array (SKA-lo) and a future lunar radio array are able to approach their thermal noise limits, they should be able to detect an HD 80606b-like planet, unless the planet's luminosity increases by substantially less than a factor of 3000.

preprint2010arXiv

Surveying the Dynamic Radio Sky with the Long Wavelength Demonstrator Array

This paper presents a search for radio transients at a frequency of 73.8 MHz (4 m wavelength) using the all-sky imaging capabilities of the Long Wavelength Demonstrator Array (LWDA). The LWDA was a 16-dipole phased array telescope, located on the site of the Very Large Array in New Mexico. The field of view of the individual dipoles was essentially the entire sky, and the number of dipoles was sufficiently small that a simple software correlator could be used to make all-sky images. From 2006 October to 2007 February, we conducted an all-sky transient search program, acquiring a total of 106 hr of data; the time sampling varied, being 5 minutes at the start of the program and improving to 2 minutes by the end of the program. We were able to detect solar flares, and in a special-purpose mode, radio reflections from ionized meteor trails during the 2006 Leonid meteor shower. We detected no transients originating outside of the solar system above a flux density limit of 500 Jy, equivalent to a limit of no more than about 10^{-2} events/yr/deg^2, having a pulse energy density >~ 1.5 x 10^{-20} J/m^2/Hz at 73.8 MHz for pulse widths of about 300 s. This event rate is comparable to that determined from previous all-sky transient searches, but at a lower frequency than most previous all-sky searches. We believe that the LWDA illustrates how an all-sky imaging mode could be a useful operational model for low-frequency instruments such as the Low Frequency Array, the Long Wavelength Array station, the low-frequency component of the Square Kilometre Array, and potentially the Lunar Radio Array.

preprint2010arXiv

The Commensal Real-time ASKAP Fast Transients (CRAFT) survey

We are developing a purely commensal survey experiment for fast (<5s) transient radio sources. Short-timescale transients are associated with the most energetic and brightest single events in the Universe. Our objective is to cover the enormous volume of transients parameter space made available by ASKAP, with an unprecedented combination of sensitivity and field of view. Fast timescale transients open new vistas on the physics of high brightness temperature emission, extreme states of matter and the physics of strong gravitational fields. In addition, the detection of extragalactic objects affords us an entirely new and extremely sensitive probe on the huge reservoir of baryons present in the IGM. We outline here our approach to the considerable challenge involved in detecting fast transients, particularly the development of hardware fast enough to dedisperse and search the ASKAP data stream at or near real-time rates. Through CRAFT, ASKAP will provide the testbed of many of the key technologies and survey modes proposed for high time resolution science with the SKA.

preprint2009arXiv

Discovery of Three Pulsars from a Galactic Center Pulsar Population

We report the discovery of three pulsars whose large dispersion measures and angular proximity to \sgr indicate the existence of a Galactic center population of neutron stars. The relatively long periods (0.98 to 1.48 s) most likely reflect strong selection against short-period pulsars from radio-wave scattering at the observation frequency of 2 GHz used in our survey with the Green Bank Telescope. One object (PSR J1746-2850I) has a characteristic spindown age of only 13 kyr along with a high surface magnetic field $\sim 4\times 10^{13}$ G. It and a second object found in the same telescope pointing, PSR J1746-2850II (which has the highest known dispersion measure among pulsars), may have originated from recent star formation in the Arches or Quintuplet clusters given their angular locations. Along with a third object, PSR J1745-2910, and two similar high-dispersion, long-period pulsars reported by Johnston et al. (2006), the five objects found so far are 10 to 15 arc min from \sgr, consistent with there being a large pulsar population in the Galactic center, most of whose members are undetectable in relatively low-frequency surveys because of pulse broadening from the same scattering volume that angularly broadens \sgr and OH/IR masers.

preprint2009arXiv

Precision Astrometry with the Very Long Baseline Array: Parallaxes and Proper Motions for 14 Pulsars

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames.

preprint2005arXiv

Pulsars as Tools for Fundamental Physics and Astrophysics

The sheer number of pulsars discovered by the SKA, in combination with the exceptional timing precision it can provide, will revolutionize the field of pulsar astrophysics. The SKA will provide a complete census of pulsars in both the Galaxy and in Galactic globular clusters that can be used to provide a detailed map of the electron density and magnetic fields, the dynamics of the systems, and their evolutionary histories. This complete census will provide examples of nearly every possible outcome of the evolution of massive stars, including the discovery of very exotic systems such as pulsar black-hole systems and sub-millisecond pulsars, if they exist. These exotic systems will allow unique tests of the strong field limit of relativistic gravity and the equation of state at extreme densities. Masses of pulsars and their binary companions -- planets, white dwarfs, other neutron stars, and black holes -- will be determined to $\sim 1$% for hundreds of objects. With the SKA we can discover and time highly-stable millisecond pulsars that comprise a pulsar-timing array for the detection of low-frequency gravitational waves. The SKA will also provide partial censuses of nearby galaxies through periodicity and single-pulse detections, yielding important information on the intergalactic medium.

preprint2005arXiv

RFI Identification and Mitigation Using Simultaneous Dual Station Observations

RFI mitigation is a critically important issue in radio astronomy using existing instruments as well as in the development of next-generation radio telescopes, such as the Square Kilometer Array (SKA). Most designs for the SKA involve multiple stations with spacings of up to a few thousands of kilometers and thus can exploit the drastically different RFI environments at different stations. As demonstrator observations and analysis for SKA-like instruments, and to develop RFI mitigation schemes that will be useful in the near term, we recently conducted simultaneous observations with Arecibo Observatory and the Green Bank Telescope (GBT). The observations were aimed at diagnosing RFI and using the mostly uncorrelated RFI between the two sites to excise RFI from several generic kinds of measurements such as giant pulses from Crab-like pulsars and weak HI emission from galaxies in bands heavily contaminated by RFI. This paper presents observations, analysis, and RFI identification and excision procedures that are effective for both time series and spectroscopy applications using multi-station data.