# Posts Tagged timing model

## Recent Postings from timing model

### Simultaneous spectral and reverberation modelling of relativistic reflection in Mrk 335

We present an X-ray spectral and timing model to investigate the broad and variable iron line seen in the high flux state of Mrk 335. The model consists of a variable X-ray source positioned along the rotation axis of the black hole that illuminates the accretion disc producing a back-scattered, ionized reflection spectrum. We compute time lags including full dilution effects and perform simultaneous fitting of the 2-10 keV spectrum and the frequency-dependent time lags of 2.5-4 vs. 4-6.5 keV bands. The best-fitting parameters are consistent with a black hole mass of approximately 1.3 x 10^7 M_sun, disc inclination of 45 degrees and the photon index of the direct continuum of 2.4. The iron abundance is 0.5 and the ionization parameter is 10^3 erg cm / s at the innermost part of the disc and decreases further out. The X-ray source height is very small, approximately 2 r_g. Furthermore, we fit the Fe L lags simultaneously with the 0.3-10 keV spectrum. The key parameters are comparable to those previously obtained. We also report the differences below 2 keV using the xillver and reflionx models which could affect the interpretation of the soft excess. While simultaneously fitting spectroscopic and timing data can break the degeneracy between the source height and the black hole mass, we find that the measurements of the source height and the central mass significantly depend on the ionization state of the disc and are possibly model-dependent.

### 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 time-correlated "red" signals in 10 of the pulsars. 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.

### Generative pulsar timing analysis

A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using shapelets’ – a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA).

### Generative pulsar timing analysis [Replacement]

A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using shapelets’ – a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA).

### Bayesian estimation of non-Gaussianity in pulsar timing analysis [Replacement]

We introduce a method for performing a robust Bayesian analysis of non-Gaussianity present in pulsar timing data, simultaneously with the pulsar timing model, and additional stochastic parameters such as those describing red spin noise and dispersion measure variations. The parameters used to define the presence of non-Gaussianity are zero for Gaussian processes, giving a simple method of defining the strength of non-Gaussian behaviour. We use simulations to show that assuming Gaussian statistics when the noise in the data is drawn from a non-Gaussian distribution can significantly increase the uncertainties associated with the pulsar timing model parameters. We then apply the method to the publicly available 15 year Parkes Pulsar Timing Array data release 1 dataset for the binary pulsar J0437$-$4715. In this analysis we present a significant detection of non-Gaussianity in the uncorrelated non-thermal noise, but we find that it does not yet impact the timing model or stochastic parameter estimates significantly compared to analysis performed assuming Gaussian statistics. The methods presented are, however, shown to be of immediate practical use for current European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) datasets.

### TempoNest: A Bayesian approach to pulsar timing analysis

A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis.

### VLBI astrometry of PSR J2222-0137: a pulsar distance measured to 0.4% accuracy

The binary pulsar J2222-0137 is an enigmatic system containing a partially recycled millisecond pulsar and a companion of unknown nature. Whilst the low eccentricity of the system favors a white dwarf companion, an unusual double neutron star system is also a possibility, and optical observations will be able to distinguish between these possibilities. In order to allow the absolute luminosity (or upper limit) of the companion object to be properly calibrated, we undertook astrometric observations with the Very Long Baseline Array to constrain the system distance via a measurement of annual geometric parallax. With these observations, we measure the parallax of the J2222-0137 system to be 3.742 +0.013 -0.016 milliarcseconds, yielding a distance of 267.3 +1.2 -0.9 pc, and measure the transverse velocity to be 57.1 +0.3 -0.2 km/s. Fixing these parameters in the pulsar timing model made it possible to obtain a measurement of Shapiro delay and hence the system inclination, which shows that the system is nearly edge-on (sin i = 0.9985 +/- 0.0005). Furthermore, we were able to detect the orbital motion of J2222-0137 in our VLBI observations and measure the longitude of ascending node. The VLBI astrometry yields the most accurate distance obtained for a radio pulsar to date, and is furthermore the most accurate parallax for any radio source obtained at "low" radio frequencies (below ~5 GHz, where the ionosphere dominates the error budget). Using the astrometric results, we show the companion to J2222-0137 will be easily detectable in deep optical observations if it is a white dwarf. Finally, we discuss the implications of this measurement for future ultra-high-precision astrometry, in particular in support of pulsar timing arrays.

### The Einstein@Home search for radio pulsars and PSR J2007+2722 discovery

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 192 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using PALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 min. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722′s pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~500 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.

### The Einstein@Home search for radio pulsars and PSR J2007+2722 discovery [Replacement]

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 193 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using Pulsar ALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 minutes. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722′s pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~ 100 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.

### The Einstein@Home search for radio pulsars and PSR J2007+2722

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 192 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using PALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 min. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722′s pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~500 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.

### Coherently dedispersed gated imaging of millisecond pulsars

Motivated by the need for rapid localisation of newly discovered faint millisecond pulsars (MSPs) we have developed a coherently dedispersed gating correlator. This gating correlator accounts for the orbital motions of MSPs in binaries while folding the visibilities with best-fit topocentric rotational model derived from periodicity search in simultaneously generated beamformer output. Unique applications of the gating correlator for sensitive interferometric studies of MSPs are illustrated using the Giant Metrewave Radio Telescope (GMRT) interferometric array. We could unambiguously localise five newly discovered Fermi MSPs in the on-off gated image plane with an accuracy of +-1". Immediate knowledge of such precise position allows the use of sensitive coherent beams of array telescopes for follow-up timing observations, which substantially reduces the use of telescope time (~ 20X for the GMRT). In addition, precise a-priori astrometric position reduces the effect of large covariances in timing fit (with discovery position, pulsar period derivative and unknown binary model), which in-turn accelerates the convergence to initial timing model. For example, while fitting with precise a-priori position (+-1"), timing model converges in about 100 days, accounting the effect of covariance between position and pulsar period derivative. Moreover, such accurate positions allows for rapid identification of pulsar counterpart at other wave-bands. We also report a new methodology of in-beam phase calibration using the on-off gated image of the target pulsar, which provides the optimal sensitivity of the coherent array removing the possible temporal and spacial decoherences.

### The Benefits of VLBI Astrometry to Pulsar Timing Array Searches for Gravitational Radiation [Replacement]

Precision astrometry is an integral component of successful pulsar timing campaigns. Astrometric parameters are commonly derived by fitting them as parameters of a timing model to a series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over several-year spans can yield position measurements with sub-milliarcsecond precision. However, timing-based astrometry can become biased if a pulsar displays any red spin noise, which can be compared to the red noise signal produced by the stochastic gravitational wave background. We investigate how noise of different spectral types is absorbed by timing models, leading to significant estimation errors in the astrometric parameters. We find that commonly used techniques for fitting timing models in the presence of red noise (Cholesky whitening) prevent the absorption of noise into the timing model remarkably well if the time baseline of observations exceeds several years, but are inadequate for dealing with shorter pulsar data sets. Independent of timing, pulsar-optimized very long baseline interferometry (VLBI) is capable of providing position estimates precise to the sub-milliarcsecond levels needed for high-precision timing. We compute a necessary transformation between the International Celestial Reference Frame and pulsar timing frames and quantitatively discuss how the transformation will improve in coming years. We find that incorporating VLBI astrometry into the timing models of pulsars for which only a couple of years of timing data exist will lead to more realistic assessments of red spin noise and could enhance the amplitude of gravitational wave signatures in post-fit timing residuals by factors of 20 or more.

### The Benefits of VLBI Astrometry to Pulsar Timing Array Searches for Gravitational Radiation [Replacement]

Precision astrometry is an integral component of successful pulsar timing campaigns. Astrometric parameters are commonly derived by fitting them as parameters of a timing model to a series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over several-year spans can yield position measurements with sub-milliarcsecond precision. However, timing-based astrometry can become biased if a pulsar displays any red spin noise, which can be compared to the red noise signal produced by the stochastic gravitational wave background. We investigate how noise of different spectral types is absorbed by timing models, leading to significant estimation errors in the astrometric parameters. We find that commonly used techniques for fitting timing models in the presence of red noise (Cholesky whitening) prevent the absorption of noise into the timing model remarkably well if the time baseline of observations exceeds several years, but are inadequate for dealing with shorter pulsar data sets. Independent of timing, pulsar-optimized very long baseline interferometry (VLBI) is capable of providing position estimates precise to the sub-milliarcsecond levels needed for high-precision timing. We compute a necessary transformation between the International Celestial Reference Frame and pulsar timing frames and quantitatively discuss how the transformation will improve in coming years. We find that incorporating VLBI astrometry into the timing models of pulsars for which only a couple of years of timing data exist will lead to more realistic assessments of red spin noise and could enhance the amplitude of gravitational wave signatures in post-fit timing residuals by factors of 20 or more.

### The Benefits of VLBI Astrometry to Pulsar Timing Array Searches for Gravitational Radiation

Precisely measured astrometric parameters are integral to successful pulsar timing campaigns. They are commonly measured by fitting the astrometric parameters of a deterministic timing model to a series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over several-year spans can in this way provide astrometric parameters precise to sub-milliarcsecond levels. However, pulsars do not pulsate in a deterministic fashion. Many display significant amounts of red spin noise. Furthermore, a stochastic background of gravitational waves can lead to red noise-like structure in TOAs. We investigate how noise of different spectral types is absorbed by timing models and leads to significant estimation errors in the astrometric parameters. Independent of timing, very long baseline interferometry (VLBI) is capable of providing sub-milliarcsecond astrometric parameters for pulsars. We find that incorporating VLBI astrometric measurements into the timing models of pulsars for which only a couple of years of timing data exist will lead to more realistic assessments of red spin noise, yield more accurate astrometric parameters, and could enhance the amplitude of certain gravitational wave signatures in post-fit timing residuals by factors of 20 or more.

### Applying Bayesian Inference to the first International Pulsar Timing Array data challenge [Cross-Listing]

This is a very brief summary of the techniques I used to analyze the IPTA challenge 1 data sets. I tried many things, and more failed than succeeded, but in the end I found two approaches that appear to work based on tests done using the open data sets. One approach works directly with the time domain data, and the other works with a specially constructed Fourier transform of the data. The raw data was run through TEMPO2 to produce reduced timing residuals for the analysis. Standard Markov Chain Monte Carlo techniques were used to produce samples from the posterior distribution function for the model parameters. The model parameters include the gravitational wave amplitude and spectral slope, and the white noise amplitude for each pulsar in the array. While red timing noise was only included in Dataset 3, I found that it was necessary to include effective red noise in all the analyses to account for some of the spurious effects introduced by the TEMPO2 timing fit. This added an additional amplitude and slope parameter for each pulsar, so my overall model for the 36 pulsars residuals has 110 parameters. As an alternative to using an effective red noise model, I also tried to simultaneously re-fit the timing model model while looking for the gravitational wave signal, but for reasons that are not yet clear, this approach was not very successful. I comment briefly on ways in which the algorithms could be improved. My best estimates for the gravitational wave amplitudes in the three closed (blind) data sets are: (1) $A=(7.3\pm 1.0)\times 10^{-15}$; (2) $A=(5.7\pm 0.6)\times 10^{-14}$; and (3) $A=(4.6\pm 1.3)\times 10^{-15}$.

### Optimal strategies for continuous gravitational wave detection in pulsar timing arrays

Supermassive black hole binaries (SMBHBs) are expected to emit continuous gravitational waves in the pulsar timing array (PTA) frequency band ($10^{-9}$–$10^{-7}$ Hz). The development of data analysis techniques aimed at efficient detection and characterization of these signals is critical to the gravitational wave detection effort. In this paper we leverage methods developed for LIGO continuous wave gravitational searches, and explore the use of the $\mathcal{F}$-statistic for such searches in pulsar timing data. Babak & Sesana 2012 have already used this approach in the context of PTAs to show that one can resolve multiple SMBHB sources in the sky. Our work improves on several aspects of prior continuous wave search methods developed for PTA data analysis. The algorithm is implemented fully in the time domain, which naturally deals with the irregular sampling typical of PTA data and avoids spectral leakage problems associated with frequency domain methods. We take into account the fitting of the timing model, and have generalized our approach to deal with both correlated and uncorrelated colored noise sources. We also develop an incoherent detection statistic that maximizes over all pulsar dependent contributions to the likelihood. To test the effectiveness and sensitivity of our detection statistics, we perform a number of monte-carlo simulations. We produce sensitivity curves for PTAs of various configurations, and outline an implementation of a fully functional data analysis pipeline. Finally, we present a derivation of the likelihood maximized over the gravitational wave phases at the pulsar locations, which results in a vast reduction of the search parameter space.

### Understanding and analysing time-correlated stochastic signals in pulsar timing

Although it is widely understood that pulsar timing observations generally contain time-correlated stochastic signals (TCSSs; red timing noise is of this type), most data analysis techniques that have been developed make an assumption that the stochastic uncertainties in the data are uncorrelated, i.e. “white”. Recent work has pointed out that this can introduce severe bias in determination of timing-model parameters, and that better analysis methods should be used. This paper presents a detailed investigation of timing-model fitting in the presence of TCSSs, and gives closed expressions for the post-fit signals in the data. This results in a Bayesian technique to obtain timing-model parameter estimates in the presence of TCSSs, as well as computationally more efficient expressions of their marginalised posterior distribution. A new method to analyse hundreds of mock dataset realisations simultaneously without significant computational overhead is presented, as well as a statistically rigorous method to check the internal consistency of the results. As a by-product of the analysis, closed expressions of the rms introduced by a stochastic background of gravitational-waves in timing-residuals are obtained. Using $T$ as the length of the dataset, and $h_c(1\rm{yr}^{-1})$ as the characteristic strain, this is: $\sigma_{\rm GWB}^2 = h_{c}(1\rm{yr}^{-1})^2 (9\sqrt[3]{2\pi^4}\Gamma(-10/3) / 8008) \rm{yr}^{-4/3} T^{10/3}$.

### Understanding and analysing time-correlated stochastic signals in pulsar timing [Replacement]

Although it is widely understood that pulsar timing observations generally contain time-correlated stochastic signals (TCSSs; red timing noise is of this type), most data analysis techniques that have been developed make an assumption that the stochastic uncertainties in the data are uncorrelated, i.e. "white". Recent work has pointed out that this can introduce severe bias in determination of timing-model parameters, and that better analysis methods should be used. This paper presents a detailed investigation of timing-model fitting in the presence of TCSSs, and gives closed expressions for the post-fit signals in the data. This results in a Bayesian technique to obtain timing-model parameter estimates in the presence of TCSSs, as well as computationally more efficient expressions of their marginalised posterior distribution. A new method to analyse hundreds of mock dataset realisations simultaneously without significant computational overhead is presented, as well as a statistically rigorous method to check the internal consistency of the results. As a by-product of the analysis, closed expressions of the rms introduced by a stochastic background of gravitational-waves in timing-residuals are obtained. Using $T$ as the length of the dataset, and $h_c(1\rm{yr}^{-1})$ as the characteristic strain, this is: $\sigma_{\rm GWB}^2 = h_{c}(1\rm{yr}^{-1})^2 (9\sqrt[3]{2\pi^4}\Gamma(-10/3) / 8008) \rm{yr}^{-4/3} T^{10/3}$.

### Limits on the Stochastic Gravitational Wave Background from the North American Nanohertz Observatory for Gravitational Waves

We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7×10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.

### Pulsar timing analysis in the presence of correlated noise

Pulsar timing observations are usually analysed with least-square-fitting procedures under the assumption that the timing residuals are uncorrelated (statistically “white”). Pulsar observers are well aware that this assumption often breaks down and causes severe errors in estimating the parameters of the timing model and their uncertainties. Ad hoc methods for minimizing these errors have been developed, but we show that they are far from optimal. Compensation for temporal correlation can be done optimally if the covariance matrix of the residuals is known using a linear transformation that whitens both the residuals and the timing model. We adopt a transformation based on the Cholesky decomposition of the covariance matrix, but the transformation is not unique. We show how to estimate the covariance matrix with sufficient accuracy to optimize the pulsar timing analysis. We also show how to apply this procedure to estimate the spectrum of any time series with a steep red power-law spectrum, including those with irregular sampling and variable error bars, which are otherwise very difficult to analyse.

### Discovery of Eclipses from the Accreting Millisecond X-ray Pulsar SWIFT J1749.4-2807

We report the discovery of X-ray eclipses in the recently discovered accreting millisecond X-ray pulsar Swift J1749.4-2807. This is the first detection of X-ray eclipses in a system of this type and should enable a precise neutron star mass measurement once the companion star is identified and studied. We present a combined pulse and eclipse timing solution that enables tight constraints on the orbital parameters and inclination and shows that the companion mass is in the range 0.6-0.8 M_sun for a likely range of neutron star masses, and that it is larger than a main sequence star of the same mass. We observed two individual eclipse egresses and a single ingress. Our timing model shows that the eclipse features are symmetric about the time of 90 deg longitude from the ascending node, as expected. Our eclipse timing solution gives an eclipse duration (from the mid-points of ingress to egress) of 2172 +/- 13 s. This represents 6.85% of the 8.82 hr orbital period. This system also presents a potential measurement of "Shapiro" delay due to General Relativity; through this technique alone, we set an upper limit to the companion mass of 2.2 M_sun.

### Discovery of Eclipses from the Accreting Millisecond X-ray Pulsar SWIFT J1749.4-2807 [Replacement]

We report the discovery of X-ray eclipses in the recently discovered accreting millisecond X-ray pulsar Swift J1749.4-2807. This is the first detection of X-ray eclipses in a system of this type and should enable a precise neutron star mass measurement once the companion star is identified and studied. We present a combined pulse and eclipse timing solution that enables tight constraints on the orbital parameters and inclination and shows that the companion mass is in the range 0.6-0.8 M_sun for a likely range of neutron star masses, and that it is larger than a main sequence star of the same mass. We observed two individual eclipse egresses and a single ingress. Our timing model shows that the eclipse features are symmetric about the time of 90 deg longitude from the ascending node, as expected. Our eclipse timing solution gives an eclipse duration (from the mid-points of ingress to egress) of 2172 +/- 13 s. This represents 6.85% of the 8.82 hr orbital period. This system also presents a potential measurement of "Shapiro" delay due to General Relativity; through this technique alone, we set an upper limit to the companion mass of 2.2 M_sun.

### The Sensitivity of the Parkes Pulsar Timing Array to Individual Sources of Gravitational Waves [Replacement]

We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed.

### The Sensitivity of the Parkes Pulsar Timing Array to Individual Sources of Gravitational Waves

We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed.

### The Vela Pulsar: Results from the First Year of Fermi LAT Observations

We report on analysis of timing and spectroscopy of the Vela pulsar using eleven months of observations with the Large Area Telescope on the Fermi Gamma-Ray Space Telescope. The intrinsic brightness of Vela at GeV energies combined with the angular resolution and sensitivity of the LAT allow us to make the most detailed study to date of the energy-dependent light curves and phase-resolved spectra, using a LAT-derived timing model. The light curve consists of two peaks (P1 and P2) connected by bridge emission containing a third peak (P3). We have confirmed the strong decrease of the P1/P2 ratio with increasing energy seen with EGRET and previous Fermi LAT data, and observe that P1 disappears above 20 GeV. The increase with energy of the mean phase of the P3 component can be followed with much greater detail, showing that P3 and P2 are present up to the highest energies of pulsation. We find significant pulsed emission at phases outside the main profile, indicating that magnetospheric emission exists over 80% of the pulsar period. With increased high-energy counts the phase-averaged spectrum is seen to depart from a power- law with simple exponential cutoff, and is better fit with a more gradual cutoff. The spectra in fixed-count phase bins are well fit with power-laws with exponential cutoffs, revealing a strong and complex phase dependence of the cutoff energy, especially in the peaks. By combining these results with predictions of the outer magnetosphere models that map emission characteristics to phase, it will be possible to probe the particle acceleration and the structure of the pulsar magnetosphere with unprecedented detail.

### Ruling out Kozai resonance in highly eccentric galactic binary millisecond pulsar PSR J1903+0327 [Replacement]

We investigate the observational signatures associated with one of the proposed formation scenario for the recently discovered highly eccentric binary millisecond pulsar (MSP) PSR J1903+0327 in the galactic plane. The scenario requires that the MSP to be part of a hierarchical triple (HT), consisting of inner and outer binaries, experiencing the Kozai resonance. Numerical modeling of a bound point mass HT, while incorporating the effects due to the quadrupolar interactions between the binary orbits and dominant contributions to the general relativistic periastron precession in the inner binary, reveals that, at the present epoch, the orbital eccentricity of the binary MSP should decrease for reasonable ranges in the HT parameters. The estimated decrements in the orbital eccentricity of the inner binary are few parts in $10^{5}$, substantially higher than the reported accuracies in the estimation of the orbital eccentricity of the binary MSP, while employing various general relativistic timing models for isolated binary pulsars. For wide ranges in the allowed orbital parameters, the estimated rate of change in the eccentricity of the inner binary is orders of magnitude higher than the value recently measured by the pulsar timing analysis. Therefore, we rule out the scenario that the MSP is part of a HT undergoing the Kozai oscillations. The origin of this system in a typical globular cluster is also shown to be less likely than inferred in the discovery paper.

### Spin-down rate and inferred dipole magnetic field of the soft gamma-ray repeater SGR 1627-41

Using Chandra data taken on 2008 June, we detected pulsations at 2.59439(4) s in the soft gamma-ray repeater SGR 1627-41. This is the second measurement of the source spin period and allows us to derive for the first time a long-term spin-down rate of (1.9 +/- 0.4)E-11 s/s. From this value we infer for SGR 1627-41 a characteristic age of 2.2 kyr, a spin-down luminosity of 4E+34 erg/s (one of the highest among sources of the same class), and a surface dipole magnetic field strength of 2E+14 G. These properties confirm the magnetar nature of SGR 1627-41; however, they should be considered with caution since they were derived on the basis of a period derivative measurement made using two epochs only and magnetar spin-down rates are generally highly variable. The pulse profile, double-peaked and with a pulsed fraction of (13 +/- 2)% in the 2-10 keV range, closely resembles that observed by XMM-Newton in 2008 September. Having for the first time a timing model for this SGR, we also searched for a pulsed signal in archival radio data collected with the Parkes radio telescope nine months after the previous X-ray outburst. No evidence for radio pulsations was found, down to a luminosity level 10-20 times fainter (for a 10% duty cycle and a distance of 11 kpc) than the peak luminosity shown by the known radio magnetars.

### You need to log in to vote

The blog owner requires users to be logged in to be able to vote for this post.

Alternatively, if you do not have an account yet you can create one here.