Researcher profile

M. A. McLaughlin

M. A. McLaughlin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

19 published item(s)

preprint2023arXiv

Single-pulse studies of three millisecond pulsars

Single-pulse studies are important to understand the pulsar emission mechanism and the noise floor in precision timing. We study total intensity and polarimetry properties of three bright millisecond pulsars - PSRs J1022+1001, J1713+0747, and B1855+09 - that have detectable single pulses at multiple frequencies. We report for the first time the detection of single pulses from PSRs J1022+1001 and J1713+0747 at 4.5 GHz. In addition, for those two pulsars the fraction of linear polarization in the average profile is significantly reduced at 4.5 GHz, compared to 1.38 GHz, which could support the expected deviation from a dipolar field closer to the pulsar surface. There is a hint of orthogonal modes in the single pulses of PSR J1713+0747. More sensitive multi-frequency observations may be useful to confirm these findings. The jitter noise contributions at 1.38 GHz, scaled to one hour, for PSRs J1022+1001, J1713+0747 and B1855+09 are ~135 ns, ~45 ns, and ~60 ns respectively and are consistent with previous studies. We also show that selective bright-pulse timing of PSR J1022+1001 yields improved root-mean-square residuals of ~22 $μ$s, which is a factor of ~3 better than timing using single pulses alone.

preprint2022arXiv

News and views regarding PSR J1757-1854, a highly-relativistic binary pulsar

We provide an update on the ongoing monitoring and study of the highly-relativistic double neutron star binary system PSR J1757-1854, a 21.5-ms pulsar in a highly eccentric, 4.4-hour orbit. The extreme nature of this pulsar's orbit allows it to probe a parameter space largely unexplored by other relativistic binary pulsars. For example, it displays one of the highest gravitational wave (GW) luminosities of any known binary pulsar, as well as the highest rate or orbital decay due to GW damping. PSR J1757-1854 is also notable in that it is an excellent candidate for exploring new tests of General Relativity and other gravitational theories, with possible measurements of both Lense-Thirring precession and relativistic orbital deformation (through the post-Keplerian parameter $δ_θ$) anticipated within the next 3-5 years. Here we present a summary of the latest interim results from the ongoing monitoring of this pulsar as part of an international, multi-telescope campaign. This includes an update of the pulsar's long-term timing and post-Keplerian parameters, new constraints on the pulsar's proper motion and corresponding Shklovskii kinematic correction, and new limits on the pulsar's geodetic precession as determined by monitoring for secular changes in the pulse profile. We also highlight prospects for future work, including an updated timeline on new relativistic tests following the introduction of MeerKAT observations.

preprint2022arXiv

Single Pulse Dispersion Measure of the Crab Pulsar

We investigate the use of bright single pulses from the Crab pulsar to determine separately the dispersion measure (DM) for the Main Pulse and Interpulse components. We develop two approaches using cross correlation functions (CCFs). The first method computes the CCF of the total intensity of each of 64 frequency channels with a reference channel and converts the time lag of maximum correlation into a DM. The second method separately computes the CCF between every pair of channels for each individual bright pulse and extracts an average DM from the distribution of all channel-pair DMs. Both methods allow the determination of the DM with a relative uncertainty of better than 10^-5 and provide robust estimates for the uncertainty of the best-fit value. We find differences in DM between the Main Pulse, the Low Frequency Interpulse, and the High Frequency Interpulse using both methods in a frequency range from 4 to 6 GHz. Earlier observations of the High Frequency Interpulse carried out by Hankins et al. (2016) resulted in DM_HFIP-DM_MP of 0.010 +- 0.016 pc cm^-3. Our results indicate a DM_HFIP-DM_MP of 0.0127 +- 0.0011 pc cm^-3 (with DM_comp being the DM value of the respective emission component), confirming earlier results with an independent method. During our studies we also find a relation between the brightness of single pulses in the High Frequency Interpulse and their DM. We also discuss the application of the developed methods on the identification of substructures in the case of Fast Radio Bursts.

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.

preprint2021arXiv

Commissioning the HI Observing Mode of the Beamformer for the Cryogenically Cooled Focal L-band Array for the GBT (FLAG)

We present the results of commissioning observations for a new digital beamforming back end for the Focal plane L-band Array for the Robert C. Byrd Green Bank Telescope (FLAG), a cryogenically cooled Phased Array Feed (PAF) with the lowest measured T_sys/eta of any PAF outfitted on a radio telescope to date. We describe the custom software used to apply beamforming weights to the raw element covariances to create research quality spectral line images for the new fine-channel mode, study the stability of the beam weights over time, characterize FLAG's sensitivity over a frequency range of 150 MHz, and compare the measured noise properties and observed distribution of neutral hydrogen emission from several extragalactic and Galactic sources with data obtained with the current single-pixel L-band receiver. These commissioning runs establish FLAG as the preeminent PAF receiver currently available for spectral line observations on the world's major radio telescopes.

preprint2021arXiv

Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary

We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 12 nulling and two intermittent pulsars. We detected eight glitches in five pulsars. Among them is PSR J1939+2609, an apparently old pulsar (characteristic age $\sim$1 Gy), and PSR J1954+2529, which likely belongs to a newly-emerging class of binary pulsars. The latter is the only pulsar among the 72 that is clearly not isolated: a non-recycled neutron star with a 931-ms spin period in an eccentric ($e\,=\,0.114$) wide ($P_b\,=\,82.7\,$d) orbit with a companion of undetermined nature having a minimum mass of $\sim0.6\,M_{\odot}$. Since operations at Arecibo ceased in 2020 August, we give a final tally of PALFA sky coverage, and compare its 207 pulsar discoveries to the known population. On average, they are 50% more distant than other Galactic plane radio pulsars; PALFA millisecond pulsars (MSP) have twice the dispersion measure per unit spin period than the known population of MSP in the Plane. The four intermittent pulsars discovered by PALFA more than double the population of such objects, which should help to improve our understanding of pulsar magnetosphere physics. The statistics for these, RRATS, and nulling pulsars suggest that there are many more of these objects in the Galaxy than was previously thought.

preprint2020arXiv

Asymmetric mass ratios for bright double neutron-star mergers

The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentially explained by a progenitor system that is asymmetric in terms of the stellar component masses, i.e. with a mass ratio $q$ of 0.7-0.8, the known Galactic population of merging double neutron star (DNS) systems (i.e. those that will coalesce within billions of years or less) has, until now, only consisted of nearly equal-mass ($q > 0.9$) binaries. PSR J1913+1102 is a DNS system in a 5-hour, low-eccentricity ($e = 0.09$) orbit, implying an orbital separation of 1.8 solar radii, with the two neutron stars predicted to coalesce in 470 million years due to gravitational-wave emission. Here we report that the masses of the two neutron stars, as measured by a dedicated pulsar timing campaign, are $1.62 \pm 0.03$ and $1.27 \pm 0.03$ solar masses for the pulsar and companion neutron star, respectively; with a measured mass ratio $q = 0.78 \pm 0.03$, it is the most asymmetric DNS among known merging systems. Based on this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30% (90% confidence) of the total population of merging DNS binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

preprint2020arXiv

Initial results from a realtime FRB search with the GBT

We present the data analysis pipeline, commissioning observations and initial results from the GREENBURST fast radio burst (FRB) detection system on the Robert C. Byrd Green Bank Telescope (GBT) previously described by Surnis et al. which uses the 21~cm receiver observing commensally with other projects. The pipeline makes use of a state-of-the-art deep learning classifier to winnow down the very large number of false positive single-pulse candidates that mostly result from radio frequency interference. In our observations totalling 156.5 days so far, we have detected individual pulses from 20 known radio pulsars which provide an excellent verification of the system performance. We also demonstrate, through blind injection analyses, that our pipeline is complete down to a signal-to-noise threshold of 12. Depending on the observing mode, this translates to peak flux sensitivities in the range 0.14--0.89~Jy. Although no FRBs have been detected to date, we have used our results to update the analysis of Lawrence et al. to constrain the FRB all-sky rate to be $1140^{+200}_{-180}$ per day above a peak flux density of 1~Jy. We also constrain the source count index $α=0.83\pm0.06$ which indicates that the source count distribution is substantially flatter than expected from a Euclidean distribution of standard candles (where $α=1.5$). We discuss this result in the context of the FRB redshift and luminosity distributions. Finally, we make predictions for detection rates with GREENBURST, as well as other ongoing and planned FRB experiments.

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

The MeerKAT Telescope as a Pulsar Facility: System verification and early science results from MeerTime

We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

preprint2020arXiv

The NANOGrav 11-year Data Set: Constraints on Planetary Masses Around 45 Millisecond Pulsars

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3).

preprint2019arXiv

A Deep Targeted Search for Fast Radio Bursts from the Sites of Low-Redshift Short Gamma-Ray Bursts

Some short gamma-ray bursts (SGRBs) are thought to be caused by the mergers of binary neutron stars which may sometimes produce massive neutron star remnants capable of producing extragalactic fast radio bursts (FRBs). We conducted a deep search for FRBs from the sites of six low-redshift SGRBs. We collected high time- and frequency-resolution data from each of the sites for 10 hours using the 2 GHz receiver of the Green Bank Telescope. Two of the SGRB sites we targeted were visible with the Arecibo Radio Telescope with which we conducted an additional 10 hours of 1.4 GHz observations for each. We searched our data for FRBs using the GPU-optimized dedispersion algorithm $\texttt{heimdall}$ and the machine-learning-based package $\texttt{FETCH}$ (Fast Extragalactic Transient Candidate Hunter). We did not discover any FRBs, but would have detected any with peak flux densities in excess of 87 mJy at the Green Bank Telescope or 21 mJy at Arecibo with a signal-to-noise ratio of at least 10. The isotropic-equivalent energy of any FRBs emitted from these sites in our bands during our observations must not have exceeded a few times $10^{38}$ erg, comparable to some of the lowest energy bursts yet seen from the first known repeating FRB 121102.

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

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 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 &#34;time slicing&#34; 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 &#34;false positive&#34; 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&#39;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.

preprint2010arXiv

Constraining Pulsar Emission Physics through Radio/Gamma-Ray Correlation of Crab Giant Pulses

To constrain the giant pulse (GP) emission mechanism and test the model of Lyutikov (2007) of GP emission, we are carrying out a campaign of simultaneous observations of the Crab pulsar between gamma-rays (Fermi) and radio wavelengths. The correlation between times of arrival of radio GPs and high-energy photons, whether it exists or not, will allow us to choose between different origins of GP emission and further constrain the emission physics. Our foremost goal was testing whether radio GPs are due to changes in the coherence of the radio emission mechanism, variations in the pair creation rate in the pulsar magnetosphere, or changes in the beaming direction. Accomplishing this goal requires an enormous number of simultaneous radio GPs and gamma-photons. Thus, we organized a radio observations campaign using the 42-ft telescope at the Jodrell Bank Observatory (UK), the 140-ft telescope, and the 100-m Robert C. Byrd Green Bank Telescope (GBT) at the Green Bank Observatory (WV). While the observations with the two first ones are ongoing, here we present the preliminary results of 20 hrs of observations with the GBT at the high frequency of 8.9 GHz. These particular observations were aimed to probe the model of GP emission by Lyutikov (2007) which predicts that GPs at frequencies > 4 GHz should be accompanied by gamma-ray photons of energies of 1-100 GeV.

preprint2009arXiv

New Limits on Radio Emission from X-ray Dim Isolated Neutron Stars

We have carried out a search for radio emission at 820 MHz from six X-ray dim isolated neutron stars with the Robert C. Byrd Green Bank Radio Telescope. No transient or pulsed emission was found using fast folding, fast Fourier transform, and single-pulse searches. The corresponding flux limits are about 0.01 mJy for pulsed emission, depending on the integration time for the particular source and assuming a duty cycle of 2%, and 20 mJy for single dispersed pulses. These are the most sensitive limits to date on radio emission from X-ray dim isolated neutron stars. There is no evidence for isolated radio pulses, as seen in a class of neutron stars known as rotating radio transients. Our results imply that either the radio luminosities of these objects are lower than those of any known radio pulsars, or they could simply be long-period nearby radio pulsars with high magnetic fields beaming away from the Earth. To test the latter possibility, we would need around 40 similar sources to provide a 1 sigma probability of at least one of them beaming toward us. We also give a detailed description of our implementation of the Fast Folding Algorithm.