Engineering topics

Harding, A. K.

Publications and source records attributed to Harding, A. K..

At least 19 records

DOE OSTI2025

A Systematic Search for MeV–GeV Pulsar Wind Nebulae without Gamma-Ray Detected Pulsars

An increasing number of pulsar wind nebulae (PWNe) are being identified in the TeV band by ground-based Imaging Air Cherenkov Telescopes such that they constitute the dominant source class of Galactic TeV emitters. However, MeV–GeV PWN counterparts are still largely lacking. To date, only a dozen PWNe are identified by the Fermi–Large Area Telescope (LAT) in the MeV–GeV band. Most PWNe are located along the Galactic plane embedded within the prominent, diffuse Galactic γ-ray emission, which makes these sources difficult to disentangle from the bright diffuse background. We present a systematic search for γ-ray counterparts to known PWNe in the 300 MeV–2 TeV energy band using the Fermi–LAT. We target the locations of previously identified PWNe that lack detected Fermi–LAT pulsars to minimize associated pulsar contamination. The sample includes six previously identified Fermi PWNe and eight Fermi–LAT sources associated with PWNe. We report the analysis of 58 regions of interest and classify Fermi–LAT detected sources as either a likely PWN or a candidate PWN counterpart based on their morphological and spectral characteristics across the broadband spectrum. There are nine unidentified Fermi–LAT sources that we consider as likely PWN counterparts, which, if confirmed to be PWNe, would greatly increase the PWN population detected by the Fermi–LAT from 12 to 21. The remaining Fermi–LAT detected sources are considered weaker PWN candidates. A second approach in the systematic search for γ-ray emitting PWNe will involve studying the off-pulse phases of Fermi–LAT detected pulsars for the presence of an obscured PWN and will be reported in a subsequent paper.

Acharyya, A. [Univ. of Southern Denmark, Odense (D↗

DOE OSTI2023

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Broadband X-ray burst spectroscopy of the fast-radio-burst-emitting Galactic magnetar

Magnetars are young, magnetically powered neutron stars that possess the strongest magnetic fields in the Universe. Fast radio bursts (FRBs) are extremely intense millisecond-long radio pulses of primarily extragalactic origin, and a leading attribution for their genesis focuses on magnetars. A hallmark signature of magnetars is their emission of bright, hard X-ray bursts of sub-second duration. On 27 April 2020, the Galactic magnetar SGR J1935+2154 emitted hundreds of X-ray bursts within a few hours. One of these temporally coincided with an FRB, the first known detection of an FRB from the Milky Way. Here, we present spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with the Neutron Star Interior Composition Explorer (NICER) mission of the National Aeronautics and Space Administration and with the Fermi Gamma-ray Burst Monitor (GBM) mission in their combined energy range of 0.2 keV to 30 MeV. These broadband spectra permit direct comparison with the spectrum of the FRB-associated X-ray burst (FRB-X). We demonstrate that all 24 NICER and GBM bursts are very similar temporally to the FRB-X, but strikingly different spectrally. The singularity of the FRB-X burst is perhaps indicative of an uncommon locale for its origin. We suggest that this event originated in quasi-polar open or closed magnetic field lines that extend to high altitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation

Here, we report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius R eq are, respectively, ${1.34}_{-0.16}^{+0.15}\,{M}_{\odot }$ and ${12.71}_{-1.19}^{+1.14}\,\mathrm{km}$, while the compactness ${GM}/{R}_{\mathrm{eq}}{c}^{2}={0.156}_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state (EoS) of this dense matter is to measure both a star's equatorial circumferential radius R e and its gravitational mass M. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. In this work, our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are ${R}_{e}={13.02}_{-1.06}^{+1.24}$ km and $M={1.44}_{-0.14}^{+0.15}\,{M}_{\odot }$ (68%). The independent analysis reported in the companion paper by Riley et al. explores different emitting spot models, but finds spot shapes and locations and estimates of R e and M that are consistent with those found in this work. We show that our measurements of R e and M for PSR J0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2019

On the Maximum Pair Multiplicity of Pulsar Cascades

We study electron–positron pair production in polar caps of energetic pulsars to determine the maximum multiplicity of pair plasma a pulsar can produce under the most favorable conditions. This paper complements and updates our study of pair cascades presented in Timokhin & Harding (2015) with a more accurate treatment of the effects of ultrastrong B approximately greater than 3 x 10 (exp 12) G magnetic fields and emission processes of primary and secondary particles. We include pairs produced by curvature and synchrotron radiation photons as well as resonant Compton-scattered photons. We develop a semianalytical model of electron–positron cascades that can efficiently simulate pair cascades with an arbitrary number of microphysical processes and use it to explore cascade properties for a wide range of pulsar parameters. We argue that the maximum cascade multiplicity cannot exceed approximately a few x 10 (exp 5) and that the multiplicity has a rather weak dependence on pulsar period. The highest multiplicity is achieved in pulsars with magnetic field 4 x 10 (exp 12) is approximately greater than B is approximately greater than 10 (exp 13) G and hot surfaces, with T is approximately greater than 10 (exp 6)K. We also derive analytical expressions for several physical quantities relevant for electromagnetic cascade in pulsars, which may be useful in future works on pulsar cascades, including the upper limit on cascade multiplicity and various approximations for the parameter χ, the exponential factor in the expression for photon attenuation in strong magnetic fields.

Timokhin, A. N.↗

DOE OSTI2018

Fermi-LAT Observations of LIGO/Virgo Event GW170817

We present the Fermi Large Area Telescope (LAT) observations of the binary neutron star merger event GW170817 and the associated short gamma-ray burst (SGRB) GRB 170817A detected by the Fermi Gamma-ray Burst Monitor. The LAT was entering the South Atlantic Anomaly at the time of the LIGO/Virgo trigger (t GW ) and therefore cannot place constraints on the existence of high-energy (E > 100 MeV) emission associated with the moment of binary coalescence. We focus instead on constraining high-energy emission on longer timescales. No candidate electromagnetic counterpart was detected by the LAT on timescales of minutes, hours, or days after the LIGO/Virgo detection. The resulting flux upper bound (at 95% C.L.) from the LAT is 4.5 × 10 -10 erg cm -2 s -1 in the 0.1–1 GeV range covering a period from t GW + 1153 s to t GW + 2027 s. At the distance of GRB 170817A, this flux upper bound corresponds to a luminosity upper bound of 9.7 × 10 43 erg s -1 , which is five orders of magnitude less luminous than the only other LAT SGRB with known redshift, GRB 090510. We also discuss the prospects for LAT detection of electromagnetic counterparts to future gravitational-wave events from Advanced LIGO/Virgo in the context of GW170817/GRB 170817A.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2018

NICER Discovers the Ultracompact Orbit of the Accreting Millisecond Pulsar IGR J17062-6143

We present results of recent Neutron Star Interior Composition Explorer (NICER) observations of the accreting millisecond X-ray pulsar (AMXP) IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38-minute orbital period. NICER observed the source for 26 kiloseconds over a 5.3-day span in 2017 August, and again for 14 and 11 kiloseconds in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hertz in Rossi X-ray Timing Explorer (RXTE) data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z squared method finds a best-fitting circular orbit with a period of 2278.21 seconds (37.97 minutes), a projected semimajor axis of 0.00390 lt-s (Localization Test Statistic), and a barycentric pulsar frequency of 163.6561105 Hertz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 times 10 (sup minus 8) solar mass, presently the smallest known for a stellar binary. The minimum donor mass ranges from approximately 0.005 to 0.007 times the solar mass for a neutron star mass from 1.2 to 2 times the solar mass. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175-0.0155 times the solar mass and 19 degrees less than i less than 27.5 degrees, where the lower and upper bounds correspond to 1.4 and 2 times the solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 plus or minus 0.11 percent (0.3-3.2 kiloelectronvolts).

Strohmayer, T. E.↗

DOE OSTI2017

3FHL: The Third Catalog of Hard Fermi -LAT Sources

We present a catalog of sources detected above 10 GeV by the Fermi Large Area Telescope (LAT) in the first 7 years of data using the Pass 8 event-level analysis. This is the Third Catalog of Hard Fermi -LAT Sources (3FHL), containing 1556 objects characterized in the 10 GeV–2 TeV energy range. The sensitivity and angular resolution are improved by factors of 3 and 2 relative to the previous LAT catalog at the same energies (1FHL). The vast majority of detected sources (79%) are associated with extragalactic counterparts at other wavelengths, including 16 sources located at very high redshift ( z > 2). Of the sources, 8% have Galactic counterparts and 13% are unassociated (or associated with a source of unknown nature). The high-latitude sky and the Galactic plane are observed with a flux sensitivity of 4.4 to 9.5 × 10{sup −11} ph cm{sup −2} s{sup −1}, respectively (this is approximately 0.5% and 1% of the Crab Nebula flux above 10 GeV). The catalog includes 214 new γ -ray sources. The substantial increase in the number of photons (more than 4 times relative to 1FHL and 10 times to 2FHL) also allows us to measure significant spectral curvature for 32 sources and find flux variability for 163 of them. Furthermore, we estimate that for the same flux limit of 10{sup −12} erg cm{sup −2} s{sup −1}, the energy range above 10 GeV has twice as many sources as the range above 50 GeV, highlighting the importance, for future Cherenkov telescopes, of lowering the energy threshold as much as possible.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Search for Extended Sources in the Galactic Plane Using Six Years of Fermi -Large Area Telescope Pass 8 Data above 10 GeV

The spatial extension of a γ -ray source is an essential ingredient to determine its spectral properties, as well as its potential multiwavelength counterpart. The capability to spatially resolve γ -ray sources is greatly improved by the newly delivered Fermi -Large Area Telescope (LAT) Pass 8 event-level analysis, which provides a greater acceptance and an improved point-spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7° from the Galactic plane, using 6 yr of Fermi -LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard Fermi -LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Cosmic-ray electron-positron spectrum from 7 GeV to 2 TeV with the Fermi Large Area Telescope

Here, we present a measurement of the cosmic-ray electron+positron spectrum between 7 GeV and 2 TeV performed with almost seven years of data collected with the Fermi Large Area Telescope. We find that the spectrum is well fit by a broken power law with a break energy at about 50 GeV. Above 50 GeV, the spectrum is well described by a single power law with a spectral index of 3.07 ± 0.02(stat + syst) ± 0.04(energy measurement). An exponential cutoff lower than 1.8 TeV is excluded at 95% CL.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Search for Cosmic-Ray Electron and Positron Anisotropies with Seven Years of Fermi Large Area Telescope Data

We present the Large Area Telescope on board the Fermi Gamma-ray Space Telescope that has collected the largest ever sample of high-energy cosmic-ray electron and positron events since the beginning of its operation. Potential anisotropies in the arrival directions of cosmic-ray electrons or positrons could be a signature of the presence of nearby sources. We use almost seven years of data with energies above 42 GeV processed with the Pass 8 reconstruction. The present data sample can probe dipole anisotropies down to a level of 10 -3 . We take into account systematic effects that could mimic true anisotropies at this level. We present a detailed study of the event selection optimization of the cosmic-ray electrons and positrons to be used for anisotropy searches. Since no significant anisotropies have been detected on any angular scale, we present upper limits on the dipole anisotropy. Lastly, the present constraints are among the strongest to date probing the presence of nearby young and middle-aged sources.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

SEARCHING THE GAMMA-RAY SKY FOR COUNTERPARTS TO GRAVITATIONAL WAVE SOURCES: FERMI GAMMA-RAY BURST MONITO R AND LARGE AREA TELESCOPE OBSERVATIONS OF LVT151012 AND GW151226

We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candidate LVT151012. At the time of the LIGO triggers on LVT151012 and GW151226, GBM was observing 68% and 83% of the localization regions, and LAT was observing 47% and 32%, respectively. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for characterizing the flux upper bounds across large areas of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, differences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

The effect of an offset polar cap dipolar magnetic field on the modeling of the vela pulsar's γ-ray light curves

We performed geometric pulsar light curve modeling using static, retarded vacuum, and offset polar cap (PC) dipole B-fields (the latter is characterized by a parameter ϵ), in conjunction with standard two-pole caustic (TPC) and outer gap (OG) emission geometries. The offset-PC dipole B-field mimics deviations from the static dipole (which corresponds to ϵ=0). In addition to constant-emissivity geometric models, we also considered a slot gap (SG) E-field associated with the offset-PC dipole B-field and found that its inclusion leads to qualitatively different light curves. Solving the particle transport equation shows that the particle energy only becomes large enough to yield significant curvature radiation at large altitudes above the stellar surface, given this relatively low E-field. Therefore, particles do not always attain the radiation-reaction limit. Our overall optimal light curve fit is for the retarded vacuum dipole field and OG model, at an inclination angle α=78{sub −1}{sup +1{sup {sup ∘}}} and observer angle ζ=69{sub −1}{sup +2{sup {sup ∘}}}. For this B-field, the TPC model is statistically disfavored compared to the OG model. For the static dipole field, neither model is significantly preferred. We found that smaller values of ϵ are favored for the offset-PC dipole field when assuming constant emissivity, and larger ϵ values favored for variable emissivity, but not significantly so. When multiplying the SG E-field by a factor of 100, we found improved light curve fits, with α and ζ being closer to best fits from independent studies, as well as curvature radiation reaction at lower altitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

A NEW, LOW BRAKING INDEX FOR THE LMC PULSAR B0540–69

We report the results of a 16 month monitoring campaign using the Swift satellite of PSR B0540–69, a young pulsar in the Large Magellanic Cloud. Phase connection was maintained throughout the campaign so that a reliable ephemeris could be determined, and the length of the campaign is adequate to accurately determine the spin frequency ν and its first and second derivatives. The braking index n is 0.031 ± 0.013 (90% confidence), a value much lower than previously reported for B0540–69 and almost all other young pulsars. We use data from the extensive monitoring campaign with Rossi X-ray Timing Explorer to show that timing noise is unlikely to significantly affect the measurement. This is the first measurement of the braking index in the pulsar's recently discovered high spin-down state. We discuss possible mechanisms for producing the low braking index.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2016

A New, Low Braking Index for the LMC Pulsar B0540-69

We report the results of a 16-month monitoring campaign using the Swift satellite of PSR 0540, a young pulsar in the Large Magellanic Cloud. Phase connection was maintained throughout the campaign so that a reliable ephemeris could be determined, and the length of the campaign is adequate to accurately determine the spin frequency and its first and second derivatives. The braking index n is 0.031 +/- 0.013 (90% confidence), a value much lower than previously reported for 0 540 and almost all other young pulsars. We use data from the extensive monitoring campaign with RXTE to showt hat timing noise is unlikely to significantly affect the measurement. This is the first measurement of the braking index in the pulsars recently discovered high spin-down state. We discuss possible mechanisms for producing the low braking index.

pulsars individual (PSR B0540âˆ’69)↗

DOE OSTI2016

3FGL demographics outside the galactic plane using supervised machine learning: pulsar and dark matter subhalo interpretations

Nearly one-third of the sources listed in the Third Fermi Large Area Telescope (LAT) catalog (3FGL) remain unassociated. It is possible that predicted and even unanticipated gamma-ray source classes are present in these data waiting to be discovered. Taking advantage of the excellent spectral capabilities achieved by the Fermi LAT, we use machine-learning classifiers (Random Forest and XGBoost) to pinpoint potentially novel source classes in the unassociated 3FGL sample outside the Galactic plane. Here we report a total of 34 high-confidence Galactic candidates at |b|⩾5{sup ∘}. The currently favored standard astrophysical interpretations for these objects are pulsars or low-luminosity globular clusters hosting millisecond pulsars (MSPs). Yet these objects could also be interpreted as dark matter annihilation taking place in ultra-faint dwarf galaxies or dark matter subhalos. Unfortunately, Fermi LAT spectra are not sufficient to break degeneracies between the different scenarios. Careful visual inspection of archival optical images reveals no obvious evidence for low-luminosity globular clusters or ultra-faint dwarf galaxies inside the 95% error ellipses. If these are pulsars, this would bring the total number of MSPs at |b|⩾5{sup ∘} to 106, down to an energy flux ≈4.0 × 10{sup −12} erg cm{sup −2} s{sup −1} between 100 MeV and 100 GeV. We find this number to be in excellent agreement with predictions from a new population synthesis of MSPs that predicts 100–126 high-latitude 3FGL MSPs depending on the choice of high-energy emission model. If, however, these are dark matter substructures, we can place upper limits on the number of Galactic subhalos surviving today and on dark matter annihilation cross sections. These limits are beginning to approach the canonical thermal relic cross section for dark matter particle masses below ∼100 GeV in the bottom quark (b b-bar ) annihilation channel.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2016

Gamma-Ray Pulsar Light Curves as Probes of Magnetospheric Structure

The large number of gamma-ray pulsars discovered by the Fermi Gamma-Ray Space Telescope since its launch in 2008 dwarfs the handful that were previously known. The variety of observed light curves makes possible a tomography of both the ensemble-averaged field structure and the high-energy emission regions of a pulsar magnetosphere. Fitting the gamma-ray pulsar light curves with model magnetospheres and emission models has revealed that most of the high-energy emission, and the particles acceleration, takes place near or beyond the light cylinder, near the current sheet. As pulsar magnetosphere models become more sophisticated, it is possible to probe magnetic field structure and emission that are self-consistently determined. Light curve modeling will continue to be a powerful tool for constraining the pulsar magnetosphere physics.

Harding, A. K.↗