In recent years, surprise discoveries of pulsed emission from the Crab and Vela pulsars above 100 GeV have drawn renewed attention to this largely unexplored region of the energy range. In this paper, we discuss example light curves due to curvature emission, with good resolution in the different energy bands. Continued light curve modelling may help to discriminate between different emission mechanisms, as well as constrain the location where emission is produced within the pulsar magnetosphere, including regions beyond the light cylinder.

RX J1713.7-3946 is a prototype in the γ-ray-bright supernova remnants (SNRs) and is in continuing debates on its hadronic versus leptonic origin of the γ-ray emission. We explore the role played by the diffusive relativistic protons that escape from the SNR shock wave in the γ-ray emission, apart from the emission of high energy particles from the inside of the SNR. In the scenario that the SNR shock propagates in a clumpy molecular cavity, we consider that the γ-ray emission from the inside of the SNR may either arise from the IC scattering or from the interaction between the trapped energetic protons and the shocked clumps. The dominant origin between them depends on the electron-to-proton number ratio. The surrounding molecular cavity wall is considered to also produce γ-ray emission due to the “illumination” by the diffusive protons that escaped from the shock wave during the expansion history. The broad-band spectrum can be well explained by this two-zone model, in which the γ-ray emission from the inside governs the TeV band, while the outer emission component substantially contributes to the GeV γ-rays. The two-zone model can also explain the TeV γ-ray radial brightness profile that significantly stretches beyond the nonthermal X-ray emitting region.

Recent very high energy observations of the galactic centre region performed by H.E.S.S. revealed the presence of a powerful PeVatron. This is the first of such objects detected, and its most plausible counterpart seems to be associated to Sgr A*, the supermassive black hole in the centre of our galaxy. The implications of this discovery will be discussed, in particular in the context of the problem of the origin of galactic cosmic rays.

In 2012 November MAGIC detected a bright flare from IC 310. The flare consisted of two sharp peaks with a typical duration of ~ 5 min. The energy released during that event has been estimated to be at the level of 2 × 1044 erg s−1.

In this work we derive an upper limit on the possible luminosity of flares generated in black hole (BH) magnetosphere, which depends very weakly on the mass of BH and is determined by disk magnetisation, viewing angle, and pair multiplicity. Since all these parameters are smaller than a unit, the luminosity 2 × 1043 erg s−1 can be taken as a strict upper limit for flare luminosity for several minutes variability time. This upper limit appears to be approximately an order of magnitude below the value measured with MAGIC. Thus, we conclude that it seems very unfeasible that the magnetospheric processes can be indeed behind the bright flaring activity recorded from IC 310.

It has been indicated that low-luminosity active galactic nuclei (LLAGNs) are accelerating high-energy cosmic-ray (CR) protons in their radiatively inefficient accretion flows (RIAFs). If this is the case, Sagittarius A* (Sgr A*) should also be generating CR protons, because Sgr A* is a LLAGN. Based on this scenario, we calculate a production rate of CR protons in Sgr A* and their diffusion in the central molecular zone (CMZ) around Sgr A*. The CR protons diffusing in the CMZ create gamma-rays through pp interaction. We show that the gamma-ray luminosity and spectrum are consistent with observations if Sgr A* was active in the past.

The origin of the IceCube astrophysical neutrinos is an outstanding question. Star-forming activities which can accelerate particles to very high energies have been suggested as possible origin of these neutrinos. I will present a scenario where a subset of the neutrino events originate from the Galactic center region and Fermi Bubbles, resulting from star-forming activities. Multi-messenger signal in high energy gamma rays and neutrinos can probe this scenario. I will also present an analysis of the statistical association of the star-forming sources in our Galaxy and outside, with astrophysical neutrinos, as well as expected neutrino signal from these sources by fitting gamma-ray data.

The Milky Way Galaxy glows with the soft gamma ray emission resulting from the annihilation of ~5 × 1043 electron-positron pairs every second. The origin of this vast quantity of antimatter and the peculiar morphology of the 511keV gamma ray line resulting from this annihilation have been the subject of debate for almost half a century. Most obvious positron sources are associated with star forming regions and cannot explain the rate of positron annihilation in the Galactic bulge, which last saw star formation some 10 Gyr ago, or else violate stringent constraints on the positron injection energy. Radioactive decay of elements formed in core collapse supernovae (CCSNe) and normal Type Ia supernovae (SNe Ia) could supply positrons matching the injection energy constraints but the distribution of such potential sources does not replicate the required morphology. We show that a single class of peculiar thermonuclear supernova - SN1991bg-like supernovae (SNe 91bg) - can supply the number and distribution of positrons we see annihilating in the Galaxy through the decay of 44Ti synthesised in these events. Such 44Ti production simultaneously addresses the observed abundance of 44Ca, the 44Ti decay product, in solar system material.

Massive stars at low metallicity are strong candidates for two of the most energetic explosions in the Universe: long duration gamma-ray bursts and superluminous supernovae. But what is the reason these explosions prefer low metallicity environments? To answer this question, we investigate how massive stellar evolution proceeds in low metallicity environments.

The high-energy emission from blazars and other relativistic jet sources indicates that electrons are accelerated to ultra-relativistic (GeV - TeV) energies in these systems. This paper summarizes recent results from numerical studies of two fundamentally different particle acceleration mechanisms potentially at work in relativistic jets: Magnetic-field generation and relativistic particle acceleration in relativistic shear layers, which are likely to be present in relativistic jets, is studied via Particle-in-Cell (PIC) simulations. Diffusive shock acceleration at relativistic shocks is investigated using Monte-Carlo simulations. The resulting magnetic-field configurations and thermal + non-thermal particle distributions are then used to predict multi-wavelength radiative (synchrotron + Compton) signatures of both acceleration scenarios. In particular, we address how anisotropic shear-layer acceleration may be able to circumvent the well-known Lorentz-factor crisis, and how the self-consistent evaluation of thermal + non-thermal particle populations in diffusive shock acceleration simulations provides tests of the bulk Comptonization model for the Big Blue Bump observed in the SEDs of several blazars.

We present an inhomogeneous synchrotron self-Compton (SSC) model for active parts of jets in low mass microquasars. We compare results of our model computed, for typical parameters of microquasars, with the sensitivity of Fermi-LAT.

We consider the shock structures that can arise in blazar jets as a consequence of variations in the jet flow velocity. There are two possible cases: (1) a double shock system consisting of both a forward and reverse shock, and (2) a single shock (either forward or reverse) together with a rarefaction wave. These possibilities depend upon the relative velocity of the two different sections of jet. Using previously calculated spherical models for estimates of the magnetic field and electron number density of the emission region in the TeV blazar Mkn 501, we show that this region is in the form of a thin disk in the plasma rest frame. It is possible to reconcile spectral and pair opacity constraints for Mkn 501 for Doppler factors in the range of 10–20. This is easiest if the corrections for TeV absorption by the infrared background are not as large as implied by recent models.

Pulsars may be born with a short rotation period of milliseconds with the magnetic field amplified through dynamo processes up to ∼1015–1016 G. Such millisecond magnetars spin down rapidly, emitting bursts of high-energy neutrinos and gamma rays. Specifically, acceleration of ions in both the polar gap (as in a normal pulsar) and the relativistic magnetar wind is considered. In both cases ions can be accelerated to ultra-high energies and these energetic ions can lead to production of high-energy neutrinos and gamma rays through interaction with thermal radiation from the hot neutron star or the heated inner boundary region of the stellar envelope as the result of the deposition of energy by the magnetar wind. The detectability of the neutrino flux by a kilometre-scale neutrino detector such as the planned IceCube neutrino observatory is discussed.

TeV gamma rays have been observed from blazars as well as from radio galaxies like M 87 and Centaurus A. In leptonic models, gamma rays above the pair production threshold can escape from the ultrarelativistic jet, because large Lorentz factors reduce the background photon densities compared to those required for isotropic emission. Here we discuss an alternative scenario, where very high energy photons are generated as secondaries from ultrahigh energy cosmic rays interactions in the cores of active galactic nuclei. We show that TeV gamma-rays can escape from the core despite large infrared and ultraviolet backgrounds. For the special case of Centaurus A, we study whether the various existing observations from the far infrared to the ultrahigh energy range can be reconciled within this picture.

Current models of pulsar magnetospheres typically assume either a complete absence of plasma or abundant ideal plasma filling the magnetosphere in order to compute the field structure. The latter condition is thought to be closer to reality, but we know of a number of pulsars in which the ideal conditions break down, resulting in dissipation and high-energy emission. In this work we formulate a resistive force-free scheme that allows us to consider the effects of resistive plasma and accelerating fields on the magnetospheric structure. We run numerical simulations to construct a family of resistive solutions that smoothly bridges the gap between the vacuum and the force-free magnetosphere solutions. We further provide a self-consistent model for the spin-down of intermittent pulsars, pulsars which appear to transition between radio-loud and radio-quiet states with different spin-down rates. Finally, we present models for high-energy emission from reconnecting current sheets in Gamma-ray pulsars.

The Galactic centre – as the closest galactic nucleus – holds both intrinsic interest and possibly represents a useful analogue to starburst nuclei which we can observe with orders of magnitude finer detail than these external systems. The environmental conditions in the GC – here taken to mean the inner 200 pc in diameter of the Milky Way – are extreme with respect to those typically encountered in the Galactic disk. The energy densities of the various GC ISM components are typically ~two orders of magnitude larger than those found locally and the star-formation rate density ~three orders of magnitude larger. Unusually within the Galaxy, the Galactic centre exhibits hard-spectrum, diffuse TeV (=1012 eV) gamma-ray emission spatially coincident with the region's molecular gas. Recently the nuclei of local starburst galaxies NGC 253 and M82 have also been detected in gamma-rays of such energies. We have embarked on an extended campaign of modelling the broadband (radio continuum to TeV gamma-ray), non- thermal signals received from the inner 200 pc of the Galaxy. On the basis of this modelling we find that star-formation and associated supernova activity is the ultimate driver of the region's non-thermal activity. This activity drives a large-scale wind of hot plasma and cosmic rays out of the GC. The wind advects the locally-accelerated cosmic rays quickly, before they can lose much energy in situ or penetrate into the densest molecular gas cores where star-formation occurs. The cosmic rays can, however, heat/ionize the lower density/warm H2 phase enveloping the cores. On very large scales (~10 kpc) the non-thermal signature of the escaping GC cosmic rays has probably been detected recently as the spectacular ‘Fermi bubbles’ and corresponding ‘YWMAP haze’.

The nearby starburst galaxies M82 and NGC 253 are now detected in GeV and TeV γ-rays, allowing us to directly study cosmic rays (CRs) in starburst galaxies. Combined with radio observations, the detections constrain the propagation and density of CRs in these starbursts. We discuss the implications for “proton calorimetry”, whether CR protons cool through pion losses before escaping these galaxies. The ratio of γ-ray and radio luminosities constrains how much of the CR electron cooling is due to synchrotron losses. As for leptonic emission, we predict that synchrotron and Inverse Compton emission make up ~1–10% of the unresolved hard X-ray emission from M82, and a few percent or less of the total X-ray emission from starbursts. A detection of these components would inform us of the magnetic field strength and 10 – 100 TeV electron spectrum. We conclude by discussing the prospects for detecting leptonic MeV γ-rays from starbursts and the cosmic γ-ray background.

I present a model for the non-thermal production of electromagnetic radiation in the jets of radio galaxies. The model goes far beyond the simple one-zone models usually applied to these sources. The transport equation is solved in the co-moving frame of the jet, taken into account the inhomogeneous structure of the outflow. Energy distributions for all types of particles are then obtained in a self-consistent way, including protons, electrons, and secondaries. The spectral energy distribution resulting from all relevant radiative processes is computed, including synchrotron radiation, relativistic Bremsstrahlung, proton-proton collisions and subsequent decays, photo-meson production, radiation from pairs formed by photon absorption and injection from decays, as well as direct pair production. Absorbing fields in the host galaxy are considered when computing the final SED. The model is applied to Centaurus A and compared with the available multi-wavelength data.

Starting from first principles, we construct a simple model for the evolution of energetic particles produced by supernovae in the starburst galaxy M82. The supernova rate, geometry, and properties of the interstellar medium are all well observed in this nearby galaxy. Assuming a uniform interstellar medium and constant cosmic-ray injection rate, we estimate the cosmic-ray proton and primary & secondary electron/positron populations. From these particle spectra, we predict the gamma ray flux and the radio synchrotron spectrum. The model is then compared to the observed radio and gamma-ray spectra of M82 as well as previous models by Torres (2004), Persic et al. (2008), and de Cea del Pozo et al. (2009). Through this project, we aim to build a better understanding of the calorimeter model, in which energetic particle fluxes reflect supernova rates, and a better understanding of the radio-FIR correlation in galaxies.

The impact of non-thermal processes on the spectral energy distributions of galaxies can be dramatic, but such processes are often neglected in considerations of their structure and evolution. Particle acceleration associated with high mass star formation and AGN activity not only leads to very broad band (radio-γ-ray) emission, but may also produce very significant feedback effects on galaxies and their environment. The recent detections of starburst galaxies at GeV and TeV energies suggest that γ-ray instruments have now reached the critical level of sensitivity to probe the connection between particle acceleration and star-formation in galaxies. In this paper we will try to summarise this recent progress, put it into a multi-wavelength context and also discuss the prospects for more precise and sensitive γ-ray measurements with the upcoming CTA observatory.

We present a lepto-hadronic model for the VHE emission from the relativistic jets of FR I radiogalaxies. We assume that protons and electrons are accelerated in a compact region near the base of the jet, and they cool emitting multiwavelength radiation as they propagate along the jet. The particle distributions are obtained using an inhomogeneous steady-state transport equation that accounts for the cooling processes as well as the convection of particles in the jet. The dominant processes that contribute to the photon SED are electron and proton synchrotron radiation, inverse Compton interactions, and the inelastic collisions pp and pγ. The accompanying neutrino output is obtained and the possibility of detection with Km3Net and IceCube is discussed for the cases of Cen A and M87.