TeV emission from BL Lacertae (BL) objects is commonly modeled as Synchrotron-Self Compton (SSC) radiation from relativistically moving homogeneous plasma blobs. In the context of these models, the blob Lorentz factors needed to reproduce the corrected for absorption by the diffuse IR background (DIRB) TeV emission are large ($\delta \gtrsim 50$). The main reason for this is that stronger beaming eases the problem of the lack of $\sim$ IR-UV synchrotron seed photons needed to produce the de-absorbed $\sim $ few TeV peak of the spectral energy distribution (SED). However, such high Doppler factors are in strong disagreement with the unified scheme, according to which BLs are FR I radio galaxies with their jets closely aligned to the line of sight. Here, motivated by the detection of sub-luminal velocities in the sub-pc scale jets of the best studied TeV blazars, MKN 421 and MKN 501, we examine the possibility that the relativistic flows in the TeV BLs decelerate. In this case, the problem of the missing seed photons is solved because of Upstream Compton (UC) scattering, a process in which the upstream energetic electrons from the fast base of the flow `see' the synchrotron seed photons produced in the slow part of the flow relativistically beamed. Modest Lorentz factors ($\Gamma \sim 15$), decelerating down to values compatible with the recent radio interferometric observations, reproduce the $\sim $ few TeV peak energy of these sources. Furthermore, such decelerating flows are shown to be in agreement with the BL - FR I unification, naturally reproducing the observed BL/FR I broad band luminosity ratios.Comment: Accepted in ApJ

We present three-dimensional (3D) global kinetic pulsar magnetosphere models, where the charged particle trajectories and the corresponding electromagnetic fields are treated self-consistently. For our study, we have developed a Cartesian 3D relativistic particle-in-cell code that incorporates radiation reaction forces. We describe our code and discuss the related technical issues, treatments, and assumptions. Injecting particles up to large distances in the magnetosphere, we apply arbitrarily low to high particle injection rates, and obtain an entire spectrum of solutions from close to the vacuum-retarded dipole to close to the force-free (FF) solution, respectively. For high particle injection rates (close to FF solutions), significant accelerating electric field components are confined only near the equatorial current sheet outside the light cylinder. A judicious interpretation of our models allows the particle emission to be calculated, and consequently, the corresponding realistic high-energy sky maps and spectra to be derived. Using model parameters that cover the entire range of spin-down powers of Fermi young and millisecond pulsars, we compare the corresponding model γ-ray light curves, cutoff energies, and total γ-ray luminosities with those observed by Fermi to discover a dependence of the particle injection rate, , on the spin-down power, , indicating an increase of with . Our models, guided by Fermi observations, provide field structures and particle distributions that are not only consistent with each other but also able to reproduce a broad range of the observed γ-ray phenomenologies of both young and millisecond pulsars.

We discuss a simple non-supersymmetric model based on the electroweak gauge group SU(2) L × SU(2) × U(1) B-L where the lightest of the right-handed neutrinos, which are part of the leptonic doublet of SU(2) , play the role of a long-lived unstable dark matter with mass in the multi-PeV range. We use a resonant s-channel annihilation to obtain the correct thermal relic density and relax the unitarity bound on dark matter mass. In this model, there exists a 3-body dark matter decay mode producing tau leptons and neutrinos, which could be the source for the PeV cascade events observed in the IceCube experiment. The model can be tested with more precise flavor information of the highest-energy neutrino events in future data.

We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh & Timokhin (WT19). 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 that directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates that 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, ≲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 ∼10³⁷–10⁴⁰ erg with a low-energy cutoff at ∼10³⁷ 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 a flattening of the global FRB fluence distribution and saturation in FRB rates.

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