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↗

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↗

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↗

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.↗

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↗

NASA NTRS2014

The Neutron Star Zoo

Neutron stars are a very diverse population, both in their observational and their physical properties. They prefer to radiate most of their energy at X-ray and gamma-ray wavelengths. But whether their emission is powered by rotation, accretion, heat, magnetic fields or nuclear reactions, they are all different species of the same animal whose magnetic field evolution and interior composition remain a mystery. This article will broadly review the properties of inhabitants of the neutron star zoo, with emphasis on their high-energy emission. XXX Neutron stars are found in a wide variety of sources, displaying an amazing array of behavior. They can be isolated or in binary systems, accreting, heating, cooling, spinning down, spinning up, pulsing, flaring and bursting. The one property that seems to determine their behavior most strongly is their magnetic field strength, structure and evolution. The hot polar caps, bursts and flares of magnetars are likely due to the rapid decay and twisting of their superstrong magnetic fields, whose very existence requires some kind of early dynamo activity. The intermediate-strength magnetic fields of RPPs determines their spin-down behavior and radiation properties. However, the overlap of the magnetar and RPP populations is not understood at present. Why don't high-field RPPs burst or flare? Why don't lower-field magnetars sometimes behave more like RPPs? INS may be old magnetars whose high fields have decayed, but they do not account for the existence of younger RPPs with magnetar-strength fields. Not only the strength of the magnetic field but also its configuration may be important in making a NS a magnetar or a RPP. Magnetic field decay is a critical link between other NS populations as well. "Decay" of the magnetic field is necessary for normal RPPs to evolve into MSPs through accretion and spin up in LMXBs. Some kind of accretion-driven field reduction is the most likely mechanism, but it is controversial since it is not clear how effective it is or on what timescale a buried field might re-emerge. One piece of evidence in favor of accretion-driven field reduction is the fact that NSs in LMXBs, which are older systems (> 108 yr), have mostly low fields and NSs in HMXBs, which are younger systems (107 − 108 yr), have higher fields. This may be an indication that accretion-driven field reduction or decay has not had enough time to operate in HMXBs but has in LMXBs. However, there does not seem to be any evidence of decaying fields in either the LMXB or HMXB populations; e.g. smaller magnetic fields in older systems. On the other hand, CCOs are very young so if they acquired their low fields through mass fallback accretion, the field submergence would have had to operate on much faster timescales than it apparently does in LMXBs. But as we continue to find new species in the NS zoo, one of these may someday be the "Rosetta Stone" that will give us the clues for solving these puzzles.

Zoo↗

NASA NTRS2013

The VELA-X-Pulsar Wind Nebula Revisited with Four Years of Fermi Large Area Telescope Observations

The Vela supernova remnant (SNR) is the closest SNR to Earth containing an active pulsar, the Vela pulsar (PSR B0833-45). This pulsar is an archetype of the middle-aged pulsar class and powers a bright pulsar wind nebula (PWN), Vela-X, spanning a region of 2deg × 3deg south of the pulsar and observed in the radio, X-ray, and very high energy γ-ray domains. The detection of the Vela-X PWN by the Fermi Large Area Telescope (LAT) was reported in the first year of the mission. Subsequently, we have reinvestigated this complex region and performed a detailed morphological and spectral analysis of this source using 4 yr of Fermi-LAT observations. This study lowers the threshold for morphological analysis of the nebula from 0.8 GeV to 0.3 GeV, allowing for the inspection of distinct energy bands by the LAT for the first time. We describe the recent results obtained on this PWN and discuss the origin of the newly detected spatial features.

Vela-X↗

NASA NTRS2012

Pulsar Polar Cap and Slot Gap Models: Confronting Fermi Data

Rotation-powered pulsars are excellent laboratories for studying particle acceleration as well as fundamental physics of strong gravity, strong magnetic fields and relativity. I will review acceleration and gamma-ray emission from the pulsar polar cap and slot gap. Predictions of these models can be tested with the data set on pulsars collected by the Large Area Telescope on the Fermi Gamma-Ray Telescope over the last four years, using both detailed light curve fitting and population synthesis.

Harding, Alice K.↗

NASA NTRS2012

Crab Flares and Magnetic Reconnection in Pulsar Winds

The striped winds of rotation-powered pulsars are ideal sites for magnetic reconnection. The magnetic fields of the wind near the current sheet outside the light cylinder alternate polarity every pulsar period and eventually encounter a termination shock. Magnetic reconnection in the wind has been proposed as a mechanism for transferring energy from electromagnetic fields to particles upstream of the shock (the "sigma" problem), but it is not clear if, where and how this occurs. Fermi and AGILE have recently observed powerful gamma-ray flares from the Crab nebula, which challenge traditional models of acceleration at the termination shock. New simulations are revealing that magnetic reconnection may be instrumental in understanding the Crab flares and in resolving the "sigma" problem in pulsar wind nebulae.

Harding, Alice K.↗

NASA NTRS2012

Pulsar Emission Geometry and Accelerating Field Strength

The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems, The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA I pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models. within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle. maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius. under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry

DeCesar, Megan E.↗

NASA NTRS2012

Pulsar Pair Cascades in Magnetic Fields with Offset Polar Caps

Neutron star magnetic fields may have polar caps (PC) that are offset from the dipole axis, through field-line sweepback near the light cylinder or non-symmetric currents within the star. The effects of such offsets on electron-positron pair cascades are investigated, using simple models of dipole magnetic fields with small distortions that shift the PCs by different amounts or directions. Using a Monte Carlo pair cascade simulation, we explore the changes in the pair spectrum, multiplicity and energy flux across the PC, as well as the trends in pair flux and pair energy flux with spin-down luminosity, L(sub sd). We also give an estimate of the distribution of heating flux from returning positrons on the PC for different offsets. We find that even modest offsets can produce significant increases in pair multiplicity, especially for pulsars that are near or beyond the pair death lines for centered PCs, primarily because of higher accelerating fields. Pair spectra cover several decades in energy, with the spectral range of millisecond pulsars (MSPs) two orders of magnitude higher than for normal pulsars, and PC offsets allow significant extension of all spectra to lower pair energies. We find that the total PC pair luminosity L(sub pair) is proportional to L(sub sd), with L(sub pair) approximates 10(exp -3) L(sub sd) for normal pulsars and L(sub pair) approximates 10(exp -2) L(sub sd) for MSPs. Remarkably, the total PC heating luminosity for even large offsets increases by less than a factor of two, even though the PC area increases by much larger factors, because most of the heating occurs near the magnetic axis.

Harding, Alice K.↗

NASA NTRS2012

Pair Cascades and Deathlines in Magnetic Fields with Offset Polar Caps

We present results of electron-positron pair cascade simulations in a dipole magnetic field whose polar cap is offset from the dipole axis. In such a field geometry, the polar cap is displaced a small fraction of the neutron star radius from the star symmetry axis and the field line radius of curvature is modified. Using the modified parallel electric field near the offset polar cap, we simulate pair cascades to determine the pair deathlines and pair multiplicities as a function of the offset. We find that the pair multiplicity can change dr;unatically with a modest offset, with a significant increase on one side of the polar cap. Lower pair deathlines allow a larger fraction of the pulsar population, that include old and millisecond pulsars, to produce cascades with high multiplicity. The results have some important implications for pulsar particle production, high-energy emission and cosmic-ray contribution.

Harding, Alice K.↗

NASA NTRS2011

Gamma-Ray Pulsar Light Curves in Vacuum and Force-Free Geometry

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profiles. We find that, compared to the profiles derived from symmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines, increases significantly relative to the off-peak emission, formed along leading field lines. The enhanced contrast produces improved slot gap model fits to Fermi pulsar light curves like Vela, with vacuum dipole fits being more favorable.

Harding, Alice K.↗