Abstract
Neutron Star Interior Composition Explorer
has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science ...target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. Two model parameters track different types of background events, while the third is used to predict a low-energy excess tied to observations conducted in sunlight. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4–12 keV; 50 detectors) of 0.87 c s
−1
, but in 5% (1%) of cases, the rate exceeds 10 (300) c s
−1
. Model residuals persist at 20%–30% of the initial rate for the brightest GTIs, implying one or more missing model parameters. Filtering criteria are given to flag GTIs likely to have unsatisfactory background predictions. With such filtering, we estimate a detection limit, 1.20 c s
−1
(3
σ
, single GTI) at 0.4–12 keV, equivalent to 3.6 × 10
−12
erg cm
−2
s
−1
for a Crab-like spectrum. The corresponding limit for soft X-ray sources is 0.51 c s
−1
at 0.3–2.0 keV, or 4.3 × 10
−13
erg cm
−2
s
−1
for a 100 eV blackbody. These limits would be four times lower if exploratory GTIs accumulate 10 ks of data after filtering at the level prescribed for faint sources. Such filtering selects background GTIs 85% of the time. An application of the model to a 1 s timescale makes it possible to distinguish source flares from possible surges in the background.
Abstract We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data ...by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot regions. We show that there is a significant offset to an antipodal spot configuration, mainly due to the non-half-cycle azimuthal separation of the two emitting spots. The offset angle from the antipode is inferred to be above 25° with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620.
Abstract
Over the past few decades, the measurement precision of some pulsar timing experiments has advanced from ∼10
μ
s to ∼10 ns, revealing many subtle phenomena. Such high precision demands both ...careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present
PINT
(
P
INT
I
s
N
ot
T
empo3
), a high-precision
Python
pulsar timing data analysis package, which is hosted on GitHub and available on the
Python
Package Index (PyPI) as
pint-pulsar
.
PINT
is well tested, validated, object oriented, and modular, enabling interactive data analysis and providing an extensible and flexible development platform for timing applications. It utilizes well-debugged public
Python
packages (e.g., the N
um
P
y
and A
stropy
libraries) and modern software development schemes (e.g., version control and efficient development with
git
and GitHub) and a continually expanding test suite for improved reliability, accuracy, and reproducibility.
PINT
is developed and implemented without referring to, copying, or transcribing the code from other traditional pulsar timing software packages (e.g.,
Tempo
/
Tempo2
) and therefore provides a robust tool for cross-checking timing analyses and simulating pulse arrival times. In this paper, we describe the design, use, and validation of
PINT
, and we compare timing results between it and
Tempo
and
Tempo2
.
In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the ...Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, allowing for pulsar timing based astrometry of this accreting millisecond pulsar.
Since the classical molecular dynamics simulator LAMMPS was released as an open source code in 2004, it has become a widely-used tool for particle-based modeling of materials at length scales ranging ...from atomic to mesoscale to continuum. Reasons for its popularity are that it provides a wide variety of particle interaction models for different materials, that it runs on any platform from a single CPU core to the largest supercomputers with accelerators, and that it gives users control over simulation details, either via the input script or by adding code for new interatomic potentials, constraints, diagnostics, or other features needed for their models. As a result, hundreds of people have contributed new capabilities to LAMMPS and it has grown from fifty thousand lines of code in 2004 to a million lines today. In this paper several of the fundamental algorithms used in LAMMPS are described along with the design strategies which have made it flexible for both users and developers. We also highlight some capabilities recently added to the code which were enabled by this flexibility, including dynamic load balancing, on-the-fly visualization, magnetic spin dynamics models, and quantum-accuracy machine learning interatomic potentials.
Program Title: Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS)
CPC Library link to program files:https://doi.org/10.17632/cxbxs9btsv.1
Developer's repository link:https://github.com/lammps/lammps
Licensing provisions: GPLv2
Programming language: C++, Python, C, Fortran
Supplementary material:https://www.lammps.org
Nature of problem: Many science applications in physics, chemistry, materials science, and related fields require parallel, scalable, and efficient generation of long, stable classical particle dynamics trajectories. Within this common problem definition, there lies a great diversity of use cases, distinguished by different particle interaction models, external constraints, as well as timescales and lengthscales ranging from atomic to mesoscale to macroscopic.
Solution method: The LAMMPS code uses parallel spatial decomposition, distributed neighbor lists, and parallel FFTs for long-range Coulombic interactions 1. The time integration algorithm is based on the Størmer-Verlet symplectic integrator 2, which provides better stability than higher-order non-symplectic methods. In addition, LAMMPS supports a wide range of interatomic potentials, constraints, diagnostics, software interfaces, and pre- and post-processing features.
Additional comments including restrictions and unusual features: This paper serves as the definitive reference for the LAMMPS code.
1S. Plimpton, Fast parallel algorithms for short-range molecular dynamics. J. Comp. Phys. 117 (1995) 1–19.2L. Verlet, Computer experiments on classical fluids: I. Thermodynamical properties of Lennard–Jones molecules, Phys. Rev. 159 (1967) 98–103.
Abstract
We present the discovery of a new optical/X-ray source likely associated with the Fermi
γ
-ray source 4FGL J1408.6–2917. Its high-amplitude periodic optical variability, large spectroscopic ...radial-velocity semiamplitude, evidence for optical emission lines and flaring, and X-ray properties together imply the source is probably a new black widow millisecond pulsar binary. We compile the properties of the 41 confirmed and suspected field black widows, finding a median secondary mass of 0.027 ± 0.003
M
⊙
. Considered jointly with the more massive redback millisecond pulsar binaries, we find that the “spider” companion mass distribution remains strongly bimodal, with essentially zero systems having companion masses of between ∼0.07 and 0.1
M
⊙
. X-ray emission from black widows is typically softer and less luminous than in redbacks, consistent with less efficient particle acceleration in the intrabinary shock in black widows, excepting a few systems that appear to have more efficient “redback-like” shocks. Together black widows and redbacks dominate the census of the fastest spinning field millisecond pulsars in binaries with known companion types, making up ≳80% of systems with
P
spin
< 2 ms. Similar to redbacks, the neutron star masses in black widows appear on average significantly larger than the canonical 1.4
M
⊙
, and many of the highest-mass neutron stars claimed to date are black widows with
M
NS
≳ 2.1
M
⊙
. Both of these observations are consistent with an evolutionary picture where spider millisecond pulsars emerge from short orbital period progenitors that had a lengthy period of mass transfer initiated while the companion was on the main sequence, leading to fast spins and high masses.
Abstract
Following the discovery of radio pulsars at the position of Fermi-LAT unassociated sources by the TRAPUM group, we conduct Swift-XRT observations of six of those 4FGL sources to determine if ...any pulsar-like X-ray sources are present and to confirm the reported detection of an X-ray counterpart via eROSITA at 4FGL J1803.1−6708. At two of the six targets, we detect no X-ray sources at the TRAPUM radio position, placing an upper limit on the 0.3–10.0 keV flux. At 4FGL J1803.1−6708 we find an X-ray source at the TRAPUM and eROSITA position. At 4FGL J1858.3−5424 we find a new X-ray counterpart at the TRAPUM position with signal-to-noise ratio S/N = 4.17, but also detect a distinct and separate X-ray source. At 4FGL J1823.8−3544 and 4FGL J1906.4−1757 we detect no X-ray flux at the TRAPUM positions, but we do detect separate X-ray sources elsewhere in the Fermi error ellipse. At these last two targets, our newly detected Swift sources are possible alternatives to the radio pulsar associations proposed by TRAPUM. Our findings confirm several of the discoveries reported by the TRAPUM group but suggest that further observations and investigations are necessary to confirm the low-energy counterpart of several unassociated sources.
Abstract
We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior ...Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be
12.39
−
0.98
+
1.30
km and
2.072
−
0.066
+
0.067
M
⊙
respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is
log
10
(
T
K
)
=
5.99
−
0.06
+
0.05
for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.
Abstract
Gamma-ray binary systems, a subclass of high-mass X-ray binaries, show nonthermal emissions from radio to TeV. While efficient electron acceleration is considered to take place in them, the ...nature of the acceleration mechanism and the physical environments in these systems have been a long-standing question. In this work, we report on long-term recurrent patterns in the short-term variability of the soft X-ray emission of LS 5039, one of the brightest gamma-ray binary systems. The Neutron star Interior Composition Explorer (NICER) observed LS 5039 four times from 2018 to 2021. By comparing them with the previous Suzaku and NuSTAR long-exposure observations, we studied the long-term evolution of the orbital light curve in the soft X-ray band. Although the observations by NICER and Suzaku are separated by ∼14 yr, i.e., more than 10
3
orbits, the orbital light curves show remarkable consistency after calculating their running averages with a window width ≳70 ks. Furthermore, all of the light curves show short-term variability with a timescale of ∼10 ks. Since the column density did not vary when the flux changed abruptly, such a short-term variability seems to be an intrinsic feature of the X-ray emission. We propose that the short-term variability is caused by clumps (or inhomogeneities) of the companion star wind impacting the X-ray production site. The observed timescale matches well with the lifetime of the clumps interacting with the pulsar wind and the dynamical timescale of the relativistic intrabinary shock in the pulsar wind scenario.
PSR J0537−6910, also known as the Big Glitcher, is the most prolific glitching pulsar known, and its spin-induced pulsations are only detectable in X-ray. We present results from analysis of 2.7 yr ...of NICER timing observations, from 2017 August to 2020 April. We obtain a rotation phase-connected timing model for the entire time span, which overlaps with the third observing run of LIGO/Virgo, thus enabling the most sensitive gravitational wave searches of this potentially strong gravitational wave-emitting pulsar. We find that the short-term braking index between glitches decreases towards a value of 7 or lower at longer times since the preceding glitch. By combining NICER and RXTE data, we measure a long-term braking index n = −1.25 ± 0.01. Our analysis reveals eight new glitches, the first detected since 2011, near the end of RXTE, with a total NICER and RXTE glitch activity of 8.88 × 10−7 yr−1. The new glitches follow the seemingly unique time-to-next-glitch–glitch-size correlation established previously using RXTE data, with a slope of 5 d μHz−1. For one glitch around which NICER observes 2 d on either side, we search for but do not see clear evidence of spectral nor pulse profile changes that may be associated with the glitch.