ABSTRACT Observations by the Fermi Large Area Telescope of γ-ray millisecond pulsar (MSP) light curves imply copious pair production in their magnetospheres, and not exclusively in those of younger ...pulsars. Such pair cascades may be a primary source of Galactic electrons and positrons, contributing to the observed enhancement in positron flux above ∼10 GeV. Fermi has also uncovered many new MSPs, impacting Galactic stellar population models. We investigate the contribution of Galactic MSPs to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day MSPs. We simulate their pair spectra invoking an offset-dipole magnetic field. We also consider positrons and electrons that have been further accelerated to energies of several TeV by strong intrabinary shocks in black widow (BW) and redback (RB) systems. Since MSPs are not surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs freely escape and undergo losses only in the intergalactic medium. We compute the transported pair spectra at Earth, following their diffusion and energy loss through the Galaxy. The predicted particle flux increases for non-zero offsets of the magnetic polar caps. Pair cascades from the magnetospheres of MSPs are only modest contributors around a few tens of GeV to the lepton fluxes measured by the Alpha Magnetic Spectrometer, PAMELA, and Fermi, after which this component cuts off. The contribution by BWs and RBs may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.
ABSTRACT Nearly one-third of the sources listed in the Third Fermi Large Area Telescope (LAT) catalog (3FGL) remain unassociated. It is possible that predicted and even unanticipated gamma-ray source ...classes are present in these data waiting to be discovered. Taking advantage of the excellent spectral capabilities achieved by the Fermi LAT, we use machine-learning classifiers (Random Forest and XGBoost) to pinpoint potentially novel source classes in the unassociated 3FGL sample outside the Galactic plane. Here we report a total of 34 high-confidence Galactic candidates at . The currently favored standard astrophysical interpretations for these objects are pulsars or low-luminosity globular clusters hosting millisecond pulsars (MSPs). Yet these objects could also be interpreted as dark matter annihilation taking place in ultra-faint dwarf galaxies or dark matter subhalos. Unfortunately, Fermi LAT spectra are not sufficient to break degeneracies between the different scenarios. Careful visual inspection of archival optical images reveals no obvious evidence for low-luminosity globular clusters or ultra-faint dwarf galaxies inside the 95% error ellipses. If these are pulsars, this would bring the total number of MSPs at to 106, down to an energy flux 4.0 × 10−12 erg cm−2 s−1 between 100 MeV and 100 GeV. We find this number to be in excellent agreement with predictions from a new population synthesis of MSPs that predicts 100-126 high-latitude 3FGL MSPs depending on the choice of high-energy emission model. If, however, these are dark matter substructures, we can place upper limits on the number of Galactic subhalos surviving today and on dark matter annihilation cross sections. These limits are beginning to approach the canonical thermal relic cross section for dark matter particle masses below ∼100 GeV in the bottom quark ( ) annihilation channel.
Context. The radio-quiet γ -ray pulsar PSR J2021+4026 is a peculiar Fermi -LAT pulsar showing repeated and quasi-periodic mode changes. Its γ -ray flux shows repeated variations between two states at ...intervals of ∼3.5 years. These events occur over timescales < 100 days and are correlated with sudden changes in the spin-down rate. Multiwavelength observations also revealed an X-ray phase shift relative to the γ -ray profile for one of the events. PSR J2021+4026 is currently the only known isolated γ -ray pulsar showing significant variability, and thus it has been the object of thorough investigations. Aims. The goal of our work is to study the mode changes of PSR J2021+4026 with improved detail. By accurately characterizing variations in the γ -ray spectrum and pulse profile, we aim to relate the Fermi -LAT observations to theoretical models. We also aim to interpret the mode changes in terms of variations in the structure of a multipolar dissipative magnetosphere. Methods. We continually monitored the rotational evolution and the γ -ray flux of PSR J2021+4026 using more than 13 years of Fermi -LAT data with a binned likelihood approach. We investigated the features of the phase-resolved spectrum and pulse profile, and from these we inferred the macroscopic conductivity, the electric field parallel to the magnetic field, and the curvature radiation cutoff energy. These physical quantities are related to the spin-down rate and the γ -ray flux and therefore are relevant to the theoretical interpretation of the mode changes. We introduced a simple magnetosphere model that combines a dipole field with a strong quadrupole component. We simulated magnetic field configurations to determine the positions of the polar caps for different sets of parameters. Results. We clearly detect the previous mode changes and confirm a more recent mode change that occurred around June 2020. We provide a full set of best-fit parameters for the phase-resolved γ -ray spectrum and the pulse profile obtained in five distinct time intervals. We computed the relative variations in the best-fit parameters, finding typical flux changes between 13% and 20%. Correlations appear between the γ -ray flux and the spectral parameters, as the peak of the spectrum shifts by ∼10% toward lower energies when the flux decreases. The analysis of the pulse profile reveals that the pulsed fraction of the light curve is larger when the flux is low. Finally, the magnetosphere simulations show that some configurations could explain the observed multiwavelength variability. However, self-consistent models are required to reproduce the observed magnitudes of the mode changes.
In more than four years of observation the Large Area Telescope on board the Fermi satellite has identified pulsed γ-ray emission from more than 80 young or middle-aged pulsars, in most cases ...providing light curves with high statistics. Fitting the observed profiles with geometrical models can provide estimates of the magnetic obliquity α and of the line of sight angle ζ, yielding estimates of the radiation beaming factor and radiated luminosity. Using different γ-ray emission geometries (Polar Cap, Slot Gap, Outer Gap, One Pole Caustic) and core plus cone geometries for the radio emission, we fit γ-ray light curves for 76 young or middle-aged pulsars and we jointly fit their γ-ray plus radio light curves when possible. We find that a joint radio plus γ-ray fit strategy is important to obtain (α,ζ) estimates that can explain simultaneously detectable radio and γ-ray emission: when the radio emission is available, the inclusion of the radio light curve in the fit leads to important changes in the (α,ζ) solutions. The most pronounced changes are observed for Outer Gap and One Pole Caustic models for which the γ-ray only fit leads to underestimated α or ζ when the solution is found to the left or to the right of the main α-ζ plane diagonal respectively. The intermediate-to-high altitude magnetosphere models, Slot Gap, Outer Gap, and One pole Caustic, are favoured in explaining the observations. We find no apparent evolution of α on a time scale of 106 years. For all emission geometries our derived γ-ray beaming factors are generally less than one and do not significantly evolve with the spin-down power. A more pronounced beaming factor vs. spin-down power correlation is observed for Slot Gap model and radio-quiet pulsars and for the Outer Gap model and radio-loud pulsars. The beaming factor distributions exhibit a large dispersion that is less pronounced for the Slot Gap case and that decreases from radio-quiet to radio-loud solutions. For all models, the correlation between γ-ray luminosity and spin-down power is consistent with a square root dependence. The γ-ray luminosities obtained by using the beaming factors estimated in the framework of each model do not exceed the spin-down power. This suggests that assuming a beaming factor of one for all objects, as done in other studies, likely overestimates the real values. The data show a relation between the pulsar spectral characteristics and the width of the accelerator gap. The relation obtained in the case of the Slot Gap model is consistent with the theoretical prediction.
Energy-transport effects can alter the structure that develops as a supernova evolves into a supernova remnant. The Rayleigh-Taylor instability is thought to produce structure at the interface ...between the stellar ejecta and the circumstellar matter, based on simple models and hydrodynamic simulations. Here we report experimental results from the National Ignition Facility to explore how large energy fluxes, which are present in supernovae, affect this structure. We observed a reduction in Rayleigh-Taylor growth. In analyzing the comparison with supernova SN1993J, a Type II supernova, we found that the energy fluxes produced by heat conduction appear to be larger than the radiative energy fluxes, and large enough to have dramatic consequences. No reported astrophysical simulations have included radiation and heat conduction self-consistently in modeling supernova remnants and these dynamics should be noted in the understanding of young supernova remnants.
ABSTRACT We study the efficiency of pair production in polar caps of young pulsars under a variety of conditions to estimate the maximum possible multiplicity of pair plasma in pulsar magnetospheres. ...We develop a semi-analytic model for calculation of cascade multiplicity which allows efficient exploration of the parameter space and corroborate it with direct numerical simulations. Pair creation processes are considered separately from particle acceleration in order to assess different factors affecting cascade efficiency, with acceleration of primary particles described by recent self-consistent non-stationary model of pair cascades. We argue that the most efficient cascades operate in the curvature radiation/synchrotron regime, the maximum multiplicity of pair plasma in pulsar magnetospheres is ∼few × 105. The multiplicity of pair plasma in magnetospheres of young energetic pulsars weakly depends on the strength of the magnetic field and the radius of curvature of magnetic field lines and has a stronger dependence on pulsar inclination angle. This result questions assumptions about very high pair plasma multiplicity in theories of pulsar wind nebulae.
Thanks to the huge amount of γ-ray pulsar photons collected by the Fermi Large Area Telescope since its launch in June 2008, it is now possible to constrain γ-ray geometrical models by comparing ...simulated and observed light-curve morphological characteristics. We assumed vacuum-retarded dipole (VRD) pulsar magnetic field and tested simulated and observed morphological light-curve characteristics in the framework of two pole emission geometries, Polar Cap (PC) and Slot Gap (SG), and one pole emission geometries, traditional Outer Gap (OG) and One Pole Caustic (OPC). Radio core plus cone emission was assumed for the pulsars of the simulated sample. We compared simulated and observed recurrence of class shapes and peak multiplicity, peak separation, radio-lag distributions, and trends of peak separation and radio lag as a function of observable and non-observable pulsar parameters. We studied how pulsar morphological characteristics change in multi-dimensional observable and non-observable pulsar parameter space. The PC model gives the poorest description of the LAT pulsar light-curve morphology. The OPC best explains both the observed γ-ray peak multiplicity and shape classes. The OPC and SG models describe the observed γ-ray peak-separation distribution for low- and high-peak separations, respectively. This suggests that the OPC geometry best explains the single-peak structure but does not manage to describe the widely separated peaks predicted in the framework of the SG model as the emission from the two magnetic hemispheres. The OPC radio-lag distribution shows higher agreement with observations suggesting that assuming polar radio emission, the γ-ray emission regions are likely to be located in the outer magnetosphere. Alternatively, the radio emission altitude could be higher that we assumed. We compared simulated non-observable parameters with the same parameters estimated for LAT pulsars in the framework of the same models. The larger agreement between simulated and LAT estimations in the framework of the OPC suggests that the OPC model best predicts the observed variety of profile shapes. The larger agreement obtained between observations and the OPC model predictions jointly with the need to explain the abundant 0.5 separated peaks with two-pole emission geometries, calls for thin OPC gaps to explain the single-peak geometry but highlights the need of two-pole caustic emission geometry to explain widely separated peaks.
Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio ...beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase, we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II), or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of the magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.
Here, we present precise phase-connected pulse timing solutions for 16 γ-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one ...very faint radio pulsar (PSR J1124–5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multiwavelength follow-up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124–5916, and J1813–1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard power-law component. Another Chandra ACIS image of PSR J1732–3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the γ-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124–5916 to establish the γ-ray to radio phase offset.
Sixteen pulsars have been discovered so far in blind searches of photons collected with the Large Area Telescope on the Fermi Gamma-ray Space Telescope. We here report the discovery of radio ...pulsations from two of them. PSR J1741-2054, with period P = 413 ms, was detected in archival Parkes telescope data and subsequently has been detected at the Green Bank Telescope (GBT). Its received flux varies greatly due to interstellar scintillation and it has a very small dispersion measure of DM = 4.7 pc cm-3, implying a distance of 0.4 kpc and possibly the smallest luminosity of any known radio pulsar. At this distance, for isotropic emission, its gamma-ray luminosity above 0.1 GeV corresponds to 28% of the spin-down luminosity of erg s-1. The gamma-ray profile occupies 1/3 of pulse phase and has three closely spaced peaks with the first peak lagging the radio pulse by Delta *d = 0.29 P. We have also identified a soft Swift source that is the likely X-ray counterpart. In many respects PSR J1741-2054 resembles the Geminga pulsar. The second source, PSR J2032+4127, was detected at the GBT. It has P = 143 ms, and its DM = 115 pc cm-3 suggests a distance of 3.6 kpc, but we consider it likely that it is located within the Cyg OB2 stellar association at half that distance. The radio emission is nearly 100% linearly polarized, and the main radio peak precedes by Delta *d = 0.15 P the first of two narrow gamma-ray peaks that are separated by Delta *D = 0.50 P. The second peak has a harder spectrum than the first one, following a trend observed in young gamma-ray pulsars. Faint, diffuse X-ray emission in a Chandra image is possibly its pulsar wind nebula. The wind of PSR J2032+4127 is responsible for the formerly unidentified HEGRA source TeV J2032+4130. PSR J2032+4127 is coincident in projection with MT91 213, a Be star in Cyg OB2, although apparently not a binary companion of it.