-
Vacuum polarization corrections to hyperfine structure in many-electron atoms
Authors:
J. C. Hasted,
C. J. Fairhall,
O. R. Smits,
B. M. Roberts,
J. S. M. Ginges
Abstract:
We perform a theoretical study of vacuum polarization corrections to the hyperfine structure in many-electron atoms. Calculations are performed for systems of interest for precision atomic tests of fundamental physics belonging to the alkali-metal atoms and singly-ionized alkaline earths. The vacuum polarization is considered in the Uehling approximation, and we study the many-body effects core re…
▽ More
We perform a theoretical study of vacuum polarization corrections to the hyperfine structure in many-electron atoms. Calculations are performed for systems of interest for precision atomic tests of fundamental physics belonging to the alkali-metal atoms and singly-ionized alkaline earths. The vacuum polarization is considered in the Uehling approximation, and we study the many-body effects core relaxation, core polarization, and valence-core correlations in the relativistic framework. We find that for s states, the relative vacuum polarization correction may be well-approximated by that for hydrogen-like ions, though for all other states account of many-body effects -- in particular, the polarization of the core -- is needed to obtain the correct sign and magnitude of the effect.
△ Less
Submitted 26 September, 2024;
originally announced September 2024.
-
Smallness of the nuclear polarization effect in the hyperfine structure of heavy muonic atoms as a stimulus for next-generation experiments
Authors:
J. Vandeleur,
G. Sanamyan,
O. R. Smits,
I. A. Valuev,
N. S. Oreshkina,
J. S. M. Ginges
Abstract:
There is renewed interest in studies of muonic atoms, which may provide detailed information on nuclear structure. A major limiting factor in the interpretation of measurements is the nuclear polarization contribution. We propose a method to determine this contribution to the hyperfine structure in muonic atoms from a combination of theory and experiment for hydrogenlike ions and muonic atoms. App…
▽ More
There is renewed interest in studies of muonic atoms, which may provide detailed information on nuclear structure. A major limiting factor in the interpretation of measurements is the nuclear polarization contribution. We propose a method to determine this contribution to the hyperfine structure in muonic atoms from a combination of theory and experiment for hydrogenlike ions and muonic atoms. Applying the method to $^{203,205}$Tl and $^{209}$Bi, for which there are H-like ion and muonic atom hyperfine experimental data, we find that the nuclear polarization contribution for these systems is small, and place a limit on its size of less than $10\%$ the total hyperfine splitting. We have also performed direct calculations of the nuclear polarization contribution using a semi-analytical model, which indicate that it may be as much as two orders of magnitude smaller. Therefore, we conclude that the nuclear polarization correction to the hyperfine structure of muonic atoms does not represent a limiting factor for next-generation experiments.
△ Less
Submitted 29 August, 2024;
originally announced August 2024.
-
On the First Eigenvalue of the $p$-Laplace Operator with Robin Boundary Conditions in the Complement of a Compact Set
Authors:
Lukas Bundrock,
Tiziana Giorgi,
Robert Smits
Abstract:
We consider the first eigenvalue $λ_1$ of the $p$-Laplace operator subject to Robin boundary conditions in the exterior of a compact set. We discuss the conditions for the existence of a variational $λ_1$, depending on the boundary parameter, the space dimension, and $p$. Our analysis involves the first $p$-harmonic Steklov eigenvalue in exterior domains. We establish properties of $λ_1$ for the e…
▽ More
We consider the first eigenvalue $λ_1$ of the $p$-Laplace operator subject to Robin boundary conditions in the exterior of a compact set. We discuss the conditions for the existence of a variational $λ_1$, depending on the boundary parameter, the space dimension, and $p$. Our analysis involves the first $p$-harmonic Steklov eigenvalue in exterior domains. We establish properties of $λ_1$ for the exterior of a ball, including general inequalities, the asymptotic behavior as the boundary parameter approaches zero, and a monotonicity result with respect to a special type of domain inclusion. In two dimensions, we generalized to $p\neq 2$ some known shape optimization results.
△ Less
Submitted 12 August, 2024;
originally announced August 2024.
-
Brownian Motion with a Singular Drift
Authors:
Dante DeBlassie,
Adina Oprisan,
Robert G. Smits
Abstract:
We study the effect of a power law drift on Brownian motion in the positive half-line, where the order of the drift at 0 and infinity is different.
We study the effect of a power law drift on Brownian motion in the positive half-line, where the order of the drift at 0 and infinity is different.
△ Less
Submitted 9 March, 2024;
originally announced March 2024.
-
Statistical and Algebraic Properties for the 4-Laplacian via Averaging
Authors:
Ishraq Al-Awamleh,
Robert Smits
Abstract:
We find discrete analogs to continuous mean value principles that are used in the numerical analysis of the normalized p-Laplacian for particular values of p, specifically when p is 4.
We find discrete analogs to continuous mean value principles that are used in the numerical analysis of the normalized p-Laplacian for particular values of p, specifically when p is 4.
△ Less
Submitted 28 February, 2023;
originally announced March 2023.
-
Pushing the Limits of the Periodic Table -- A Review on Atomic Relativistic Electronic Structure Theory and Calculations for the Superheavy Elements
Authors:
O. R. Smits,
P. Indelicato,
W. Nazarewicz,
M. Piibeleht,
P. Schwerdtfege
Abstract:
We review the progress in atomic structure theory with a focus on superheavy elements and the aim to predict their ground state configuration and element's placement in the periodic table. To understand the electronic structure and correlations in the regime of large atomic numbers, it is important to correctly solve the Dirac equation in strong Coulomb fields, and also to take into account quantu…
▽ More
We review the progress in atomic structure theory with a focus on superheavy elements and the aim to predict their ground state configuration and element's placement in the periodic table. To understand the electronic structure and correlations in the regime of large atomic numbers, it is important to correctly solve the Dirac equation in strong Coulomb fields, and also to take into account quantum electrodynamic effects. We specifically focus on the fundamental difficulties encountered when dealing with the many-particle Dirac equation. We further discuss the possibility for future many-electron atomic structure calculations going beyond the critical nuclear charge \(Z_{\rm crit}\approx 170\), where levels such as the \(1s\) shell dive into the negative energy continuum (\(E_{nκ}<-m_ec^2\)). The nature of the resulting Gamow states within a rigged Hilbert space formalism is highlighted.
△ Less
Submitted 6 January, 2023;
originally announced January 2023.
-
First release of Apertif imaging survey data
Authors:
Elizabeth A. K. Adams,
B. Adebahr,
W. J. G. de Blok,
H. Denes,
K. M. Hess,
J. M. van der Hulst,
A. Kutkin,
D. M. Lucero,
R. Morganti,
V. A. Moss,
T. A. Oosterloo,
E. Orru,
R. Schulz,
A. S. van Amesfoort,
A. Berger,
O. M. Boersma,
M. Bouwhuis,
R. van den Brink,
W. A. van Cappellen,
L. Connor,
A. H. W. M. Coolen,
S. Damstra,
G. N. J. van Diepen,
T. J. Dijkema,
N. Ebbendorf
, et al. (34 additional authors not shown)
Abstract:
(Abridged) Apertif is a phased-array feed system for WSRT, providing forty instantaneous beams over 300 MHz of bandwidth. A dedicated survey program started on 1 July 2019, with the last observations taken on 28 February 2022. We describe the release of data products from the first year of survey operations, through 30 June 2020. We focus on defining quality control metrics for the processed data…
▽ More
(Abridged) Apertif is a phased-array feed system for WSRT, providing forty instantaneous beams over 300 MHz of bandwidth. A dedicated survey program started on 1 July 2019, with the last observations taken on 28 February 2022. We describe the release of data products from the first year of survey operations, through 30 June 2020. We focus on defining quality control metrics for the processed data products. The Apertif imaging pipeline, Apercal, automatically produces non-primary beam corrected continuum images, polarization images and cubes, and uncleaned spectral line and dirty beam cubes for each beam of an Apertif imaging observation. For this release, processed data products are considered on a beam-by-beam basis within an observation. We validate the continuum images by using metrics that identify deviations from Gaussian noise in the residual images. If the continuum image passes validation, we release all processed data products for a given beam. We apply further validation to the polarization and line data products. We release all raw observational data from the first year of survey observations, for a total of 221 observations of 160 independent target fields, covering approximately one thousand square degrees of sky. Images and cubes are released on a per beam basis, and 3374 beams are released. The median noise in the continuum images is 41.4 uJy/bm, with a slightly lower median noise of 36.9 uJy/bm in the Stokes V polarization image. The median angular resolution is 11.6"/sin(Dec). The median noise for all line cubes, with a spectral resolution of 36.6 kHz, is 1.6 mJy/bm, corresponding to a 3-sigma HI column density sensitivity of 1.8 x 10^20 atoms cm^-2 over 20 km/s (for a median angular resolution of 24" x 15"). We also provide primary beam images for each individual Apertif compound beam. The data are made accessible using a Virtual Observatory interface.
△ Less
Submitted 22 November, 2022; v1 submitted 10 August, 2022;
originally announced August 2022.
-
Characterising the Apertif primary beam response
Authors:
H. Dénes,
K. M. Hess,
E. A. K. Adams,
A. Kutkin,
R. Morganti,
J. M. van der Hulst,
T. A. Oosterloo,
V. A. Moss,
B. Adebahr,
W. J. G. de Blok,
M. V. Ivashina,
A. H. W. M. Coolen,
S. Damstra,
B. Hut,
G. M. Loose,
D. M. Lucero,
Y. Maan,
Á. Mika,
M. J. Norden,
L. C. Oostrum,
D. J. Pisano,
R. Smits,
W. A. van Cappellen,
R. van den Brink,
D. van der Schuur
, et al. (5 additional authors not shown)
Abstract:
Context. Phased Array Feeds (PAFs) are multi element receivers in the focal plane of a telescope that make it possible to form simultaneously multiple beams on the sky by combining the complex gains of the individual antenna elements. Recently the Westerbork Synthesis Radio Telescope (WSRT) was upgraded with PAF receivers and carried out several observing programs including two imaging surveys and…
▽ More
Context. Phased Array Feeds (PAFs) are multi element receivers in the focal plane of a telescope that make it possible to form simultaneously multiple beams on the sky by combining the complex gains of the individual antenna elements. Recently the Westerbork Synthesis Radio Telescope (WSRT) was upgraded with PAF receivers and carried out several observing programs including two imaging surveys and a time domain survey. The Apertif imaging surveys use a configuration, where 40 partially overlapping compound beams (CBs) are simultaneously formed on the sky and arranged in an approximately rectangular shape. Aims. This manuscript aims to characterise the response of the 40 Apertif CBs to create frequency-resolved, I, XX and YY polarization empirical beam shapes. The measured CB maps can be used for image deconvolution, primary beam correction and mosaicing of Apertif imaging data. Methods. We use drift scan measurements to measure the response of each of the 40 CBs of Apertif. We derive beam maps for all individual beams in I, XX and YY polarisation in 10 or 18 frequency bins over the same bandwidth as the Apertif imaging surveys. We sample the main lobe of the beams and the side lobes up to a radius of 0.6 degrees from the beam centres. In addition, we derive beam maps for each individual WSRT dish as well. Results. We present the frequency and time dependence of the beam shapes and sizes. We compare the compound beam shapes derived with the drift scan method to beam shapes derived with an independent method using a Gaussian Process Regression comparison between the Apertif continuum images and the NRAO VLA Sky Survey (NVSS) catalogue. We find a good agreement between the beam shapes derived with the two independent methods.
△ Less
Submitted 2 August, 2022; v1 submitted 19 May, 2022;
originally announced May 2022.
-
The Apertif science verification campaign - Characteristics of polarised radio sources
Authors:
B. Adebahr,
A. Berger,
E. A. K. Adams,
K. M. Hess,
W. J. G. de Blok,
H. Dénes,
V. A. Moss,
R. Schulz,
J. M. van der Hulst,
L. Connor,
S. Damstra,
B. Hut,
M. V. Ivashina,
G. M. Loose,
Y. Maan,
A. Mika,
H. Mulder,
M. J. Norden,
L. C. Oostrum,
E. Orrú,
M. Ruiter,
R. Smits,
W. A. van Cappellen,
J. van Leeuwen,
N. J. Vermaas
, et al. (2 additional authors not shown)
Abstract:
We analyse five early science datasets from the APERture Tile in Focus (Apertif) phased array feed system to verify the polarisation capabilities of Apertif in view of future larger data releases. We aim to characterise the source population of the polarised sky in the L-Band using polarised source information in combination with IR and optical data. We use automatic routines to generate full fiel…
▽ More
We analyse five early science datasets from the APERture Tile in Focus (Apertif) phased array feed system to verify the polarisation capabilities of Apertif in view of future larger data releases. We aim to characterise the source population of the polarised sky in the L-Band using polarised source information in combination with IR and optical data. We use automatic routines to generate full field-of-view Q- and U-cubes and perform RM-Synthesis, source finding, and cross-matching with published radio, optical, and IR data to generate polarised source catalogues. SED-fitting routines were used to determine photometric redshifts, star-formation rates, and galaxy masses. IR colour information was used to classify sources as AGN or star-forming-dominated and early- or late-type. We surveyed an area of 56deg$^2$ and detected 1357 polarised source components in 1170 sources. The fraction of polarised sources is 10.57% with a median fractional polarisation of 4.70$\pm$0.14%. We confirmed the reliability of the Apertif measurements by comparing them with polarised cross-identified NVSS sources. Average RMs of the individual fields lie within the error of the best Milky Way foreground measurements. All of our polarised sources were found to be dominated by AGN activity in the radio regime with most of them being radio-loud (79%) and of the FRII class (87%). The host galaxies of our polarised source sample are dominated by intermediate disc and star-forming disc galaxies. The contribution of star formation to the radio emission is on the order of a few percent for $\approx$10% of the polarised sources while for $\approx$90% it is completely dominated by the AGN. We do not see any change in fractional polarisation for different star-formation rates of the AGN host galaxies.
△ Less
Submitted 31 March, 2022;
originally announced March 2022.
-
A fast radio burst with sub-millisecond quasi-periodic structure
Authors:
Inés Pastor-Marazuela,
Joeri van Leeuwen,
Anna Bilous,
Liam Connor,
Yogesh Maan,
Leon Oostrum,
Emily Petroff,
Samayra Straal,
Dany Vohl,
E. A. K. Adams,
B. Adebahr,
Jisk Attema,
Oliver M. Boersma,
R. van den Brink,
W. A. van Cappellen,
A. H. W. M. Coolen,
S. Damstra,
H. Dénes,
K. M. Hess,
J. M. van der Hulst,
B. Hut,
A. Kutkin,
G. Marcel Loose,
D. M. Lucero,
Á. Mika
, et al. (9 additional authors not shown)
Abstract:
Fast radio bursts (FRBs) are extragalactic radio transients of extraordinary luminosity. Studying the diverse temporal and spectral behaviour recently observed in a number of FRBs may help determine the nature of the entire class. For example, a fast spinning or highly magnetised neutron star might generate the rotation-powered acceleration required to explain the bright emission. Periodic, sub-se…
▽ More
Fast radio bursts (FRBs) are extragalactic radio transients of extraordinary luminosity. Studying the diverse temporal and spectral behaviour recently observed in a number of FRBs may help determine the nature of the entire class. For example, a fast spinning or highly magnetised neutron star might generate the rotation-powered acceleration required to explain the bright emission. Periodic, sub-second components, suggesting such rotation, were recently reported in one FRB, and potentially in two more. Here we report the discovery of FRB 20201020A with Apertif, an FRB showing five components regularly spaced by 0.415 ms. This sub-millisecond structure in FRB 20201020A carries important clues about the progenitor of this FRB specifically, and potentially about that of FRBs in general. We thus contrast its features to the predictions of the main FRB source models. We perform a timing analysis of the FRB 20201020A components to determine the significance of the periodicity. We compare these against the timing properties of the previously reported CHIME FRBs with sub-second quasi-periodic components, and against two Apertif bursts from repeating FRB 20180916B that show complex time-frequency structure. We find the periodicity of FRB 20201020A to be marginally significant at 2.5$σ$. Its repeating subcomponents cannot be explained as a pulsar rotation since the required spin rate of over 2 kHz exceeds the limits set by typical neutron star equations of state and observations. The fast periodicity is also in conflict with a compact object merger scenario. These quasi-periodic components could, however, be caused by equidistant emitting regions in the magnetosphere of a magnetar. The sub-millisecond spacing of the components in FRB 20201020A, the smallest observed so far in a one-off FRB, may rule out both neutron-star rotation and binary mergers as the direct source of quasi-periodic FRBs.
△ Less
Submitted 16 February, 2022;
originally announced February 2022.
-
Apercal -- The Apertif Calibration Pipeline
Authors:
B. Adebahr,
R. Schulz,
T. J. Dijkema,
V. A. Moss,
A. R. Offringa,
A. Kutkin,
J. M. van der Hulst,
B. S. Frank,
N. P. E. Vilchez,
J. Verstappen,
E. K. Adams,
W. J. G. de Blok,
H. Denes,
K. M. Hess,
D. Lucero,
R. Morganti,
T. Oosterloo,
D. -J. Pisano,
M. V. Ivashina,
W. A. van Cappellen,
L. D. Connor,
A. H. W. M. Coolen,
S. Damstra,
G. M. Loose,
Y. Maan
, et al. (11 additional authors not shown)
Abstract:
Apertif (APERture Tile In Focus) is one of the Square Kilometre Array (SKA) pathfinder facilities. The Apertif project is an upgrade to the 50-year-old Westerbork Synthesis Radio Telescope (WSRT) using phased-array feed technology. The new receivers create 40 individual beams on the sky, achieving an instantaneous sky coverage of 6.5 square degrees. The primary goal of the Apertif Imaging Survey i…
▽ More
Apertif (APERture Tile In Focus) is one of the Square Kilometre Array (SKA) pathfinder facilities. The Apertif project is an upgrade to the 50-year-old Westerbork Synthesis Radio Telescope (WSRT) using phased-array feed technology. The new receivers create 40 individual beams on the sky, achieving an instantaneous sky coverage of 6.5 square degrees. The primary goal of the Apertif Imaging Survey is to perform a wide survey of 3500 square degrees (AWES) and a medium deep survey of 350 square degrees (AMES) of neutral atomic hydrogen (up to a redshift of 0.26), radio continuum emission and polarisation. Each survey pointing yields 4.6 TB of correlated data. The goal of Apercal is to process this data and fully automatically generate science ready data products for the astronomical community while keeping up with the survey observations. We make use of common astronomical software packages in combination with Python based routines and parallelisation. We use an object oriented module-based approach to ensure easy adaptation of the pipeline. A Jupyter notebook based framework allows user interaction and execution of individual modules as well as a full automatic processing of a complete survey observation. If nothing interrupts processing, we are able to reduce a single pointing survey observation on our five node cluster with 24 physical cores and 256 GB of memory each within 24h keeping up with the speed of the surveys. The quality of the generated images is sufficient for scientific usage for 44 % of the recorded data products with single images reaching dynamic ranges of several thousands. Future improvements will increase this percentage to over 80 %. Our design allowed development of the pipeline in parallel to the commissioning of the Apertif system.
△ Less
Submitted 7 December, 2021;
originally announced December 2021.
-
Removal and replacement of interference in tied-array radio pulsar observations using the spectral kurtosis estimator
Authors:
M. Purver,
C. G. Bassa,
I. Cognard,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. J. Lee,
K. Liu,
J. W. McKee,
D. Perrodin,
S. Sanidas,
R. Smits,
B. W. Stappers
Abstract:
We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is re…
▽ More
We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is rejected at a specified confidence level. We derive signal-to-noise ratios (SNRs) as a function of amount of zapping for folded radio pulsar observations consisting of a sum of signals from multiple telescopes in independent radio-frequency interference (RFI) environments, comparing four methods to compensate for lost data with coherent (tied-array) and incoherent summation. For coherently summed amplitudes, scaling amplitudes from non-zapped telescopes achieves a higher SNR than replacing zapped amplitudes with artificial noise. For incoherently summed power values, the highest SNR is given by scaling power from non-zapped telescopes to maintain a constant mean. We use spectral kurtosis to clean a tied-array radio pulsar observation by the Large European Array for Pulsars (LEAP): the signal from one telescope is zapped with time and frequency resolutions of 6.25 $μ$s and 0.16 MHz, removing interference along with 0.27 per cent of `good' data, giving an uncertainty of 0.25 $μ$s in pulse time of arrival (TOA) for PSR J1022+1001. We use a single-telescope observation to demonstrate recovery of the pulse profile shape, with 0.6 per cent of data zapped and a reduction from 1.22 to 0.70 $μ$s in TOA uncertainty.
△ Less
Submitted 24 February, 2023; v1 submitted 28 November, 2021;
originally announced November 2021.
-
Apertif, Phased Array Feeds for the Westerbork Synthesis Radio Telescope
Authors:
W. A. van Cappellen,
T. A. Oosterloo,
M. A. W. Verheijen,
E. A. K. Adams,
B. Adebahr,
R. Braun,
K. M. Hess,
H. Holties,
J. M. van der Hulst,
B. Hut,
E. Kooistra,
J. van Leeuwen,
G. M. Loose,
R. Morganti,
V. A. Moss,
E. Orrú,
M. Ruiter,
A. P. Schoenmakers,
N. J. Vermaas,
S. J. Wijnholds,
A. S. van Amesfoort,
M. J. Arts,
J. J. Attema,
L. Bakker,
C. G. Bassa
, et al. (65 additional authors not shown)
Abstract:
We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the s…
▽ More
We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the survey speed of the telescope. With this upgraded instrument, an imaging survey covering an area of 2300 deg2 is being performed which will deliver both continuum and spectral line data sets, of which the first data has been publicly released. In addition, a time domain transient and pulsar survey covering 15,000 deg2 is in progress. An overview of the Apertif science drivers, hardware and software of the upgraded telescope is presented, along with its key performance characteristics.
△ Less
Submitted 30 September, 2021; v1 submitted 29 September, 2021;
originally announced September 2021.
-
Dual-frequency single-pulse study of PSR B0950+08
Authors:
A. V. Bilous,
J. M. Griessmeier,
T. Pennucci,
Z. Wu,
L. Bondonneau,
V. Kondratiev,
J. van Leeuwen,
Y. Maan,
L. Connor,
L. C. Oostrum,
E. Petroff,
J. P. W. Verbiest,
D. Vohl,
J. W. McKee,
G. Shaifullah,
G. Theureau,
O. M. Ulyanov,
B. Cecconi,
A. H. Coolen,
S. Corbel,
S. Damstra,
H. Denes,
J. N. Girard,
B. Hut,
M. Ivashina
, et al. (11 additional authors not shown)
Abstract:
PSR B0950+08 is a bright non-recycled pulsar whose single-pulse fluence variability is reportedly large. Based on observations at two widely separated frequencies, 55 MHz (NenuFAR) and 1.4 GHz (Westerbork Synthesis Radio Telescope), we review the properties of these single pulses. We conclude that they are more similar to ordinary pulses of radio emission than to a special kind of short and bright…
▽ More
PSR B0950+08 is a bright non-recycled pulsar whose single-pulse fluence variability is reportedly large. Based on observations at two widely separated frequencies, 55 MHz (NenuFAR) and 1.4 GHz (Westerbork Synthesis Radio Telescope), we review the properties of these single pulses. We conclude that they are more similar to ordinary pulses of radio emission than to a special kind of short and bright Giant Pulses, observed from only a handful of pulsars. We argue that temporal variation of properties of interstellar medium along the line of sight to this nearby pulsar, namely the fluctuating size of decorrelation bandwidth of diffractive scintillation makes important contribution to observed single-pulse fluence variability. We further present interesting structures in the low-frequency single-pulse spectra that resemble the "sad trombones" seen in Fast Radio Bursts (FRBs); although for PSR B0950+08 the upward frequency drift is also routinely present. We explain these spectral features with radius-to-frequency mapping, similar to the model developed by Wang et al. (2019) for FRBs. Finally, we speculate that microsecond-scale fluence variability of the general pulsar population remains poorly known, and that its further study may bring important clues about the nature of FRBs.
△ Less
Submitted 26 November, 2021; v1 submitted 17 September, 2021;
originally announced September 2021.
-
A search for radio emission from double-neutron star merger GW190425 using Apertif
Authors:
Olivér Boersma,
Joeri van Leeuwen,
Elizabeth A. K. Adams,
Björn Adebahr,
Alexander Kutkin,
Tom Oosterloo,
W. J. G. de Blok,
R. van den Brink,
A. H. W. M. Coolen,
L. Connor,
S. Damstra,
H. Dénes,
K. M. Hess,
J. M. van der Hulst,
B. Hut,
M. Ivashina,
G. M. Loose,
D. M. Lucero,
Y. Maan,
Á. Mika,
V. A. Moss,
H. Mulder,
L. C. Oostrum,
M. Ruiter,
D. van der Schuur
, et al. (4 additional authors not shown)
Abstract:
Detection of the electromagnetic emission from coalescing binary neutron stars (BNS) is important for understanding the merger and afterglow. We present a search for a radio counterpart to the gravitational-wave source GW190425, a BNS merger, using Apertif on the Westerbork Synthesis Radio Telescope (WSRT). We observe a field of high probability in the associated localisation region for 3 epochs a…
▽ More
Detection of the electromagnetic emission from coalescing binary neutron stars (BNS) is important for understanding the merger and afterglow. We present a search for a radio counterpart to the gravitational-wave source GW190425, a BNS merger, using Apertif on the Westerbork Synthesis Radio Telescope (WSRT). We observe a field of high probability in the associated localisation region for 3 epochs at 68, 90 and 109 days post merger. We identify all sources that exhibit flux variations consistent with the expected afterglow emission of GW190425. We also look for possible transients. These are sources which are only present in one epoch. In addition, we quantify our ability to search for radio afterglows in fourth and future observing runs of the gravitational-wave detector network using Monte Carlo simulations. We found 25 afterglow candidates based on their variability. None of these could be associated with a possible host galaxy at the luminosity distance of GW190425. We also found 55 transient afterglow candidates that were only detected in one epoch. All turned out to be image artefacts. In the fourth observing run, we predict that up to three afterglows will be detectable by Apertif. While we did not find a source related to the afterglow emission of GW190425, the search validates our methods for future searches of radio afterglows.
△ Less
Submitted 9 April, 2021;
originally announced April 2021.
-
Chromatic periodic activity down to 120 MHz in a Fast Radio Burst
Authors:
Inés Pastor-Marazuela,
Liam Connor,
Joeri van Leeuwen,
Yogesh Maan,
Sander ter Veen,
Anna Bilous,
Leon Oostrum,
Emily Petroff,
Samayra Straal,
Dany Vohl,
Jisk Attema,
Oliver M. Boersma,
Eric Kooistra,
Daniel van der Schuur,
Alessio Sclocco,
Roy Smits,
Elizabeth A. K. Adams,
Björn Adebahr,
Willem J. G. de Blok,
Arthur H. W. M. Coolen,
Sieds Damstra,
Helga Dénes,
Kelley M. Hess,
Thijs van der Hulst,
Boudewijn Hut
, et al. (12 additional authors not shown)
Abstract:
Fast radio bursts (FRBs) are extragalactic astrophysical transients whose brightness requires emitters that are highly energetic, yet compact enough to produce the short, millisecond-duration bursts. FRBs have thus far been detected between 300 MHz and 8 GHz, but lower-frequency emission has remained elusive. A subset of FRBs is known to repeat, and one of those sources, FRB 20180916B, does so wit…
▽ More
Fast radio bursts (FRBs) are extragalactic astrophysical transients whose brightness requires emitters that are highly energetic, yet compact enough to produce the short, millisecond-duration bursts. FRBs have thus far been detected between 300 MHz and 8 GHz, but lower-frequency emission has remained elusive. A subset of FRBs is known to repeat, and one of those sources, FRB 20180916B, does so with a 16.3 day activity period. Using simultaneous Apertif and LOFAR data, we show that FRB 20180916B emits down to 120 MHz, and that its activity window is both narrower and earlier at higher frequencies. Binary wind interaction models predict a narrower periodic activity window at lower frequencies, which is the opposite of our observations. Our detections establish that low-frequency FRB emission can escape the local medium. For bursts of the same fluence, FRB 20180916B is more active below 200 MHz than at 1.4 GHz. Combining our results with previous upper-limits on the all-sky FRB rate at 150 MHz, we find that there are 3-450 FRBs/sky/day above 50 Jy ms at 90% confidence. We are able to rule out the scenario in which companion winds cause FRB periodicity. We also demonstrate that some FRBs live in clean environments that do not absorb or scatter low-frequency radiation.
△ Less
Submitted 15 December, 2020;
originally announced December 2020.
-
The Lennard Jones Potential Revisited -- Analytical Expressions for Vibrational Effects in Cubic and Hexagonal Close-Packed Lattices
Authors:
Peter Schwerdtfeger,
Antony Burrows,
Odile R. Smits
Abstract:
Analytical formulae are derived for the zero-point vibrational energy and anharmonicity corrections of the cohesive energy and the mode Grüneisen parameter within the Einstein model for the cubic lattices (sc, bcc and fcc) and for the hexagonal close-packed structure. This extends the work done by Lennard Jones and Ingham in 1924, Corner in 1939 and Wallace in 1965. The formulae are based on the d…
▽ More
Analytical formulae are derived for the zero-point vibrational energy and anharmonicity corrections of the cohesive energy and the mode Grüneisen parameter within the Einstein model for the cubic lattices (sc, bcc and fcc) and for the hexagonal close-packed structure. This extends the work done by Lennard Jones and Ingham in 1924, Corner in 1939 and Wallace in 1965. The formulae are based on the description of two-body energy contributions by an inverse power expansion (extended Lennard-Jones potential). These make use of three-dimensional lattice sums, which can be transformed to fast converging series and accurately determined by various expansion techniques. We apply these new lattice sum expressions to the rare gas solids and discuss associated critical points. The derived formulae give qualitative but nevertheless deep insight into vibrational effects in solids from the lightest (helium) to the heaviest rare gas element (oganesson), both presenting special cases because of strong quantum effects for the former and strong relativistic effects for the latter.
△ Less
Submitted 9 December, 2020;
originally announced December 2020.
-
Grey-Box Learning of Register Automata
Authors:
Bharat Garhewal,
Frits Vaandrager,
Falk Howar,
Timo Schrijvers,
Toon Lenaerts,
Rob Smits
Abstract:
Model learning (a.k.a. active automata learning) is a highly effective technique for obtaining black-box finite state models of software components. Thus far, generalisation to infinite state systems with inputs/outputs that carry data parameters has been challenging. Existing model learning tools for infinite state systems face scalability problems and can only be applied to restricted classes of…
▽ More
Model learning (a.k.a. active automata learning) is a highly effective technique for obtaining black-box finite state models of software components. Thus far, generalisation to infinite state systems with inputs/outputs that carry data parameters has been challenging. Existing model learning tools for infinite state systems face scalability problems and can only be applied to restricted classes of systems (register automata with equality/inequality). In this article, we show how we can boost the performance of model learning techniques by extracting the constraints on input and output parameters from a run, and making this grey-box information available to the learner. More specifically, we provide new implementations of the tree oracle and equivalence oracle from RALib, which use the derived constraints. We extract the constraints from runs of Python programs using an existing tainting library for Python, and compare our grey-box version of RALib with the existing black-box version on several benchmarks, including some data structures from Python's standard library. Our proof-of-principle implementation results in almost two orders of magnitude improvement in terms of numbers of inputs sent to the software system. Our approach, which can be generalised to richer model classes, also enables RALib to learn models that are out of reach of black-box techniques, such as combination locks.
△ Less
Submitted 21 September, 2020;
originally announced September 2020.
-
A bright, high rotation-measure FRB that skewers the M33 halo
Authors:
Liam Connor,
Joeri van Leeuwen,
L. C. Oostrum,
E. Petroff,
Yogesh Maan,
E. A. K. Adams,
J. J. Attema,
J. E. Bast,
O. M. Boersma,
H. Dénes,
D. W. Gardenier,
J. E. Hargreaves,
E. Kooistra,
I. Pastor-Marazuela,
R. Schulz,
A. Sclocco,
R. Smits,
S. M. Straal,
D. van der Schuur,
Dany Vohl,
B. Adebahr,
W. J. G. de Blok,
W. A. van Cappellen,
A. H. W. M. Coolen,
S. Damstra
, et al. (15 additional authors not shown)
Abstract:
We report the detection of a bright fast radio burst, FRB\,191108, with Apertif on the Westerbork Synthesis Radio Telescope (WSRT). The interferometer allows us to localise the FRB to a narrow $5\arcsec\times7\arcmin$ ellipse by employing both multibeam information within the Apertif phased-array feed (PAF) beam pattern, and across different tied-array beams. The resulting sight line passes close…
▽ More
We report the detection of a bright fast radio burst, FRB\,191108, with Apertif on the Westerbork Synthesis Radio Telescope (WSRT). The interferometer allows us to localise the FRB to a narrow $5\arcsec\times7\arcmin$ ellipse by employing both multibeam information within the Apertif phased-array feed (PAF) beam pattern, and across different tied-array beams. The resulting sight line passes close to Local Group galaxy M33, with an impact parameter of only 18\,kpc with respect to the core. It also traverses the much larger circumgalactic medium of M31, the Andromeda Galaxy. We find that the shared plasma of the Local Group galaxies could contribute $\sim$10\% of its dispersion measure of 588\,pc\,cm$^{-3}$. FRB\,191108 has a Faraday rotation measure of +474\,$\pm\,3$\,rad\,m$^{-2}$, which is too large to be explained by either the Milky Way or the intergalactic medium. Based on the more moderate RMs of other extragalactic sources that traverse the halo of M33, we conclude that the dense magnetised plasma resides in the host galaxy. The FRB exhibits frequency structure on two scales, one that is consistent with quenched Galactic scintillation and broader spectral structure with $Δν\approx40$\,MHz. If the latter is due to scattering in the shared M33/M31 CGM, our results constrain the Local Group plasma environment. We found no accompanying persistent radio sources in the Apertif imaging survey data.
△ Less
Submitted 22 September, 2020; v1 submitted 4 February, 2020;
originally announced February 2020.
-
Repeating fast radio bursts with WSRT/Apertif
Authors:
L. C. Oostrum,
Y. Maan,
J. van Leeuwen,
L. Connor,
E. Petroff,
J. J. Attema,
J. E. Bast,
D. W. Gardenier,
J. E. Hargreaves,
E. Kooistra,
D. van der Schuur,
A. Sclocco,
R. Smits,
S. M. Straal,
S. ter Veen,
D. Vohl,
E. A. K. Adams,
B. Adebahr,
W. J. G. de Blok,
R. H. van den Brink,
W. A. van Cappellen,
A. H. W. M. Coolen,
S. Damstra,
G. N. J. van Diepen,
B. S. Frank
, et al. (18 additional authors not shown)
Abstract:
Repeating fast radio bursts (FRBs) present excellent opportunities to identify FRB progenitors and host environments, as well as decipher the underlying emission mechanism. Detailed studies of repeating FRBs might also hold clues to the origin of FRBs as a population. We aim to detect the first two repeating FRBs: FRB 121102 (R1) and FRB 180814.J0422+73 (R2), and characterise their repeat statisti…
▽ More
Repeating fast radio bursts (FRBs) present excellent opportunities to identify FRB progenitors and host environments, as well as decipher the underlying emission mechanism. Detailed studies of repeating FRBs might also hold clues to the origin of FRBs as a population. We aim to detect the first two repeating FRBs: FRB 121102 (R1) and FRB 180814.J0422+73 (R2), and characterise their repeat statistics. We also want to significantly improve the sky localisation of R2. We use the Westerbork Synthesis Radio Telescope to conduct extensive follow-up of these two repeating FRBs. The new phased-array feed system, Apertif, allows covering the entire sky position uncertainty of R2 with fine spatial resolution in one pointing. We characterise the energy distribution and the clustering of detected R1 bursts. We detected 30 bursts from R1. Our measurements indicate a dispersion measure of 563.5(2) pc cm$^{-3}$, suggesting a significant increase in DM over the past few years. We place an upper limit of 8% on the linear polarisation fraction of the brightest burst. We did not detect any bursts from R2. A single power-law might not fit the R1 burst energy distribution across the full energy range or widely separated detections. Our observations provide improved constraints on the clustering of R1 bursts. Our stringent upper limits on the linear polarisation fraction imply a significant depolarisation, either intrinsic to the emission mechanism or caused by the intervening medium, at 1400 MHz that is not observed at higher frequencies. The non-detection of any bursts from R2 implies either a highly clustered nature of the bursts, a steep spectral index, or a combination of both. Alternatively, R2 has turned off completely, either permanently or for an extended period of time.
△ Less
Submitted 28 January, 2020; v1 submitted 27 December, 2019;
originally announced December 2019.
-
Statistics of VHE gamma-Rays in Temporal Association with Radio Giant Pulses from the Crab Pulsar
Authors:
MAGIC Collaboration,
M. L. Ahnen,
S. Ansoldi,
L. A. Antonelli,
C. Arcaro,
A. Babić,
B. Banerjee,
P. Bangale,
U. Barres de Almeida,
J. A. Barrio,
J. Becerra González,
W. Bednarek,
E. Bernardini,
A. Berti,
W. Bhattacharyya,
B. Biasuzzi,
A. Biland,
O. Blanch,
G. Bonnoli,
R. Carosi,
A. Carosi,
A. Chatterjee,
S. M. Colak,
P. Colin,
E. Colombo
, et al. (124 additional authors not shown)
Abstract:
The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E>100 GeV) gamma-rays from the Crab pulsar. 16 hours of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV with the Major Atmospheric Gamma-r…
▽ More
The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E>100 GeV) gamma-rays from the Crab pulsar. 16 hours of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes were performed. We searched for a statistical correlation between the radio and VHE gamma-ray emission with search windows of different lengths and different time lags to the arrival times of a radio GP. A dedicated search for an enhancement in the number of VHE gamma-rays correlated with the occurrence of radio GPs was carried out separately for the P1 and P2 phase ranges respectively. 99444 radio GPs have been detected in the radio data sample. We find no significant correlation between the GPs and VHE photons in any of the search windows. Depending on phase cuts and the chosen search windows we find upper limits at 95% confidence level on an increase in VHE gamma-ray events correlated with radio GPs between 7% and 61% of the average Crab pulsar VHE flux for the P1 and P2 phase ranges respectively. This puts upper limits on the flux increase during a radio GP of 12% to 2900% (depending on search window duration and phase cuts) of the pulsed VHE flux. This is the most stringent upper limit on a correlation between gamma-ray emission and radio GPs reported so far.}
△ Less
Submitted 31 January, 2020; v1 submitted 1 November, 2019;
originally announced November 2019.
-
Improved Bounds for Hermite-Hadamard Inequalities in Higher Dimensions
Authors:
Thomas Beck,
Barbara Brandolini,
Krzysztof Burdzy,
Antoine Henrot,
Jeffrey J. Langford,
Simon Larson,
Robert G. Smits,
Stefan Steinerberger
Abstract:
Let $Ω\subset \mathbb{R}^n$ be a convex domain and let $f:Ω\rightarrow \mathbb{R}$ be a positive, subharmonic function (i.e. $Δf \geq 0$). Then $$ \frac{1}{|Ω|} \int_Ω{f dx} \leq \frac{c_n}{ |\partial Ω| } \int_{\partial Ω}{ f dσ},$$ where $c_n \leq 2n^{3/2}$. This inequality was previously only known for convex functions with a much larger constant. We also show that the optimal constant satisfie…
▽ More
Let $Ω\subset \mathbb{R}^n$ be a convex domain and let $f:Ω\rightarrow \mathbb{R}$ be a positive, subharmonic function (i.e. $Δf \geq 0$). Then $$ \frac{1}{|Ω|} \int_Ω{f dx} \leq \frac{c_n}{ |\partial Ω| } \int_{\partial Ω}{ f dσ},$$ where $c_n \leq 2n^{3/2}$. This inequality was previously only known for convex functions with a much larger constant. We also show that the optimal constant satisfies $c_n \geq n-1$. As a byproduct, we establish a sharp geometric inequality for two convex domains where one contains the other $ Ω_2 \subset Ω_1 \subset \mathbb{R}^n$:
$$ \frac{|\partial Ω_1|}{|Ω_1|} \frac{| Ω_2|}{|\partial Ω_2|} \leq n.$$
△ Less
Submitted 13 July, 2019;
originally announced July 2019.
-
A detailed study of giant pulses from PSR B1937+21 using the Large European Array for Pulsars
Authors:
J. W. McKee,
B. W. Stappers,
C. G. Bassa,
S. Chen,
I. Cognard,
M. Gaikwad,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. J. Lee,
K. Liu,
D. Perrodin,
S. A. Sanidas,
R. Smits,
L. Wang,
W. W. Zhu
Abstract:
We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $α= -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken…
▽ More
We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $α= -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken power law, with the break occurring at $\sim7$ Jy $μ$s, with power law indices $α_{\text{low}} = -3.48 \pm 0.04$ and $α_{\text{high}} = -2.10 \pm 0.09$. The modulation indices of the GP emission are measured, which are found to vary by $\sim0.5$ at pulse phases close to the centre of the GP phase distributions. We find the frequency-resolved structure of GPs to vary significantly, and in a manner that cannot be attributed to the interstellar medium influence on the observed pulses. We examine the distribution of polarisation fractions of the GPs and find no correlation between GP emission phase and fractional polarisation. We use the GPs to time PSR B1937+21 and although the achievable time of arrival precision of the GPs is approximately a factor of two greater than that of the average pulse profile, there is a negligible difference in the precision of the overall timing solution when using the GPs.
△ Less
Submitted 7 November, 2018;
originally announced November 2018.
-
Tests of Gravitational Symmetries with Pulsar Binary J1713+0747
Authors:
W. W. Zhu,
G. Desvignes,
N. Wex,
R. N. Caballero,
D. J. Champion,
P. B. Demorest,
J. A. Ellis,
G. H. Janssen,
M. Kramer,
A. Krieger,
L. Lentati,
D. J. Nice,
S. M. Ransom,
I. H. Stairs,
B. W. Stappers,
J. P. W. Verbiest,
Z. Arzoumanian,
C. G. Bassa,
M. Burgay,
I. Cognard,
K. Crowter,
T. Dolch,
R. D. Ferdman,
E. Fonseca,
M. E. Gonzalez
, et al. (28 additional authors not shown)
Abstract:
Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pu…
▽ More
Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hatα_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|Δ|< 0.002$ for the UFF test and $-3\times10^{-20} < \hatα_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hatα_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.
△ Less
Submitted 26 February, 2018;
originally announced February 2018.
-
Pulsar science with data from the Large European Array for Pulsars
Authors:
J. W. McKee,
C. G. Bassa,
S. Chen,
M. Gaikwad,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. J. Lee,
K. Liu,
D. Perrodin,
S. A. Sanidas,
R. Smits,
B. W. Stappers,
L. Wang,
W. W. Zhu
Abstract:
The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nançay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz…
▽ More
The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nançay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz and a bandwidth of 128 MHz at each telescope, and are correlated offline on a cluster at Jodrell Bank Observatory using a purpose-built correlator, detailed in Smits et al. (2017). LEAP offers a substantial increase in sensitivity over that of the individual telescopes, and can operate in timing and imaging modes (notably in observations of the galactic centre radio magnetar; Wucknitz 2015). To date, 4 years of observations have been reduced. Here, we report on the scientific projects which have made use of LEAP data.
△ Less
Submitted 6 December, 2017;
originally announced December 2017.
-
The beamformer and correlator for the Large European Array for Pulsars
Authors:
R. Smits,
C. G. Bassa,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. J. Lee,
K. Liu,
J. McKee,
D. Perrodin,
M. Purver,
S. Sanidas,
B. W. Stappers,
W. W. Zhu
Abstract:
The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tie…
▽ More
The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tied-array beam from the raw voltages from each of telescopes. We explain the concepts and techniques involved in the process of adding the raw voltages coherently. We further present the software processing pipeline that is specifically designed to deal with data from widely spaced, inhomogeneous radio telescopes and describe the steps involved in preparing, correlating and creating the tied-array beam. This includes polarization calibration, bandpass correction, frequency dependent phase correction, interference mitigation and pulsar gating. A link is provided where the software can be obtained.
△ Less
Submitted 19 March, 2017;
originally announced March 2017.
-
Variability, polarimetry, and timing properties of single pulses from PSR J1713+0747 using the Large European Array for Pulsars
Authors:
K. Liu,
C. G. Bassa,
G. H. Janssen,
R. Karuppusamy,
J. McKee,
M. Kramer,
K. J. Lee,
D. Perrodin,
M. Purver,
S. Sanidas,
R. Smits,
B. W. Stappers,
P. Weltevrede,
W. W. Zhu
Abstract:
Single pulses preserve information about the pulsar radio emission and propagation in the pulsar magnetosphere, and understanding the behaviour of their variability is essential for estimating the fundamental limit on the achievable pulsar timing precision. Here we report the findings of our analysis of single pulses from PSR J1713+0747 with data collected by the Large European Array for Pulsars (…
▽ More
Single pulses preserve information about the pulsar radio emission and propagation in the pulsar magnetosphere, and understanding the behaviour of their variability is essential for estimating the fundamental limit on the achievable pulsar timing precision. Here we report the findings of our analysis of single pulses from PSR J1713+0747 with data collected by the Large European Array for Pulsars (LEAP). We present statistical studies of the pulse properties that include distributions of their energy, phase and width. Two modes of systematic sub-pulse drifting have been detected, with a periodicity of 7 and 3 pulse periods. The two modes appear at different ranges of pulse longitude but overlap under the main peak of the integrated profile. No evidence for pulse micro-structure is seen with a time resolution down to 140 ns. In addition, we show that the fractional polarisation of single pulses increases with their pulse peak flux density. By mapping the probability density of linear polarisation position angle with pulse longitude, we reveal the existence of two orthogonal polarisation modes. Finally, we find that the resulting phase jitter of integrated profiles caused by single pulse variability can be described by a Gaussian probability distribution only when at least 100 pulses are used for integration. Pulses of different flux densities and widths contribute approximately equally to the phase jitter, and no improvement on timing precision is achieved by using a sub-set of pulses with a specific range of flux density or width.
△ Less
Submitted 1 September, 2016;
originally announced September 2016.
-
Pulsar observations with European telescopes for testing gravity and detecting gravitational waves
Authors:
D. Perrodin,
C. G. Bassa,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. Lee,
K. Liu,
J. McKee,
M. Purver,
S. Sanidas,
R. Smits,
B. W. Stappers,
W. Zhu,
R. Concu,
A. Melis,
M. Burgay,
S. Casu,
A. Corongiu,
E. Egron,
N. Iacolina,
A. Pellizzoni,
M. Pilia,
A. Trois
Abstract:
A background of nanohertz gravitational waves from supermassive black hole binaries could soon be detected by pulsar timing arrays, which measure the times-of-arrival of radio pulses from millisecond pulsars with very high precision. The European Pulsar Timing Array uses five large European radio telescopes to monitor high-precision millisecond pulsars, imposing in this way strong constraints on a…
▽ More
A background of nanohertz gravitational waves from supermassive black hole binaries could soon be detected by pulsar timing arrays, which measure the times-of-arrival of radio pulses from millisecond pulsars with very high precision. The European Pulsar Timing Array uses five large European radio telescopes to monitor high-precision millisecond pulsars, imposing in this way strong constraints on a gravitational wave background. To achieve the necessary precision needed to detect gravitational waves, the Large European Array for Pulsars (LEAP) performs simultaneous observations of pulsars with all five telescopes, which allows us to coherently add the radio pulses, maximize the signal-to-noise of pulsar signals and increase the precision of times-of-arrival. We report on the progress made and results obtained by the LEAP collaboration, and in particular on the addition of the Sardinia Radio Telescope to the LEAP observations during its scientific validation phase. In addition, we discuss how LEAP can be used to monitor strong-gravity systems such as double neutron star systems and impose strong constraints on post-keplerian parameters.
△ Less
Submitted 5 August, 2016;
originally announced August 2016.
-
A millisecond pulsar in an extremely wide binary system
Authors:
C. G. Bassa,
G. H. Janssen,
B. W. Stappers,
T. M. Tauris,
T. Wevers,
P. G. Jonker,
L. Lentati,
J. P. W. Verbiest,
G. Desvignes,
E. Graikou,
L. Guillemot,
P. C. C. Freire,
P. Lazarus,
R. N. Caballero,
D. J. Champion,
I. Cognard,
A. Jessner,
C. Jordan,
R. Karuppusamy,
M. Kramer,
K. Lazaridis,
K. J. Lee,
K. Liu,
A. G. Lyne,
J. McKee
, et al. (8 additional authors not shown)
Abstract:
We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photome…
▽ More
We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.
△ Less
Submitted 1 April, 2016;
originally announced April 2016.
-
High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array
Authors:
G. Desvignes,
R. N. Caballero,
L. Lentati,
J. P. W. Verbiest,
D. J. Champion,
B. W. Stappers,
G. H. Janssen,
P. Lazarus,
S. Osłowski,
S. Babak,
C. G. Bassa,
P. Brem,
M. Burgay,
I. Cognard,
J. R. Gair,
E. Graikou,
L. Guillemot,
J. W. T. Hessels,
A. Jessner,
C. Jordan,
R. Karuppusamy,
M. Kramer,
A. Lassus,
K. Lazaridis,
K. J. Lee
, et al. (18 additional authors not shown)
Abstract:
We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing…
▽ More
We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.
△ Less
Submitted 26 February, 2016;
originally announced February 2016.
-
From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release
Authors:
L. Lentati,
R. M. Shannon,
W. A. Coles,
J. P. W. Verbiest,
R. van Haasteren,
J. A. Ellis,
R. N. Caballero,
R. N. Manchester,
Z. Arzoumanian,
S. Babak,
C. G. Bassa,
N. D. R. Bhat,
P. Brem,
M. Burgay,
S. Burke-Spolaor,
D. Champion,
S. Chatterjee,
I. Cognard,
J. M. Cordes,
S. Dai,
P. Demorest,
G. Desvignes,
T. Dolch,
R. D. Ferdman,
E. Fonseca
, et al. (58 additional authors not shown)
Abstract:
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single…
▽ More
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.
△ Less
Submitted 16 February, 2016;
originally announced February 2016.
-
The International Pulsar Timing Array: First Data Release
Authors:
J. P. W. Verbiest,
L. Lentati,
G. Hobbs,
R. van Haasteren,
P. B. Demorest,
G. H. Janssen,
J. -B. Wang,
G. Desvignes,
R. N. Caballero,
M. J. Keith,
D. J. Champion,
Z. Arzoumanian,
S. Babak,
C. G. Bassa,
N. D. R. Bhat,
A. Brazier,
P. Brem,
M. Burgay,
S. Burke-Spolaor,
S. J. Chamberlin,
S. Chatterjee,
B. Christy,
I. Cognard,
J. M. Cordes,
S. Dai
, et al. (67 additional authors not shown)
Abstract:
The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the…
▽ More
The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.
△ Less
Submitted 11 February, 2016;
originally announced February 2016.
-
A discrete in continuous mathematical model of cardiac progenitor cells formation and growth as spheroid clusters (Cardiospheres)
Authors:
Ezio Di Costanzo,
Alessandro Giacomello,
Elisa Messina,
Roberto Natalini,
Giuseppe Pontrelli,
Fabrizio Rossi,
Robert Smits,
Monika Twarogowska
Abstract:
We propose a discrete in continuous mathematical model describing the in vitro growth process of biophsy-derived mammalian cardiac progenitor cells growing as clusters in the form of spheres (Cardiospheres). The approach is hybrid: discrete at cellular scale and continuous at molecular level. In the present model cells are subject to the self-organizing collective dynamics mechanism and, additiona…
▽ More
We propose a discrete in continuous mathematical model describing the in vitro growth process of biophsy-derived mammalian cardiac progenitor cells growing as clusters in the form of spheres (Cardiospheres). The approach is hybrid: discrete at cellular scale and continuous at molecular level. In the present model cells are subject to the self-organizing collective dynamics mechanism and, additionally, they can proliferate and differentiate, also depending on stochastic processes. The two latter processes are triggered and regulated by chemical signals present in the environment. Numerical simulations show the structure and the development of the clustered progenitors and are in a good agreement with the results obtained from in vitro experiments.
△ Less
Submitted 22 December, 2015;
originally announced December 2015.
-
LEAP: the large European array for pulsars
Authors:
C. G. Bassa,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. J. Lee,
K. Liu,
J. McKee,
D. Perrodin,
M. Purver,
S. Sanidas,
R. Smits,
B. W. Stappers
Abstract:
The Large European Array for Pulsars (LEAP) is an experiment that harvests the collective power of Europe's largest radio telescopes in order to increase the sensitivity of high-precision pulsar timing. As part of the ongoing effort of the European Pulsar Timing Array (EPTA), LEAP aims to go beyond the sensitivity threshold needed to deliver the first direct detection of gravitational waves. The f…
▽ More
The Large European Array for Pulsars (LEAP) is an experiment that harvests the collective power of Europe's largest radio telescopes in order to increase the sensitivity of high-precision pulsar timing. As part of the ongoing effort of the European Pulsar Timing Array (EPTA), LEAP aims to go beyond the sensitivity threshold needed to deliver the first direct detection of gravitational waves. The five telescopes presently included in LEAP are: the Effelsberg telescope, the Lovell telescope at Jodrell Bank, the Nan\c cay radio telescope, the Sardinia Radio Telescope and the Westerbork Synthesis Radio Telescope. Dual polarization, Nyquist-sampled time-series of the incoming radio waves are recorded and processed offline to form the coherent sum, resulting in a tied-array telescope with an effective aperture equivalent to a 195-m diameter circular dish. All observations are performed using a bandwidth of 128 MHz centered at a frequency of 1396 MHz. In this paper, we present the design of the LEAP experiment, the instrumentation, the storage and transfer of data, and the processing hardware and software. In particular, we present the software pipeline that was designed to process the Nyquist-sampled time-series, measure the phase and time delays between each individual telescope and a reference telescope and apply these delays to form the tied-array coherent addition. The pipeline includes polarization calibration and interference mitigation. We also present the first results from LEAP and demonstrate the resulting increase in sensitivity, which leads to an improvement in the pulse arrival times.
△ Less
Submitted 20 November, 2015;
originally announced November 2015.
-
The noise properties of 42 millisecond pulsars from the European Pulsar Timing Array and their impact on gravitational wave searches
Authors:
R. N. Caballero,
K. J. Lee,
L. Lentati,
G. Desvignes,
D. J. Champion,
J. P. W. Verbiest,
G. H. Janssen,
B. W. Stappers,
M. Kramer,
P. Lazarus,
A. Possenti,
C. Tiburzi,
D. Perrodin,
S. Osłowski,
S. Babak,
C. G. Bassa,
P. Brem,
M. Burgay,
I. Cognard,
J. R. Gair,
E. Graikou,
L. Guillemot,
J. W. T. Hessels,
R. Karuppusamy,
A. Lassus
, et al. (14 additional authors not shown)
Abstract:
The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispe…
▽ More
The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.
△ Less
Submitted 19 January, 2016; v1 submitted 30 October, 2015;
originally announced October 2015.
-
Single-Source Gravitational Wave Limits from the J1713+0747 24-hr Global Campaign
Authors:
T. Dolch,
J. A. Ellis,
S. Chatterjee,
J. M. Cordes,
M. T. Lam,
C. Bassa,
B. Bhattacharyya,
D. J. Champion,
I. Cognard,
K. Crowter,
P. B. Demorest,
J. W. T. Hessels,
G. Janssen,
F. A. Jenet,
G. Jones,
C. Jordan,
R. Karuppusamy,
M. Keith,
V. I. Kondratiev,
M. Kramer,
P. Lazarus,
T. J. W. Lazio,
D. R. Lorimer,
D. R. Madison,
M. A. McLaughlin
, et al. (12 additional authors not shown)
Abstract:
Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all u…
▽ More
Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.
△ Less
Submitted 20 August, 2016; v1 submitted 17 September, 2015;
originally announced September 2015.
-
European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries
Authors:
Stanislav Babak,
Antoine Petiteau,
Alberto Sesana,
Patrick Brem,
Pablo A. Rosado,
Stephen R. Taylor,
Antoine Lassus,
Jason W. T. Hessels,
Cees G. Bassa,
Marta Burgay,
R. Nicolas Caballero,
David J. Champion,
Ismael Cognard,
Gregory Desvignes,
Jonathan R. Gair,
Lucas Guillemot,
Gemma H. Janssen,
Ramesh Karuppusamy,
Michael Kramer,
Patrick Lazarus,
K. J. Lee,
Lindley Lentati,
Kuo Liu,
Chiara M. F. Mingarelli,
Stefan Oslowski
, et al. (11 additional authors not shown)
Abstract:
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algori…
▽ More
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.
△ Less
Submitted 11 September, 2015; v1 submitted 7 September, 2015;
originally announced September 2015.
-
Limits on anisotropy in the nanohertz stochastic gravitational-wave background
Authors:
S. R. Taylor,
C. M. F. Mingarelli,
J. R. Gair,
A. Sesana,
G. Theureau,
S. Babak,
C. G. Bassa,
P. Brem,
M. Burgay,
R. N. Caballero,
D. J. Champion,
I. Cognard,
G. Desvignes,
L. Guillemot,
J. W. T. Hessels,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
A. Lassus,
P. Lazarus,
L. Lentati,
K. Liu,
S. Osłowski,
D. Perrodin,
A. Petiteau
, et al. (10 additional authors not shown)
Abstract:
The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in…
▽ More
The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.
△ Less
Submitted 29 June, 2015;
originally announced June 2015.
-
European Pulsar Timing Array Limits On An Isotropic Stochastic Gravitational-Wave Background
Authors:
Lindley Lentati,
Stephen R. Taylor,
Chiara M. F. Mingarelli,
Alberto Sesana,
Sotiris A. Sanidas,
Alberto Vecchio,
R. Nicolas Caballero,
K. J. Lee,
Rutger van Haasteren,
Stanislav Babak,
Cees G. Bassa,
Patrick Brem,
Marta Burgay,
David J. Champion,
Ismael Cognard,
Gregory Desvignes,
Jonathan R. Gair,
Lucas Guillemot,
Jason W. T. Hessels,
Gemma H. Janssen,
Ramesh Karuppusamy,
Michael Kramer,
Antoine Lassus,
Patrick Lazarus,
Kuo Liu
, et al. (11 additional authors not shown)
Abstract:
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtai…
▽ More
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $Ω_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $Gμ/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $Gμ/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $Ω^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.
△ Less
Submitted 9 September, 2015; v1 submitted 14 April, 2015;
originally announced April 2015.
-
Single pulse and profile variability study of PSR J1022+1001
Authors:
K. Liu,
R. Karuppusamy,
K. J. Lee,
B. W. Stappers,
M. Kramer,
R. Smits,
M. B. Purver,
G. H. Janssen,
D. Perrodin
Abstract:
Millisecond pulsars (MSPs) are known as highly stable celestial clocks. Nevertheless, recent studies have revealed the unstable nature of their integrated pulse profiles, which may limit the achievable pulsar timing precision. In this paper, we present a case study on the pulse profile variability of PSR J1022+1001. We have detected approximately 14,000 sub-pulses (components of single pulses) in…
▽ More
Millisecond pulsars (MSPs) are known as highly stable celestial clocks. Nevertheless, recent studies have revealed the unstable nature of their integrated pulse profiles, which may limit the achievable pulsar timing precision. In this paper, we present a case study on the pulse profile variability of PSR J1022+1001. We have detected approximately 14,000 sub-pulses (components of single pulses) in 35-hr long observations, mostly located at the trailing component of the integrated profile. Their flux densities and fractional polarisation suggest that they represent the bright end of the energy distribution in ordinary emission mode and are not giant pulses. The occurrence of sub-pulses from the leading and trailing components of the integrated profile is shown to be correlated. For sub-pulses from the latter, a preferred pulse width of approximately 0.25 ms has been found. Using simultaneous observations from the Effelsberg 100-m telescope and the Westerbork Synthesis Radio Telescope, we have found that the integrated profile varies on a timescale of a few tens of minutes. We show that improper polarisation calibration and diffractive scintillation cannot be the sole reason for the observed instability. In addition, we demonstrate that timing residuals generated from averages of the detected sub-pulses are dominated by phase jitter, and place an upper limit of ~700 ns for jitter noise based on continuous 1-min integrations.
△ Less
Submitted 24 February, 2015;
originally announced February 2015.
-
Magnetoresistance from broken spin helicity
Authors:
D. P. Leusink,
R. G. J. Smits,
P. Ngabonziza,
X. L. Wang,
S. Wiedmann,
U. Zeitler,
A. Brinkman
Abstract:
The propensity of some materials and multilayers to have a magnetic field dependent resistance, called magnetoresistance, has found commercial applications such as giant magnetoresistance harddisk read heads. But magnetoresistance can also be a powerful probe of electronic and magnetic interactions in matter. For example, magnetoresistance can be used to analyze multiband conductivity, conduction…
▽ More
The propensity of some materials and multilayers to have a magnetic field dependent resistance, called magnetoresistance, has found commercial applications such as giant magnetoresistance harddisk read heads. But magnetoresistance can also be a powerful probe of electronic and magnetic interactions in matter. For example, magnetoresistance can be used to analyze multiband conductivity, conduction inhomogeneity, localized magnetic moments, and (fractional) Landau level structure. For materials with strong spin-orbit interaction, magnetoresistance can be used as a probe for weak antilocalization or a nontrivial Berry phase, such as in topological insulator surface states. For the three dimensional topological insulators a large and linear magnetoresistance is often used as indication for underlying non-trivial topology, although the origin of this effect has not yet been established. Here, we observe a large magnetoresistance in the conducting bulk state of Bi$_2$Te$_3$. We show that this type of large magnetoresistance is due to the competition between helical spin-momentum locking (i.e. spin rotates with momentum direction) and the unidirectional spin alignment by an applied magnetic field. Warping effects are found to provide the (quasi) linear dependence on magnetic field. We provide a quantitative model for the helicity breaking induced magnetoresistance that can be applied to a vast range of materials, surfaces or interfaces with weak to strong spin-orbit interactions, such as the contemporary oxide interfaces, bulk Rashba systems, and topological insulator surface states.
△ Less
Submitted 12 December, 2014;
originally announced December 2014.
-
Frequency independent quenching of pulsed emission
Authors:
Vishal Gajjar,
Bhal Chandra Joshi,
Michael Kramer,
Ramesh Karuppusamy,
Roy Smits
Abstract:
Simultaneous observations at four different frequencies viz. 313, 607, 1380 and 4850 MHz, for three pulsars, PSRs B0031-07, B0809+74 and B2319+60, are reported in this paper. Identified null and burst pulses are highly concurrent across more than decade of frequency. Small fraction of non-concurrent pulses (less than or equal to 3%) are observed, most of which occur at the transition instances. We…
▽ More
Simultaneous observations at four different frequencies viz. 313, 607, 1380 and 4850 MHz, for three pulsars, PSRs B0031-07, B0809+74 and B2319+60, are reported in this paper. Identified null and burst pulses are highly concurrent across more than decade of frequency. Small fraction of non-concurrent pulses (less than or equal to 3%) are observed, most of which occur at the transition instances. We report, with very high significance for the first time, full broadband nature of the nulling phenomenon in these three pulsars. These results suggest that nulling invokes changes on the global magnetospheric scale.
△ Less
Submitted 22 September, 2014; v1 submitted 19 September, 2014;
originally announced September 2014.
-
A 24-Hour Global Campaign To Assess Precision Timing of the Millisecond Pulsar J1713+0747
Authors:
T. Dolch,
M. T. Lam,
J. M. Cordes,
S. Chatterjee,
C. Bassa,
B. Bhattacharyya,
D. J. Champion,
I. Cognard,
K. Crowter,
P. B. Demorest,
J. W. T. Hessels,
G. H. Janssen,
F. A. Jenet,
G. Jones,
C. Jordan,
R. Karuppusamy,
M. Keith,
V. I. Kondratiev,
M. Kramer,
P. Lazarus,
T. J. W. Lazio,
K. J. Lee,
M. A. McLaughlin,
J. Roy,
R. M. Shannon
, et al. (18 additional authors not shown)
Abstract:
The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) ti…
▽ More
The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.
△ Less
Submitted 16 December, 2014; v1 submitted 7 August, 2014;
originally announced August 2014.
-
Model-based asymptotically optimal dispersion measure correction for pulsar timing
Authors:
K. J. Lee,
C. G. Bassa,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. Liu,
D. Perrodin,
R. Smits,
B. W. Stappers,
R. van Haasteren,
L. Lentati
Abstract:
In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists o…
▽ More
In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists of two major steps. First, the spectral index and amplitude of dispersion measure variations are measured via a time-domain spectral analysis. Second, the linear optimal filter is constructed based on the model parameters found in the first step, and is used to extract the dispersion measure variation waveforms. Compared to current existing methods, this method has better time resolution for the study of short timescale dispersion variations, and generally produces smaller errors in waveform estimations. This method can process irregularly sampled data without any interpolation because of its time-domain nature. Furthermore, it offers the possibility to interpolate or extrapolate the waveform estimation to regions where no data is available. Examples using simulated data sets are included for demonstration.
△ Less
Submitted 8 April, 2014;
originally announced April 2014.
-
A semi-Lagrangian scheme for the game $p$-Laplacian via $p$-averaging
Authors:
M. Falcone,
S. Finzi Vita,
T. Giorgi,
R. G. Smits
Abstract:
We present and analyze an approximation scheme for the two-dimensional game $p$-Laplacian in the framework of viscosity solutions. The approximation is based on a semi-Lagrangian scheme which exploits the idea of $p$-averages. We study the properties of the scheme and prove that it converges, in particular cases, to the viscosity solution of the game $p$-Laplacian. We also present a numerical impl…
▽ More
We present and analyze an approximation scheme for the two-dimensional game $p$-Laplacian in the framework of viscosity solutions. The approximation is based on a semi-Lagrangian scheme which exploits the idea of $p$-averages. We study the properties of the scheme and prove that it converges, in particular cases, to the viscosity solution of the game $p$-Laplacian. We also present a numerical implementation of the scheme for different values of $p$; the numerical tests show that the scheme is accurate.
△ Less
Submitted 31 January, 2013;
originally announced January 2013.
-
Supra-oscillatory critical temperature dependence of Nb-Ho bilayers
Authors:
F. Chiodi,
J. D. S. Witt,
R. G. J. Smits,
L. Qu,
G. B. Halász,
C. -T. Wu,
O. T. Valls,
K. Halterman,
J. W. A. Robinson,
M. G. Blamire
Abstract:
We investigate the critical temperature Tc of a thin s-wave superconductor (Nb) proximity coupled to a helical rare earth ferromagnet (Ho). As a function of the Ho layer thickness, we observe multiple oscillations of Tc superimposed on a slow decay, that we attribute to the influence of the Ho on the Nb proximity effect. Because of Ho inhomogeneous magnetization, singlet and triplet pair correlati…
▽ More
We investigate the critical temperature Tc of a thin s-wave superconductor (Nb) proximity coupled to a helical rare earth ferromagnet (Ho). As a function of the Ho layer thickness, we observe multiple oscillations of Tc superimposed on a slow decay, that we attribute to the influence of the Ho on the Nb proximity effect. Because of Ho inhomogeneous magnetization, singlet and triplet pair correlations are present in the bilayers. We take both into consideration when solving the self consistent Bogoliubov-de Gennes equations, and we observe a reasonable agreement. We also observe non-trivial transitions into the superconducting state, the zero resistance state being attained after two successive transitions which appear to be associated with the magnetic structure of Ho.
△ Less
Submitted 6 November, 2012;
originally announced November 2012.
-
The optimal schedule for pulsar timing array observations
Authors:
K. J. Lee,
C. G. Bassa,
R. Karuppusamy,
M. Kramer,
R. Smits,
B. W. Stappers
Abstract:
In order to maximize the sensitivity of pulsar timing arrays to a stochastic gravitational wave background, we present computational techniques to optimize observing schedules. The techniques are applicable to both single and multi-telescope experiments. The observing schedule is optimized for each telescope by adjusting the observing time allocated to each pulsar while keeping the total amount of…
▽ More
In order to maximize the sensitivity of pulsar timing arrays to a stochastic gravitational wave background, we present computational techniques to optimize observing schedules. The techniques are applicable to both single and multi-telescope experiments. The observing schedule is optimized for each telescope by adjusting the observing time allocated to each pulsar while keeping the total amount of observing time constant. The optimized schedule depends on the timing noise characteristics of each individual pulsar as well as the performance of instrumentation. Several examples are given to illustrate the effects of different types of noise. A method to select the most suitable pulsars to be included in a pulsar timing array project is also presented.
△ Less
Submitted 19 April, 2012;
originally announced April 2012.
-
Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data
Authors:
R. van Haasteren,
Y. Levin,
G. H. Janssen,
K. Lazaridis,
M. Kramer B. W. Stappers,
G. Desvignes,
M. B. Purver,
A. G. Lyne,
R. D. Ferdman,
A. Jessner,
I. Cognard,
G. Theureau,
N. D'Amico,
A. Possenti,
M. Burgay,
A. Corongiu,
J. W. T. Hessels,
R. Smits,
J. P. W. Verbiest
Abstract:
Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^α$, where $f$ is the gravitational-wave fr…
▽ More
Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^α$, where $f$ is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude $A$ as a function of the unknown spectral slope $α$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $α=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.
△ Less
Submitted 2 March, 2011;
originally announced March 2011.
-
Gravitational wave astronomy of single sources with a pulsar timing array
Authors:
K. J. Lee,
N. Wex,
M. Kramer,
B. W. Stappers,
C. G. Bassa,
G. H. Janssen,
R. Karuppusamy,
R. Smits
Abstract:
Abbreviated:
We investigate the potential of detecting the gravitational wave from individual binary black hole systems using pulsar timing arrays (PTAs) and calculate the accuracy for determining the GW properties. This is done in a consistent analysis, which at the same time accounts for the measurement of the pulsar distances via the timing parallax.
We find that, at low redshift, a PTA is…
▽ More
Abbreviated:
We investigate the potential of detecting the gravitational wave from individual binary black hole systems using pulsar timing arrays (PTAs) and calculate the accuracy for determining the GW properties. This is done in a consistent analysis, which at the same time accounts for the measurement of the pulsar distances via the timing parallax.
We find that, at low redshift, a PTA is able to detect the nano-Hertz GW from super massive black hole binary systems with masses of $\sim10^8 - 10^{10}\,M_{\sun}$ less than $\sim10^5$\,years before the final merger, and those with less than $\sim10^3 - 10^4$ years before merger may allow us to detect the evolution of binaries.
We derive an analytical expression to describe the accuracy of a pulsar distance measurement via timing parallax. We consider five years of bi-weekly observations at a precision of 15\,ns for close-by ($\sim 0.5 - 1$\,kpc) pulsars. Timing twenty pulsars would allow us to detect a GW source with an amplitude larger than $5\times 10^{-17}$. We calculate the corresponding GW and binary orbital parameters and their measurement precision. The accuracy of measuring the binary orbital inclination angle, the sky position, and the GW frequency are calculated as functions of the GW amplitude. We note that the "pulsar term", which is commonly regarded as noise, is essential for obtaining an accurate measurement for the GW source location.
We also show that utilizing the information encoded in the GW signal passing the Earth also increases the accuracy of pulsar distance measurements. If the gravitational wave is strong enough, one can achieve sub-parsec distance measurements for nearby pulsars with distance less than $\sim 0.5 - 1$\,kpc.
△ Less
Submitted 1 March, 2011;
originally announced March 2011.
-
Evidence for gravitational quadrupole moment variations in the companion of PSR J2051-0827
Authors:
K. Lazaridis,
J. P. W. Verbiest,
T. M. Tauris,
B. W. Stappers,
M. Kramer,
N. Wex,
A. Jessner,
I. Cognard,
G. Desvignes,
G. H. Janssen,
M. B. Purver,
G. Theureau,
C. G. Bassa,
R. Smits
Abstract:
We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, μ_t = 7.3 {\pm} 0.4 mas/yr. Our tim…
▽ More
We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, μ_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.
△ Less
Submitted 28 February, 2011;
originally announced February 2011.