-
An X-ray and radio view of the 2022 reactivation of the magnetar SGRJ1935+2154
Authors:
A. Y. Ibrahim,
A. Borghese,
F. Coti Zelati,
E. Parent,
A. Marino,
O. S. Ould-Boukattine,
N. Rea,
S. Ascenzi,
D. P. Pacholski,
S. Mereghetti,
G. L. Israel,
A. Tiengo,
A. Possenti,
M. Burgay,
R. Turolla,
S. Zane,
P. Esposito,
D. Gotz,
S. Campana,
F. Kirsten,
M. P. Gawronski,
J. W. T. Hessels
Abstract:
Recently, the Galactic magnetar SGR J1935+2154 has garnered attention due to its emission of an extremely luminous radio burst, reminiscent of Fast Radio Bursts (FRBs). SGR J1935+2154 is one of the most active magnetars, displaying flaring events nearly every year, including outbursts as well as short and intermediate bursts. Here, we present our results on the properties of the persistent and bur…
▽ More
Recently, the Galactic magnetar SGR J1935+2154 has garnered attention due to its emission of an extremely luminous radio burst, reminiscent of Fast Radio Bursts (FRBs). SGR J1935+2154 is one of the most active magnetars, displaying flaring events nearly every year, including outbursts as well as short and intermediate bursts. Here, we present our results on the properties of the persistent and bursting X-ray emission from SGR J1935+2154, during the initial weeks following its outburst on October 10, 2022. The source was observed with XMM-Newton and NuSTAR (quasi-)simultaneously during two epochs, separated by $\sim$5 days. The persistent emission spectrum is well described by an absorbed blackbody plus power-law model up to an energy of $\sim$25 keV. No significant changes were observed in the blackbody temperature ($kT_{\rm BB}\sim$ 0.4 keV) and emitting radius ($R_{\rm BB}\sim$ 1.9 km) between the two epochs. However, we observed a slight variation in the power-law parameters. Moreover, we detected X-ray pulsations in all the datasets and derived a spin period derivative of $\dot{P} = 5.52(5) \times 10^{-11}$ ss. This is 3.8 times larger than the value measured after the first recorded outburst in 2014. Additionally, we performed quasi-simultaneous radio observations using three 25--32-m class radio telescopes for a total of 92.5 hr to search for FRB-like radio bursts and pulsed emission. However, our analysis did not reveal any radio bursts or periodic emission.
△ Less
Submitted 13 February, 2024;
originally announced February 2024.
-
Long-term study of the 2020 magnetar-like outburst of the young pulsar PSRJ1846-0258 in Kes 75
Authors:
Rajath Sathyaprakash,
Emilie Parent,
Nanda Rea,
Francesco Coti Zelati,
Alice Borghese,
Maura Pilia,
Matteo Trudu,
Marta Burgay,
Roberto Turolla,
Silvia Zane,
Paolo Esposito,
Sandro Mereghetti,
Sergio Campana,
Diego Götz,
Abubakr Ibrahim,
GianLuca Israel,
Andrea Possenti,
Andrea Tiengo
Abstract:
Magnetar-like activity has been observed in a large variety of neutron stars. PSR J1846-0258 is a young 327 ms radio-quiet pulsar with a large rotational power ($\sim 8 \times 10^{36}$ erg s$^{-1}$), and resides at the center of the supernova remnant Kes 75. It is one of the rare examples of a high magnetic field pulsar showing characteristics both of magnetars and radio pulsars, and can thus prov…
▽ More
Magnetar-like activity has been observed in a large variety of neutron stars. PSR J1846-0258 is a young 327 ms radio-quiet pulsar with a large rotational power ($\sim 8 \times 10^{36}$ erg s$^{-1}$), and resides at the center of the supernova remnant Kes 75. It is one of the rare examples of a high magnetic field pulsar showing characteristics both of magnetars and radio pulsars, and can thus provide important clues on the differences in the emission mechanisms between these two classes. In 2006, PSR J1846-0258 was detected to undergo an outburst for the first time, accompanied by a large flux increase, millisecond X-ray bursts, significant spectral changes and a large timing glitch. On August 1st 2020, after nineteen years of quiescent stable emission, the source underwent a second magnetar-like outburst, which was followed up with several observations by NICER, XMM-Newton, NuSTAR and Swift. In this work, we report on the long-term timing and X-ray spectral properties of the source following the 2020 outburst, and place upper limits on any radio activity. We demonstrate that the pulsed flux increased by a factor $> 6$ during the outburst, followed by non-trivial variability in the spin-down rate. Our timing analysis shows that the spin frequency and its derivative are clearly affected by magnetospheric activity due to the outburst. We find hints for an oscillation in the frequency derivative with a timescale of 50-60 days, recovering later on to stable quiescence.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
The Green Bank North Celestial Cap Survey IX: Timing Follow-up for 128 Pulsars
Authors:
A. E. McEwen,
J. K. Swiggum,
D. L. Kaplan,
C. M. Tan,
B. W. Meyers,
E. Fonseca,
G. Y. Agazie,
P. Chawla,
K. Crowter,
M. E. DeCesar,
T. Dolch,
F. A. Dong,
W. Fiore,
E. Fonseca,
D. C. Good,
A. G. Istrate,
V. M. Kaspi,
V. I. Kondratiev,
J. van Leeuwen,
L. Levin,
E. F. Lewis,
R. S. Lynch,
K. W. Masui,
J. W. McKee,
M. A. McLaughlin
, et al. (6 additional authors not shown)
Abstract:
The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations take…
▽ More
The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including 9 previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6)\Msun and 1.412(6)\Msun, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves: $\dot{P}_{\rm B} = -1.37(7)\times10^{-12}$, which is in agreement with the rate predicted by general relativity for these masses.
△ Less
Submitted 26 July, 2024; v1 submitted 12 December, 2023;
originally announced December 2023.
-
The Green Bank North Celestial Cap Survey. VIII. 21 New Pulsar Timing Solutions
Authors:
William Fiore,
Lina Levin,
Maura A. McLaughlin,
Akash Anumarlapudi,
David L. Kaplan,
Joseph K. Swiggum,
Gabriella Y. Agazie,
Robert Bavisotto,
Pragya Chawla,
Megan E. DeCesar,
Timothy Dolch,
Emmanuel Fonseca,
Victoria M. Kaspi,
Zachary Komassa,
Vlad I. Kondratiev,
Joeri van Leeuwen,
Evan F. Lewis,
Ryan S. Lynch,
Alexander E. McEwen,
Rusty Mundorf,
Hind Al Noori,
Emilie Parent,
Ziggy Pleunis,
Scott M. Ransom,
Xavier Siemens
, et al. (4 additional authors not shown)
Abstract:
We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957-0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528-day…
▽ More
We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957-0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528-day orbit, PSR J0032+6946 joins a small group of five other mildly recycled wide binary pulsars, for which the duration of recycling through accretion is limited by the length of the companion's giant phase. PSRs J0141+6303 and J1327+3423 are new disrupted recycled pulsars. We incorporate Arecibo observations from the NANOGrav pulsar timing array into our analysis of the latter. We also observed PSR J1327+3423 with the Long Wavelength Array, and our data suggest a frequency-dependent dispersion measure. PSR J0957-0619 was discovered as a rotating radio transient, but is a nulling pulsar at 820 MHz. PSR J1239+3239 is a new millisecond pulsar (MSP) in a 4-day orbit with a low-mass companion. Four of our pulsars already have published timing solutions, which we update in this work: the recycled wide binary PSR J0214+5222, the non-eclipsing black widow PSR J0636+5128, the disrupted recycled pulsar J1434+7257, and the eclipsing binary MSP J1816+4510, which is in an 8.7 hr orbit with a redback-mass companion.
△ Less
Submitted 22 May, 2023;
originally announced May 2023.
-
PODIUM:A Pulsar Navigation Unit for Science Missions
Authors:
Francesco Cacciatore,
Víctor Gómez Ruiz,
Gonzalo Taubmann,
Jacinto Muñoz,
Pablo Hermosín,
Marcello Sciarra,
Martiño Saco,
Nanda Rea,
Margarita Hernanz,
Emilie Parent,
Jeroen Vandersteen
Abstract:
PODIUM is a compact spacecraft navigation unit, currently being designed to provide interplanetary missions with autonomous position and velocity estimations. The unit will make use of Pulsar X-ray observations to measure the distance and distance rate from the host spacecraft to the Solar System Barycenter. Such measurements will then be used by the onboard orbit determination function to estimat…
▽ More
PODIUM is a compact spacecraft navigation unit, currently being designed to provide interplanetary missions with autonomous position and velocity estimations. The unit will make use of Pulsar X-ray observations to measure the distance and distance rate from the host spacecraft to the Solar System Barycenter. Such measurements will then be used by the onboard orbit determination function to estimate the complete orbital elements of the spacecraft. The design aims at 6 kg of mass and 20 W of power, in a volume of 150 mm by 240 mm by 600 mm. PODIUM is designed to minimize the impact on the mission operational and accommodation constraints. The architecture is based on a grazing incidence X-ray telescope with focal distance limited to 50 cm. The effective area shall be in the range 25 to 50 cm2 for photon energies in the range 0.2-10 keV, requiring nesting of several mirrors in the Wolter-1 geometry. Grazing incidence angles will be very small, below 2 deg. The current target FOV is 0.25 deg. The pulsars photon arrivals are detected with a single pixel Silicon Drift Detector (SDD) sensor with timing accuracy below 1usec. The unit has no gimbaling to meet the applicable power, size and mass requirements. Instead, the host spacecraft shall slew and point to allow pulsar observation. The avionics architecture is based on a radiation hardened LEON4 processor, to allow a synchronous propagation task and measurement generation and orbit determination step in an asynchronous task. PODIUM will enable higher autonomy and lower cost for interplanetary missions. L2 space observatories and planetary flybys are the current reference use cases. Onboard autonomous state estimation can reduce the ground support effort required for navigation and orbit correction/maintenance computation, and reduce the turnaround time, thus enabling more accurate maneuvers, reducing the orbit maintenance mass budget.
△ Less
Submitted 20 January, 2023;
originally announced January 2023.
-
The Green Bank North Celestial Cap Survey. VII. 12 New Pulsar Timing Solutions
Authors:
Joseph K. Swiggum,
Ziggy Pleunis,
Emilie Parent,
David L. Kaplan,
Maura A. McLaughlin,
Ingrid H. Stairs,
Renée Spiewak,
Gabriella Y. Agazie,
Pragya Chawla,
Megan E. DeCesar,
Timothy Dolch,
William Fiore,
Emmanuel Fonseca,
Alina G. Istrate,
Victoria M. Kaspi,
Vlad I. Kondratiev,
Joeri van Leeuwen,
Lina Levin,
Evan F. Lewis,
Ryan S. Lynch,
Alex E. McEwen,
Hind Al Noori,
Scott M. Ransom,
Xavier Siemens,
Mayuresh Surnis
Abstract:
We present timing solutions for 12 pulsars discovered in the Green Bank North Celestial Cap (GBNCC) 350 MHz pulsar survey, including six millisecond pulsars (MSPs), a double neutron star (DNS) system, and a pulsar orbiting a massive white dwarf companion. Timing solutions presented here include 350 and 820 MHz Green Bank Telescope data from initial confirmation and follow-up as well as a dedicated…
▽ More
We present timing solutions for 12 pulsars discovered in the Green Bank North Celestial Cap (GBNCC) 350 MHz pulsar survey, including six millisecond pulsars (MSPs), a double neutron star (DNS) system, and a pulsar orbiting a massive white dwarf companion. Timing solutions presented here include 350 and 820 MHz Green Bank Telescope data from initial confirmation and follow-up as well as a dedicated timing campaign spanning one year. PSR J1122$-$3546 is an isolated MSP, PSRs J1221$-$0633 and J1317$-$0157 are MSPs in black widow systems and regularly exhibit eclipses, and PSRs J2022+2534 and J2039$-$3616 are MSPs that can be timed with high precision and have been included in pulsar timing array experiments seeking to detect low-frequency gravitational waves. PSRs J1221$-$0633 and J2039$-$3616 have Fermi Large Area Telescope $γ$-ray counterparts and also exhibit significant $γ$-ray pulsations. We measure proper motion for three of the MSPs in this sample and estimate their space velocities, which are typical compared to those of other MSPs. We have detected the advance of periastron for PSR J1018$-$1523 and therefore measure the total mass of the double neutron star system, $m_{\rm tot}=2.3\pm0.3$ M$_{\odot}$. Long-term pulsar timing with data spanning more than one year is critical for classifying recycled pulsars, carrying out detailed astrometry studies, and shedding light on the wealth of information in these systems post-discovery.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Deep X-ray and radio observations of the first outburst of the young magnetar Swift J1818.0-1607
Authors:
A. Y. Ibrahim,
A. Borghese,
N. Rea,
F. Coti Zelati,
E. Parent,
T. D. Russell,
S. Ascenzi,
R. Sathyaprakash,
D. Gotz,
S. Mereghetti,
M. Topinka,
M. Rigoselli,
V. Savchenko,
S. Campana,
G. L. Israel,
A. Tiengo,
R. Perna,
R. Turolla,
S. Zane,
P. Esposito,
G. A. Rodrıguez Castillo,
V. Graber,
A. Possenti,
C. Dehman,
M. Ronchi
, et al. (1 additional authors not shown)
Abstract:
Swift J1818.0-1607 is a radio-loud magnetar with a spin period of 1.36 s and a dipolar magnetic field strength of B~3E14 G, which is very young compared to the Galactic pulsar population. We report here on the long-term X-ray monitoring campaign of this young magnetar using XMM-Newton, NuSTAR, and Swift from the activation of its first outburst in March 2020 until October 2021, as well as INTEGRAL…
▽ More
Swift J1818.0-1607 is a radio-loud magnetar with a spin period of 1.36 s and a dipolar magnetic field strength of B~3E14 G, which is very young compared to the Galactic pulsar population. We report here on the long-term X-ray monitoring campaign of this young magnetar using XMM-Newton, NuSTAR, and Swift from the activation of its first outburst in March 2020 until October 2021, as well as INTEGRAL upper limits on its hard X-ray emission. The 1-10 keV magnetar spectrum is well modeled by an absorbed blackbody with a temperature of kT_BB~1.1 keV, and apparent reduction in the radius of the emitting region from ~0.6 to ~0.2 km. We also confirm the bright diffuse X-ray emission around the source extending between ~50'' and ~110''. A timing analysis revealed large torque variability, with an average spin-down rate nudot~-2.3E-11 Hz^2 that appears to decrease in magnitude over time. We also observed Swift J1818.0-1607 with the Karl G. Jansky Very Large Array (VLA) on 2021 March 22. We detected the radio counterpart to Swift J1818.0-1607 measuring a flux density of S_v = 4.38+/-0.05 mJy at 3 GHz, and a half ring-like structure of bright diffuse radio emission located at ~90'' to the west of the magnetar. We tentatively suggest that the diffuse X-ray emission is due to a dust scattering halo and that the radio structure may be associated with the supernova remnant of this young pulsar, based on its morphology.
△ Less
Submitted 22 November, 2022;
originally announced November 2022.
-
Relativistic effects in a mildly recycled pulsar binary: PSR J1952+2630
Authors:
T. Gautam,
P. C. C. Freire,
A. Batrakov,
M. Kramer,
C. C. Miao,
E. Parent,
W. W. Zhu
Abstract:
We report the results of timing observations of PSR J1952+2630, a 20.7 ms pulsar in orbit with a massive white dwarf companion. With the increased timing baseline, we obtain improved estimates for astrometric, spin, and binary parameters for this system. We get an improvement of an order of magnitude on the proper motion, and, for the first time, we detect three post-Keplerian parameters in this s…
▽ More
We report the results of timing observations of PSR J1952+2630, a 20.7 ms pulsar in orbit with a massive white dwarf companion. With the increased timing baseline, we obtain improved estimates for astrometric, spin, and binary parameters for this system. We get an improvement of an order of magnitude on the proper motion, and, for the first time, we detect three post-Keplerian parameters in this system: the advance of periastron, the orbital decay, and the Shapiro delay. We constrain the pulsar mass to 1.20$^{+0.28}_{-0.29}\rm M_{\odot}$ and the mass of its companion to 0.97$^{+0.16}_{-0.13}\rm M_{\odot}$. The current value of $\dot{P}_{\rm b}$ is consistent with GR expectation for the masses obtained using $\dotω$ and $h_3$. The excess represents a limit on the emission of dipolar GWs from this system. This results in a limit on the difference in effective scalar couplings for the pulsar and companion (predicted by scalar-tensor theories of gravity; STTs) of $|α_{\rm p}-α_{\rm c}| < 4.8 \times 10^{-3}$, which does not yield a competitive test for STTs. However, our simulations of future campaigns of this system show that by 2032, the precision of $\dot{P}_{\rm b}$ and $\dotω$ will allow for much more precise masses and much tighter constraints on the orbital decay contribution from dipolar GWs, resulting in $|α_{\rm p}-α_{\rm c}|<1.3 \times 10^{-3}$. We also present the constraints this system will place on the $\{α_0,β_0\}$ parameters of DEF gravity by 2032. They are comparable to those of PSR J1738+0333. Unlike PSR J1738+0333, PSR J1952+2630 will not be limited in its mass measurement and has the potential to place even more restrictive limits on DEF gravity in the future. Further improvements to this test will likely be limited by uncertainties in the kinematic contributions to $\dot{P}_{\rm b}$ due to lack of precise distance measurements.
△ Less
Submitted 7 October, 2022;
originally announced October 2022.
-
Constraining the nature of the 18-min periodic radio transient GLEAM-X J162759.5-523504.3 via multi-wavelength observations and magneto-thermal simulations
Authors:
N. Rea,
F. Coti Zelati,
C. Dehman,
N. Hurley-Walker,
D. De Martino,
A. Bahramian,
D. A. H. Buckley,
J. Brink,
A. Kawka,
J. A. Pons,
D. Vigano,
V. Graber,
M. Ronchi,
C. Pardo,
A. Borghese,
E. Parent
Abstract:
We observed the periodic radio transient GLEAM-X J162759.5-523504.3 (GLEAM-X J1627) using the Chandra X-ray Observatory for about 30-ks on January 22-23, 2022, simultaneously with radio observations from MWA, MeerKAT and ATCA. Its radio emission and 18-min periodicity led the source to be tentatively interpreted as an extreme magnetar or a peculiar highly magnetic white dwarf. The source was not d…
▽ More
We observed the periodic radio transient GLEAM-X J162759.5-523504.3 (GLEAM-X J1627) using the Chandra X-ray Observatory for about 30-ks on January 22-23, 2022, simultaneously with radio observations from MWA, MeerKAT and ATCA. Its radio emission and 18-min periodicity led the source to be tentatively interpreted as an extreme magnetar or a peculiar highly magnetic white dwarf. The source was not detected in the 0.3-8 keV energy range with a 3-sigma upper-limit on the count rate of 3x10^{-4} counts/s. No radio emission was detected during our X-ray observations either. Furthermore, we studied the field around GLEAM-X J1627 using archival ESO and DECam data, as well as recent SALT observations. Many sources are present close to the position of GLEAM-X J1627, but only two within the 2" radio position uncertainty. Depending on the assumed spectral distribution, the upper limits converted to an X-ray luminosity of L_{X}<6.5x10^{29} erg/s for a blackbody with temperature kT=0.3 keV, or L_{X}<9x10^{29} erg/s for a power-law with photon index Gamma = 2 (assuming a 1.3 kpc distance). Furthermore, we performed magneto-thermal simulations for neutron stars considering crust- and core-dominated field configurations. Based on our multi-band limits, we conclude that: i) in the magnetar scenario, the X-ray upper limits suggest that GLEAM-X J1627 should be older than ~1 Myr, unless it has a core-dominated magnetic field or has experienced fast-cooling; ii) in the white dwarf scenario, we can rule out most binary systems, a hot sub-dwarf and a hot magnetic isolated white dwarf (T>10.000 K), while a cold isolated white dwarf is still compatible with our limits.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary
Authors:
E. Parent,
H. Sewalls,
P. C. C. Freire,
T. Matheny,
A. G. Lyne,
B. B. P. Perera,
F. Cardoso,
M. A. McLaughlin,
B. Allen,
A. Brazier,
F. Camilo,
S. Chatterjee,
J. M. Cordes,
F. Crawford,
J. S. Deneva,
F. A. Dong,
R. D. Ferdman,
E. Fonseca,
J. W. T. Hessels,
V. M. Kaspi,
B. Knispel,
J. van Leeuwen,
R. S. Lynch,
B. M. Meyers,
J. W. McKee
, et al. (9 additional authors not shown)
Abstract:
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, 1…
▽ More
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.
△ Less
Submitted 17 November, 2021; v1 submitted 4 August, 2021;
originally announced August 2021.
-
The Green Bank Northern Celestial Cap Pulsar Survey. VI. Timing and Discovery of PSR J1759+5036: A Double Neutron Star Binary Pulsar
Authors:
Gabriella Agazie,
Michael Mingyar,
Maura McLaughlin,
Joseph Swiggum,
David Kaplan,
Harsha Blumer,
Pragya Chawla,
Megan DeCesar,
Paul Demorest,
William Fiore,
Emmanuel Fonseca,
Joseph Gelfand,
Victoria Kaspi,
Vladislav Kondratiev,
Malcolm LaRose,
Joeri van Leeuwen,
Lina Levin,
Evan Lewis,
Ryan Lynch,
Alexander McEwen,
Hind Al Noori,
Emilie Parent,
Scott Ransom,
Mallory Roberts,
Ann Schmiedekamp
, et al. (5 additional authors not shown)
Abstract:
The Green Bank North Celestial Cap (GBNCC) survey is a 350-MHz all-sky survey for pulsars and fast radio transients using the Robert C. Byrd Green Bank Telescope. To date, the survey has discovered over 190 pulsars, including 33 millisecond pulsars (MSPs) and 24 rotating radio transients(RRATs). Several exotic pulsars have been discovered in the survey, including PSR J1759+5036, a binary pulsar wi…
▽ More
The Green Bank North Celestial Cap (GBNCC) survey is a 350-MHz all-sky survey for pulsars and fast radio transients using the Robert C. Byrd Green Bank Telescope. To date, the survey has discovered over 190 pulsars, including 33 millisecond pulsars (MSPs) and 24 rotating radio transients(RRATs). Several exotic pulsars have been discovered in the survey, including PSR J1759+5036, a binary pulsar with a 176-ms spin period in an orbit with a period of 2.04 days, an eccentricity of 0.3,and a projected semi-major axis of 6.8 light seconds. Using seven years of timing data, we are able to measure one post-Keplerian parameter, advance of periastron, which has allowed us to constrain the total system mass to 2.62(3) solar masses. This constraint, along with the spin period and orbital parameters, suggests that this is a double neutron star system, although we cannot entirely rule out a pulsar-white dwarf binary. This pulsar is only detectable in roughly 45% of observations, most likely due to scintillation. However, additional observations are required to determine whether there may be other contributing effects.
△ Less
Submitted 12 July, 2021; v1 submitted 19 February, 2021;
originally announced February 2021.
-
First Discovery of a Fast Radio Burst at 350 MHz by the GBNCC Survey
Authors:
E. Parent,
P. Chawla,
V. M. Kaspi,
G. Y. Agazie,
H. Blumer,
M. DeCesar,
W. Fiore,
E. Fonseca,
J. W. T. Hessels,
D. L. Kaplan,
V. I. Kondratiev,
M. LaRose,
L. Levin,
E. F. Lewis,
R. S. Lynch,
A. E. McEwen,
M. A. McLaughlin,
M. Mingyar,
H. Al Noori,
S. M. Ransom,
M. S. E. Roberts,
A. Schmiedekamp,
C. Schmiedekamp,
X. Siemens,
R. Spiewak
, et al. (4 additional authors not shown)
Abstract:
We report the first discovery of a fast radio burst (FRB), FRB 20200125A, by the Green Bank Northern Celestial Cap (GBNCC) Pulsar Survey conducted with the Green Bank Telescope at 350 MHz. FRB 20200125A was detected at a Galactic latitude of 58.43 degrees with a dispersion measure of 179 pc cm$^{-3}$, while electron density models predict a maximum Galactic contribution of 25 pc cm$^{-3}$ along th…
▽ More
We report the first discovery of a fast radio burst (FRB), FRB 20200125A, by the Green Bank Northern Celestial Cap (GBNCC) Pulsar Survey conducted with the Green Bank Telescope at 350 MHz. FRB 20200125A was detected at a Galactic latitude of 58.43 degrees with a dispersion measure of 179 pc cm$^{-3}$, while electron density models predict a maximum Galactic contribution of 25 pc cm$^{-3}$ along this line of sight. Moreover, no apparent Galactic foreground sources of ionized gas that could account for the excess DM are visible in multi-wavelength surveys of this region. This argues that the source is extragalactic. The maximum redshift for the host galaxy is $z_{max}=0.17$, corresponding to a maximum comoving distance of approximately 750 Mpc. The measured peak flux density for FRB 20200125A is 0.37 Jy, and we measure a pulse width of 3.7 ms, consistent with the distribution of FRB widths observed at higher frequencies. Based on this detection and assuming an Euclidean flux density distribution of FRBs, we calculate an all-sky rate at 350 MHz of $3.4^{+15.4}_{-3.3} \times 10^3$ FRBs sky$^{-1}$ day$^{-1}$ above a peak flux density of 0.42 Jy for an unscattered pulse having an intrinsic width of 5 ms, consistent with rates reported at higher frequencies. Given the recent improvements in our single-pulse search pipeline, we also revisit the GBNCC survey sensitivity to various burst properties. Finally, we find no evidence of interstellar scattering in FRB 20200125A, adding to the growing evidence that some FRBs have circumburst environments where free-free absorption and scattering are not significant.
△ Less
Submitted 10 August, 2020;
originally announced August 2020.
-
Asymmetric mass ratios for bright double neutron-star mergers
Authors:
R. D. Ferdman,
P. C. C. Freire,
B. B. P. Perera,
N. Pol,
F. Camilo,
S. Chatterjee,
J. M. Cordes,
F. Crawford,
J. W. T. Hessels,
V. M. Kaspi,
M. A. McLaughlin,
E. Parent,
I. H. Stairs,
J. van Leeuwen
Abstract:
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 potentiall…
▽ More
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.
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
Eight Millisecond Pulsars Discovered in the Arecibo PALFA Survey
Authors:
E. Parent,
V. M. Kaspi,
S. M. Ransom,
P. C. C. Freire,
A. Brazier,
F. Camilo,
S. Chatterjee,
J. M. Cordes,
F. Crawford,
J. S. Deneva,
R. D. Ferdman,
J. W. T. Hessels,
J. van Leeuwen,
A. G. Lyne,
E. C. Madsen,
M. A. McLaughlin,
C. Patel,
P. Scholz,
I. H. Stairs,
B. W. Stappers,
W. W. Zhu
Abstract:
We report on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo PALFA survey. Phase-coherent timing solutions derived from 2.5 to 5 years of observations carried out at Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65-ms pulsar in a 39.6-day orbit for which we detect $γ$-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mas…
▽ More
We report on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo PALFA survey. Phase-coherent timing solutions derived from 2.5 to 5 years of observations carried out at Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65-ms pulsar in a 39.6-day orbit for which we detect $γ$-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mass-function system with an orbital period of 3.3 hours that belongs to the non-eclipsing black widow population. We also present PSR J1932+1756, the longest-orbital-period (41.5 days) intermediate-mass binary pulsar known to date. In light of the numerous discoveries of binary MSPs over the past years, we characterize the Galactic distribution of known MSP binaries in terms of binary class. Our results support and strengthen previous claims that the scatter in the Galactic scale height distribution correlates inversely with the binary mass function. We provide evidence of observational biases against detecting the most recycled pulsars near the Galactic plane, which overestimates the scale height of lighter systems. A possible bimodality in the mass function of MSPs with massive white dwarfs is also reported.
△ Less
Submitted 26 August, 2019;
originally announced August 2019.
-
Computer code validation via mixture model estimation
Authors:
Kaniav Kamary,
Merlin Keller,
Pierre Barbillon,
Cédric Gœury,
Éric Parent
Abstract:
When computer codes are used for modeling complex physical systems, their unknown parameters are tuned by calibration techniques. A discrepancy function may be added to the computer code in order to capture its discrepancy with the real physical process. By considering the validation question of a computer code as a Bayesian selection model problem, Damblin et al. (2016) have highlighted a possibl…
▽ More
When computer codes are used for modeling complex physical systems, their unknown parameters are tuned by calibration techniques. A discrepancy function may be added to the computer code in order to capture its discrepancy with the real physical process. By considering the validation question of a computer code as a Bayesian selection model problem, Damblin et al. (2016) have highlighted a possible confounding effect in certain configurations between the code discrepancy and a linear computer code by using a Bayesian testing procedure based on the intrinsic Bayes factor. In this paper, we investigate the issue of code error identifiability by applying another Bayesian model selection technique which has been recently developed by Kamary et al. (2014). By embedding the competing models within an encompassing mixture model, Kamary et al. (2014)'s method allows each observation to belong to a different mixing component, providing a more flexible inference, while remaining competitive in terms of computational cost with the intrinsic Bayesian approach. By using the technique of sharing parameters mentioned in Kamary et al. (2014), an improper non-informative prior can be used for some computer code parameters and we demonstrate that the resulting posterior distribution is proper. We then check the sensitivity of our posterior estimates to the choice of the parameter prior distributions. We illustrate that the value of the correlation length of the discrepancy Gaussian process prior impacts the Bayesian inference of the mixture model parameters and that the model discrepancy can be identified by applying the Kamary et al. (2014) method when the correlation length is not too small. Eventually, the proposed method is applied on a hydraulic code in an industrial context.
△ Less
Submitted 8 March, 2019;
originally announced March 2019.
-
Probing type Ia supernova properties using bolometric light curves from the Carnegie Supernova Project and the CfA Supernova Group
Authors:
Richard A. Scalzo,
Emilie Parent,
Christopher Burns,
Michael J. Childress,
Brad E. Tucker,
Peter J. Brown,
Carlos Contreras,
Eric Hsiao,
Kevin Krisciunas,
Nidia Morrell,
Mark M. Phillips,
Anthony L. Piro,
Maximilian Stritzinger,
Nicholas Suntzeff
Abstract:
We present bolometric light curves constructed from multi-wavelength photometry of Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project and the CfA Supernova Group, using near-infrared observations to provide robust constraints on host galaxy dust extinction. This set of light curves form a well-measured reference set for comparison with theoretical models. Ejected mass and synthesized…
▽ More
We present bolometric light curves constructed from multi-wavelength photometry of Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project and the CfA Supernova Group, using near-infrared observations to provide robust constraints on host galaxy dust extinction. This set of light curves form a well-measured reference set for comparison with theoretical models. Ejected mass and synthesized $^{56}$Ni mass are inferred for each SN Ia from its bolometric light curve using a semi-analytic Bayesian light curve model, and fitting formulae provided in terms of light curve width parameters from the SALT2 and SNooPy light curve fitters. A weak bolometric width-luminosity relation is confirmed, along with a correlation between ejected mass and the bolometric light curve width. SNe Ia likely to have sub-Chandrasekhar ejected masses belong preferentially to the broad-line and cool-photosphere spectroscopic subtypes, and have higher photospheric velocities and populate older, higher-mass host galaxies than SNe Ia consistent with Chandrasekhar-mass explosions. Two peculiar events, SN 2006bt and SN 2006ot, have normal peak luminosities but appear to have super-Chandrasekhar ejected masses.
△ Less
Submitted 21 November, 2018;
originally announced November 2018.
-
Validation of a computer code for the energy consumption of a building, with application to optimal electric bill pricing
Authors:
M. Keller,
G. Damblin,
A. Pasanisi,
M. Schuman,
P. Barbillon,
F. Ruggeri,
E. Parent
Abstract:
In this paper, we propose a practical Bayesian framework for the calibration and validation of a computer code, and apply it to a case study concerning the energy consumption forecasting of a building. Validation allows to quantify forecasting uncertainties in view of the code's final use. Here we explore the situation where an energy provider promotes new energy contracts for residential building…
▽ More
In this paper, we propose a practical Bayesian framework for the calibration and validation of a computer code, and apply it to a case study concerning the energy consumption forecasting of a building. Validation allows to quantify forecasting uncertainties in view of the code's final use. Here we explore the situation where an energy provider promotes new energy contracts for residential buildings, tailored to each customer's needs, and including a guarantee of energy performance.
Based on power field measurements, collected from an experimental building cell over a certain time period, the code is calibrated, effectively reducing the epistemic uncertainty affecting some code parameters (here albedo, thermal bridge factor and convective coefficient). Validation is conducted by testing the goodness of fit of the code with respect to field measures, and then by propagating the a posteriori parametric uncertainty through the code, yielding probabilistic forecasts of the average electric power delivered inside the cell over a given time period.
To illustrate the benefits of the proposed Bayesian validation framework, we address the decision problem for an energy supplier offering a new type of contract, wherein the customer pays a fixed fee chosen in advance, based on an overall energy consumption forecast. According to Bayesian decision theory, we show how to choose such a fee optimally from the point of view of the supplier, in order to balance short-terms benefits with customer loyalty.
△ Less
Submitted 21 September, 2018;
originally announced October 2018.
-
The Carnegie Supernova Project: Absolute Calibration and the Hubble Constant
Authors:
Christopher. R. Burns,
Emilie Parent,
M. M. Phillips,
Maximillian Stritzinger,
Kevin Krisciunas,
Nicholas B. Suntzeff,
Eric Y. Hsiao,
Carlos Contreras,
Jorge Anais,
Luis Boldt,
Luis Busta,
Abdo Campillay,
Sergio Castellon,
Gaston Folatelli,
Wendy L. Freedman,
Consuelo Gonzalez,
Mario Hamuy,
Peter Heoflich,
Wojtek Krzeminski,
Barry F. Madore,
Nidia Morrell,
S. E. Persson,
Miguel Roth,
Francisco Salgado,
Jacqueline Seron
, et al. (1 additional authors not shown)
Abstract:
We present an analysis of the final data release of the Carnegie Supernova Project I, focusing on the absolute calibration of the luminosity-decline-rate relation for Type Ia supernovae (SNeIa) using new intrinsic color relations with respect to the color-stretch parameter, $s_{BV}$, enabling improved dust extinction corrections. We investigate to what degree the so-called fast-declining SNeIa can…
▽ More
We present an analysis of the final data release of the Carnegie Supernova Project I, focusing on the absolute calibration of the luminosity-decline-rate relation for Type Ia supernovae (SNeIa) using new intrinsic color relations with respect to the color-stretch parameter, $s_{BV}$, enabling improved dust extinction corrections. We investigate to what degree the so-called fast-declining SNeIa can be used to determine accurate extragalactic distances. We estimate the intrinsic scatter in the luminosity-decline-rate relation, and find it ranges from $\pm 0.13$ mag to $\pm 0.18$ mag with no obvious dependence on wavelength. Using the Cepheid variable star data from the SH0ES project (Riess et al., 2016), the SNIa distance scale is calibrated and the Hubble constant is estimated using our optical and near-infrared sample, and these results are compared to those determined exclusively from a near-infrared sub-sample. The systematic effect of the supernova's host galaxy mass is investigated as a function of wavelength and is found to decrease toward redder wavelengths, suggesting this effect may be due to dust properties of the host. Using estimates of the dust extinction derived from optical and NIR wavelengths, and applying these to H band, we derive a Hubble constant $H_0 = 73.2 \pm 2.3$ km/s/Mpc, whereas using a simple $B-V$ color-correction applied to B band yields $H_0 = 72.7 \pm 2.1$ km/s/Mpc. Photometry of two calibrating SNeIa from the CSP-II sample, SN2012ht and SN2015F, is presented and used to improve the calibration of the SNIa distance ladder.
△ Less
Submitted 19 September, 2018; v1 submitted 17 September, 2018;
originally announced September 2018.
-
PALFA Single-Pulse Pipeline: New Pulsars, Rotating Radio Transients and a Candidate Fast Radio Burst
Authors:
C. Patel,
D. Agarwal,
M. Bhardwaj,
M. M. Boyce,
A. Brazier,
S. Chatterjee,
P. Chawla,
V. M. Kaspi,
D. R. Lorimer,
M. M. McLaughlin,
E. Parent,
Z. Pleunis,
S. M. Ransom,
P. Scholz,
R. S. Wharton,
W. W. Zhu,
M. Alam,
K. Caballero Valdez,
F. Camilo,
J. M. Cordes,
F. Crawford,
J. S. Deneva,
R. D. Ferdman,
P. C. C. Freire,
J. W. T. Hessels
, et al. (4 additional authors not shown)
Abstract:
We present a newly implemented single-pulse pipeline for the PALFA survey to efficiently identify single radio pulses from pulsars, Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs). We have conducted a sensitivity analysis of this new pipeline in which multiple single pulses with a wide range of parameters were injected into PALFA data sets and run through the pipeline. Based on the…
▽ More
We present a newly implemented single-pulse pipeline for the PALFA survey to efficiently identify single radio pulses from pulsars, Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs). We have conducted a sensitivity analysis of this new pipeline in which multiple single pulses with a wide range of parameters were injected into PALFA data sets and run through the pipeline. Based on the recovered pulses, we find that for pulse widths $\rm < 5\ ms$ the sensitivity of the PALFA pipeline is at most a factor of $\rm \sim 2$ less sensitive to single pulses than our theoretical predictions. For pulse widths $\rm > 10\ ms$, as the $\rm DM$ decreases, the degradation in sensitivity gets worse and can increase up to a factor of $\rm \sim 4.5$. Using this pipeline, we have thus far discovered 7 pulsars and 2 RRATs and identified 3 candidate RRATs and 1 candidate FRB. The confirmed pulsars and RRATs have DMs ranging from 133 to 386 pc cm$^{-3}$ and flux densities ranging from 20 to 160 mJy. The pulsar periods range from 0.4 to 2.1 s. We report on candidate FRB 141113, which we argue is likely astrophysical and extragalactic, having $\rm DM \simeq 400\ pc~cm^{-3}$, which represents an excess over the Galactic maximum along this line of sight of $\rm \sim$ 100 - 200 pc cm$^{-3}$. We consider implications for the FRB population and show via simulations that if FRB 141113 is real and extragalactic, the slope $α$ of the distribution of integral source counts as a function of flux density ($N (>S) \propto S^{-α}$) is $1.4 \pm 0.5$ (95% confidence range). However this conclusion is dependent on several assumptions that require verification.
△ Less
Submitted 10 August, 2018;
originally announced August 2018.
-
CaliCo: a R package for Bayesian calibration
Authors:
Mathieu Carmassi,
Pierre Barbillon,
Matthieu Chiodetti,
Merlin Keller,
Eric Parent
Abstract:
In this article, we present a recently released R package for Bayesian calibration. Many industrial fields are facing unfeasible or costly field experiments. These experiments are replaced with numerical/computer experiments which are realized by running a numerical code. Bayesian calibration intends to estimate, through a posterior distribution, input parameters of the code in order to make the c…
▽ More
In this article, we present a recently released R package for Bayesian calibration. Many industrial fields are facing unfeasible or costly field experiments. These experiments are replaced with numerical/computer experiments which are realized by running a numerical code. Bayesian calibration intends to estimate, through a posterior distribution, input parameters of the code in order to make the code outputs close to the available experimental data. The code can be time consuming while the Bayesian calibration implies a lot of code calls which makes studies too burdensome. A discrepancy might also appear between the numerical code and the physical system when facing incompatibility between experimental data and numerical code outputs. The package CaliCo deals with these issues through four statistical models which deal with a time consuming code or not and with discrepancy or not. A guideline for users is provided in order to illustrate the main functions and their arguments. Eventually, a toy example is detailed using CaliCo. This example (based on a real physical system) is in five dimensions and uses simulated data.
△ Less
Submitted 29 August, 2018; v1 submitted 3 August, 2018;
originally announced August 2018.
-
Dating and localizing an invasion from post-introduction data and a coupled reaction-diffusion-absorption model
Authors:
Candy Abboud,
Olivier Bonnefon,
Eric Parent,
Samuel Soubeyrand
Abstract:
Invasion of new territories by alien organisms is of primary concern for environmental and health agencies and has been a core topic in mathematical modeling, in particular in the intents of reconstructing the past dynamics of the alien organisms and predicting their future spatial extents. Partial differential equations offer a rich and flexible modeling framework that has been applied to a large…
▽ More
Invasion of new territories by alien organisms is of primary concern for environmental and health agencies and has been a core topic in mathematical modeling, in particular in the intents of reconstructing the past dynamics of the alien organisms and predicting their future spatial extents. Partial differential equations offer a rich and flexible modeling framework that has been applied to a large number of invasions. In this article, we are specifically interested in dating and localizing the introduction that led to an invasion using mathematical modeling, post-introduction data and an adequate statistical inference procedure. We adopt a mechanistic-statistical approach grounded on a coupled reaction-diffusion-absorption model representing the dynamics of an organism in an heterogeneous domain with respect to growth. Initial conditions (including the date and site of the introduction) and model parameters related to diffusion, reproduction and mortality are jointly estimated in the Bayesian framework by using an adaptive importance sampling algorithm. This framework is applied to the invasion of \textit{Xylella fastidiosa}, a phytopathogenic bacterium detected in South Corsica in 2015, France.
△ Less
Submitted 1 August, 2018;
originally announced August 2018.
-
The implementation of a Fast-Folding pipeline for long-period pulsar searching in the PALFA survey
Authors:
E. Parent,
V. M. Kaspi,
S. M. Ransom,
M. Krasteva,
C. Patel,
P. Scholz,
A. Brazier,
M. A. McLaughlin,
M. Boyce,
W. W. Zhu,
Z. Pleunis,
B. Allen,
S. Bogdanov,
K. Caballero,
F. Camilo,
R. Camuccio,
S. Chatterjee,
J. M. Cordes,
F. Crawford,
J. S. Deneva,
R. Ferdman,
P. C. C. Freire,
J. W. T. Hessels,
F. A. Jenet,
B. Knispel
, et al. (9 additional authors not shown)
Abstract:
The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fo…
▽ More
The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fourier analyses (typically $\sim$16), or due to the strong effect of red noise at low modulation frequencies. We address this reduction in sensitivity by using a time-domain search technique: the Fast-Folding Algorithm (FFA). We designed a program that implements a FFA-based search in the PALFA processing pipeline, and tested the efficiency of the algorithm by performing tests under both ideal, white noise conditions, as well as with real PALFA observational data. In the two scenarios, we show that the time-domain algorithm has the ability to outperform the FFT-based periodicity search implemented in the survey. We perform simulations to compare the previously reported PALFA sensitivity with that obtained using our new FFA implementation. These simulations show that for a pulsar having a pulse duty cycle of roughly 3%, the performance of our FFA pipeline exceeds that of our FFT pipeline for pulses with DM $\lesssim$ 40 pc cm$^{-3}$ and for periods as short as $\sim$500 ms, and that the survey sensitivity is improved by at least a factor of two for periods $\gtrsim$ 6 sec. Discoveries from the implementation of the algorithm in PALFA are also presented in this paper.
△ Less
Submitted 21 May, 2018;
originally announced May 2018.
-
Post-processing multi-ensemble temperature and precipitation forecasts through an Exchangeable Gamma Normal model and its Tobit extension
Authors:
Marie Courbariaux,
Pierre Barbillon,
Luc Perreault,
Éric Parent
Abstract:
Meteorological ensembles are a collection of scenarios for future weather delivered by a meteorological center. Such ensembles form the main source of valuable information for probabilistic forecasting which aims at producing a predictive probability distribution of the quantity of interest instead of a single best guess estimate. Unfortunately, ensembles cannot generally be considered as a sample…
▽ More
Meteorological ensembles are a collection of scenarios for future weather delivered by a meteorological center. Such ensembles form the main source of valuable information for probabilistic forecasting which aims at producing a predictive probability distribution of the quantity of interest instead of a single best guess estimate. Unfortunately, ensembles cannot generally be considered as a sample from such a predictive probability distribution without a preliminary post-processing treatment to calibrate the ensemble. Two main families of post-processing methods, either competing such as BMA or collaborative such as EMOS, can be found in the literature. This paper proposes a mixed effect model belonging to the collaborative family. The structure of the model is based on the hypothesis of invariance under the relabelling of the ensemble members. Its interesting specificities are as follows: 1) exchangeability, which contributes to parsimony, with a latent pivot variable synthesizing the essential meteorological features of the ensembles, 2) a multi-ensemble implementation, allowing to take advantage of various information so as to increase the sharpness of the forecasting procedure. Focus is cast onto Normal statistical structures, first with a direct application for temperatures, then with its Tobit extension for precipitation. Inference is performed by EM algorithms with recourse made to stochastic conditional simulations in the precipitation case. After checking its good behavior on artificial data, the proposed post-processing technique is applied to temperature and precipitation ensemble forecasts produced over five river basins managed by Hydro-Qu$é$bec. These ensemble forecasts were extracted from the THORPEX Interactive Grand Global Ensemble (TIGGE) database. The results indicate that post-processed ensemble are calibrated and generally sharper than the raw ensembles.
△ Less
Submitted 5 March, 2019; v1 submitted 24 April, 2018;
originally announced April 2018.
-
PALFA Discovery of a Highly Relativistic Double Neutron Star Binary
Authors:
K. Stovall,
P. C. C. Freire,
S. Chatterjee,
P. B. Demorest,
D. R. Lorimer,
M. A. McLaughlin,
N. Pol,
J. van Leeuwen,
R. S. Wharton,
B. Allen,
M. Boyce,
A. Brazier,
K. Caballero,
F. Camilo,
R. Camuccio,
J. M. Cordes,
F. Crawford,
J. S. Deneva,
R. D. Ferdman,
J. W. T. Hessels,
F. A. Jenet,
V. M. Kaspi,
B. Knispel,
P. Lazarus,
R. Lynch
, et al. (10 additional authors not shown)
Abstract:
We report the discovery and initial follow-up of a double neutron star (DNS) system, PSR J1946$+$2052, with the Arecibo L-Band Feed Array pulsar (PALFA) survey. PSR J1946$+$2052 is a 17-ms pulsar in a 1.88-hour, eccentric ($e \, =\, 0.06$) orbit with a $\gtrsim 1.2 \, M_\odot$ companion. We have used the Jansky Very Large Array to localize PSR J1946$+$2052 to a precision of 0.09 arcseconds using a…
▽ More
We report the discovery and initial follow-up of a double neutron star (DNS) system, PSR J1946$+$2052, with the Arecibo L-Band Feed Array pulsar (PALFA) survey. PSR J1946$+$2052 is a 17-ms pulsar in a 1.88-hour, eccentric ($e \, =\, 0.06$) orbit with a $\gtrsim 1.2 \, M_\odot$ companion. We have used the Jansky Very Large Array to localize PSR J1946$+$2052 to a precision of 0.09 arcseconds using a new phase binning mode. We have searched multiwavelength catalogs for coincident sources but did not find any counterparts. The improved position enabled a measurement of the spin period derivative of the pulsar ($\dot{P} \, = \, 9\,\pm \, 2 \,\times 10^{-19}$); the small inferred magnetic field strength at the surface ($B_S \, = \, 4 \, \times \, 10^9 \, \rm G$) indicates that this pulsar has been recycled. This and the orbital eccentricity lead to the conclusion that PSR J1946$+$2052 is in a DNS system. Among all known radio pulsars in DNS systems, PSR J1946$+$2052 has the shortest orbital period and the shortest estimated merger timescale, 46 Myr; at that time it will display the largest spin effects on gravitational wave waveforms of any such system discovered to date. We have measured the advance of periastron passage for this system, $\dotω \, = \, 25.6 \, \pm \, 0.3\, °\rm yr^{-1}$, implying a total system mass of only 2.50 $\pm$ 0.04 $M_\odot$, so it is among the lowest mass DNS systems. This total mass measurement combined with the minimum companion mass constrains the pulsar mass to $\lesssim 1.3 \, M_\odot$.
△ Less
Submitted 5 February, 2018;
originally announced February 2018.
-
Bayesian calibration of a numerical code for prediction
Authors:
Mathieu Carmassi,
Pierre Barbillon,
Merlin Keller,
Eric Parent,
Matthieu Chiodetti
Abstract:
Field experiments are often difficult and expensive to make. To bypass these issues, industrial companies have developed computational codes. These codes intend to be representative of the physical system, but come with a certain amount of problems. The code intends to be as close as possible to the physical system. It turns out that, despite continuous code development, the difference between the…
▽ More
Field experiments are often difficult and expensive to make. To bypass these issues, industrial companies have developed computational codes. These codes intend to be representative of the physical system, but come with a certain amount of problems. The code intends to be as close as possible to the physical system. It turns out that, despite continuous code development, the difference between the code outputs and experiments can remain significant. Two kinds of uncertainties are observed. The first one comes from the difference between the physical phenomenon and the values recorded experimentally. The second concerns the gap between the code and the physical system. To reduce this difference, often named model bias, discrepancy, or model error, computer codes are generally complexified in order to make them more realistic. These improvements lead to time consuming codes. Moreover, a code often depends on parameters to be set by the user to make the code as close as possible to field data. This estimation task is called calibration. This paper proposes a review of Bayesian calibration methods and is based on an application case which makes it possible to discuss the various methodological choices and to illustrate their divergences. This example is based on a code used to predict the power of a photovoltaic plant.
△ Less
Submitted 25 March, 2019; v1 submitted 5 January, 2018;
originally announced January 2018.
-
Two long-term intermittent pulsars discovered in the PALFA Survey
Authors:
A. G. Lyne,
B. W. Stappers,
P. C. C. Freire,
J. W. T. Hessels,
V. M. Kaspi,
B. Allen,
S. Bogdanov,
A. Brazier,
F. Camilo,
F. Cardoso,
S. Chatterjee,
J. M. Cordes,
F. Crawford,
J. S. Deneva,
R. D. Ferdman,
F. A. Jenet,
B. Knispel,
P. Lazarus,
J. van Leeuwen,
R. Lynch,
E. Madsen,
M. A. McLaughlin,
E. Parent,
C. Patel,
S. M. Ransom
, et al. (9 additional authors not shown)
Abstract:
We report the discovery of two long-term intermittent radio pulsars in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, extended observations of these pulsars over several years at Jodrell Bank Observatory have revealed the details of their rotation and radiation properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bi-modal intermittency, switch…
▽ More
We report the discovery of two long-term intermittent radio pulsars in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, extended observations of these pulsars over several years at Jodrell Bank Observatory have revealed the details of their rotation and radiation properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bi-modal intermittency, switching between active (ON) and inactive (OFF) emission states and indicating the presence of a large, hitherto unrecognised, underlying population of such objects. For PSR J1929+1357, the initial duty cycle was fON=0.008, but two years later this changed, quite abruptly, to fON=0.16. This is the first time that a significant evolution in the activity of an intermittent pulsar has been seen and we show that the spin-down rate of the pulsar is proportional to the activity. The spin-down rate of PSR J1929+1357 is increased by a factor of 1.8 when it is in active mode, similar to the increase seen in the other three known long-term intermittent pulsars.
△ Less
Submitted 8 December, 2016; v1 submitted 31 August, 2016;
originally announced August 2016.
-
Timing of 29 Pulsars Discovered in the PALFA Survey
Authors:
A. G. Lyne,
B. W. Stappers,
S. Bogdanov,
R. Ferdman,
P. C. C. Freire,
V. M. Kaspi,
B. Knispel,
R. Lynch,
B. Allen,
A. Brazier,
F. Camilo,
F. Cardoso,
S. Chatterjee,
J. M. Cordes,
F. Crawford,
J. S. Deneva,
J. W. T. Hessels,
F. A. Jenet,
P. Lazarus,
J. van Leeuwen,
D. R. Lorimer,
E. Madsen,
J. McKee,
M. A. McLaughlin,
E. Parent
, et al. (11 additional authors not shown)
Abstract:
We report on the discovery and timing observations of 29 distant long-period pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, confirmation and timing observations of these pulsars over several years at Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation and radiation properties. We have used multi…
▽ More
We report on the discovery and timing observations of 29 distant long-period pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, confirmation and timing observations of these pulsars over several years at Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation and radiation properties. We have used multi-frequency data to measure the interstellar scattering properties of some of these pulsars. Most of the pulsars have properties that mirror those of the previously known pulsar population, although four show some notable characteristics. PSRs J1907+0631 and J1925+1720 are young and are associated with supernova remnants or plerionic nebulae: J1907+0631 lies close to the center of SNR G40.5-0.5, while J1925+1720 is coincident with a high-energy Fermi gamma-ray source. One pulsar, J1932+1500, is in a surprisingly eccentric, 199-day binary orbit with a companion having a minimum mass of 0.33 solar masses. Several of the sources exhibit timing noise, and two, PSRs J0611+1436 and J1907+0631, have both suffered large glitches, but with very different post-glitch rotation properties. In particular, the rotational period of PSR J0611+1436 will not recover to its pre-glitch value for about 12 years, a far greater recovery timescale than seen following any other large glitches.
△ Less
Submitted 14 December, 2016; v1 submitted 31 August, 2016;
originally announced August 2016.
-
The repeating Fast Radio Burst FRB 121102: Multi-wavelength observations and additional bursts
Authors:
P. Scholz,
L. G. Spitler,
J. W. T. Hessels,
S. Chatterjee,
J. M. Cordes,
V. M. Kaspi,
R. S. Wharton,
C. G. Bassa,
S. Bogdanov,
F. Camilo,
F. Crawford,
J. Deneva,
J. van Leeuwen,
R. Lynch,
E. C. Madsen,
M. A. McLaughlin,
M. Mickaliger,
E. Parent,
C. Patel,
S. M. Ransom,
A. Seymour,
I. H. Stairs,
B. W. Stappers,
S. P. Tendulkar
Abstract:
We report on radio and X-ray observations of the only known repeating Fast Radio Burst (FRB) source, FRB 121102. We have detected six additional radio bursts from this source: five with the Green Bank Telescope at 2 GHz, and one at 1.4 GHz at the Arecibo Observatory for a total of 17 bursts from this source. All have dispersion measures consistent with a single value ($\sim559$ pc cm$^{-3}$) that…
▽ More
We report on radio and X-ray observations of the only known repeating Fast Radio Burst (FRB) source, FRB 121102. We have detected six additional radio bursts from this source: five with the Green Bank Telescope at 2 GHz, and one at 1.4 GHz at the Arecibo Observatory for a total of 17 bursts from this source. All have dispersion measures consistent with a single value ($\sim559$ pc cm$^{-3}$) that is three times the predicted maximum Galactic value. The 2-GHz bursts have highly variable spectra like those at 1.4 GHz, indicating that the frequency structure seen across the individual 1.4 and 2-GHz bandpasses is part of a wideband process. X-ray observations of the FRB 121102 field with the Swift and Chandra observatories show at least one possible counterpart; however, the probability of chance superposition is high. A radio imaging observation of the field with the Jansky Very Large Array at 1.6 GHz yields a 5$σ$ upper limit of 0.3 mJy on any point-source continuum emission. This upper limit, combined with archival WISE 22-$μ$m and IPHAS H$α$ surveys, rules out the presence of an intervening Galactic HII region. We update our estimate of the FRB detection rate in the PALFA survey to be 1.1$^{+3.7}_{-1.0} \times 10^4$ FRBs sky$^{-1}$ day$^{-1}$ (95% confidence) for peak flux density at 1.4 GHz above 300 mJy. We find that the intrinsic widths of the 12 FRB 121102 bursts from Arecibo are, on average, significantly longer than the intrinsic widths of the 13 single-component FRBs detected with the Parkes telescope.
△ Less
Submitted 18 October, 2016; v1 submitted 29 March, 2016;
originally announced March 2016.
-
Strong near-infrared carbon in the Type Ia supernova iPTF13ebh
Authors:
E. Y. Hsiao,
C. R. Burns,
C. Contreras,
P. Höflich,
D. Sand,
G. H. Marion,
M. M. Phillips,
M. Stritzinger,
S. González-Gaitán,
R. E. Mason,
G. Folatelli,
E. Parent,
C. Gall,
R. Amanullah,
G. C. Anupama,
I. Arcavi,
D. P. K. Banerjee,
Y. Beletsky,
G. A. Blanc,
J. S. Bloom,
P. J. Brown,
A. Campillay,
Y. Cao,
A. De Cia,
T. Diamond
, et al. (29 additional authors not shown)
Abstract:
We present near-infrared (NIR) time-series spectroscopy, as well as complementary ultraviolet (UV), optical, and NIR data, of the Type Ia supernova (SN Ia) iPTF13ebh, which was discovered within two days from the estimated time of explosion. The first NIR spectrum was taken merely 2.3 days after explosion and may be the earliest NIR spectrum yet obtained of a SN Ia. The most striking features in t…
▽ More
We present near-infrared (NIR) time-series spectroscopy, as well as complementary ultraviolet (UV), optical, and NIR data, of the Type Ia supernova (SN Ia) iPTF13ebh, which was discovered within two days from the estimated time of explosion. The first NIR spectrum was taken merely 2.3 days after explosion and may be the earliest NIR spectrum yet obtained of a SN Ia. The most striking features in the spectrum are several NIR C I lines, and the C I λ1.0693 μm line is the strongest ever observed in a SN Ia. Interestingly, no strong optical C II counterparts were found, even though the optical spectroscopic time series began early and is densely-cadenced. Except at the very early epochs, within a few days from the time of explosion, we show that the strong NIR C I compared to the weaker optical C II appears to be general in SNe Ia. iPTF13ebh is a fast decliner with Δm15(B) = 1.79 $\pm$ 0.01, and its absolute magnitude obeys the linear part of the width-luminosity relation. It is therefore categorized as a "transitional" event, on the fast-declining end of normal SNe Ia as opposed to subluminous/91bg-like objects. iPTF13ebh shows NIR spectroscopic properties that are distinct from both the normal and subluminous/91bg-like classes, bridging the observed characteristics of the two classes. These NIR observations suggest composition and density of the inner core similar to that of 91bg-like events, and a deep reaching carbon burning layer not observed in slower declining SNe Ia. There is also a substantial difference between the explosion times inferred from the early-time light curve and the velocity evolution of the Si II λ0.6355 μm line, implying a long dark phase of ~ 4 days.
△ Less
Submitted 8 March, 2015;
originally announced March 2015.
-
Adaptive numerical designs for the calibration of computer codes
Authors:
Guillaume Damblin,
Pierre Barbillon,
Merlin Keller,
Alberto Pasanisi,
Eric Parent
Abstract:
Making good predictions of a physical system using a computer code requires the inputs to be carefully specified. Some of these inputs called control variables have to reproduce physical conditions whereas other inputs, called parameters, are specific to the computer code and most often uncertain. The goal of statistical calibration consists in estimating these parameters with the help of a statis…
▽ More
Making good predictions of a physical system using a computer code requires the inputs to be carefully specified. Some of these inputs called control variables have to reproduce physical conditions whereas other inputs, called parameters, are specific to the computer code and most often uncertain. The goal of statistical calibration consists in estimating these parameters with the help of a statistical model which links the code outputs with the field measurements. In a Bayesian setting, the posterior distribution of these parameters is normally sampled using MCMC methods. However, they are impractical when the code runs are high time-consuming. A way to circumvent this issue consists of replacing the computer code with a Gaussian process emulator, then sampling a cheap-to-evaluate posterior distribution based on it. Doing so, calibration is subject to an error which strongly depends on the numerical design of experiments used to fit the emulator. We aim at reducing this error by building a proper sequential design by means of the Expected Improvement criterion. Numerical illustrations in several dimensions assess the efficiency of such sequential strategies.
△ Less
Submitted 3 April, 2018; v1 submitted 25 February, 2015;
originally announced February 2015.
-
On the Role of Decision Theory in Uncertainty Analysis
Authors:
Merlin Keller,
Eric Parent,
Alberto Pasanisi
Abstract:
Maximum likelihood estimation (MLE) and heuristic predictive estimation (HPE) are two widely used approaches in industrial uncertainty analysis. We review them from the point of view of decision theory, using Bayesian inference as a gold standard for comparison. The main drawback of MLE is that it may fail to properly account for the uncertainty on the physical process generating the data, especia…
▽ More
Maximum likelihood estimation (MLE) and heuristic predictive estimation (HPE) are two widely used approaches in industrial uncertainty analysis. We review them from the point of view of decision theory, using Bayesian inference as a gold standard for comparison. The main drawback of MLE is that it may fail to properly account for the uncertainty on the physical process generating the data, especially when only a small amount of data are available. HPE offers an improvement in that it takes this uncertainty into account. However, we show that this approach is actually equivalent to Bayes estimation for a particular cost function that is not explicitly chosen by the decision maker. This may produce results that are suboptimal from a decisional perspective. These results plead for a systematic use of Bayes estimators based on carefully defined cost functions.
△ Less
Submitted 22 September, 2010;
originally announced September 2010.
-
Random effects compound Poisson model to represent data with extra zeros
Authors:
Marie-Pierre Etienne,
Eric Parent,
Benoit Hugues,
Bernier Jacques
Abstract:
This paper describes a compound Poisson-based random effects structure for modeling zero-inflated data. Data with large proportion of zeros are found in many fields of applied statistics, for example in ecology when trying to model and predict species counts (discrete data) or abundance distributions (continuous data). Standard methods for modeling such data include mixture and two-part conditio…
▽ More
This paper describes a compound Poisson-based random effects structure for modeling zero-inflated data. Data with large proportion of zeros are found in many fields of applied statistics, for example in ecology when trying to model and predict species counts (discrete data) or abundance distributions (continuous data). Standard methods for modeling such data include mixture and two-part conditional models. Conversely to these methods, the stochastic models proposed here behave coherently with regards to a change of scale, since they mimic the harvesting of a marked Poisson process in the modeling steps. Random effects are used to account for inhomogeneity. In this paper, model design and inference both rely on conditional thinking to understand the links between various layers of quantities : parameters, latent variables including random effects and zero-inflated observations. The potential of these parsimonious hierarchical models for zero-inflated data is exemplified using two marine macroinvertebrate abundance datasets from a large scale scientific bottom-trawl survey. The EM algorithm with a Monte Carlo step based on importance sampling is checked for this model structure on a simulated dataset : it proves to work well for parameter estimation but parameter values matter when re-assessing the actual coverage level of the confidence regions far from the asymptotic conditions.
△ Less
Submitted 28 July, 2009;
originally announced July 2009.
-
A model structure a la Thomason on 2-Cat
Authors:
K. Worytkiewicz,
K. Hess,
P. E. Parent,
A. Tonks
Abstract:
We exhibit a model structure on 2-Cat, obtained by transfer from sSet across the adjunction C_2 o Sd^2 -| Ex^2 o N_2.
We exhibit a model structure on 2-Cat, obtained by transfer from sSet across the adjunction C_2 o Sd^2 -| Ex^2 o N_2.
△ Less
Submitted 7 November, 2004;
originally announced November 2004.