Researcher profile

Alice K. Harding

Alice K. Harding contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

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 on and improves past and currently proposed missions at such energies. We outline the motivations for combining a coded mask and a Compton telescope and we define the scientific goals of 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 enable 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.

preprint2022arXiv

Incremental Fermi Large Area Telescope Fourth Source Catalog

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

preprint2022arXiv

The Fundamental Plane Relation for Gamma-Ray Pulsars Implied by 4FGL

We explore the validity of the recently reported fundamental plane (FP) relation of gamma-ray pulsars using 190 pulsars included in the latest 4FGL-DR3 catalog. This sample number is more than twice as large as that of the original study. The FP relation incorporates 4 parameters, i.e., the spin-down power, $\dot{\mathcal{E}}$, the surface magnetic field, $B_{\star}$, the total gamma-ray luminosity, $L_γ$, and a spectral cutoff energy, $ε_{\rm cut}$. The derivation of $ε_{\rm cut}$ is the most intriguing one because $ε_{\rm cut}$ depends on the proper interpretation of the available phase-averaged spectra. We construct synthetic phase-averaged spectra, guided by the few existing phase-resolved ones, to find that the best fit cutoff energy, $ε_{\rm c1}$, corresponding to a purely exponential cutoff (plus a power law) spectral form, is the parameter that optimally probes the maximum cutoff energy of the emission that originates from the core of the dissipative region, i.e., the equatorial current sheet. Computing this parameter for the 190 4FGL pulsars, we find that the resulting FP relation, i.e. the gamma-ray luminosity in terms of the other observables, reads $L_γ=10^{14.3\pm 1.3}(ε_{\rm c1}/{\rm MeV})^{1.39\pm0.17}(B_{\star}/{\rm G})^{0.12\pm 0.03}(\dot{\mathcal{E}}/{\rm erg\;s^{-1}})^{0.39\pm 0.05}{\rm ~erg\;s^{-1}}$; this is in good agreement with both the empirical relation reported by Kalapotharakos et al. (2019) and the theoretically predicted relation for curvature radiation. Finally, we revisit the radiation reaction limited condition, to find it is a sufficient but not necessary condition for the theoretical derivation of the FP relation. However, the assumption of the radiation reaction limited acceleration reveals the underlying accelerating electric field component and its scaling with $\dot{\mathcal{E}}$.

preprint2021arXiv

Probing the High-energy $γ$-ray Emission Mechanism in the Vela Pulsar via Phase-resolved Spectral and Energy-dependent Light Curve Modeling

Recent kinetic simulations sparked a debate regarding the emission mechanism responsible for pulsed GeV $γ$-ray emission from pulsars. Some models invoke curvature radiation, while other models assume synchrotron radiation in the current-sheet. We interpret the curved spectrum of the Vela pulsar as seen by H.E.S.S. II (up to $\sim$100 GeV) and the $Fermi$ Large Area Telescope (LAT) to be the result of curvature radiation due to primary particles in the pulsar magnetosphere and current sheet. We present phase-resolved spectra and energy-dependent light curves using an extended slot gap and current sheet model, invoking a step function for the accelerating electric field as motivated by kinetic simulations. We include a refined calculation of the curvature radius of particle trajectories in the lab frame, impacting the particle transport, predicted light curves, and spectra. Our model reproduces the decrease of the flux of the first peak relative to the second one, evolution of the bridge emission, near-constant phase positions of peaks, and narrowing of pulses with increasing energy. We can explain the first of these trends because we find that the curvature radii of the particle trajectories in regions where the second $γ$-ray light curve peak originates are systematically larger than those associated with the first peak, implying that the spectral cutoff of the second peak is correspondingly larger. However, an unknown azimuthal dependence of the $E$-field as well as uncertainty in the precise spatial origin of the GeV emission, precludes a simplistic discrimination of emission mechanisms.

preprint2021arXiv

The Emission Physics of Millisecond Pulsars

Understanding the physics of rotation-powered millisecond pulsars (MSPs) presents a number of challenges compared to that of the non-recycled pulsar population. Even though their fast rotation rates can produce high spin-down power and accelerating electric fields, their relatively low surface magnetic fields make the production of electron-positron pairs required for radio emission difficult. The Fermi Gamma-Ray Space Telescope has discovered pulsed W-rays from a large fraction of the MSP population that have light curves surprisingly similar to those of young pulsars. However, their very compact magnetospheres enable magnetic fields at the light cylinder that are comparable to those of the most energetic pulsars. This fact and recent global magnetosphere models showing that particle acceleration takes place near and beyond the light cylinder, now makes the W-rays from MSPs plausible. The large increase in binary systems harboring MSPs has revitalized the study of shock acceleration and high-energy emission in such systems, with many showing orbitally-modulated X-rays. This Chapter will review the history and our current studies of the mechanisms for multiwavelength emission from MSPs.

preprint2021arXiv

Very-High-Energy Emission From Pulsars

Air-Cherenkov telescopes have detected pulsations at energies above 50 GeV from a growing number of Fermi pulsars. These include the Crab, Vela, PSR B1706-44 and Geminga, with the first two having pulsed detections above 1 TeV. In some cases, there appears to be very-high-energy (VHE) emission that is an extension of the Fermi spectra to high energies, while in other cases, additional higher-energy spectral components that require a separate emission mechanism may be present. We present results of broad-band spectral modeling using global magnetosphere fields and multiple emission mechanisms that include synchro-curvature (SC) and inverse Compton scattered (ICS) radiation from accelerated particles (primaries) and synchrotron-self Compton (SSC) emission from lower-energy pairs. Our models predict three distinct VHE components: SC from primaries whose high-energy tail can extend to 100 GeV, SSC from pairs that can extend to several TeV and ICS from primary particles accelerated in the current sheet, scattering pair synchrotron radiation, that appears beyond 10 TeV. Our models suggest that H.E.S.S.-II and MAGIC have detected the high-energy tail of the primary SC component that produces the Fermi spectrum in Vela, Geminga and PSR B1706-44. We argue that the ICS component peaking above 10 TeV from Vela has been seen by H.E.S.S. Detection of this emission component from the Crab and other pulsars is possible with HAWC and CTA, and directly measures the maximum particle energy in pulsars.

preprint2020arXiv

Simultaneous magnetic polar cap heating during a flaring episode from the magnetar 1RXS J170849.0-400910

During a pointed 2018 NuSTAR observation, we detected a flare with a 2.2 hour duration from the magnetar 1RXS J170849.0$-$400910. The flare, which rose in $\sim25$ seconds to a maximum flux 6 times larger than the persistent emission, is highly pulsed with an rms pulsed fraction of $53\%$. The pulse profile shape consists of two peaks separated by half a rotational cycle, with a peak flux ratio of $\sim$2. The flare spectrum is thermal with an average temperature of 2.1 keV. Phase resolved spectroscopy show that the two peaks possess the same temperature, but differ in size. These observational results along with simple light curve modeling indicate that two identical antipodal spots, likely the magnetic poles, are heated simultaneously at the onset of the flare and for its full duration. Hence, the origin of the flare has to be connected to the global dipolar structure of the magnetar. This might best be achieved externally, via twists to closed magnetospheric dipolar field lines seeding bombardment of polar footpoint locales with energetic pairs. Approximately 1.86 hours following the onset of the flare, a short burst with its own 3-minute thermal tail occurred. The burst tail is also pulsating at the spin period of the source and phase-aligned with the flare profile, implying an intimate connection between the two phenomena. The burst may have been caused by a magnetic reconnection event in the same twisted dipolar field lines anchored to the surface hot spots, with subsequent return currents supplying extra heat to these polar caps.

preprint2020arXiv

The Fast Radio Burst Luminosity Function and Death Line in the Low-Twist Magnetar Model

We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). Motivated by the power-law fluence distributions of FRB 121102, we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol which directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar short bursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts (required in WT19) for individual repeaters implies the predicted burst fluence distribution is narrow, $\lesssim 3$ decades for yielding strains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and charge starvation, we obtain a death line for FRBs which segregates magnetars from the normal pulsar population, suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individual magnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-law luminosity function's low energy character depends on the population model, while the high energy index traces that of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalent energy/luminosity function peaks at $\sim10^{37}-10^{40}$ erg with a low-energy cutoff at $\sim 10^{37}$ erg. Lastly, we consider the local fluence distribution of FRBs, and find that it can constrain the subset of FRB-producing magnetar progenitors. Our model suggests that improvements in sensitivity may reveal flattening of the global FRB fluence distribution and saturation in FRB rates.

preprint2020arXiv

The Multipolar Magnetic Field of Millisecond Pulsar PSR J0030+0451

Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed to constrain the neutron star mass and radius, has inferred surface hot-spots (the magnetic polar caps) that imply significantly non-dipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole plus quadrupole components using static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. (2019) and Riley et al. (2019), 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 to the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous gamma-ray light curves following Kalapotharakos et al. (2014) these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the gamma-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 gamma-ray light curves. We conclude with a discussion of the broader implications of our study.

preprint2020arXiv

VTXO: The Virtual Telescope for X-ray Observations

The Virtual Telescope for X-ray Observations (VTXO) will use lightweight Phase Fresnel Lenses (PFLs) in a virtual X-ray telescope with $\sim$1 km focal length and with $\sim$50 milli-arcsecond angular resolution. VTXO is formed by using precision formation flying of two SmallSats: a smaller OpticsSat that houses the PFLs and navigation beacons while a larger DetectorSat contains an X-ray camera, a precision start tracker, and the propulsion for the formation flying. The baseline flight dynamics uses a highly elliptical supersynchronous orbit allow the formation to hold in an inertial frame around the 90,000 km apogee for 10 hours of the 32.5 hour orbit with nearly a year mission lifetime. VTXO's fine angular resolution enables measuring the environments close to the central engines of bright compact X-ray sources. This X-ray imaging capability allows for the study of the effects of dust scattering near to the central objects such as Cyg X-3 and GX 5-1, for the search for jet structure near to the compact object in X-ray novae such as Cyg X-1 and GRS 1915+105, and for the search for structure in the termination shock of in the Crab pulsar wind nebula. The VTXO SmallSat and instrument designs, mission parameters, and science performance are described. VTXO development was supported as one of the selected 2018 NASA Astrophysics SmallSat Study (AS$^3$) missions.

preprint2019arXiv

Constraining the Emission Geometry and Mass of the White Dwarf Pulsar AR Sco using the Rotating Vector Model

We apply the standard radio pulsar rotating vector model to the white dwarf pulsar AR Sco's optical polarization position angle swings folded at the white dwarf's spin period as obtained by Buckley et al. (2017). Owing to the long duty cycle of spin pulsations with a good signal-to-noise ratio over the entire spin phase, in contrast to neutron star radio pulsars, we find well-constrained values for the magnetic obliquity $α$ and observer viewing direction $ζ$ with respect to the spin axis. We find $\cosα=0.060^{+0.050}_{-0.053}$ and $\cosζ=0.49^{+0.09}_{-0.08}$, implying an orthogonal rotator with an observer angle $ζ={60.4^\circ}^{+5.3^\circ}_{-6.0^\circ}$. This orthogonal nature of the rotator is consistent with the optical light curve consisting of two pulses per spin period, separated by $180^\circ$ in phase. Under the assumption that $ζ\approx i$, where $i$ is the orbital inclination, and a Roche-lobe-filling companion M star, we obtain $m_{\rm WD} = 1.00^{+0.16}_{-0.10} M_\odot$ for the white dwarf mass. These polarization modeling results suggest the that non-thermal emission arises from a dipolar white dwarf magnetosphere and close to the star, with synchrotron radiation (if non-zero pitch angles can be maintained) being the plausible loss mechanism, marking AR Sco as an exceptional system for future theoretical and observational study.