Showing 1–50 of 175 results for author: Flynn, C

MeerKAT Pulsar Timing Array parallaxes and proper motions

Authors: Mohsen Shamohammadi, Matthew Bailes, Christopher Flynn, Daniel J. Reardon, Ryan M. Shannon, Sarah Buchner, Andrew D. Cameron, Fernando Camilo, Alessandro Coronigu, Marisa Geyer, Michael Kramer, Matthew Miles, Renee Spiewak

Abstract: We have determined positions, proper motions, and parallaxes of $77$ millisecond pulsars (MSPs) from $\sim3$ years of MeerKAT radio telescope observations. Our timing and noise analyses enable us to measure $35$ significant parallaxes ($12$ of them for the first time) and $69$ significant proper motions. Eight pulsars near the ecliptic have an accurate proper motion in ecliptic longitude only. PSR… ▽ More We have determined positions, proper motions, and parallaxes of $77$ millisecond pulsars (MSPs) from $\sim3$ years of MeerKAT radio telescope observations. Our timing and noise analyses enable us to measure $35$ significant parallaxes ($12$ of them for the first time) and $69$ significant proper motions. Eight pulsars near the ecliptic have an accurate proper motion in ecliptic longitude only. PSR~J0955$-$6150 has a good upper limit on its very small proper motion ($<$0.4 mas yr$^{-1}$). We used pulsars with accurate parallaxes to study the MSP velocities. This yields $39$ MSP transverse velocities, and combined with MSPs in the literature (excluding those in Globular Clusters) we analyse $66$ MSPs in total. We find that MSPs have, on average, much lower velocities than normal pulsars, with a mean transverse velocity of only $78(8)$ km s$^{-1}$ (MSPs) compared with $246(21)$ km s$^{-1}$ (normal pulsars). We found no statistical differences between the velocity distributions of isolated and binary millisecond pulsars. From Galactocentric cylindrical velocities of the MSPs, we derive 3-D velocity dispersions of $σ_ρ$, $σ_φ$, $σ_{z}$ = $63(11)$, $48(8)$, $19(3)$ km s$^{-1}$. We measure a mean asymmetric drift with amplitude $38(11)$ km s$^{-1}$, consistent with expectation for MSPs, given their velocity dispersions and ages. The MSP velocity distribution is consistent with binary evolution models that predict very few MSPs with velocities $>300$ km s$^{-1}$ and a mild anticorrelation of transverse velocity with orbital period. △ Less

Submitted 12 January, 2024; originally announced January 2024.

North-South asymmetries in the Galactic thin disk associated with the vertical phase spiral as seen using LAMOST-Gaia stars

Authors: Jun Lin, Rui Guo, Sarah A. Bird, Haijun Tian, Chao Liu, Chris Flynn, Gaochao Liu, Sheng Cui

Abstract: We select 1,052,469 (754,635) thin disk stars from {\it Gaia} eDR3 and LAMOST DR7 in the range of Galactocentric radius $R$ (guiding center radius $R_\mathrm{g}$) from 8 to 11\,kpc to investigate the asymmetries between the North and South of the disk midplane. More specifically we analyze the vertical velocity dispersion profiles ($σ_{v_{z}}(z$)) in different bins of $R$ ($R_\mathrm{g}$) and… ▽ More We select 1,052,469 (754,635) thin disk stars from {\it Gaia} eDR3 and LAMOST DR7 in the range of Galactocentric radius $R$ (guiding center radius $R_\mathrm{g}$) from 8 to 11\,kpc to investigate the asymmetries between the North and South of the disk midplane. More specifically we analyze the vertical velocity dispersion profiles ($σ_{v_{z}}(z$)) in different bins of $R$ ($R_\mathrm{g}$) and $[\mathrm{Fe/H}]$. We find troughs in the profiles of $σ_{v_{z}}(z)$ located in both the North ($z \sim 0.7$\,kpc) and South ($z \sim -0.5$\,kpc) of the disk at all radial and chemical bins studied. The difference between the Northern and Southern vertical velocity dispersion profiles ($Δσ_{v_{z}}(|z|)$) shows a shift between curves of different $R$ and $R_\mathrm{g}$. A similar shift exists in these NS asymmetry profiles further divided into different $[\mathrm{Fe/H}]$ ranges. The sample binned with $R_\mathrm{g}$ more clearly displays the features in the velocity dispersion profiles. The shift in the peaks of the $Δσ_{v_{z}}$ profiles and the variation in the phase spiral shape binned by metallicity indicate the variation of the vertical potential profiles and the radial metallicity gradient. The wave-like signal in NS asymmetry of $σ_{v_{z}}(z)$ largely originates from phase spiral; while the NS asymmetry profiles of [Fe/H] only display a weak wave-like feature near solar radius. We perform a test particle simulation to qualitatively reproduce the observed results. A quantitative explanation of the NS asymmetry in the metallicity profile needs careful consideration of the spiral shape and the perturbation model, and we leave this for future work. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 13 pages, 12 figures, accepted for publication in MNRAS

The High Time Resolution Universe Pulsar Survey -- XVIII. The reprocessing of the HTRU-S Low Lat survey around the Galactic centre using a Fast Folding Algorithm pipeline for accelerated pulsars

Authors: J. Wongphechauxsorn, D. J. Champion, M. Bailes, V. Balakrishnan, E. D. Barr, M. C. i Bernadich, N. D. R. Bhat, M. Burgay, A. D. Cameron, W. Chen, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, C. Ng, A. Possenti, R. Sengar, R. M. Shannon, B. Stappers, W. van Straten

Abstract: The HTRU-S Low Latitude survey data within 1$^{\circ}$of the Galactic Centre (GC) were searched for pulsars using the Fast Folding Algorithm (FFA). Unlike traditional Fast Fourier Transform (FFT) pipelines, the FFA optimally folds the data for all possible periods over a given range, which is particularly advantageous for pulsars with low-duty cycle. For the first time, a search over acceleration… ▽ More The HTRU-S Low Latitude survey data within 1$^{\circ}$of the Galactic Centre (GC) were searched for pulsars using the Fast Folding Algorithm (FFA). Unlike traditional Fast Fourier Transform (FFT) pipelines, the FFA optimally folds the data for all possible periods over a given range, which is particularly advantageous for pulsars with low-duty cycle. For the first time, a search over acceleration was included in the FFA to improve its sensitivity to binary pulsars. The steps in dispersion measure (DM) and acceleration were optimised, resulting in a reduction of the number of trials by 86 per cent. This was achieved over a search period range from 0.6-s to 432-s, i.e. 10 per cent of the observation time (4320s), with a maximum DM of 4000 pc cm$^{-3}$ and an acceleration range of $\pm 128$m s$^{-2}$. The search resulted in the re-detections of four known pulsars, including a pulsar which was missed in previous FFT processing of this survey. This result indicates that the FFA pipeline is more sensitive than the FFT pipeline used in the previous processing of the survey within our parameter range. Additionally, we discovered a 1.89-s pulsar, PSR J1746-2829, with a large DM, located~0.5 from the GC. Follow-up observations revealed that this pulsar has a relatively flat spectrum($α=-0.9\pm0.1$) and has a period derivative of $\sim1.3\times10^{-12}$ s s$^{-1}$, implying a surface magnetic field of $\sim5.2\times10^{13}$ G and a characteristic age of $\sim23000$ yr. While the period, spectral index, and surface magnetic field strength are similar to many radio magnetars, other characteristics such as high linear polarization are absent. △ Less

Submitted 21 October, 2023; originally announced October 2023.

Comments: 12 pages, 6 figures, 4 tables, Accepted for publication on Monthly Notices of the Royal Astronomical Society

Flux density monitoring of 89 millisecond pulsars with MeerKAT

Authors: P. Gitika, M. Bailes, R. M. Shannon, D. J. Reardon, A. D. Cameron, M. Shamohammadi, M. T. Miles, C. M. L. Flynn, A. Corongiu, M. Kramer

Abstract: We present a flux density study of 89 millisecond pulsars (MSPs) regularly monitored as part of the MeerKAT Pulsar Timing Array (MPTA) using the L-Band receiver with an approximately two week cadence between 2019-2022. For each pulsar, we have determined the mean flux densities at each epoch in eight $\sim$97 MHz sub-bands ranging from 944 to 1625 MHz. From these we have derived their modulation i… ▽ More We present a flux density study of 89 millisecond pulsars (MSPs) regularly monitored as part of the MeerKAT Pulsar Timing Array (MPTA) using the L-Band receiver with an approximately two week cadence between 2019-2022. For each pulsar, we have determined the mean flux densities at each epoch in eight $\sim$97 MHz sub-bands ranging from 944 to 1625 MHz. From these we have derived their modulation indices, their average and peak-to-median flux densities in each sub-band, as well as their mean spectral indices across the entire frequency range. We find that the vast majority of the MSPs have spectra that are well described by a simple power law, with a mean spectral index of -1.86(6). Using the temporal variation of the flux densities we measured the structure functions and determined the refractive scintillation timescale for seven. The structure functions provide strong evidence that the intrinsic radio luminosities of MSPs are stable. As a population, the average modulation index at 20 cm wavelengths peaks near unity at dispersion measures (DMs) of $\sim$20 pc cm$^{-3}$ and by a DM of 100 pc cm$^{-3}$ are closer to 0.2, due to refractive scintillation. We find that timing arrays can improve their observing efficiency by reacting to scintillation maxima, and that 20 cm FRB surveys should prioritise highly scintillating mid-latitude regions of the Galactic sky where they will find $\sim$30% more events and bursts at greater distances. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 16 pages, 8 figures, accepted for publication in MNRAS

Rotational and radio emission properties of PSR J0738-4042 over half a century

Authors: M. E. Lower, S. Johnston, A. Karastergiou, P. R. Brook, M. Bailes, S. Buchner, A. T. Deller, L. Dunn, C. Flynn, M. Kerr, R. N. Manchester, A. Mandlik, L. S. Oswald, A. Parthasarathy, R. M. Shannon, C. Sobey, P. Weltevrede

Abstract: We present a comprehensive study of the rotational and emission properties of PSR J0738$-$4042 using a combination of observations taken by the Deep Space Network, Hartebeesthoek, Parkes (Murriyang) and Molonglo observatories between 1972 and 2023. Our timing of the pulsar is motivated by previously reported profile/spin-down events that occurred in September 2005 and December 2015, which result i… ▽ More We present a comprehensive study of the rotational and emission properties of PSR J0738$-$4042 using a combination of observations taken by the Deep Space Network, Hartebeesthoek, Parkes (Murriyang) and Molonglo observatories between 1972 and 2023. Our timing of the pulsar is motivated by previously reported profile/spin-down events that occurred in September 2005 and December 2015, which result in an anomalously large braking index of $n = 23300 \pm 1800$. Using a Gaussian process regression framework, we develop continuous models for the evolution of the pulsar spin-down rate ($\dotν$) and profile shape. We find that the pulse profile variations are similar regardless of radio observing frequency and polarisation. Small-scale differences can be ascribed to changes in the interstellar medium along the line of sight and frequency-dependent changes in magnetospheric radio emission height. No new correlated spin-down or profile events were identified in our extended dataset. However, we found that the disappearance of a bright emission component in the leading edge of archival profiles between 1981-1988 was not associated with a substantial change in $\dotν$. This marks a notable departure from the previous profile/spin-down events in this pulsar. We discuss the challenges these observations pose for physical models and conclude that interactions between the pulsar and in-falling asteroids or a form of magnetospheric state-switching with a long periodicity are plausible explanations. △ Less

Submitted 24 July, 2023; v1 submitted 21 July, 2023; originally announced July 2023.

Comments: 14 pages, 8 figures. Fixed typo on page 5. Accepted for publication in MNRAS

The Third Fermi Large Area Telescope Catalog of Gamma-ray Pulsars

Authors: David A. Smith, Philippe Bruel, Colin J. Clark, Lucas Guillemot, Matthew T. Kerr, Paul Ray, Soheila Abdollahi, Marco Ajello, Luca Baldini, Jean Ballet, Matthew Baring, Cees Bassa, Josefa Becerra Gonzalez, Ronaldo Bellazzini, Alessandra Berretta, Bhaswati Bhattacharyya, Elisabetta Bissaldi, Raffaella Bonino, Eugenio Bottacini, Johan Bregeon, Marta Burgay, Toby Burnett, Rob Cameron, Fernando Camilo, Regina Caputo , et al. (134 additional authors not shown)

Abstract: We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray M… ▽ More We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 142 pages. Accepted by the Astrophysical Journal Supplement

A Geometric Calibration of the Tip of the Red Giant Branch in the Milky Way using Gaia DR3

Authors: M. Dixon, J. Mould, C. Flynn, E. N. Taylor, C. Lidman, A. R. Duffy

Abstract: We use the latest parallaxes measurements from Gaia DR3 to obtain a geometric calibration of the tip of the red giant branch (TRGB) in Cousins $I$ magnitudes as a standard candle for cosmology. We utilise the following surveys: SkyMapper DR3, APASS DR9, ATLAS Refcat2, and Gaia DR3 synthetic photometry to obtain multiple zero-point calibrations of the TRGB magnitude, $M_{I}^{TRGB}$. Our sample cont… ▽ More We use the latest parallaxes measurements from Gaia DR3 to obtain a geometric calibration of the tip of the red giant branch (TRGB) in Cousins $I$ magnitudes as a standard candle for cosmology. We utilise the following surveys: SkyMapper DR3, APASS DR9, ATLAS Refcat2, and Gaia DR3 synthetic photometry to obtain multiple zero-point calibrations of the TRGB magnitude, $M_{I}^{TRGB}$. Our sample contains Milky Way halo stars at high galactic latitudes ($|b| > 36$) where the impact of metallicity, dust, and crowding are minimised. The magnitude of the TRGB is identified using Sobel edge detection, but this approach introduced a systematic offset. To address this issue, we utilised simulations with PARSEC isochrones and showed how to calibrate and remove this bias. Applying our method within the colour range where the slope of the TRGB is relatively flat for metal-poor halo stars (1.55 $<$ $(BP-RP)$ $<$ 2.25), we find a weighted average $M_{I}^{TRGB} = -4.042 \pm 0.041$ (stat) $\pm0.031$ (sys) mag. A geometric calibration of the Milky Way TRGB has the benefit of being independent of other distance indicators and will help probe systematics in the local distance ladder, leading to improved measurements of the Hubble constant. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 14 pages, 13 figures. Accepted for publication in MNRAS

Spectropolarimetric variability in the repeating fast radio burst source FRB 20180301A

Authors: Pravir Kumar, Rui Luo, Danny C. Price, Ryan M. Shannon, Adam T. Deller, Shivani Bhandari, Yi Feng, Chris Flynn, Jinchen Jiang, Pavan A. Uttarkar, Shuangqiang Wang, Songbo Zhang

Abstract: As the sample size of repeating fast radio bursts (FRBs) has grown, an increasing diversity of phenomenology has emerged. Through long-term multi-epoch studies of repeating FRBs, it is possible to assess which phenomena are common to the population and which are unique to individual sources. We present a multi-epoch monitoring campaign of the repeating FRB source 20180301A using the ultra-wideband… ▽ More As the sample size of repeating fast radio bursts (FRBs) has grown, an increasing diversity of phenomenology has emerged. Through long-term multi-epoch studies of repeating FRBs, it is possible to assess which phenomena are common to the population and which are unique to individual sources. We present a multi-epoch monitoring campaign of the repeating FRB source 20180301A using the ultra-wideband low (UWL) receiver observations with Murriyang, the Parkes 64-m radio telescope. The observations covered a wide frequency band spanning approximately 0.7--4 GHz, and yielded the detection of 46 bursts. None of the repeat bursts displayed radio emission in the range of 1.8--4 GHz, while the burst emission peaked at 1.1 GHz. We discover evidence for secular trends in the burst dispersion measure, indicating a decline at a rate of $-2.7\pm0.2\,{\rm pc\,cm^{-3}\,yr^{-1}}$. We also found significant variation in the Faraday rotation measure of the bursts across the follow-up period, including evidence of a sign reversal. While a majority of bursts did not exhibit any polarization, those that did show a decrease in the linear polarization fraction as a function of frequency, consistent with spectral depolarization due to scattering, as observed in other repeating FRB sources. Surprisingly, no significant variation in the polarization position angles was found, which is in contrast with earlier measurements reported for the FRB source. We measure the burst rate and sub-pulse drift rate variation and compare them with the previous results. These novel observations, along with the extreme polarization properties observed in other repeating FRBs, suggest that a sub-sample of FRB progenitors possess highly dynamic magneto-ionic environments. △ Less

Submitted 15 October, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 21 pages, 14 figures; accepted for publication in MNRAS

Survey for Distant Solar Twins (SDST) -- III. Identification of new solar twin and solar analogue stars

Authors: Christian Lehmann, Michael T. Murphy, Fan Liu, Chris Flynn, Daniel Smith, Daniel A. Berke

Abstract: The Survey for Distant Solar Twins (SDST) aims to find stars very similar to the Sun at distances 1-4 kpc, several times more distant than any currently known solar twins and analogues. The goal is to identify the best stars with which to test whether the fine-structure constant, alpha, varies with dark matter density in our Galaxy. Here we use EPIC, our line-by-line differential technique, to mea… ▽ More The Survey for Distant Solar Twins (SDST) aims to find stars very similar to the Sun at distances 1-4 kpc, several times more distant than any currently known solar twins and analogues. The goal is to identify the best stars with which to test whether the fine-structure constant, alpha, varies with dark matter density in our Galaxy. Here we use EPIC, our line-by-line differential technique, to measure the stellar parameters - effective temperature Teff, surface gravity log g, metallicity [Fe/H] - from moderate resolution (R < 32,000) spectra of 877 solar twin and analogue candidates (547 at 1-4 kpc) observed with the HERMES spectrograph on the Anglo-Australian Telescope. These are consistent with expectations for Teff and log g from photometry, and for [Fe/H] from the Besancon stellar population model. EPIC provides small enough uncertainties (~90 K, 0.08 dex, 0.05 dex, respectively), even at the low signal-to-noise ratios available (S/N >~ 25 per pixel), to identify 299 new solar analogues (> 90% confidence), and 20 solar twins (>50% confidence), 206 and 12 of which are at 1-4 kpc. By extending EPIC to measure line broadening and lithium abundance from HERMES spectra, and with ages derived from isochrone fitting with our stellar parameters, we identify 174 solar analogues at 1-4 kpc which are relatively inactive, slowly rotating, and with no evidence of spectroscopic binarity. These are the preferred targets for follow-up spectroscopy to measure alpha. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 12 pages, 12 figures. Accepted by MNRAS. The full table of stars used in this work can be found at https://doi.org/10.5281/zenodo.7332606

Discovery of 37 new pulsars through GPU-accelerated reprocessing of archival data of the Parkes Multibeam Pulsar Survey

Authors: R. Sengar, M. Bailes, V. Balakrishnan, M. C. i Bernadich, M. Burgay, E. D. Barr, C. M. L. Flynn, R. Shannon, S. Stevenson, J. Wongphechauxsorn

Abstract: We present the discovery of 37 pulsars from $\sim$ 20 years old archival data of the Parkes Multibeam Pulsar Survey using a new FFT-based search pipeline optimised for discovering narrow-duty cycle pulsars. When developing our pulsar search pipeline, we noticed that the signal-to-noise ratios of folded and optimised pulsars often exceeded that achieved in the spectral domain by a factor of two or… ▽ More We present the discovery of 37 pulsars from $\sim$ 20 years old archival data of the Parkes Multibeam Pulsar Survey using a new FFT-based search pipeline optimised for discovering narrow-duty cycle pulsars. When developing our pulsar search pipeline, we noticed that the signal-to-noise ratios of folded and optimised pulsars often exceeded that achieved in the spectral domain by a factor of two or greater, in particular for narrow duty cycle ones. Based on simulations, we verified that this is a feature of search codes that sum harmonics incoherently and found that many promising pulsar candidates are revealed when hundreds of candidates per beam with even with modest spectral signal-to-noise ratios of S/N$\sim$5--6 in higher-harmonic folds (up to 32 harmonics) are folded. Of these candidates, 37 were confirmed as new pulsars and a further 37 would have been new discoveries if our search strategies had been used at the time of their initial analysis. While 19 of these newly discovered pulsars have also been independently discovered in more recent pulsar surveys, 18 are exclusive to only the Parkes Multibeam Pulsar Survey data. Some of the notable discoveries include: PSRs J1635$-$47 and J1739$-$31, which show pronounced high-frequency emission; PSRs J1655$-$40 and J1843$-$08, which belong to the nulling/intermittent class of pulsars; and PSR J1636$-$51, which is an interesting binary system in a $\sim$0.75 d orbit and shows hints of eclipsing behaviour -- unusual given the 340 ms rotation period of the pulsar. Our results highlight the importance of reprocessing archival pulsar surveys and using refined search techniques to increase the normal pulsar population. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 21 pages, 13 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Searches for Shapiro delay in seven binary pulsars using the MeerKAT telescope

Authors: Mohsen Shamohammadi, Matthew Bailes, Paulo C. C. Freire, Aditya Parthasarathy, Daniel J. Reardon, Ryan M. Shannon, Vivek Venkatraman Krishnan, Miquel C. i. Bernadich, Andrew D. Cameron, David J. Champion, Alessandro Corongiu, Christopher Flynn, Marisa Geyer, Michael Kramer, Matthew T. Miles, Andrea Possenti, Renee Spiewak

Abstract: Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope th… ▽ More Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope that show evidence of Shapiro delay: PSRs~J0101$-$6422, J1101$-$6424, J1125$-$6014, J1514$-$4946, J1614$-$2230, J1732$-$5049, and J1909$-$3744. Evidence for Shapiro delay was found in all of the systems, and for three the orientations and data quality enabled strong constraints on their orbital inclinations and component masses. For PSRs~J1125$-$6014, J1614$-$2230 and J1909$-$3744, we determined pulsar masses to be $M_{\rm p} = 1.68\pm 0.17 \, {\rm M_{\odot}} $, $1.94\pm 0.03 \, {\rm M_{\odot}} $ and $1.45 \pm 0.03 \, {\rm M_{\odot}}$, and companion masses to be $M_{\rm c} = 0.33\pm 0.02 \, {\rm M_{\odot}} $, $0.495\pm 0.005 \, {\rm M_{\odot}} $ and $0.205 \pm 0.003 \, {\rm M_{\odot}}$, respectively. This provides the first independent confirmation of PSR~J1614$-$2230's mass, one of the highest known. The Shapiro delays measured for PSRs~J0101$-$6422, J1101$-$6424, J1514$-$4946, and J1732$-$5049 were only weak, and could not provide interesting component mass limits. Despite a large number of millisecond pulsars being routinely timed, relatively few have accurate masses via Shapiro delays. We use simulations to show that this is expected, and provide a formula for observers to assess how accurately a pulsar mass can be determined. We also discuss the observed correlation between pulsar companion masses and spin period, and the anti-correlation between recycled pulsar mass and their companion masses. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: 19 pages, 8 figures

A limit on variations in the fine-structure constant from spectra of nearby Sun-like stars

Authors: Michael T. Murphy, Daniel A. Berke, Fan Liu, Chris Flynn, Christian Lehmann, Vladimir A. Dzuba, Victor V. Flambaum

Abstract: The fine structure constant, $α$, sets the strength of the electromagnetic force. The Standard Model of particle physics provides no explanation for its value, which could potentially vary. The wavelengths of stellar absorption lines depend on $α$, but are subject to systematic effects owing to astrophysical processes in stellar atmospheres. We measured precise line wavelengths using 17 stars, sel… ▽ More The fine structure constant, $α$, sets the strength of the electromagnetic force. The Standard Model of particle physics provides no explanation for its value, which could potentially vary. The wavelengths of stellar absorption lines depend on $α$, but are subject to systematic effects owing to astrophysical processes in stellar atmospheres. We measured precise line wavelengths using 17 stars, selected to have almost identical atmospheric properties to those of the Sun (solar twins), which reduces those systematic effects. We found that $α$ varies by $\lesssim$50 parts-per-billion (ppb) within 50 parsecs from Earth. Combining the results from all 17 stars provides an empirical, local reference for stellar measurements of $α$ with an ensemble precision of 12 ppb. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 33 pages, 6 figures. Published in Science (11 November 2022). This is the accepted version which includes 20 pages of Supplementary Materials

Survey for Distant Solar Twins (SDST) -- II. Design, observations and data

Authors: Fan Liu, Michael T. Murphy, Christian Lehmann, Chris Flynn, Daniel Smith, Janez Kos, Daniel A. Berke, Sarah L. Martell

Abstract: Studies of solar twins have key impacts on the astronomical community, but only $\sim$100--200 nearby solar twins ($<$ 1 kpc) have been reliably identified over the last few decades. The aim of our survey (SDST) is to identify $\sim$150--200 distant solar twins and analogues (up to $\lesssim$ 4 kpc) closer to the Galactic Centre. We took advantage of the precise Gaia and Skymapper surveys to selec… ▽ More Studies of solar twins have key impacts on the astronomical community, but only $\sim$100--200 nearby solar twins ($<$ 1 kpc) have been reliably identified over the last few decades. The aim of our survey (SDST) is to identify $\sim$150--200 distant solar twins and analogues (up to $\lesssim$ 4 kpc) closer to the Galactic Centre. We took advantage of the precise Gaia and Skymapper surveys to select Sun-like candidates in a 2-degree field, which were observed with the HERMES spectrograph on the Anglo-Australian Telescope. We successfully built up the required signal-to-noise ratio (25-per-pixel in the HERMES red band) for most targets as faint as Gaia G of 17.4 mag. The stellar photometric/astrometric parameters (e.g., \teff, \logg, mass) of our candidates are derived in this paper, while the spectroscopic parameters will be presented in the third paper in this SDST series. The selection success rate - the fraction of targets which belong to solar twins or analogues - was estimated from simulated survey data and the Besançon stellar population model, and compared with the actual success rate of the survey. We find that expected and actual success rates agree well, indicating that the numbers of solar twins and analogues we discover in SDST are consistent with expectations, affirming the survey approach. These distant solar analogues are prime targets for testing for any variation in the strength of electromagnetism in regions of higher dark matter density, and can make additional contributions to our understanding of, e.g., Galactic chemical evolution in the inner Milky Way. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 15 pages, 13 figures. Accepted by MNRAS

Probing Galactic variations in the fine-structure constant using solar twin stars: systematic errors

Authors: Daniel A. Berke, Michael T. Murphy, Chris Flynn, Fan Liu

Abstract: Sun-like stars are a new probe of variations in the fine-structure constant, $α$, via the solar twins approach: velocity separations of close pairs of absorption lines are compared between stars with very similar stellar parameters, i.e. effective temperature, metallicity and surface gravity within 100K, 0.1 dex and 0.2 dex of the Sun's values. Here we assess possible systematic errors in this app… ▽ More Sun-like stars are a new probe of variations in the fine-structure constant, $α$, via the solar twins approach: velocity separations of close pairs of absorption lines are compared between stars with very similar stellar parameters, i.e. effective temperature, metallicity and surface gravity within 100K, 0.1 dex and 0.2 dex of the Sun's values. Here we assess possible systematic errors in this approach by analysing $\gtrsim$10,000 archival exposures from the High-Accuracy Radial velocity Planetary Searcher (HARPS) of 130 stars covering a much broader range of stellar parameters. We find that each transition pair's separation shows broad, low-order variations with stellar parameters which can be accurately modelled, leaving only a small residual, intrinsic star-to-star scatter of 0-33 m/s (average $\approx$7 m/s, $\approx$10$^{-4}$Å at 5000Å). This limits the precision available from a single pair in one star. We consider potential systematic errors from a range of instrumental and astrophysical sources (e.g. wavelength calibration, charge transfer inefficiency, stellar magnetic activity, line blending) and conclude that variations in elemental abundances, isotope ratios and stellar rotational velocities may explain this star-to-star scatter. Finally, we find that the solar twins approach can be extended to solar analogues - within 300K, 0.3 dex and 0.4 dex of the Sun's parameters - without significant additional systematic errors, allowing a much larger number of stars to be used as probes of variation in $α$, including at much larger distances. △ Less

Submitted 15 October, 2022; originally announced October 2022.

Comments: 18 pages, 12 figures. Accepted by MNRAS. The VarConLib (Varying Constants Library) software used in this work is available at https://github.com/DBerke/varconlib. The measurements and results in this work are available at https://github.com/DBerke/Berke_et_alia_2022_supplemental_data

Probing Galactic variations in the fine-structure constant using solar twin stars: methodology and results

Authors: Daniel A. Berke, Michael T. Murphy, Chris Flynn, Fan Liu

Abstract: The rich absorption spectra of Sun-like stars are enticing probes for variations in the fine-structure constant, $α$, which gauges the strength of electromagnetism. While individual line wavelengths are sensitive to $α$, they are also sensitive to physical processes in the stellar atmospheres, which has precluded their use so far. Here we demonstrate a new, differential approach using solar twins:… ▽ More The rich absorption spectra of Sun-like stars are enticing probes for variations in the fine-structure constant, $α$, which gauges the strength of electromagnetism. While individual line wavelengths are sensitive to $α$, they are also sensitive to physical processes in the stellar atmospheres, which has precluded their use so far. Here we demonstrate a new, differential approach using solar twins: velocity separations between close pairs of transitions are compared across stars with very similar physical properties, strongly suppressing astrophysical and instrumental systematic errors. We utilise 423 archival exposures of 18 solar twins from the High-Accuracy Radial velocity Planetary Searcher (HARPS), in which calibration errors can be reduced to $\lesssim$3 m/s. For stars with $\approx$10 high signal-to-noise ratio spectra ($\ge$200 per pixel), velocity separations between pairs are measured with $\approx$10 m/s statistical precision. A companion paper assesses a range of systematic error sources using 130 stars, with a greater range of stellar parameters, providing accurate corrections for astrophysical effects and a residual, intrinsic star-to-star scatter of 0-13 m/s. Within these uncertainties, we find no evidence for velocity separation differences in 17 transition pairs between solar twins. In a second companion paper, this is found to limit local ($\lesssim$50 pc) variations in $α$ to $\approx$50 parts per billion, $\sim$2 orders of magnitude less than other Galactic constraints. △ Less

Submitted 15 October, 2022; originally announced October 2022.

Comments: 15 pages, 5 figures. Accepted by MNRAS. The VarConLib (Varying Constants Library) software used in this work is available at https://github.com/DBerke/varconlib. The measurements and results in this work are available at https://github.com/DBerke/Berke_et_alia_2022_supplemental_data

The ultra narrow FRB20191107B, and the origins of FRB scattering

Authors: Vivek Gupta, Chris Flynn, Wael Farah, Matthew Bailes, Adam T. Deller, Cherie K. Day, Marcus E. Lower

Abstract: We report the detection of FRB20191107B with the UTMOST radio telescope at a dispersion measure (DM) of 714.9 ${\rm pc~cm^{-3}}$. The burst consists of three components, the brightest of which has an intrinsic width of only 11.3 $μ$s and a scattering tail with an exponentially decaying time-scale of 21.4 $μ$s measured at 835 MHz. We model the sensitivity of UTMOST and other major FRB surveys to su… ▽ More We report the detection of FRB20191107B with the UTMOST radio telescope at a dispersion measure (DM) of 714.9 ${\rm pc~cm^{-3}}$. The burst consists of three components, the brightest of which has an intrinsic width of only 11.3 $μ$s and a scattering tail with an exponentially decaying time-scale of 21.4 $μ$s measured at 835 MHz. We model the sensitivity of UTMOST and other major FRB surveys to such narrow events. We find that $>60\%$ of FRBs like FRB20191107B are being missed, and that a significant population of very narrow FRBs probably exists and remains underrepresented in these surveys. The high DM and small scattering timescale of FRB20191107B allows us to place an upper limit on the strength of turbulence in the Intergalactic Medium (IGM), quantified as scattering measure (SM), of ${\rm SM_{IGM} < 8.4 \times 10^{-7} ~kpc~m^{-20/3}}$. Almost all UTMOST FRBs have full phase information due to real-time voltage capture which provides us with the largest sample of coherently dedispersed single burst FRBs. Our 10.24 $μ$s time resolution data yields accurately measured FRB scattering timescales. We combine the UTMOST FRBs with 10 FRBs from the literature and find no obvious evidence for a DM-scattering relation, suggesting that IGM is not the dominant source of scattering in FRBs. We support the results of previous studies and identify the local environment of the source in the host galaxy as the most likely region which dominates the observed scattering of our FRBs. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Journal ref: MNRAS 514, 5866-5878 (2022)

Milky Way Mass with K Giants and BHB Stars Using LAMOST, SDSS/SEGUE, and Gaia: 3D Spherical Jeans Equation and Tracer Mass Estimator

Authors: Sarah A. Bird, Xiang-Xiang Xue, Chao Liu, Chris Flynn, Juntai Shen, Jie Wang, Chengqun Yang, Meng Zhai, Ling Zhu, Gang Zhao, Hai-Jun Tian

Abstract: We measure the enclosed Milky Way mass profile to Galactocentric distances of $\sim70$ and $\sim50$ kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are LAMOST and SDSS/SEGUE K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and {\it… ▽ More We measure the enclosed Milky Way mass profile to Galactocentric distances of $\sim70$ and $\sim50$ kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are LAMOST and SDSS/SEGUE K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and {\it Gaia} DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to $\sim70$ and $\sim50$ kpc, the Jeans method yields total enclosed masses of $4.3\pm0.95$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ and $4.1\pm1.2$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ for the KG and BHB stars, respectively. For the KG and BHB samples we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) $\pm0.083$ (systematic) $\times10^{12}$ M$_\odot$ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) $\pm0.15$ (systematic) $\times10^{12}$ M$_\odot$, respectively. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS. Comments welcomed

Using Host Galaxy Spectroscopy to Explore Systematics in the Standardisation of Type Ia Supernovae

Authors: M. Dixon, C. Lidman, J. Mould, L. Kelsey, D. Brout, A. Möller, P. Wiseman, M. Sullivan, L. Galbany, T. M. Davis, M. Vincenzi, D. Scolnic, G. F. Lewis, M. Smith, R. Kessler, A. Duffy, E. Taylor, C. Flynn, T. M. C. Abbott, M. Aguena, S. Allam, F. Andrade-Oliveir, J. Annis, J. Asorey, E. Bertin , et al. (53 additional authors not shown)

Abstract: We use stacked spectra of the host galaxies of photometrically identified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilising full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and prope… ▽ More We use stacked spectra of the host galaxies of photometrically identified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilising full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and properties of the host galaxies that rely on spectral absorption features ($< 1.3σ$), such as stellar population age, metallicity, and mass-to-light ratio. However, we find significant trends between the Hubble residuals and the strengths of [OII] ($4.4σ$) and the Balmer emission lines ($3σ$). These trends are weaker than the well known trend between Hubble residuals and host galaxy stellar mass ($7.2σ$) that is derived from broad band photometry. After light curve corrections, we see fainter SNe Ia residing in galaxies with larger line strengths. We also find a trend (3$σ$) between Hubble residual and the Balmer decrement (a measure of reddening by dust) using H$β$ and H$γ$. The trend, quantified by correlation coefficients, is slightly more significant in the redder SNe Ia, suggesting that bluer SNe Ia are relatively unaffected by dust in the interstellar medium of the host and that dust contributes to current Hubble diagram scatter impacting the measurement of cosmological parameters. △ Less

Submitted 24 October, 2022; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 15 pages, 10 figures. Accepted for publication in MNRAS

Report number: DES-2018-0379 FERMILAB-PUB-22-428-PPD

The High Time Resolution Universe Pulsar Survey -- XVII. PSR J1325-6253, a low eccentricity double neutron star system from an ultra-stripped supernova

Authors: R. Sengar, V. Balakrishnan, S. Stevenson, M. Bailes, E. D. Barr, N. D. R. Bhat, M. Burgay, M. C. i Bernadich, A. D. Cameron, D. J. Champion, W. Chen, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, V. Morello, C. Ng, A. Possenti, B. Stappers, R. M. Shannon, W. van Straten, J. Wongphechauxsorn

Abstract: The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surpr… ▽ More The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surprisingly low eccentricity of just $e = 0.064$. Through 1.4 yr of dedicated timing with the Parkes radio telescope we have been able to measure its rate of advance of periastron, $\dotω=0.138 \pm 0.002$ $\rm deg$ $\rm yr^{-1}$. If this induced $\dotω$ is solely due to general relativity then the total mass of the system is, $M_{\rm sys} = 2.57 \pm 0.06$ M$_{\odot}$. Assuming an edge-on orbit the minimum companion mass is constrained to be $M_\mathrm{c,min}>0.98$ M$_{\odot}$ which implies the pulsar mass is $M_\mathrm{p,max}<1.59 $ M$_{\odot}$. Its location in the $P$-$\dot{P}$ diagram suggests that, like other DNS systems, PSR J1325$-$6253 is a recycled pulsar and if its mass is similar to the known examples ($>1.3$ M$_\odot$), then the companion neutron star is probably less than $\sim1.25$ M$_\odot$ and the system is inclined at about $50^{\circ}$-$60^{\circ}$. The low eccentricity along with the wide orbit of the system strongly favours a formation scenario involving an ultra-stripped supernova explosion. △ Less

Submitted 26 May, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 12 pages, 7 figures, 1 table, accepted for publication in MNRAS

Systematic upper limits on the size of missing pulsar glitches in the first UTMOST open data release

Authors: L. Dunn, A. Melatos, S. Suvorova, W. Moran, R. J. Evans, S. Osłowski, M. E. Lower, M. Bailes, C. Flynn, V. Gupta

Abstract: A systematic, semi-automated search for pulsar glitches in the first UTMOST public data release is presented. The search is carried out using a hidden Markov model which incorporates both glitches and timing noise into the model of the assumed phase evolution of the pulsar. Glitches are detected through Bayesian model selection between models with and without glitches present with minimal human in… ▽ More A systematic, semi-automated search for pulsar glitches in the first UTMOST public data release is presented. The search is carried out using a hidden Markov model which incorporates both glitches and timing noise into the model of the assumed phase evolution of the pulsar. Glitches are detected through Bayesian model selection between models with and without glitches present with minimal human intervention. Nine glitches are detected among seven objects, all of which have been previously reported. No new glitches were detected. Injection studies are used to place 90\% frequentist upper limits on the size of undetected glitches in each of the 282 objects searched. The mean upper limit obtained is $Δf^{90\%}/f = 1.9 \times 10^{-8}$, with a range of $4.1 \times 10^{-11} \leq Δf^{90\%}/f \leq 2.7 \times 10^{-7}$, assuming step events with no post-glitch recoveries. It is demonstrated that including glitch recovery has a mild effect, in most cases increasing the upper limit by a factor of $\lesssim 5$ conservatively assuming complete recovery on a timescale of $100\,\mathrm{d}$. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 16 pages, accepted for publication in MNRAS

Survey for Distant Solar Twins (SDST) -- I. EPIC method for stellar parameter measurement

Authors: Christian Lehmann, Michael T. Murphy, Fan Liu, Chris Flynn, Daniel A. Berke

Abstract: Solar twins are stars of key importance to the field of astronomy and offer a multitude of science applications. Only a small number ($\lesssim200$) of solar twins are known today, all of which are relatively close to our Sun ($\lesssim800\,pc$). The goal of our Survey for Distant Solar Twins (SDST) is to identify many more solar twin and solar analogue stars out to much larger distances (… ▽ More Solar twins are stars of key importance to the field of astronomy and offer a multitude of science applications. Only a small number ($\lesssim200$) of solar twins are known today, all of which are relatively close to our Sun ($\lesssim800\,pc$). The goal of our Survey for Distant Solar Twins (SDST) is to identify many more solar twin and solar analogue stars out to much larger distances ($\sim4\,kpc$). In this paper, we present a new method to identify solar twins using relatively low $S/N$, medium resolving power ($R\sim 28{,}000$) spectra that will be typical of such distant targets observed with HERMES on the $3.9\,m$ Anglo-Australian Telescope (AAT). We developed a novel approach, namely EPIC, to measure stellar parameters which we use to identify stars similar to our Sun. EPIC determines the stellar atmospheric parameters (effective temperature $T_{\mathrm{eff}}$, surface gravity $\log g$ and metallicity [Fe/H]) using differential equivalent width (EW) measurements of selected spectroscopic absorption features and a simple model, trained on previously analysed spectra, that connects these EWs to the stellar parameters. The reference for the EW measurements is a high $S/N$ solar spectrum which is used to minimise several systematic effects. EPIC is fast, optimised for Sun-like stars and yields stellar parameter measurements with small enough uncertainties to enable spectroscopic identification of solar twin and analogue stars up to $\sim4\,kpc$ away using AAT/HERMES, i.e.\ $σ\left(T_{\mathrm{eff}}, \log g, \textrm{[Fe/H]}\right) = \left(50\,K, 0.08\,dex, 0.03\,dex\right)$ on average at $S/N=25$. △ Less

Submitted 13 February, 2022; originally announced February 2022.

Comments: Accepted by MNRAS. The current release of the EPIC algorithm is available at https://dx.doi.org/10.5281/zenodo.6029325

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato , et al. (1636 additional authors not shown)

Abstract: Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational… ▽ More Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets. △ Less

Submitted 27 June, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 37 pages, 9 figures, submitted to ApJ

Report number: LIGO-P2100267

Journal ref: ApJ, 932, 133 (2022)

Detecting Fast Radio Bursts in the Milky Way

Authors: Nayab Gohar, Chris Flynn

Abstract: Fast Radio Bursts (FRBs) are highly energetic transient events with duration of order of microseconds to milliseconds and of unknown origin. They are known to lie at cosmological distances, through localisation to host galaxies. Recently, an FRB-like event was seen from the Milky Way magnetar SGR 1935+2154 by the CHIME and STARE2 telescopes. This is the only magnetar that has produced FRB events i… ▽ More Fast Radio Bursts (FRBs) are highly energetic transient events with duration of order of microseconds to milliseconds and of unknown origin. They are known to lie at cosmological distances, through localisation to host galaxies. Recently, an FRB-like event was seen from the Milky Way magnetar SGR 1935+2154 by the CHIME and STARE2 telescopes. This is the only magnetar that has produced FRB events in our galaxy. Finding similar events in the Milky Way is of great interest to understanding FRB progenitors. Such events will be strongly affected by the turbulent interstellar medium in the Milky Way, their intrinsic energy distribution and their spatial locations within the plane of the Milky Way. We examine these effects using models for the distribution of electrons in the ISM to estimate the dispersion measure and pulse scattering of mock events, and a range of models for the spatial distribution and luminosity functions, including models motivated by the spatial distribution of the Milky Way's magnetars. We evaluate the fraction of FRB events in the Milky Way that are detectable by STARE2 for a range of ISM models, spatial distributions and burst luminosity functions. In all the models examined, only a fraction of burst events are detectable, mainly due to the scattering effects of the ISM. We find that GReX, a proposed all-sky experiment, could increase the detection rate of Milky Way FRB events by an order of magnitude, depending on assumptions made about the luminosity function and scale-height of the FRBs. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: 11 pages, 14 figures, accepted for publication in MNRAS

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1672 additional authors not shown)

Abstract: We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both… ▽ More We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory. △ Less

Submitted 20 July, 2022; v1 submitted 25 November, 2021; originally announced November 2021.

Comments: 37 pages

Report number: LIGO-P2100049

Fast radio bursts as probes of feedback from active galactic nuclei

Authors: Adam J. Batten, Alan R. Duffy, Chris Flynn, Vivek Gupta, Emma Ryan-Weber, Nastasha Wijers

Abstract: Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the… ▽ More Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the EAGLE simulation series. We measured the mean DM-redshift relation, and the scatter around it, using $2.56 \times 10^8$ sightlines at 131 redshift ($z$) bins between $0 \leq z \leq 3$. While the DM-redshift relation itself is highly robust against different AGN feedback models, significant differences are detected in the scatter around the mean: weaker feedback leads to more scatter. We find $\sim 10^4$ localised FRBs are needed to discriminate between the scatter in standard feedback and stronger, more intermittent feedback models. The number of FRBs required is dependent on the redshift distribution of the detected population. A log-normal redshift distribution at $z=0.5$ requires approximately 50% fewer localised FRBs than a distribution centred at $z=1$. With the Square Kilometre Array expected to detect $>10^3$ FRBs per day, in the future, FRBs will be able to provide constraints on AGN feedback. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 5 pages, 4 figures; Submitted to MNRAS Letters. Comments welcome

Circularly polarized radio emission from the repeating fast radio burst source FRB 20201124A

Authors: Pravir Kumar, Ryan M. Shannon, Marcus E. Lower, Shivani Bhandari, Adam T. Deller, Chris Flynn, Evan F. Keane

Abstract: The mechanism that produces fast radio burst (FRB) emission is poorly understood. Targeted monitoring of repeating FRB sources provides the opportunity to fully characterize the emission properties in a manner impossible with one-off bursts. Here, we report observations of the source of FRB 20201124A, with the Australian Square Kilometre Array Pathfinder (ASKAP) and the ultra-wideband low (UWL) re… ▽ More The mechanism that produces fast radio burst (FRB) emission is poorly understood. Targeted monitoring of repeating FRB sources provides the opportunity to fully characterize the emission properties in a manner impossible with one-off bursts. Here, we report observations of the source of FRB 20201124A, with the Australian Square Kilometre Array Pathfinder (ASKAP) and the ultra-wideband low (UWL) receiver at the Parkes 64-m radio telescope (Murriyang). The source entered a period of emitting bright bursts during early 2021 April. We have detected 16 bursts from this source. One of the bursts detected with ASKAP is the brightest burst ever observed from a repeating FRB source with an inferred fluence of $640\pm70$ Jy ms. Of the five bursts detected with the Parkes UWL, none display any emission in the range 1.1--4 GHz. All UWL bursts are highly polarized, with their Faraday rotation measures (RMs) showing apparent variations. We obtain an average RM of $-614$ rad m$^{-2}$ for this FRB source with a standard deviation of $16$ rad m$^{-2}$ in the UWL bursts. In one of the UWL bursts, we see evidence of significant circularly polarized emission with a fractional extent of $47\pm1$ per cent. Such a high degree of circular polarization has never been seen before in bursts from repeating FRB sources. We also see evidence for significant variation in the linear polarization position angle in the pulse profile of this UWL repeat burst. Models for repeat burst emission will need to account for the increasing diversity in the burst polarization properties. △ Less

Submitted 31 March, 2022; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: 14 pages, 8 figures; accepted for publication in MNRAS

Characterizing the FRB host galaxy population and its connection to transients in the local and extragalactic Universe

Authors: Shivani Bhandari, Kasper E. Heintz, Kshitij Aggarwal, Lachlan Marnoch, Cherie K. Day, Jessica Sydnor, Sarah Burke-Spolaor, Casey J. Law, J. Xavier Prochaska, Nicolas Tejos, Keith W. Bannister, Bryan J. Butler, Adam T. Deller, R. D. Ekers, Chris Flynn, Wen-fai Fong, Clancy W. James, T. Joseph W. Lazio, Rui Luo, Elizabeth K. Mahony, Stuart D. Ryder, Elaine M. Sadler, Ryan M. Shannon, JinLin Han, Kejia Lee , et al. (1 additional authors not shown)

Abstract: We present the localization and host galaxies of one repeating and two apparently non-repeating Fast Radio Bursts. FRB20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at $z=0.3304$. FRB20191228A, and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at $z=0.2430$ and $z=0.3688$, respec… ▽ More We present the localization and host galaxies of one repeating and two apparently non-repeating Fast Radio Bursts. FRB20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at $z=0.3304$. FRB20191228A, and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at $z=0.2430$ and $z=0.3688$, respectively. We combine these with 13 other well-localized FRBs in the literature, and analyze the host galaxy properties. We find no significant differences in the host properties of repeating and apparently non-repeating FRBs. FRB hosts are moderately star-forming, with masses slightly offset from the star-forming main-sequence. Star formation and low-ionization nuclear emission-line region (LINER) emission are major sources of ionization in FRB host galaxies, with the former dominant in repeating FRB hosts. FRB hosts do not track stellar mass and star formation as seen in field galaxies (more than 95% confidence). FRBs are rare in massive red galaxies, suggesting that progenitor formation channels are not solely dominated by delayed channels which lag star formation by Gigayears. The global properties of FRB hosts are indistinguishable from core-collapse supernovae (CCSNe) and short gamma-ray bursts (SGRBs) hosts, and the spatial offset (from galaxy centers) of FRBs is mostly inconsistent with that of the Galactic neutron star population (95% confidence). The spatial offsets of FRBs (normalized to the galaxy effective radius) also differ from those of globular clusters (GCs) in late- and early-type galaxies with 95% confidence. △ Less

Submitted 16 November, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: 27 pages, 12 figures, 8 tables. Accepted for publication in the Astrophysical Journal

A broadband radio view of transient jet ejecta in the black hole candidate X-ray binary MAXI J1535-571

Authors: J. Chauhan, J. C. A. Miller-Jones, G. E. Anderson, A. Paduano, M. Sokolowski, C. Flynn, P. J. Hancock, N. Hurley-Walker, D. L. Kaplan, T. D. Russell, A. Bahramian, S. W. Duchesne, D. Altamirano, S. Croft, H. A. Krimm, G. R. Sivakoff, R. Soria, C. M. Trott, R. B. Wayth, V. Gupta, M. Johnston-Hollitt, S. J. Tingay

Abstract: We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535-571, which underwent a prolonged outburst beginning on 2 September 2017. We monitored MAXI J1535-571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 20 September to 14 October 2017. The source was quasi-simultaneously observed over… ▽ More We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535-571, which underwent a prolonged outburst beginning on 2 September 2017. We monitored MAXI J1535-571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 20 September to 14 October 2017. The source was quasi-simultaneously observed over the frequency range 0.84-19 GHz by UTMOST (the upgraded Molonglo Observatory Synthesis Telescope), the Australian Square Kilometre Array Pathfinder, the Australia Telescope Compact Array (ATCA), and the Australian Long Baseline Array (LBA). Using the LBA observations from 23 September 2017, we measured the source size to be $34\pm1$ mas. During the brightest radio flare on 21 September 2017, the source was detected down to 119 MHz by the MWA, and the radio spectrum indicates a turnover between 250 and 500 MHz, which is most likely due to synchrotron self-absorption (SSA). By fitting the radio spectrum with a SSA model and using the LBA size measurement, we determined various physical parameters of the jet knot (identified in ATCA data), including the jet opening angle (= $4.5\pm1.2^{\circ}$) and the magnetic field strength (= $104^{+80}_{-78}$ mG). Our fitted magnetic field strength agrees reasonably well with that inferred from the standard equipartition approach, suggesting the jet knot to be close to equipartition. Our study highlights the capabilities of the Australian suite of radio telescopes to jointly probe radio jets in black hole X-ray binaries (BH-XRBs) via simultaneous observations over a broad frequency range, and with differing angular resolutions. This suite allows us to determine the physical properties of XRB jets. Finally, our study emphasizes the potential contributions that can be made by the low-frequency part of the Square Kilometre Array (SKA-Low) in the study of BH-XRBs. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: This article has been accepted for publication in PASA. This article has 16 pages, 7 figures and 1 table

A comparison of Galactic electron density models using PyGEDM

Authors: Danny C. Price, Adam Deller, Chris Flynn

Abstract: Galactic electron density distribution models are crucial tools for estimating the impact of the ionised interstellar medium on the impulsive signals from radio pulsars and fast radio bursts. The two prevailing Galactic electron density models are YMW16 (Yao et al., 2017) and NE2001 (Cordes & Lazio, 2002). Here, we introduce a software package PyGEDM which provides a unified application programmin… ▽ More Galactic electron density distribution models are crucial tools for estimating the impact of the ionised interstellar medium on the impulsive signals from radio pulsars and fast radio bursts. The two prevailing Galactic electron density models are YMW16 (Yao et al., 2017) and NE2001 (Cordes & Lazio, 2002). Here, we introduce a software package PyGEDM which provides a unified application programming interface (API) for these models and the YT20 (Yamasaki & Totani, 2020) model of the Galactic halo. We use PyGEDM to compute all-sky maps of Galactic dispersion measure (DM) for YMW16 and NE2001, and compare the large-scale differences between the two. In general, YMW16 predicts higher DM values toward the Galactic anticentre. YMW16 predicts higher DMs at low Galactic latitudes, but NE2001 predicts higher DMs in most other directions. We identify lines of sight for which the models are most discrepant, using pulsars with independent distance measurements. YMW16 performs better on average than NE2001, but both models show significant outliers. We suggest that future campaigns to determine pulsar distances should focus on targets where the models show large discrepancies, so future models can use those measurements to better estimate distances along those line of sight. We also suggest that the Galactic halo should be considered as a component in future GEDMs, to avoid overestimating the Galactic DM contribution for extragalactic sources such as FRBs. △ Less

Submitted 30 June, 2021; originally announced June 2021.

Comments: 13 pages, 7 figures. Accepted to PASA

The Deeper, Wider, Faster Program: Exploring stellar flare activity with deep, fast cadenced DECam imaging via machine learning

Authors: Sara Webb, Chris Flynn, Jeff Cooke, Jielai Zhang, Ashish Mahabal, Tim Abbott, Rebecca Allen, Igor Andreoni, Sarah Bird, Simon Goode, Michelle Lochner, Tyler Pritchard

Abstract: We present our 500 pc distance-limited study of stellar fares using the Dark Energy Camera as part of the Deeper, Wider, Faster Program. The data was collected via continuous 20-second cadence g band imaging and we identify 19,914 sources with precise distances from Gaia DR2 within twelve, ~3 square-degree, fields over a range of Galactic latitudes. An average of ~74 minutes is spent on each field… ▽ More We present our 500 pc distance-limited study of stellar fares using the Dark Energy Camera as part of the Deeper, Wider, Faster Program. The data was collected via continuous 20-second cadence g band imaging and we identify 19,914 sources with precise distances from Gaia DR2 within twelve, ~3 square-degree, fields over a range of Galactic latitudes. An average of ~74 minutes is spent on each field per visit. All light curves were accessed through a novel unsupervised machine learning technique designed for anomaly detection. We identify 96 flare events occurring across 80 stars, the majority of which are M dwarfs. Integrated are energies range from $\sim 10^{31}-10^{37}$ erg, with a proportional relationship existing between increased are energy with increased distance from the Galactic plane, representative of stellar age leading to declining yet more energetic are events. In agreement with previous studies we observe an increase in flaring fraction from M0 -> M6 spectral types. Furthermore, we find a decrease in the flaring fraction of stars as vertical distance from the galactic plane is increased, with a steep decline present around ~100 pc. We find that ~70% of identified flares occur on short timescales of ~8 minutes. Finally we present our associated are rates, finding a volumetric rate of $2.9 \pm 0.3 \times 10^{-6}$ flares pc$^{-3}$ hr$^{-1}$. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 20 pages, 14 figures, 6 tables

Constraining bright optical counterparts of Fast Radio Bursts

Authors: C. Nuñez, N. Tejos, G. Pignata, C. D. Kilpatrick, J. X. Prochaska, K. E. Heintz, K. W. Bannister, S. Bhandari, C. K. Day, A. T. Deller, C. Flynn, E. K. Mahony, D. Majewski, L. Marnoch, H. Qiu, S. D. Ryder, R. M. Shannon

Abstract: Fast Radio Bursts (FRBs) are extremely energetic pulses of millisecond duration and unknown origin. In order to understand the phenomenon that emits these pulses, targeted and untargeted searches have been performed for multi-wavelength counterparts, including the optical. The objective of this work is to search for optical transients at the position of 8 well-localized FRBs, after the arrival of… ▽ More Fast Radio Bursts (FRBs) are extremely energetic pulses of millisecond duration and unknown origin. In order to understand the phenomenon that emits these pulses, targeted and untargeted searches have been performed for multi-wavelength counterparts, including the optical. The objective of this work is to search for optical transients at the position of 8 well-localized FRBs, after the arrival of the burst on different time-scales (typically at one day, several months, and one year after FRB detection) in order to compare with known transient optical light curves. We used the Las Cumbres Observatory Global Telescope Network (LCOGT), which allows us to promptly take images owing to its network of twenty-three telescopes working around the world. We used a template subtraction technique on all the images we collected at different epochs. We have divided the subtractions into two groups, in one group we use the image of the last epoch as a template and in the other group we use the image of the first epoch as a template. We have searched for bright optical transients at the localizations of the FRBs (<1 arcsec) in the template subtracted images. We have found no optical transients, so we have set limiting magnitudes of optical counterparts. Typical limiting magnitudes in apparent (absolute) magnitudes for our LCOGT data are ~22 (-19) mag in the r-band. We have compared our limiting magnitudes with light curves of superluminous supernovae (SLSNe), type Ia supernovae (SNe), supernovae associated with gamma-ray bursts (GRB SNe), a kilonova, and tidal disruption events (TDEs). We rule out that FRBs are associated with SLSN at a confidence of ~99.9%. We can also rule out the brightest sub-types of type Ia SNe, GRB SNe and TDEs (under some conditions) at similar confidence, though we cannot exclude scenarios where FRBs are associated with the faintest sub-type of each of these transient classes. △ Less

Submitted 21 July, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: Accepted to Astronomy & Astrophysics on 27 June, 2021. 8 pages, 3 figures, 3 tables

Journal ref: A&A 653, A119 (2021)

The Cosmic Dispersion Measure in the EAGLE Simulations

Authors: Adam J. Batten, Alan R. Duffy, Nastasha Wijers, Vivek Gupta, Chris Flynn, Joop Schaye, Emma Ryan-Weber

Abstract: The dispersion measure (DM) of fast radio bursts (FRBs) provides a unique way to probe ionised baryons in the intergalactic medium (IGM). Cosmological models with different parameters lead to different DM-redshift ($\mathrm{DM}-z$) relations. Additionally, the over/under-dense regions in the IGM and the circumgalactic medium of intervening galaxies lead to scatter around the mean $\mathrm{DM}-z$ r… ▽ More The dispersion measure (DM) of fast radio bursts (FRBs) provides a unique way to probe ionised baryons in the intergalactic medium (IGM). Cosmological models with different parameters lead to different DM-redshift ($\mathrm{DM}-z$) relations. Additionally, the over/under-dense regions in the IGM and the circumgalactic medium of intervening galaxies lead to scatter around the mean $\mathrm{DM}-z$ relations. We have used the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations to measure the mean $\mathrm{DM}-z$ relation and the scatter around it using over one billion lines-of-sight between redshifts $0<z<3$. We investigated two techniques to estimate line-of-sight DM: `pixel scrambling' and `box transformations'. We find that using box transformations (a technique from the literature) causes strong correlations due to repeated replication of structure. Comparing a linear and non-linear model, we find that the non-linear model with cosmological parameters, provides a better fit to the $\mathrm{DM}-z$ relation. The differences between these models are the most significant at low redshifts ($z<0.5$). The scatter around the $\mathrm{DM}-z$ relation is highly asymmetric, especially at low redshift $\left(z<0.5\right)$, and becomes more Gaussian as redshift approaches $z\sim3$, the limit of this study. The increase in Gaussianity with redshift is indicative of the large scale structures that is better probed with longer lines-of-sight. The minimum simulation size suitable for investigations into the scatter around the $\mathrm{DM}-z$ relation is 100~comoving~Mpc. The $\mathrm{DM}-z$ relation measured in EAGLE is available with an easy-to-use python interface in the open-source FRB redshift estimation package FRUITBAT. △ Less

Submitted 29 November, 2020; originally announced November 2020.

Comments: Submitted to MNRAS. 14 pages, 15 figures, Comments Welcome. Data and model will be released upon being accepted

Estimating fast transient detection pipeline efficiencies at UTMOST via real-time injection of mock FRBs

Authors: Vivek Gupta, Chris Flynn, Wael Farah, Andrew Jameson, Vivek Venkatraman Krishnan, Matthew Bailes, Timothy Bateman, Adam T. Deller, Ayushi Mandlik, Angus Sutherland

Abstract: Dedicated surveys using different detection pipelines are being carried out at multiple observatories to find more Fast Radio Bursts (FRBs). Understanding the efficiency of detection algorithms and the survey completeness function is important to enable unbiased estimation of the underlying FRB population properties. One method to achieve end-to-end testing of the system is by injecting mock FRBs… ▽ More Dedicated surveys using different detection pipelines are being carried out at multiple observatories to find more Fast Radio Bursts (FRBs). Understanding the efficiency of detection algorithms and the survey completeness function is important to enable unbiased estimation of the underlying FRB population properties. One method to achieve end-to-end testing of the system is by injecting mock FRBs in the live data-stream and searching for them blindly. Mock FRB injection is particularly effective for machine-learning-based classifiers, for which analytic characterisation is impractical. We describe a first-of-its-kind implementation of a real-time mock FRB injection system at the upgraded Molonglo Observatory Synthesis Telescope (UTMOST) and present our results for a set of 20,000 mock FRB injections. The injections have yielded clear insight into the detection efficiencies and have provided a survey completeness function for pulse width, fluence and DM. Mock FRBs are recovered with uniform efficiency over the full range of injected DMs, however the recovery fraction is found to be a strong function of the width and Signal-to-Noise (SNR). For low widths ($\lesssim 20$ ms) and high SNR ($\gtrsim$ 9) the recovery is highly effective with recovery fractions exceeding 90%. We find that the presence of radio frequency interference causes the recovered SNR values to be systematically lower by up to 20% compared to the injected values. We find that wider FRBs become increasingly hard to recover for the machine-learning-based classifier employed at UTMOST. We encourage other observatories to implement live injection set-ups for similar testing of their surveys. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Comments: 12 pages, 9 figures. Accepted for publication in MNRAS

Extremely band-limited repetition from a fast radio burst source

Authors: Pravir Kumar, Ryan M. Shannon, Chris Flynn, Stefan Osłowski, Shivani Bhandari, Cherie K. Day, Adam T. Deller, Wael Farah, Jane F. Kaczmarek, Matthew Kerr, Chris Phillips, Danny C. Price, Hao Qiu, Nithyanandan Thyagarajan

Abstract: The fast radio burst (FRB) population is observationally divided into sources that have been observed to repeat and those that have not. There is tentative evidence that the bursts from repeating sources have different properties than the non-repeating ones. In order to determine the occurrence rate of repeating sources and characterize the nature of repeat emission, we have been conducting sensit… ▽ More The fast radio burst (FRB) population is observationally divided into sources that have been observed to repeat and those that have not. There is tentative evidence that the bursts from repeating sources have different properties than the non-repeating ones. In order to determine the occurrence rate of repeating sources and characterize the nature of repeat emission, we have been conducting sensitive searches for repetitions from bursts detected with the Australian Square Kilometre Array Pathfinder (ASKAP) with the 64-m Parkes radio telescope, using the recently commissioned Ultra-wideband Low (UWL) receiver system, over a band spanning 0.7$-$4.0 GHz. We report the detection of a repeat burst from the source of FRB 20190711A. The detected burst is 1 ms wide and has a bandwidth of just 65 MHz. We find no evidence of any emission in the remaining part of the 3.3 GHz UWL band. While the emission bandwidths of the ASKAP and UWL bursts show $ν^{-4}$ scaling consistent with a propagation effect, the spectral occupancy is inconsistent with diffractive scintillation. This detection rules out models predicting broad-band emission from the FRB 20190711A source and puts stringent constraints on the emission mechanism. The low spectral occupancy highlights the importance of sub-banded search methods in detecting FRBs. △ Less

Submitted 7 December, 2020; v1 submitted 2 September, 2020; originally announced September 2020.

Comments: 7 pages, 5 figures; accepted for publication in MNRAS

Limits on precursor and afterglow radio emission from a fast radio burst in a star-forming galaxy

Authors: Shivani Bhandari, Keith W. Bannister, Emil Lenc, Hyerin Cho, Ron Ekers, Cherie K. Day, Adam T. Deller, Chris Flynn, Clancy W. James, Jean-Pierre Macquart, Elizabeth K. Mahony, Lachlan Marnoch, Vanessa A. Moss, Chris Phillips, J. Xavier Prochaska, Hao Qiu, Stuart D. Ryder, Ryan M. Shannon, Nicolas Tejos, O. Ivy Wong

Abstract: We present a new fast radio burst (FRB) at 920 MHz discovered during commensal observations conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) as part of the Commensal Real-time ASKAP Fast Transients (CRAFT) survey. FRB 191001 was detected at a dispersion measure (DM) of 506.92(4) pc cm$^{-3}$ and its measured fluence of 143(15) Jy ms is the highest of the bursts localized to… ▽ More We present a new fast radio burst (FRB) at 920 MHz discovered during commensal observations conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) as part of the Commensal Real-time ASKAP Fast Transients (CRAFT) survey. FRB 191001 was detected at a dispersion measure (DM) of 506.92(4) pc cm$^{-3}$ and its measured fluence of 143(15) Jy ms is the highest of the bursts localized to host galaxies by ASKAP to date. The subarcsecond localization of the FRB provided by ASKAP reveals that the burst originated in the outskirts of a highly star-forming spiral in a galaxy pair at redshift $z=0.2340(1)$. Radio observations show no evidence for a compact persistent radio source associated with the FRB 191001 above a flux density of $15 μ$Jy. However, we detect diffuse synchrotron radio emission from the disk of the host galaxy that we ascribe to ongoing star formation. FRB 191001 was also detected as an image-plane transient in a single 10 s snapshot with a flux density of 19.3 mJy in the low-time-resolution visibilities obtained simultaneously with CRAFT data. The commensal observation facilitated a search for repeating and slowly varying radio emissions 8 hr before and 1 hr after the burst. We found no variable radio emission on timescales ranging from 1 ms to 1.4 hr. We report our upper limits and briefly review FRB progenitor theories in the literature that predict radio afterglows. Our data are still only weakly constraining of any afterglows at the redshift of the FRB. Future commensal observations of more nearby and bright FRBs will potentially provide stronger constraints. △ Less

Submitted 11 October, 2020; v1 submitted 28 August, 2020; originally announced August 2020.

Comments: 12 pages, 6 figures, Published in ApJ Letters. We have corrected the R.A. uncertainty in the position of FRB 191001 from 0.006s to 0.02s

Journal ref: 2020ApJ...901L..20B

A magnetar parallax

Authors: H. Ding, A. T. Deller, M. E. Lower, C. Flynn, S. Chatterjee, W. Brisken, N. Hurley-Walker, F. Camilo, J. Sarkissian, V. Gupta

Abstract: XTE J1810-197 (J1810) was the first magnetar identified to emit radio pulses, and has been extensively studied during a radio-bright phase in 2003$-$2008. It is estimated to be relatively nearby compared to other Galactic magnetars, and provides a useful prototype for the physics of high magnetic fields, magnetar velocities, and the plausible connection to extragalactic fast radio bursts. Upon the… ▽ More XTE J1810-197 (J1810) was the first magnetar identified to emit radio pulses, and has been extensively studied during a radio-bright phase in 2003$-$2008. It is estimated to be relatively nearby compared to other Galactic magnetars, and provides a useful prototype for the physics of high magnetic fields, magnetar velocities, and the plausible connection to extragalactic fast radio bursts. Upon the re-brightening of the magnetar at radio wavelengths in late 2018, we resumed an astrometric campaign on J1810 with the Very Long Baseline Array, and sampled 14 new positions of J1810 over 1.3 years. The phase calibration for the new observations was performed with two phase calibrators that are quasi-colinear on the sky with J1810, enabling substantial improvement of the resultant astrometric precision. Combining our new observations with two archival observations from 2006, we have refined the proper motion and reference position of the magnetar and have measured its annual geometric parallax, the first such measurement for a magnetar. The parallax of $0.40\pm0.05\,$mas corresponds to a most probable distance $2.5^{+0.4}_{-0.3}\,$kpc for J1810. Our new astrometric results confirm an unremarkable transverse peculiar velocity of $\approx200\,\mathrm{km~s^{-1}}$ for J1810, which is only at the average level among the pulsar population. The magnetar proper motion vector points back to the central region of a supernova remnant (SNR) at a compatible distance at $\approx70\,$kyr ago, but a direct association is disfavored by the estimated SNR age of ~3 kyr. △ Less

Submitted 14 August, 2020; originally announced August 2020.

Comments: 10 pages, 3 figures, accepted for publication in MNRAS

Unsupervised machine learning for transient discovery in Deeper, Wider, Faster light curves

Authors: Sara Webb, Michelle Lochner, Daniel Muthukrishna, Jeff Cooke, Chris Flynn, Ashish Mahabal, Simon Goode, Igor Andreoni, Tyler Pritchard, Timothy M. C. Abbott

Abstract: Identification of anomalous light curves within time-domain surveys is often challenging. In addition, with the growing number of wide-field surveys and the volume of data produced exceeding astronomers ability for manual evaluation, outlier and anomaly detection is becoming vital for transient science. We present an unsupervised method for transient discovery using a clustering technique and the… ▽ More Identification of anomalous light curves within time-domain surveys is often challenging. In addition, with the growing number of wide-field surveys and the volume of data produced exceeding astronomers ability for manual evaluation, outlier and anomaly detection is becoming vital for transient science. We present an unsupervised method for transient discovery using a clustering technique and the Astronomaly package. As proof of concept, we evaluate 85553 minute-cadenced light curves collected over two 1.5 hour periods as part of the Deeper, Wider, Faster program, using two different telescope dithering strategies. By combining the clustering technique HDBSCAN with the isolation forest anomaly detection algorithm via the visual interface of Astronomaly, we are able to rapidly isolate anomalous sources for further analysis. We successfully recover the known variable sources, across a range of catalogues from within the fields, and find a further 7 uncatalogued variables and two stellar flare events, including a rarely observed ultra fast flare (5 minute) from a likely M-dwarf. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: Accepted 7 Aug 2020, 19 pages, 8 figures,

Neutron Star Extreme Matter Observatory: A kilohertz-band gravitational-wave detector in the global network

Authors: K. Ackley, V. B. Adya, P. Agrawal, P. Altin, G. Ashton, M. Bailes, E. Baltinas, A. Barbuio, D. Beniwal, C. Blair, D. Blair, G. N. Bolingbroke, V. Bossilkov, S. Shachar Boublil, D. D. Brown, B. J. Burridge, J. Calderon Bustillo, J. Cameron, H. Tuong Cao, J. B. Carlin, S. Chang, P. Charlton, C. Chatterjee, D. Chattopadhyay, X. Chen , et al. (139 additional authors not shown)

Abstract: Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provid… ▽ More Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica. △ Less

Submitted 5 November, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: Accepted for publication in PASA

Journal ref: PASA (2020) 37, e047

A population analysis of pulse broadening in ASKAP Fast Radio Bursts

Authors: Hao Qiu, Ryan M. Shannon, Wael Farah, Jean-Pierre Macquart, Adam T. Deller, Keith W. Bannister, Clancy W. James, Chris Flynn, Cherie K. Day, Shivani Bhandari, Tara Murphy

Abstract: The pulse morphology of fast radio bursts (FRBs) provides key information in both understanding progenitor physics and the plasma medium through which the burst propagates. We present a study of the profiles of 33 bright FRBs detected by the Australian Square Kilometre Array Pathfinder. We identify seven FRBs with measureable intrinsic pulse widths, including two FRBs that have been seen to repeat… ▽ More The pulse morphology of fast radio bursts (FRBs) provides key information in both understanding progenitor physics and the plasma medium through which the burst propagates. We present a study of the profiles of 33 bright FRBs detected by the Australian Square Kilometre Array Pathfinder. We identify seven FRBs with measureable intrinsic pulse widths, including two FRBs that have been seen to repeat. In our modest sample we see no evidence for bimodality in the pulse width distribution. We also identify five FRBs with evidence of millisecond timescale pulse broadening caused by scattering in inhomogeneous plasma. We find no evidence for a relationship between pulse broadening and extragalactic dispersion measure. The scattering could be either caused by extreme turbulence in the host galaxy or chance propagation through foreground galaxies. With future high time resolution observations and detailed study of host galaxy properties we may be able to probe line-of-sight turbulence on gigaparsec scales. △ Less

Submitted 7 July, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: 10 pages, 6 figures, Accepted to MNRAS

The SUrvey for Pulsars and Extragalactic Radio Bursts V: Recent Discoveries and Full Timing Solutions

Authors: R. Spiewak, C. Flynn, S. Johnston, E. F. Keane, M. Bailes, E. D. Barr, S. Bhandari, M. Burgay, F. Jankowski, M. Kramer, V. Morello, A. Possenti, V. Venkatraman Krishnan

Abstract: The SUrvey for Pulsars and Extragalactic Radio Bursts ran from 2014 April to 2019 August, covering a large fraction of the southern hemisphere at mid- to high-galactic latitudes, and consisting of 9-minute pointings taken with the 20-cm multibeam receiver on the Parkes Radio Telescope. Data up to 2017 September 21 have been searched using standard Fourier techniques, single-pulse searches, and Fas… ▽ More The SUrvey for Pulsars and Extragalactic Radio Bursts ran from 2014 April to 2019 August, covering a large fraction of the southern hemisphere at mid- to high-galactic latitudes, and consisting of 9-minute pointings taken with the 20-cm multibeam receiver on the Parkes Radio Telescope. Data up to 2017 September 21 have been searched using standard Fourier techniques, single-pulse searches, and Fast Folding Algorithm searches. We present 19 new discoveries, bringing the total to 27 discoveries in the programme, and we report the results of follow-up timing observations at Parkes for 26 of these pulsars, including the millisecond pulsar PSR J1421-4409; the faint, highly-modulated, slow pulsar PSR J1646-1910; and the nulling pulsar PSR J1337-4441. We present new timing solutions for 23 pulsars, and we report flux densities, modulation indices, and polarization properties. △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: Accepted to MNRAS; data available at https://doi.org/10.5281/zenodo.3900980

The fates of massive stars: exploring uncertainties in stellar evolution with METISSE

Authors: Poojan Agrawal, Jarrod Hurley, Simon Stevenson, Dorottya Szécsi, Chris Flynn

Abstract: In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new… ▽ More In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new approach to the commonly used Single-Star Evolution (SSE) fitting formulae, one that is more adaptable: Method of Interpolation for Single Star Evolution (METISSE). It makes use of interpolation between sets of pre-computed stellar tracks to approximate evolution parameters for a population of stars. We have used METISSE with detailed stellar tracks computed by the Modules for Experiments in Stellar Astrophysics (MESA), Bonn Evolutionary Code (BEC) and Cambridge STARS code. METISSE better reproduces stellar tracks computed using the STARS code compared to SSE, and is on average three times faster. Using stellar tracks computed with MESA and BEC, we apply METISSE to explore the differences in the remnant masses, the maximum radial expansion, and the main-sequence lifetime of massive stars. We find that different physical ingredients used in the evolution of stars, such as the treatment of radiation dominated envelopes, can impact their evolutionary outcome. For stars in the mass range 9 to 100 M$_\odot$, the predictions of remnant masses can vary by up to 20 M$_\odot$, while the maximum radial expansion achieved by a star can differ by an order of magnitude between different stellar models. △ Less

Submitted 24 July, 2020; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS

High time resolution and polarisation properties of ASKAP-localised fast radio bursts

Authors: Cherie K. Day, A. T. Deller, R. M. Shannon, Hao Qiu, Keith W. Bannister, S. Bhandari, Ron Ekers, Chris Flynn, C. W. James, J. -P. Macquart, Chris J. Phillips, Elizabeth K. Mahony, J. Xavier Prochaska

Abstract: Combining high time and frequency resolution full-polarisation spectra of Fast Radio Bursts (FRBs) with knowledge of their host galaxy properties provides an opportunity to study both the emission mechanism generating them and the impact of their propagation through their local environment, host galaxy, and the intergalactic medium. The Australian Square Kilometre Array Pathfinder (ASKAP) telescop… ▽ More Combining high time and frequency resolution full-polarisation spectra of Fast Radio Bursts (FRBs) with knowledge of their host galaxy properties provides an opportunity to study both the emission mechanism generating them and the impact of their propagation through their local environment, host galaxy, and the intergalactic medium. The Australian Square Kilometre Array Pathfinder (ASKAP) telescope has provided the first ensemble of bursts with this information. In this paper, we present the high time and spectral resolution, full polarisation observations of five localised FRBs to complement the results published for the previously studied ASKAP FRB~181112. We find that every FRB is highly polarised, with polarisation fractions ranging from 80 -- 100\%, and that they are generally dominated by linear polarisation. While some FRBs in our sample exhibit properties associated with an emerging archetype (i.e., repeating or apparently non-repeating), others exhibit characteristic features of both, implying the existence of a continuum of FRB properties. When examined at high time resolution, we find that all FRBs in our sample have evidence for multiple sub-components and for scattering at a level greater than expected from the Milky Way. We find no correlation between the diverse range of FRB properties (e.g., scattering time, intrinsic width, and rotation measure) and any global property of their host galaxy. The most heavily scattered bursts reside in the outskirts of their host galaxies, suggesting that the source-local environment rather than the host interstellar medium is likely the dominant origin of the scattering in our sample. △ Less

Submitted 23 August, 2020; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: 18 pages, 4 figures, 5 tables, accepted by MNRAS

Journal ref: MNRAS 497 3 (2020) 3335-3350

The host galaxies and progenitors of Fast Radio Bursts localized with the Australian Square Kilometre Array Pathfinder

Authors: Shivani Bhandari, Elaine M. Sadler, J. Xavier Prochaska, Sunil Simha, Stuart D. Ryder, Lachlan Marnoch, Keith W. Bannister, Jean-Pierre Macquart, Chris Flynn, Ryan M. Shannon, Nicolas Tejos, Felipe Corro-Guerra, Cherie K. Day, Adam T. Deller, Ron Ekers, Sebastian Lopez, Elizabeth K. Mahony, Consuelo Nuñez, Chris Phillips

Abstract: The Australian SKA Pathfinder (ASKAP) telescope has started to localize Fast Radio Bursts (FRBs) to arcsecond accuracy from the detection of a single pulse, allowing their host galaxies to be reliably identified. We discuss the global properties of the host galaxies of the first four FRBs localized by ASKAP, which lie in the redshift range $0.11<z<0.48$. All four are massive galaxies (log(… ▽ More The Australian SKA Pathfinder (ASKAP) telescope has started to localize Fast Radio Bursts (FRBs) to arcsecond accuracy from the detection of a single pulse, allowing their host galaxies to be reliably identified. We discuss the global properties of the host galaxies of the first four FRBs localized by ASKAP, which lie in the redshift range $0.11<z<0.48$. All four are massive galaxies (log( $M_{*}/ M_{\odot}$) $\sim 9.4 -10.4$) with modest star-formation rates of up to $2M_{\odot}$yr$^{-1}$ -- very different to the host galaxy of the first repeating FRB 121102, which is a dwarf galaxy with a high specific star-formation rate. The FRBs localized by ASKAP typically lie in the outskirts of their host galaxies, which appears to rule out FRB progenitor models that invoke active galactic nuclei (AGN) or free-floating cosmic strings. The stellar population seen in these host galaxies also disfavors models in which all FRBs arise from young magnetars produced by superluminous supernovae (SLSNe), as proposed for the progenitor of FRB 121102. A range of other progenitor models (including compact-object mergers and magnetars arising from normal core-collapse supernovae) remain plausible. △ Less

Submitted 27 May, 2020; originally announced May 2020.

Comments: 15 pages, 5 figures, 5 tables. Accepted for publication in ApJL

Constraints on the assembly history of the Milky Way's smooth, diffuse stellar halo from the metallicity-dependent, radially-dominated velocity anisotropy profiles probed with K giants and BHB stars using LAMOST, SDSS/SEGUE, and Gaia

Authors: Sarah A. Bird, Xiang-Xiang Xue, Chao Liu, Juntai Shen, Chris Flynn, Chengqun Yang, Gang Zhao, Hai-Jun Tian

Abstract: We analyze the anisotropy profile of the Milky Way's smooth, diffuse stellar halo using SDSS/SEGUE blue horizontal branch stars and SDSS/SEGUE and LAMOST K giants. These intrinsically luminous stars allow us to probe the halo to approximately 100 kpc from the Galactic center. Line-of-sight velocities, distances, metallicities, and proper motions are available for all stars via SDSS/SEGUE, LAMOST,… ▽ More We analyze the anisotropy profile of the Milky Way's smooth, diffuse stellar halo using SDSS/SEGUE blue horizontal branch stars and SDSS/SEGUE and LAMOST K giants. These intrinsically luminous stars allow us to probe the halo to approximately 100 kpc from the Galactic center. Line-of-sight velocities, distances, metallicities, and proper motions are available for all stars via SDSS/SEGUE, LAMOST, and Gaia, and we use these data to construct a full 7D set consisting of positions, space motions, and metallicity. We remove substructure from our samples using integrals of motion based on the method of Xue et al. We find radially dominated kinematic profiles with nearly constant anisotropy within 20 kpc, beyond which the anisotropy profile gently declines although remains radially dominated to the furthest extents of our sample. Independent of star type or substructure removal, the anisotropy depends on metallicity, such that the orbits of the stars become less radial with decreasing metallicity. For $-1.7<$ [Fe/H] $<-1$, the smooth, diffuse halo anisotropy profile begins to decline at Galactocentric distances $\sim20$ kpc, from $β\sim0.9$ to 0.7 for K giants and from $β\sim0.8$ to 0.1 for blue horizontal branch stars. For [Fe/H] $<-1.7$, the smooth, diffuse halo anisotropy remains constant along all distances with $0.2<β<0.7$ depending on the metallicity range probed, although independent on star type. These samples are ideal for estimating the total Galactic mass as they represent the virialized stellar halo system. △ Less

Submitted 30 April, 2021; v1 submitted 12 May, 2020; originally announced May 2020.

Comments: 35 pages, 20 figures, 2 tables, accepted by ApJ, comments welcome

The Rotation Curve, Mass Distribution and Dark Matter Content of the Milky Way from Classical Cepheids

Authors: Iminhaji Ablimit, Gang Zhao, Chris Flynn, Sarah A. Bird

Abstract: With the increasing numbers of large stellar survey projects, the quality and quantity of excellent tracers to study the Milky Way is rapidly growing, one of which is the classical Cepheids. Classical Cepheids are high precision standard candles with very low typical uncertainties ($<$ 3%) available via the mid-infrared period-luminosity relation. About 3500 classical Cepheids identified from OGLE… ▽ More With the increasing numbers of large stellar survey projects, the quality and quantity of excellent tracers to study the Milky Way is rapidly growing, one of which is the classical Cepheids. Classical Cepheids are high precision standard candles with very low typical uncertainties ($<$ 3%) available via the mid-infrared period-luminosity relation. About 3500 classical Cepheids identified from OGLE, ASAS-SN, Gaia, WISE and ZTF survey data have been analyzed in this work, and their spatial distributions show a clear signature of Galactic warp. Two kinematical methods are adopted to measure the Galactic rotation curve in the Galactocentric distance range of $4\lesssim R_{\rm GC} \lesssim 19$ kpc. Gently declining rotation curves are derived by both the proper motion (PM) method and 3-dimensional velocity vector (3DV) method. The largest sample of classical Cepheids with most accurate 6D phase-space coordinates available to date are modeled in the 3DV method, and the resulting rotation curve is found to decline at the relatively smaller gradient of ($-1.33\pm0.1$) ${\rm km\,s^{-1}\,kpc^{-1}}$. Comparing to results from the PM method, a higher rotation velocity (($232.5\pm0.83$) ${\rm km\,s^{-1}}$) is derived at the position of Sun in the 3DV method. The virial mass and local dark matter density are estimated from the 3DV method which is the more reliable method, $M_{\rm vir} = (0.822\pm0.052)\times 10^{12}\,M_\odot$ and $ρ_{\rm DM,\odot} = 0.33\pm0.03$ GeV ${\rm cm^{-3}}$, respectively. △ Less

Submitted 30 April, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: 21 pages, 5 figures, 1 table. Accepted for publication in Astrophysical Journal Letters. Following Minor Revision

Journal ref: Astrophysical Journal Letters, 895:L12 (8pp), 2020 May 20. Open Access

The UTMOST pulsar timing programme II: Timing noise across the pulsar population

Authors: Marcus E. Lower, Matthew Bailes, Ryan M. Shannon, Simon Johnston, Chris Flynn, Stefan Osłowski, Vivek Gupta, Wael Farah, Timothy Bateman, Anne J. Green, Richard Hunstead, Andrew Jameson, Fabian Jankowski, Aditya Parthasarathy, Daniel C. Price, Angus Sutherland, David Temby, Vivek Venkatraman Krishnan

Abstract: While pulsars possess exceptional rotational stability, large scale timing studies have revealed at least two distinct types of irregularities in their rotation: red timing noise and glitches. Using modern Bayesian techniques, we investigated the timing noise properties of 300 bright southern-sky radio pulsars that have been observed over 1.0-4.8 years by the upgraded Molonglo Observatory Synthesi… ▽ More While pulsars possess exceptional rotational stability, large scale timing studies have revealed at least two distinct types of irregularities in their rotation: red timing noise and glitches. Using modern Bayesian techniques, we investigated the timing noise properties of 300 bright southern-sky radio pulsars that have been observed over 1.0-4.8 years by the upgraded Molonglo Observatory Synthesis Telescope (MOST). We reanalysed the spin and spin-down changes associated with nine previously reported pulsar glitches, report the discovery of three new glitches and four unusual glitch-like events in the rotational evolution of PSR J1825$-$0935. We develop a refined Bayesian framework for determining how red noise strength scales with pulsar spin frequency ($ν$) and spin-down frequency ($\dotν$), which we apply to a sample of 280 non-recycled pulsars. With this new method and a simple power-law scaling relation, we show that red noise strength scales across the non-recycled pulsar population as $ν^{a} |\dotν|^{b}$, where $a = -0.84^{+0.47}_{-0.49}$ and $b = 0.97^{+0.16}_{-0.19}$. This method can be easily adapted to utilise more complex, astrophysically motivated red noise models. Lastly, we highlight our timing of the double neutron star PSR J0737$-$3039, and the rediscovery of a bright radio pulsar originally found during the first Molonglo pulsar surveys with an incorrectly catalogued position. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: Accepted by MNRAS. 28 pages, 8 figures, 8 tables

Lense-Thirring frame dragging induced by a fast-rotating white dwarf in a binary pulsar system

Authors: V. Venkatraman Krishnan, M. Bailes, W. van Straten, N. Wex, P. C. C. Freire, E. F. Keane, T. M. Tauris, P. A. Rosado, N. D. R. Bhat, C. Flynn, A. Jameson, S. Osłowski

Abstract: Radio pulsars in short-period eccentric binary orbits can be used to study both gravitational dynamics and binary evolution. The binary system containing PSR J1141$-$6545 includes a massive white dwarf (WD) companion that formed before the gravitationally bound young radio pulsar. We observe a temporal evolution of the orbital inclination of this pulsar that we infer is caused by a combination of… ▽ More Radio pulsars in short-period eccentric binary orbits can be used to study both gravitational dynamics and binary evolution. The binary system containing PSR J1141$-$6545 includes a massive white dwarf (WD) companion that formed before the gravitationally bound young radio pulsar. We observe a temporal evolution of the orbital inclination of this pulsar that we infer is caused by a combination of a Newtonian quadrupole moment and Lense-Thirring precession of the orbit resulting from rapid rotation of the WD. Lense-Thirring precession, an effect of relativistic frame-dragging, is a prediction of general relativity. This detection is consistent with the evolutionary scenario in which the WD accreted matter from the pulsar progenitor, spinning up the WD to a period $< 200$ seconds. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: Accepted for publication in Science

The High Time Resolution Universe Pulsar Survey -- XVI. Discovery and timing of 40 pulsars from the southern Galactic plane

Authors: A. D. Cameron, D. J. Champion, M. Bailes, V. Balakrishnan, E. D. Barr, C. G. Bassa, S. Bates, S. Bhandari, N. D. R. Bhat, M. Burgay, S. Burke-Spolaor, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, L. Levin, A. G. Lyne, C. Ng, E. Petroff, A. Possenti, D. A. Smith, B. W. Stappers, W. van Straten, C. Tiburzi , et al. (1 additional authors not shown)

Abstract: We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration leng… ▽ More We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration lengths will allow us to discover pulsars at the lower end of the pulsar luminosity distribution. We report the discovery of 40 pulsars, including three millisecond pulsar-white dwarf binary systems (PSRs J1537-5312, J1547-5709 and J1618-4624), a black-widow binary system (PSR J1745-23) and a candidate black-widow binary system (PSR J1727-2951), a glitching pulsar (PSR J1706-4434), an eclipsing binary pulsar with a 1.5-yr orbital period (PSR J1653-45), and a pair of long spin-period binary pulsars which display either nulling or intermittent behaviour (PSRs J1812-15 and J1831-04). We show that the total population of 100 pulsars discovered in the HTRU-S LowLat survey to date represents both an older and lower-luminosity population, and indicates that we have yet to reach the bottom of the luminosity distribution function. We present evaluations of the performance of our search technique and of the overall yield of the survey, considering the 94% of the survey which we have processed to date. We show that our pulsar yield falls below earlier predictions by approximately 25% (especially in the case of millisecond pulsars), and discuss explanations for this discrepancy as well as future adaptations in RFI mitigation and searching techniques which may address these shortfalls. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Comments: 28 pages, 9 figures, 13 tables

Detection of a Glitch in PSR J0908$-$4913 by UTMOST

Authors: Marcus E. Lower, Matthew Bailes, Ryan M. Shannon, Simon Johnston, Chris Flynn, Timothy Bateman, Duncan Campbell-Wilson, Cherie K. Day, Adam Deller, Wael Farah, Anne J. Green, Vivek Gupta, Richard W. Hunstead, Andrew Jameson, Ayushi Mandlik, Stefan Osłowski, Aditya Parthasarathy, Daniel C. Price, Angus Sutherland, David Temby, Glen Torr, Glenn Urquhart, Vivek Venkatraman Krishnan

Abstract: We report the first detection of a glitch in the radio pulsar PSR J0908$-$4913 (PSR B0906$-$49) during regular timing observations by the Molonglo Observatory Synthesis Telescope (MOST) as part of the UTMOST project. We report the first detection of a glitch in the radio pulsar PSR J0908$-$4913 (PSR B0906$-$49) during regular timing observations by the Molonglo Observatory Synthesis Telescope (MOST) as part of the UTMOST project. △ Less

Submitted 17 December, 2019; originally announced December 2019.

Comments: 2 pages, 1 table

Journal ref: Res. Notes AAS 3 192 (2019)

Which bright fast radio bursts repeat?

Authors: C. W. James, S. Oslowski, C. Flynn, P. Kumar, K. Bannister, S. Bhandari, W. Farah, M. Kerr, D. R. Lorimer, J. -P. Macquart, C. Ng, C. Phillips, D. C. Price, H. Qiu, R. M. Shannon, R. Spiewak

Abstract: A handful of fast radio bursts (FRBs) are now known to repeat. However, the question remains --- do they all? We report on an extensive observational campaign with the Australian Square Kilometre Array Pathfinder (ASKAP), Parkes, and Robert C. Byrd Green Bank Telescope, searching for repeat bursts from FRBs detected by the Commensal Real-time ASKAP Fast Transients survey. In 383.2 hr of follow-up… ▽ More A handful of fast radio bursts (FRBs) are now known to repeat. However, the question remains --- do they all? We report on an extensive observational campaign with the Australian Square Kilometre Array Pathfinder (ASKAP), Parkes, and Robert C. Byrd Green Bank Telescope, searching for repeat bursts from FRBs detected by the Commensal Real-time ASKAP Fast Transients survey. In 383.2 hr of follow-up observations covering 27 FRBs initially detected as single bursts, only two repeat bursts from a single FRB, FRB 171019, were detected, which have been previously reported by Kumar et al. We use simulations of repeating FRBs that allow for clustering in burst arrival times to calculate new estimates for the repetition rate of FRB 171019, finding only slight evidence for incompatibility with the properties of FRB 121102. Our lack of repeat bursts from the remaining FRBs set limits on the model of all bursts being attributable to repeating FRBs. Assuming a reasonable range of repetition behaviour, at most 60% (90% C.L.) of these FRBs having an intrinsic burst distribution similar to FRB~121102. This result is shown to be robust against different assumptions on the nature of repeating FRB behaviour, and indicates that if indeed all FRBs repeat, the majority must do so very rarely. △ Less

Submitted 23 March, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: 12 pages, 7 figures, 5 tables. Updated after error discovered in original (v1) submission