BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Paul S. Ray

Paul S. Ray contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.IMastro-ph.COastro-ph.GAastro-ph.SRgr-qchep-phnucl-thphysics.space-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
9topics
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

15 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The Radius of PSR J0740+6620 from NICER with NICER Background Estimates

We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot regions. We show that there is a significant offset to an antipodal spot configuration, mainly due to the non-half-cycle azimuthal separation of the two emitting spots. The offset angle from the antipode is inferred to be above 25 degrees with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A New Flaring Black Widow Candidate and Demographics of Black Widow Millisecond Pulsars in the Galactic Field

We present the discovery of a new optical/X-ray source likely associated with the Fermi $γ$-ray source 4FGL J1408.6-2917. Its high-amplitude periodic optical variability, large spectroscopic radial velocity semi-amplitude, evidence for optical emission lines and flaring, and X-ray properties together imply the source is probably a new black widow millisecond pulsar binary. We compile the properties of the 41 confirmed and suspected field black widows, finding a median secondary mass of $0.027\pm0.003\,M_{\odot}$. Considered jointly with the more massive redback millisecond pulsar binaries, we find that the &#34;spider&#34; companion mass distribution remains strongly bimodal, with essentially zero systems having companion masses between $\sim0.07-0.1\,M_{\odot}$. X-ray emission from black widows is typically softer and less luminous than in redbacks, consistent with less efficient particle acceleration in the intrabinary shock in black widows, excepting a few systems that appear to have more efficient &#34;redback-like&#34; shocks. Together black widows and redbacks dominate the census of the fastest-spinning field millisecond pulsars in binaries with known companion types, making up $\gtrsim$80% of systems with $P_{\rm{spin}}<2\,\rm{ms}$. Similar to redbacks, the neutron star masses in black widows appear on average significantly larger than the canonical $1.4\,M_{\odot}$, and many of the highest-mass neutron stars claimed to date are black widows with $M_{\rm{NS}}\gtrsim2.1\,M_{\odot}$. Both of these observations are consistent with an evolutionary picture where spider millisecond pulsars emerge from short orbital period progenitors that had a lengthy period of mass transfer initiated while the companion was on the main sequence, leading to fast spins and high masses.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Discovery, Timing, and Multiwavelength Observations of the Black Widow Millisecond Pulsar PSR J1555-2908

We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the radio and gamma-ray pulsation discovery and timing, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ~10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The GMRT High Resolution Southern Sky Survey for pulsars and transients -- III: searching for long period pulsars

Searching for periodic non-accelerated signals in presence of ideal white noise using the fully phase-coherent Fast Folding Algorithm (FFA) is theoretically established as a more sensitive search method than the Fast Fourier Transform (FFT) search with incoherent harmonic summing. In this paper, we present a comparison of the performance of an FFA search implementation using RIPTIDE and an FFT search implementation using PRESTO, over a range of signal parameters with white noise and with real telescope noise from the GHRSS survey with the uGMRT. We find that FFA search with appropriate de-reddening of time series, performs better than FFT search with spectral whitening for long period pulsars in real GHRSS noise conditions. We describe an FFA search pipeline implemented for the GHRSS survey looking for pulsars over a period range of 0.1 s to 100 s and up to dispersion measure of 500 pc cm$^{-3}$. We processed GHRSS survey data covering $\sim$ 1500 degree$^2$ of the sky with this pipeline. We re-detected 43 known pulsars with better signal-to-noise in the FFA search than in the FFT search. We also report discovery of two new pulsars including a long period pulsar having a short duty-cycle with this FFA search pipeline. The population of long period pulsars with periods of several seconds or higher can help to constrain the pulsar death-line.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Detection of Red Noise in PSR J1824$-$2452A and Projections for PSR B1937+21 using NICER X-ray Timing Data

We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations (&#34;red noise&#34;) in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar medium. Our Bayesian search methodology yields strong evidence (natural log Bayes factor of $9.634 \pm 0.016$) for red noise in PSR J1824$-$2452A, but is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulate longer and higher precision X-ray datasets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5-year mission with a time-of-arrival (TOA) error of 2 microseconds and an observing cadence of 20 observations per month compared to the 5 microsecond TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that an injected stochastic gravitational wave background (GWB) with an amplitude of $A_{\rm GWB}=2\times10^{-15}$ and spectral index of $γ_{\rm GWB}=13/3$ can be detected in a pulsar with similar TOA precision to PSR B1937+21, but with no additional red noise, in a 10-year mission that observes the pulsar 15 times per month and has an average TOA error of 1 microsecond.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An Empirical Background Model for the NICER X-ray Timing Instrument

NICER has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4-12 keV; 50 detectors) of 0.87 c/s, but in 5 percent (1 percent) of cases, the rate exceeds 10 (300) c/s. Model residuals persist at 20-30 percent of the initial rate for the brightest GTIs, implying one or more missing model parameters. Filtering criteria are given to flag GTIs likely to have unsatisfactory background predictions. With such filtering, we estimate a detection limit, 1.20 c/s (3 sigma, single GTI) at 0.4-12 keV, equivalent to 3.6e-12 erg/cm^2/s for a Crab-like spectrum. The corresponding limit for soft X-ray sources is 0.51 c/s at 0.3-2.0 keV, or 4.3e-13 erg/cm^2/s for a 100 eV blackbody. Faint-source filtering selects 85 percent of the background GTIs, and higher rates are expected for targets scheduled more favorably. An application of the model to 1 s timescale makes it possible to distinguish source flares from possible surges in the background.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bayesian Solar Wind Modeling with Pulsar Timing Arrays

Using Bayesian analyses we study the solar electron density with the NANOGrav 11-year pulsar timing array (PTA) dataset. Our model of the solar wind is incorporated into a global fit starting from pulse times-of-arrival. We introduce new tools developed for this global fit, including analytic expressions for solar electron column densities and open source models for the solar wind that port into existing PTA software. We perform an ab initio recovery of various solar wind model parameters. We then demonstrate the richness of information about the solar electron density, $n_E$, that can be gleaned from PTA data, including higher order corrections to the simple $1/r^2$ model associated with a free-streaming wind (which are informative probes of coronal acceleration physics), quarterly binned measurements of $n_E$ and a continuous time-varying model for $n_E$ spanning approximately one solar cycle period. Finally, we discuss the importance of our model for chromatic noise mitigation in gravitational-wave analyses of pulsar timing data and the potential of developing synergies between sophisticated PTA solar electron density models and those developed by the solar physics community.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

NICER Discovery of Millisecond X-ray Pulsations and an Ultracompact Orbit in IGR J17494-3030

We report the detection of 376.05 Hz (2.66 ms) coherent X-ray pulsations in NICER observations of a transient outburst of the low-mass X-ray binary IGR J17494-3030 in 2020 October/November. The system is an accreting millisecond X-ray pulsar in a 75 minute ultracompact binary. The mass donor is most likely a $\simeq 0.02 M_\odot$ finite-entropy white dwarf composed of He or C/O. The fractional rms pulsed amplitude is 7.4%, and the soft (1-3 keV) X-ray pulse profile contains a significant second harmonic. The pulsed amplitude and pulse phase lag (relative to our mean timing model) are energy-dependent, each having a local maximum at 4 keV and 1.5 keV, respectively. We also recovered the X-ray pulsations in archival 2012 XMM-Newton observations, allowing us to measure a long-term pulsar spin-down rate of $\dotν= -2.1(7)\times10^{-14}$ Hz/s and to infer a pulsar surface dipole magnetic field strength of $\simeq 10^9$ G. We show that the mass transfer in the binary is likely non-conservative, and we discuss various scenarios for mass loss from the system.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A NICER View of Spectral and Profile Evolution for Three X-ray Emitting Millisecond Pulsars

We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at $\gt 3σ$ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821-24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821-24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bright X-ray and Radio Pulses from a Recently Reactivated Magnetar

Magnetars are young, rotating neutron stars that possess larger magnetic fields ($B$ $\approx$ $10^{13}$-$10^{15}$ G) and longer rotational periods ($P$ $\approx$ 1-12 s) than ordinary pulsars. In contrast to rotation-powered pulsars, magnetar emission is thought to be fueled by the evolution and decay of their powerful magnetic fields. They display highly variable radio and X-ray emission, but the processes responsible for this behavior remain a mystery. We report the discovery of bright, persistent individual X-ray pulses from XTE J1810-197, a transient radio magnetar, using the Neutron star Interior Composition Explorer (NICER) following its recent radio reactivation. Similar behavior has only been previously observed from a magnetar during short time periods following a giant flare. However, the X-ray pulses presented here were detected outside of a flaring state. They are less energetic and display temporal structure that differs from the impulsive X-ray events previously observed from the magnetar class, such as giant flares and short X-ray bursts. Our high frequency radio observations of the magnetar, carried out simultaneously with the X-ray observations, demonstrate that the relative alignment between the X-ray and radio pulses varies on rotational timescales. No correlation was found between the amplitudes or temporal structure of the X-ray and radio pulses. The magnetar&#39;s 8.3 GHz radio pulses displayed frequency structure, which was not observed in the pulses detected simultaneously at 31.9 GHz. Many of the radio pulses were also not detected simultaneously at both frequencies, which indicates that the underlying emission mechanism producing these pulses is not broadband. We find that the radio pulses from XTE J1810-197 share similar characteristics to radio bursts detected from fast radio burst (FRB) sources, some of which are now thought to be produced by active magnetars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind&#39; pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Timing the pulsations of the accreting millisecond pulsar SAX J1808.4-3658 during its 2019 outburst

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth&#39;s orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

0 citations0 reviewsNo code linkedOpen access