BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alice K. Harding

Alice K. Harding appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-phastro-ph.SRastro-ph.GAastro-ph.IMhep-ph

Catalog footprint

What is connected

35works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}}$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Aid of Optical Studies in Understanding Millisecond Pulsar Binaries

A large number of new "black widow" and "redback" energetic millisecond pulsars with irradiated stellar companions have been discovered through radio and optical searches of unidentified \textit{Fermi} sources. Synchrotron emission, from particles accelerated up to several TeV in the intrabinary shock, exhibits modulation at the binary orbital period. Our simulated double-peaked X-ray light curves modulated at the orbital period, produced by relativistic Doppler-boosting along the intrabinary shock, are found to qualitatively match those observed in many sources. In this model, redbacks and transitional pulsar systems where the double-peaked X-ray light curve is observed at inferior conjunction have intrinsically different shock geometry than other millisecond pulsar binaries where the light curve is centered at superior conjunction. We discuss, and advocate, how current and future optical observations may aid in constraining the emission geometry, intrabinary shock and the unknown physics of pulsar winds.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gamma-ray pulsars: a gold mine

The most energetic neutron stars, powered by their rotation, are capable of producing pulsed radiation from the radio up to gamma rays with nearly TeV energies. These pulsars are part of the universe of energetic and powerful particle accelerators, using their uniquely fast rotation and formidable magnetic fields to accelerate particles to ultra-relativistic speed. The extreme properties of these stars provide an excellent testing ground, beyond Earth experience, for nuclear, gravitational, and quantum-electrodynamical physics. A wealth of gamma-ray pulsars has recently been discovered with the Fermi Gamma-Ray Space Telescope. The energetic gamma rays enable us to probe the magnetospheres of neutron stars and particle acceleration in this exotic environment. We review the latest developments in this field, beginning with a brief overview of the properties and mysteries of rotation-powered pulsars, and then discussing gamma-ray observations and magnetospheric models in more detail.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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_+ = 3 \times 10^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^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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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. (2014), 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, is assumed to be negligible, while outside the light-cylinder it rescales with a finite conductivity (σ). 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 γ-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 σ values that best describe each of the pulsars in this group show an increase with the spin-down rate $(\dot{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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Discovery of PSR J1227-4853: A transition from a low-mass X-ray binary to a redback millisecond pulsar

XSS J12270-4859 is an X-ray binary associated with the Fermi LAT gamma-ray source 1FGL J1227.9-4852. In 2012 December, this source underwent a transition where the X-ray and optical luminosity dropped and the spectral signatures of an accretion disc disappeared. We report the discovery of a 1.69 millisecond pulsar (MSP), PSR J1227-4853, at a dispersion measure of 43.4 pc cm$^{-3}$ associated with this source, using the GMRT at 607 MHz. This demonstrates that, post-transition, the system hosts an active radio MSP. This is the third system after PSR J1023+0038 and PSR J1824-2452I showing evidence of state switching between radio MSP and low-mass X-ray binary (LMXB) states. We report timing observations of PSR J1227-4853 with the GMRT and Parkes, which give a precise determination of the rotational and orbital parameters of the system. The companion mass measurement of 0.17 to 0.46 M$_{sun}$ suggests that this is a redback system. PSR J1227-4853 is eclipsed for about 40% of its orbit at 607 MHz; with additional short-duration eclipses at all orbital phases. We also find that the pulsar is very energetic, with a spin-down luminosity of ~ 10$^{35}$ erg s$^{-1}$. We report simultaneous imaging and timing observations with the GMRT, which suggests that eclipses are caused by absorption, rather than dispersion smearing or scattering.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

We compute the patterns of $γ$-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 $γ$-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 $γ$-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 $γ$-ray emission on the $γ$-ray peak-separation. We perform a statistical comparison of our model radio-lag vs 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 $γ$-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, these models produce GeV photon cut-off energies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-Time-Resolution Measurements of the Polarization of the Crab Pulsar at 1.38 GHz

Using the Westerbork Synthesis Radio Telescope (WSRT), we obtained high-time-resolution measurements of the full (linear and circular) polarization of the Crab pulsar. Taken at a resolution of 1/8192 of the 34-ms pulse period (i.e., $4.1~μ{\rm s}$), the 1.38-GHz linear-polarization measurements are in general agreement with previous lower-time-resolution 1.4-GHz measurements of linear polarization in the main pulse (MP), in the interpulse (IP), and in the low-frequency component (LFC). We find the MP and IP to be linearly polarized at about $24\%$ and $21\%$, with no discernible difference in polarization position angle. However, and contrary to theoretical expectations and measurements in the visible, we find no evidence for significant variation (sweep) in polarization position angle over the MP, the IP, or the LFC. Although, the main pulse exhibits a small but statistically significant quadratic variation in the degree of linear polarization. We discuss the implications which appear to be in contradiction to theoretical expectations. In addition, we detect weak circular polarization in the main pulse and interpulse, and strong ($\approx 20\%$) circular polarization in the low-frequency component, which also exhibits very strong ($\approx 98\%$) linear polarization at a position angle about $40\degree$ from that of the MP or IP. The pulse-mean polarization properties are consistent with the LFC being a low-altitude component and the MP and IP being high-altitude caustic components. Nevertheless, current models for the MP and IP emission do not readily account for the observed absence of pronounced polarization changes across the pulse. Finally, we measure IP and LFC pulse phases relative to the MP that are consistent with recent measurements, which have shown that the phases of these pulse components are evolving with time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gamma-Ray Light Curves from Pulsar Magnetospheres with Finite Conductivity

We investigate the shapes of γ-ray pulsar light curves using 3D pulsar magnetosphere models of finite conductivity. These models, covering the entire spectrum of solutions between vacuum and force-free magnetospheres, for the first time afford mapping the GeV emission of more realistic, dissipative pulsar magnetospheres. To this end we generate model light curves following two different approaches: (a) We employ the emission patterns of the slot and outer gap models in the field geometries of magnetospheres with different conductivity σ. (b) We define realistic trajectories of radiating particles in magnetospheres of different σand compute their Lorentz factor under the influence of magnetospheric electric fields and curvature radiation-reaction; with these at hand we then calculate the emitted radiation intensity. The light curves resulting from these prescriptions are quite sensitive to the value of σ, especially in the second approach. While still not self-consistent, these results are a step forward in understanding the physics of pulsar γ-radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Population Synthesis of Normal Radio and Gamma-ray Pulsars Using Markov Chain Monte Carlo Techniques

We present preliminary results of a pulsar population synthesis of normal pulsars from the Galactic disk using a Markov Chain Monte Carlo method to better understand the parameter space of the assumed model. We use the Kuiper test, similar to the Kolmogorov-Smirnov test, to compare the cumulative distributions of chosen observables of detected radio pulsars with those simulated for various parameters. Our code simulates pulsars at birth using Monte Carlo techniques and evolves them to the present assuming initial spatial, kick velocity, magnetic field, and period distributions. Pulsars are spun down to the present, given radio and gamma-ray emission characteristics, filtered through ten selected radio surveys, and a {\it Fermi} all-sky threshold map. Each chain begins with a different random seed and searches a ten-dimensional parameter space for regions of high probability for a total of one thousand different simulations before ending. The code investigates both the "large world" as well as the "small world" of the parameter space. We apply the K-means clustering algorithm to verify if the chains reveal a single or multiple regions of significance. The outcome of the combined set of chains is the weighted average and deviation of each of the ten parameters describing the model. While the model reproduces reasonably well the detected distributions of normal radio pulsars, it does not replicate the predicted detected $\dot P - P$ distribution of {\it Fermi} pulsars. The simulations do not produce sufficient numbers of young, high-$\dot E$ pulsars in the Galactic plane.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling and Maximum Likelihood Fitting of Gamma-Ray and Radio Light Curves of Millisecond Pulsars Detected with Fermi

Pulsed gamma rays have been detected with the Fermi Large Area Telescope (LAT) from more than 20 millisecond pulsars (MSPs), some of which were discovered in radio observations of bright, unassociated LAT sources. We have fit the radio and gamma-ray light curves of 19 LAT-detected MSPs in the context of geometric, outer-magnetospheric emission models assuming the retarded vacuum dipole magnetic field using a Markov chain Monte Carlo maximum likelihood technique. We find that, in many cases, the models are able to reproduce the observed light curves well and provide constraints on the viewing geometries that are in agreement with those from radio polarization measurements. Additionally, for some MSPs we constrain the altitudes of both the gamma-ray and radio emission regions. The best-fit magnetic inclination angles are found to cover a broader range than those of non-recycled gamma-ray pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Population Synthesis of Radio and Gamma-ray Pulsars in the Fermi Era

We present results of our pulsar population synthesis of normal pulsars from the Galactic disk using our previously developed computer code. On the same footing, we use slot gap and outer gap models for gamma-ray emission from normal pulsars to obtain statistics of radio-loud and radio-quiet gamma-ray pulsars. From recently improved understanding of HII and star forming regions in the Galaxy, we develop a new surface density model of the birth location of neutron stars. We explore models of neutron star evolution with magnetic field-decay, and with different initial period and magnetic field distributions. We present preliminary results including simulated population statistics that are compared with recent detections by Fermi of normal, isolated pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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_{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_{pair} is proportional to L_{sd}, with L_{pair} ~ 10^{-3} L_{sd} for normal pulsars and L_{pair} ~ 10^{-2} L_{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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Photospheric radius expansion during magnetar bursts

On August 24th 2008 the new magnetar SGR 0501+4516 (discovered by SWIFT) emitted a bright burst with a pronounced double-peak structure in hard X-rays, reminiscent of the double-peak temporal structure seen in some bright thermonuclear bursts on accreting neutron stars. In the latter case this is due to Photospheric Radius Expansion (PRE): when the flux reaches the Eddington limit, the photosphere expands and cools so that emission becomes softer and drops temporarily out of the X-ray band, re-appearing as the photosphere settles back down. We consider the factors necessary to generate double-peaked PRE events, and show that such a mechanism could plausibly operate in magnetar bursts, despite the vastly different emission process. Identification of the magnetic Eddington limit in a magnetar would constrain magnetic field and distance and could, in principle, enable a measurement of gravitational redshift. It would also locate the emitting region at the neutron star surface, constraining the burst trigger mechanism. Conclusive confirmation of PRE events will require more detailed radiative models for bursts. However for SGR 0501+4516 the predicted critical flux (using the magnetic field strength inferred from timing and the distance suggested by its probable location in the Perseus arm of our Galaxy) is consistent with that observed in the August 24th burst.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pulsar Pair Cascades in a Distorted Magnetic Dipole Field

We investigate the effect of a distorted neutron star dipole magnetic field on pulsar pair cascade multiplicity and pair death lines. Using a simple model for a distorted dipole field that produces an offset polar cap, we derive the accelerating electric field above the polar cap in space charge limited flow. We find that even a modest azimuthally asymmetric distortion can significantly increase the accelerating electric field on one side of the polar cap and, combined with a smaller field line radius of curvature, leads to larger pair multiplicity. The death line for producing pairs by curvature radiation moves downward in the P-Pdot diagram, allowing high pair multiplicities in a larger percentage of the radio pulsar population. These results could have important implications for the radio pulsar population, high energy pulsed emission and the pulsar contribution to cosmic ray positrons.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Population Synthesis of Radio Pulsars in the Fermi Era

We present results of our pulsar population synthesis of normal pulsars from the Galactic disk using our previously developed Monte-Carlo code. From our studies of observed radio pulsars that have clearly identifiable core and cone components, in which we fit the polarization sweep as well as the pulse profiles in order to constrain the viewing geometry, we develop a model describing the luminosity and ratio of radio core-to-cone peak fluxes. In this model, short period pulsars are more cone-dominated. We explore models of neutron star evolution with and without magnetic field decay, and with different initial period distributions. We present preliminary results including simulated population statistics that are compared with the observed radio pulsar population. The evolved neutron star populations resulting from this simulation can be used to model distributions of gamma-ray pulsars for comparison to Fermi results.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Pair-Starved Pulsar Magnetospheres

We propose a simple analytic model for the innermost (within the light cylinder of canonical radius) structure of open-magnetic-field lines of a rotating neutron star (NS) with relativistic outflow of charged particles (electrons/positrons) and arbitrary angle between the NS spin and magnetic axes. We present the self-consistent solution of Maxwell's equations for the magnetic field and electric current in the pair-starved regime where the density of electron-positron plasma generated above the pulsar polar cap is not sufficient to completely screen the accelerating electric field and thus establish the E \cdot B = 0 condition above the pair-formation front up to very high altitudes within the light cylinder. The proposed model may provide a theoretical framework for developing a refined model of the global pair-starved pulsar magnetosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

1-100 GeV Emission from Millisecond Pulsars

A number of rotation-powered millisecond pulsars are powerful sources of X-ray emission. We present predictions for the spectral characteristics of these sources at gamma-ray energies, using a model for acceleration and pair cascades on open magnetic field lines above the polar caps. Since these pulsars have low surface magnetic fields, the majority do not produce sufficient pairs to completely screen the accelerating electric field allowing particle acceleration to high altitude. The resulting emission above 1 GeV comes from curvature radiation by primary electrons with radiation-reaction-limited Lorentz factors. The spectra are very hard power-laws with exponential cutoffs between 1 and 50 GeV, and the spectral power peaks near the cutoff energy. Millisecond pulsars are thus ideal targets for air-Cherenkov detectors that are able to reach energy thresholds below 50 GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Emission From Rotation-Powered Pulsars and Magnetars

I will review the latest developments in understanding the high-energy emission of rotation-powered pulsars and magnetically-powered Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs). These fields have been extremely active in the last few years, both observationally and theoretically, driven partly by new X-ray data from Chandra, XMM-Newton and RXTE. At the same time, the Parkes Multibeam Survey has discovered over 700 new radio pulsars, some of them young and coincident with EGRET sources, and others having magnetar-strength magnetic fields. These new observations are raising important questions about neutron star birth and evolution, as well as the properties of their high-energy emission.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

High-Energy Emission From Millisecond Pulsars

The X-ray and gamma-ray spectrum of rotation-powered millisecond pulsars is investigated in a model for acceleration and pair cascades on open field lines above the polar caps. Although these pulsars have low surface magnetic fields, their short periods allow them to have large magnetospheric potential drops, but the majority do not produce sufficient pairs to completely screen the accelerating electric field. The accelerating particles maintain high Lorentz factors and undergo cyclotron resonant absorption of radio emission, that produces and maintains a large pitch angle, resulting in a strong synchrotron component. The resulting spectra consist of several distinct components: curvature radiation from primary electrons dominating from 1 - 100 GeV, synchrotron radiation from primary and secondary electrons dominating up to about 100 MeV, and much weaker inverse-Compton radiation from primary electrons at 0.1 - 1 TeV. We find that the relative size of these components depends on pulsar period, period derivative, and neutron star mass and radius with the level of the synchrotron component also depending sensitively on the radio emission properties. This model is successful in describing the observed X-ray and gamma-ray spectrum of PSR J0218+4232 as synchrotron radiation, peaking around 100 MeV and extending up to a turnover around several GeV. The predicted curvature radiation components from a number of millisecond pulsars, as well as the collective emission from the millisecond pulsars in globular clusters, should be detectable with AGILE and GLAST. We also discuss a hidden population of X-ray-quiet and radio-quiet millisecond pulsars which have evolved below the pair death line, some of which may be detectable by telescopes sensitive above 1 GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Regimes of Pulsar Pair Formation and Particle Energetics

We investigate the conditions required for the production of electron-positron pairs above a pulsar polar cap (PC) and the influence of pair production on the energetics of the primary particle acceleration. Assuming space-charge limited flow acceleration including the inertial frame-dragging effect, we allow both one-photon and two-photon pair production by either curvature radiation (CR) photons or photons resulting from inverse-Compton scattering of thermal photons from the PC by primary electrons. We find that, while only the younger pulsars can produce pairs through CR, nearly all known radio pulsars are capable of producing pairs through non-resonant inverse-Compton scatterings. The effect of the neutron star equations of state on the pair death lines is explored. We show that pair production is facilitated in more compact stars and more massive stars. Therefore accretion of mass by pulsars in binary systems may allow pair production in most of the millisecond pulsar population. We also find that two-photon pair production may be important in millisecond pulsars if their surface temperatures are above about three million degrees K. Pulsars that produce pairs through CR will have their primary acceleration limited by the effect of screening of the electric field. In this regime, the high-energy luminosity should follow a $L_{HE} \propto \dot E_{rot}^{1/2}$ dependence. The acceleration voltage drop in pulsars that produce pairs only through inverse-Compton emission will not be limited by electric field screening. In this regime, the high-energy luminosity should follow a $L_{HE} \propto \dot E_{rot}$ dependence. Thus, older pulsars will have significantly lower gamma-ray luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

The Attenuation of Gamma-Ray Emission in Strongly-Magnetized Pulsars

Gamma-rays from pulsars can be efficiently attenuated in their magnetospheres via the mechanism of single-photon pair production and also the exotic QED process of photon splitting, which become prolific in fields approaching the quantum critical value of $B_{cr}=4.41\times 10^{13}$ Gauss. Recently we have published results of our modelling of strongly-magnetized $γ$-ray pulsars, which focused on the escape or attenuation of photons emitted near the pole at the neutron star surface in dipole fields, in a Schwarzschild metric. We found that pair production and splitting totally inhibit emission above around 10--30 MeV in PSR1509-58, whose surface field is inferred to be as high as $0.7B_{cr}$. Our model pulsar spectra are consistent with the EGRET upper limits for PSR1509-58 for a wide range of polar cap sizes. Here we review the principal predictions of our attenuation analysis, and identify how its powerful observational diagnostic capabilities relate to current and future gamma-ray experiments. Diagnostics include the energy of the gamma-ray turnover and the spectral polarization, which constrain the estimated polar cap size and field strength, and can determine the relative strength of splitting and pair creation.

0 citations0 reviewsNo code linkedOpen access