Engineering topics

Harding, Alice K.

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

At least 19 records

DOE OSTI2024

A Growing Braking Index and Spin-down Swings for the Pulsar PSR B0540–69

The way pulsars spin down is not understood in detail, but a number of possible physical mechanisms produce a spin-down rate that scales as a power of the rotation rate ($\dot{v}$ α-$v$ n ), with the power-law index n called the braking index. PSR B0540−69 is a pulsar that in 2011, after 16 yr of spinning down with a constant braking index of 2.1, experienced a giant spin-down change and a reduction of its braking index to nearly zero. Here, we show that, following this episode, the braking index monotonically increased during a period of at least 4 yr and stabilized at ∼1.1. We also present an alternative interpretation of a more modest rotational irregularity that occurred in 2023, which was modeled as an anomalous negative step of the rotation rate. Our analysis shows that the 2023 observations can be equally well described as a transient swing of the spin-down rate (lasting ∼65 days), and the Bayesian evidence indicates that this model is strongly preferred.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Rapid spin changes around a magnetar fast radio burst

Magnetars are neutron stars with extremely high magnetic fields (≳10 14 gauss) that exhibit various X-ray phenomena such as sporadic subsecond bursts, long-term persistent flux enhancements and variable rotation-period derivative. In 2020, a fast radio burst (FRB), akin to cosmological millisecond-duration radio bursts, was detected from the Galactic magnetar SGR 1935+2154, confirming the long-suspected association between some FRBs and magnetars. However, the mechanism for FRB generation in magnetars remains unclear. Here we report the X-ray observation of two glitches in SGR 1935+2154 within a time interval of approximately nine hours, bracketing an FRB that occurred on 14 October 2022. Each glitch involved a significant increase in the magnetar’s spin frequency, being among the largest abrupt changes in neutron-star rotation observed so far. Between the glitches, the magnetar exhibited a rapid spin-down phase, accompanied by an increase and subsequent decline in its persistent X-ray emission and burst rate. Here we postulate that a strong, ephemeral, magnetospheric wind provides the torque that rapidly slows the star’s rotation. The trigger for the first glitch couples the star’s crust to its magnetosphere, enhances the various X-ray signals and spawns the wind that alters magnetospheric conditions that might produce the FRB.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Exploring the MeV sky with a combined coded mask and Compton telescope: the Galactic Explorer with a Coded aperture mask Compton telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

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 ${\mathrm{log}}_{10}(T\,[{\rm{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. Furthermore, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The Multipolar Magnetic Field of the Millisecond Pulsar PSR J0030+0451

Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed at constraining the neutron star mass and radius, has inferred surface hot spots (the magnetic polar caps) that imply significantly nondipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole-plus-quadrupole components using a static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. and Riley et al., we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov Chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions of the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous γ-ray light curves following Kalapotharakos et al.; these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the γ-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and γ-ray light curves. We conclude with a discussion of the broader implications of our study.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2019

A Fundamental Plane for Gamma-Ray Pulsars

We show that the γ-ray pulsar observables, i.e., their total γ-ray luminosity, L(sub γ), spectral cutoff energy, ϵ (sub cut), stellar surface magnetic field, B(sub ⋆), and spin-down power ε, obey a relation of the form L(sub γ) = f (ϵ (sub cut), B(sub ⋆), ε ), which represents a 3D plane in their 4D logspace. Fitting the data of 88 pulsars of the second Fermi pulsar catalog, we show this relation to be L(sub γ) ∝ ϵ (sub cut)(sup 1.18±0.24 B(sub ⋆)(sup 0.17± 0.05) ε(sup 0.41±0.08), a pulsar fundamental plane (FP). We show that the observed FP is remarkably close to the theoretical relation L(sub γ) ∝ ϵ (sub cut)(sup 4/3) B(sub ⋆)(sup 1/6) ε(sup 5/12) obtained assuming that the pulsar γ-ray emission is due to curvature radiation by particles accelerated at the pulsar equatorial current sheet just outside the light cylinder. Interestingly, the FP seems incompatible with emission by synchrotron radiation. The corresponding scatter about the FP is ∼0.35 dex and can only partly be explained by the observational errors while the rest is probably due to the variation of the inclination and observer angles. We predict also that ϵ (sub cut) ∝ ε (sup 7/16) toward low ε for both young and millisecond pulsars implying that the observed death line of γ-ray pulsars is due to ϵ (sub cut) dropping below the Fermi band. Our results provide a comprehensive interpretation of the observations of γ-ray pulsars, setting requirements for successful theoretical modeling.

Kalapotharakos, Constantinos↗

NASA NTRS2019

Astro2020 Science White Paper: Prospects for Pulsar Studies at MeV Energies

Enabled by the Fermi Large Area Telescope, we now know young and recycled pulsars fill the gamma-ray sky, and we are beginning to understand their emission mechanism and their distribution throughout the Galaxy. However, key questions remain: Is there a large population of pulsars near the Galactic center? Why do the most energetic pulsars shine so brightly in MeV gamma rays but not always at GeV energies? What is the source and nature of the pair plasma in pulsar magnetospheres, and what role does the polar cap accelerator play? Addressing these questions calls for a sensitive, wide-field MeV telescope, which can detect the population of MeV-peaked pulsars hinted at by Fermi and hard X-ray telescopes and characterize their spectral shape and polarization.

Harding, Alice K.↗

NASA NTRS2019

The Physics of Accretion Onto Highly Magnetized Neutron Stars

Studying the physical processes occurring in the region just above the magnetic polesof strongly magnetized, accreting binary neutron stars is essential to our understanding of stellarand binary system evolution. Perhaps more importantly, it provides us with a natural laboratoryfor studying the physics of high temperature and density plasmas exposed to extreme radiation,gravitational, and magnetic fields. Observations over the past decade have shed new light on themanner in which plasma falling at near the speed of light onto a neutron star surface is halted. Recentadvances in modeling these processes have resulted in direct measurement of the magnetic fieldsand plasma properties. On the other hand, numerous physical processes have been identified thatchallenge our current picture of how the accretion process onto neutron stars works. Observationand theory are our essential tools in this regime because the extreme conditions cannot be duplicatedon Earth. This white paper gives an overview of the current theory, the outstanding theoreticaland observational challenges, and the importance of addressing them in contemporary astrophysicsresearch.

Wolff, Michael T.↗

NASA NTRS2019

Magnetars as Astrophysical Laboratories of Extreme Quantum Electrodynamics: The Case for a Compton Telescope

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

Wadiasingh, Zorawar↗

DOE OSTI2017

X-Ray and Radio Observations of the Magnetar SGR J1935+2154 during Its 2014, 2015, and 2016 Outbursts

We analyzed broadband X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a blackbody+power-law (BB+PL) or 2BB model during all three outbursts. Nuclear Spectroscopic Telescope Array observations revealed a hard X-ray tail, with a PL photon index Γ = 0.9, extending up to 50 keV, with flux comparable to the one detected <10 keV. Imaging analysis of Chandra data did not reveal small-scale extended emission around the source. Following the outbursts, the total 0.5–10 keV flux from SGR J1935+2154 increased in concordance to its bursting activity, with the flux at activation onset increasing by a factor of ∼7 following its strongest 2016 June outburst. A Swift /X-Ray Telescope observation taken 1.5 days prior to the onset of this outburst showed a flux level consistent with quiescence. We show that the flux increase is due to the PL or hot BB component, which increased by a factor of 25 compared to quiescence, while the cold BB component kT = 0.47 keV remained more or less constant. The 2014 and 2015 outbursts decayed quasi-exponentially with timescales of ∼40 days, while the stronger 2016 May and June outbursts showed a quick short-term decay with timescales of about four days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 μ Jy for the 4.6 GHz observations and 7 μ Jy for the 1.4 GHz observations. We discuss these results in the framework of the current magnetar theoretical models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Fermi Gamma-Ray Pulsars: Understanding the High-energy Emission from Dissipative Magnetospheres

Based on the Fermi observational data, we reveal meaningful constraints for the dependence of the macroscopic conductivity (σ) of dissipative pulsar magnetosphere models on the corresponding spin-down rate, E-dot . Our models are refinements of the FIDO (Force-free Inside, Dissipative Outside) models, which have dissipative regions that are restricted on the equatorial current sheet outside of the the light-cylinder. Taking into account the observed cutoff energies of all of the Fermi pulsars and assuming that (a) the corresponding γ-ray pulsed emission is due to curvature radiation at the radiation-reaction-limit regime, and (b) this emission is produced at the equatorial current sheet near the light cylinder, we show that the Fermi data provide clear indications about the corresponding accelerating electric-field components. A direct comparison between the Fermi cutoff energies and the model ones reveals that σ increases with E-dot for high E-dot -values, while it saturates for low ones. This comparison indicates also that the corresponding gap width increases toward low E-dot -values. Assuming the Goldreich–Julian flux for the emitting particles, we calculate the total γ-ray luminosity (L {sub γ}). A comparison between the dependence of the Fermi L {sub γ}-values and the model ones on E-dot indicates an increase of the emitting particle multiplicity with E-dot . Our modeling, guided by the Fermi data alone, enhances our understanding of the physical mechanisms behind the high-energy emission in pulsar magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2017

Fermi Gamma-Ray Pulsars: Understanding the High-Energy Emission from Dissipative Magnetospheres

Based on the Fermi observational data, we reveal meaningful constraints for the dependence of the macroscopic conductivity (sigma) of dissipative pulsar magnetosphere models on the corresponding spin-down rate, epsilon(sup dot). Our models are refinements of the FIDO (Force-free Inside, Dissipative Outside) models, which have dissipative regions that are restricted on the equatorial current sheet outside of the the light-cylinder. Taking into account the observed cutoff energies of all of the Fermi pulsars and assuming that (a) the corresponding gamma-ray pulsed emission is due to curvature radiation at the radiation-reaction-limit regime, and (b) this emission is produced at the equatorial current sheet near the light cylinder, we show that the Fermi data provide clear indications about the corresponding accelerating electric-field components. A direct comparison between the Fermi cutoff energies and the model ones reveals that if sigma increases with epsilon(sup dot) for high epsilon(sup dot)-values, while it saturates for low ones. This comparison indicates also that the corresponding gap width increases toward low epsilon(sup dot)-values. Assuming the Goldreich-Julian flux for the emitting particles, we calculate the total gamma-ray luminosity (L(sub gamma)). A comparison between the dependence of the Fermi L(sub gamma)-values and the model ones on epsilon(sup dot) indicates an increase of the emitting particle multiplicity with epsilon(sup dot). Our modeling, guided by the Fermi data alone, enhances our understanding of the physical mechanisms behind the high-energy emission in pulsar magnetospheres.

Kalapotharakos, Constantinos↗

DOE OSTI2017

Multiwavelength Polarization of Rotation-powered Pulsars

Polarization measurements provide strong constraints on models for emission from rotation-powered pulsars. We present multiwavelength polarization predictions showing that measurements over a range of frequencies can be particularly important for constraining the emission location, radiation mechanisms, and system geometry. The results assume a generic model for emission from the outer magnetosphere and current sheet in which optical to hard X-ray emission is produced by synchrotron radiation (SR) from electron–positron pairs and γ -ray emission is produced by curvature radiation (CR) or SR from accelerating primary electrons. The magnetic field structure of a force-free magnetosphere is assumed and the phase-resolved and phase-averaged polarization is calculated in the frame of an inertial observer. We find that large position angle (PA) swings and deep depolarization dips occur during the light-curve peaks in all energy bands. For synchrotron emission, the polarization characteristics are strongly dependent on photon emission radius with larger, nearly 180°, PA swings for emission outside the light cylinder (LC) as the line of sight crosses the current sheet. The phase-averaged polarization degree for SR is less that 10% and around 20% for emission starting inside and outside the LC, respectively, while the polarization degree for CR is much larger, up to 40%–60%. Observing a sharp increase in polarization degree and a change in PA at the transition between X-ray and γ -ray spectral components would indicate that CR is the γ -ray emission mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2017

Multiwavelength Polarization of Rotation-Powered Pulsars

Polarization measurements provide strong constraints on models for emission from rotation-powered pulsars. We present multiwavelength polarization predictions showing that measurements over a range of frequencies can be particularly important for constraining the emission location, radiation mechanisms, and system geometry. The results assume a generic model for emission from the outer magnetosphere and current sheet in which optical to hard X-ray emission is produced by synchrotron radiation (SR) from electron-positron pairs and gamma-ray emission is produced by curvature radiation (CR) or SR from accelerating primary electrons. The magnetic field structure of a force-free magnetosphere is assumed and the phase-resolved and phase-averaged polarization is calculated in the frame of an inertial observer. We find that large position angle (PA) swings and deep depolarization dips occur during the light-curve peaks in all energy bands. For synchrotron emission, the polarization characteristics are strongly dependent on photon emission radius with larger, nearly 180deg, PA swings for emission outside the light cylinder (LC)‚ as the line of sight crosses the current sheet. The phase-averaged polarization degree for SR is less that 10% and around 20% for emission starting inside and outside the LC, respectively, while the polarization degree for CR is much larger, up to 40%-60%. Observing a sharp increase in polarization degree and a change in PA at the transition between X-ray and gamma-ray spectral components would indicate that CR is the gamma-ray emission mechanism.

Harding, Alice K.↗

DOE OSTI2017

Constraining Relativistic Bow Shock Properties in Rotation-powered Millisecond Pulsar Binaries

Multiwavelength follow-up of unidentified Fermi sources has vastly expanded the number of known galactic-field “black widow” and “redback” millisecond pulsar binaries. Focusing on their rotation-powered state, we interpret the radio to X-ray phenomenology in a consistent framework. We advocate the existence of two distinct modes differing in their intrabinary shock orientation, distinguished by the phase centering of the double-peaked X-ray orbital modulation originating from mildly relativistic Doppler boosting. By constructing a geometric model for radio eclipses, we constrain the shock geometry as functions of binary inclination and shock standoff R {sub 0}. We develop synthetic X-ray synchrotron orbital light curves and explore the model parameter space allowed by radio eclipse constraints applied on archetypal systems B1957+20 and J1023+0038. For B1957+20, from radio eclipses the standoff is R {sub 0} ∼ 0.15–0.3 fraction of binary separation from the companion center, depending on the orbit inclination. Constructed X-ray light curves for B1957+20 using these values are qualitatively consistent with those observed, and we find occultation of the shock by the companion as a minor influence, demanding significant Doppler factors to yield double peaks. For J1023+0038, radio eclipses imply R {sub 0} ≲ 0.4, while X-ray light curves suggest 0.1 ≲ R {sub 0} ≲ 0.3 (from the pulsar). Degeneracies in the model parameter space encourage further development to include transport considerations. Generically, the spatial variation along the shock of the underlying electron power-law index should yield energy dependence in the shape of light curves, motivating future X-ray phase-resolved spectroscopic studies to probe the unknown physics of pulsar winds and relativistic shock acceleration therein.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

Synchrotron Self-Compton Emission from the Crab and Other Pulsars

Results of a simulation of synchrotron-self Compton (SSC) emission from a rotation-powered pulsar are presented. The radiating particles are assumed to be both accelerated primary electrons and a spectrum of electron-positron pairs produced in cascades near the polar cap. They follow trajectories in a slot gap using 3D force-free magnetic field geometry, gaining pitch angles through resonant cyclotron absorption of radio photons, radiating and scattering synchrotron emission at high altitudes out to and beyond the light cylinder. Full angular dependence of the synchrotron photon density is simulated in the scattering and all processes are treated in the inertial observer frame. Spectra for the Crab and Vela pulsars as well as two energetic millisecond pulsars, B1821-24 and B1937+21 are simulated using this model. The simulation of the Crab pulsar radiation can reproduce both the flux level and the shape of the observed optical to hard X-ray emission assuming a pair multiplicity of M+ = 3x10(exp 5), as well as the very-high- energy emission above 50 GeV detected by MAGIC and VERITAS, with both the synchrotron and SSC components reflecting the shape of the pair spectrum. Simulations of Vela, B1821-24 and B1937+21, for M+ up to 10(exp 5), do not produce pair SSC emission that is detectable by current telescopes, indicating that only Crab-like pulsars produce significant SSC components. The pair synchrotron emission matches the observed X-ray spectrum of the millisecond pulsars and the predicted peak of this emission at 1-10 MeV would be detectable with planned Compton telescopes.

Self-Compton↗

NASA NTRS2015

Testing Dissipative Magnetosphere Model Light Curves and Spectra with Fermi Pulsars

We explore the emission properties of a dissipative pulsar magnetosphere model introduced by Kalapotharakos et al. comparing its high-energy light curves and spectra, due to curvature radiation, with data collected by the Fermi LAT. The magnetosphere structure is assumed to be near the force-free solution. The accelerating electric field, inside the light cylinder (LC), is assumed to be negligible, while outside the LC it rescales with a finite conductivity (sigma). In our approach we calculate the corresponding high-energy emission by integrating the trajectories of test particles that originate from the stellar surface, taking into account both the accelerating electric field components and the radiation reaction forces. First, we explore the parameter space assuming different value sets for the stellar magnetic field, stellar period, and conductivity. We show that the general properties of the model are in a good agreement with observed emission characteristics of young gamma-ray pulsars, including features of the phase-resolved spectra. Second, we find model parameters that fit each pulsar belonging to a group of eight bright pulsars that have a published phase-resolved spectrum. The sigma values that best describe each of the pulsars in this group show an increase with the spin-down rate (E˙ ) and a decrease with the pulsar age, expected if pair cascades are providing the magnetospheric conductivity. Finally, we explore the limits of our analysis and suggest future directions for improving such models.

acceleration of particles↗

NASA NTRS2014

Gamma-Ray Emission in Dissipative Pulsar Magnetospheres: from Theory to Fermi Observations

We compute the patterns of gamma-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed gamma-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model gamma-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and gamma-ray emission on the gamma-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the gamma-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.

stars↗