The hot accretion flow around Kerr black holes is strongly magnetized. Magnetic field loops sustained by a surrounding accretion disk can close within the event horizon. We performed particle-in-cell simulations in Kerr metric to capture the dynamics of the electromagnetic field and of the ambient collisionless plasma in this coupled configuration. We find that a hybrid magnetic topology develops with a closed magnetosphere co-existing with open field lines threading the horizon reminiscent of the Blandford-Znajek solution. Further in the disk, highly inclined open magnetic field lines can launch a magnetically-driven wind. While the plasma is essentially force-free, a current sheet forms above the disk where magnetic reconnection produces macroscopic plasmoids and accelerates particles up to relativistic Lorentz factors. A highly dynamic Y-point forms on the furthest closed magnetic field line, with episodic reconnection events responsible for transient synchrotron emission and coronal heating.

The Radio Neutrino Observatory Greenland (RNO-G) is currently being deployed and it is currently gathering data. As a precursor and complementary detector to the future radio array of IceCube-Gen2 in Antarctica, it will explore mainly the Northern sky via in-ice radio detection technique. The total array configuration includes 35 radio stations and will be fully completed within three years from now. The antennas will register the radio signals produced by the Askaryan effect in cascades generated in ice by neutrinos. RNO-G’s scientific purpose is to detect UHE neutrinos at energies above 10 PeV. Due to the attenuation length of radio waves in ice (order of 1 km) the radio detection allows to address neutrino energies above several PeV. The detector will reach unprecedented sensitivity in the scale from tens of PeV up to EeV. Models predict GRBs induced by binary neutron star mergers as likely transient sources of such highly energetic neutrinos. The current study of NS-NS mergers will therefore possibly be complemented by future RNO-G detections through multimessenger temporal and spatial coincidence, including an alert system. In this presentation, we will describe the instrument capabilities and explore the possibility of detection of such sources with RNO-G.

The focus of this work is to comprehensively understand hydro-dynamical back-flows and their role in dynamics and non-thermal spectral signatures particularly during the initial phase of X-shaped radio galaxies. In this regard, we have performed axisymmetric (2D) and three dimensional (3D) simulations of relativistic magneto-hydrodynamic jet propagation from tri-axial galaxies. High-resolution dynamical modelling of axisymmetric jets has demonstrated the effect of magnetic field strengths on lobe and wing formation. Distinct X-shape formation due to back-flow and pressure gradient of ambient is also observed in our 3D dynamical run. Furthermore, the effect of radiative losses and diffusive shock acceleration on the particle spectral evolution is demonstrated, which particularly highlights how crucial their contributions are in the emission signature of these galaxies. This imparts a significant effect on the galaxy’s equipartition condition, indicating that one must be careful in extending its use in estimating other parameters, as the criterion evolves with time.

Astrophysical systems possess various sites of particle acceleration, which gives rise to the observed non-thermal spectra. Diffusive shock acceleration (DSA) and stochastic turbulent acceleration (STA) are the candidates for producing very high energy particles in weakly magnetized regions. While DSA is a systematic acceleration process, STA is a random energization process, usually modelled as a biased random walk in energy space with a Fokker-Planck equation. In astrophysical systems, different acceleration processes work in an integrated manner along with various energy losses.

Here we study the interplay of both STA and DSA in addition to various energy losses, in a simulated RMHD jet cocoon. Further, we consider a phenomenologically motivated STA timescale and discuss its effect on the emission profile of the RMHD jet. A parametric study on the turbulent acceleration timescale is also conducted to showcase the effect of turbulence damping on the emission structure of the simulated jet.

Particle acceleration induced by fast magnetic reconnection may help to solve current puzzles related to the interpretation of the very high energy (VHE) and neutrino emissions from AGNs and compact sources in general. Our general relativistic-MHD simulations of accretion disk-corona systems reveal the growth of turbulence driven by MHD instabilities that lead to the development of fast magnetic reconnection in the corona. In addition, our simulations of relativistic MHD jets reveal the formation of several sites of fast reconnection induced by current-driven kink turbulence. The injection of thousands of test particles in these regions causes acceleration up to energies of several PeVs, thus demonstrating the ability of this process to accelerate particles and produce VHE and neutrino emission, specially in blazars. Finally, we discuss how reconnection can also explain the observed VHE luminosity-black hole mass correlation, involving hundreds of non-blazar sources like Perseus A, and black hole binaries.

Diffuse radio emission from galaxy clusters in the form of radio halos and relics are tracers of the shocks and turbulence in the intra-cluster medium. The imprints of the physical processes that govern their origin and evolution can be found in their radio morphologies and spectra. The role of mildly relativistic population of electrons may be crucial for the acceleration mechanisms to work efficiently. Low frequency observations with telescopes that allow imaging of extended sources over a broad range of low frequencies (<2 GHz) offer the best tools to study these sources. I will review the Giant Metrewave Radio Telescope (GMRT) observations in the past few years that have led to: i) statistical studies of large samples of galaxy clusters, ii) opening of the discovery space in low mass clusters and iii) tracing the spectra of seed relativistic electrons using the Upgraded GMRT.

We present recent developments on numerical algorithms for computing photon and particle trajectories in the surrounding of compact objects. Strong gravity around neutron stars or black holes causes relativistic effects on the motion of massive particles and distorts light rays due to gravitational lensing. Efficient numerical methods are required for solving the equations of motion and compute i) the black hole shadow obtained by tracing light rays from the object to a distant observer, and ii) obtain information on the dynamics of the plasma at the microscopic scale. Here, we present generalized algorithms capable of simulating ensembles of photons or massive particles in any spacetime, with the option of including external forces. The coupling of these tools with GRMHD simulations is the key point for obtaining insight on the complex dynamics of accretion disks and jets and for comparing simulations with upcoming observational results from the Event Horizon Telescope.

We calculate the luminosity of the synchrotron radiation from the vicinity of the light cylinder. We find that even if the thermal emission from the entire surface is included as the seed photon, the γ-ray to X-ray flux ratio for young pulsars is significantly higher than the observations. For these pulsars, most of γ-ray photons may be absorbed in the magnetosphere.

Over the last fifty years since the discovery of pulsars, our understanding of where and how pulsars emit the radiation we observe has undergone significant revision. The location and mechanisms of high-energy radiation are intimately tied to the sites of particle acceleration. The evolution of emission models has paralleled the development of increasingly more sensitive telescopes, especially at high energies. I will review the history of pulsar emission modeling, from the early days of gaps at the polar caps, to outer gaps and slot gaps in the outer magnetosphere, to the present era of global magnetosphere simulations that locate most acceleration and high-energy emission in the current sheets.

The standard shock acceleration model of pulsar wind nebulae (PWNe) does not account for the hard spectrum in radio wavelengths. The origin of the radio-emitting particles is also important to determine the pair production efficiency in the pulsar magnetosphere. Here, we propose a possible resolution for the particle energy distribution in PWNe; the radio-emitting particles are not accelerated at the pulsar wind termination shock but are stochastically accelerated by turbulence inside PWNe. We upgrade our past one-zone spectral evolution model including the energy diffusion, i.e., the stochastic acceleration, and apply to the Crab Nebula. For a particle injection to the stochastic acceleration process, we consider the continuous injection from the supernova ejecta or the impulsive injection associated with supernova explosion. The observed broadband spectrum and the decay of the radio flux are reproduced by tuning the amount of the particle injected to the stochastic acceleration process. Our results imply that some unveiled mechanisms, such as back reaction to the turbulence, are required to make the energies of stochastically and shock accelerated particles comparable.

Follow-up of unidentified Fermi sources has expanded the number of known galactic-field “black widow” and “redback” millisecond pulsar binaries from four to nearly 30. Several systems observed by Chandra, XMM-Newton, Suzaku, and NuSTAR exhibit double-peaked X-ray orbital modulation. This is attributed to synchrotron emission from electrons accelerated in an intrabinary shock and Doppler boosting by mildly relativistic bulk flow. We briefly discuss the rich complexity of these systems, their astrophysical utility, and open questions.

Using the appropriate kinetic equation, we considered the problem of propagation of accelerated electrons into the solar corona and chromosphere. Its analytical solution was used for modelling the M7.7 class limb flare occurred on July 19, 2012. Coronal above-the-loop-top hard X-Ray source was interpreted in the thin-target approximation, the foot-point source - in the thick-target approximation with account of the reverse-current electric field. For the foot-point source we found a good accordance with the RHESSI observations. For the coronal source we also got very accurate estimate of the power-law spectral index, but significant differences between the modelled and observed hard X-ray intensities were noticed. The last discrepancy was solved by adding the coronal magnetic trap model to the thin target model. The former one implies that the trap collapses in two dimensions, locks and accelerates particles inside itself. In our report, we confirm an existence and high efficiency of the electron acceleration in collapsing magnetic traps during solar flares. Our new results represent (e.g. for RHESSI observations) the theoretical prediction of the double step particle acceleration in solar flares, when the first step is the acceleration in reconnection area and the second one – the acceleration in coronal trap.

Additional acceleration of protons and electrons passing through a 3D reconnecting current sheet (RCS) of the solar corona and heliosphere is investigated with PIC approach. The simulation confirms spatial separation of electrons and protons and generation of a polarisation electric field induced by separated particles. In the heliospheric current sheet with a weak magnetic field there are two populations of particles: transit and bounced. The transit particles (both protons and electrons) are accelerated to high energies while the bounced electrons fail to reach the midplane. Instead they form an electron cloud of horseshoe or medallion types at some distance D from its midplane, which is larger for bigger guiding field magnitudes. These energetic electrons and protons appearing near the HCS boundaries can be a great danger for the satellites crossing the sector boundaries.

Gamma-ray observations for Supernova remnant (SNR)-molecular cloud (MC) association systems play an important role in the research on the acceleration and propagation of cosmic-ray protons. Through the analysis of 5.6 years of Fermi-Large Area Telescope observation data, here we report on the detection of a gamma-ray emission source near the SNR Kesteven 41 with a significance of 24σ in 0.2–300 GeV. The best-fit location of the gamma-ray source is consistent with the MC with which the SNR interacts. Several hypotheses including both leptonic and hadronic scenarios are considered to investigate the origin of these gamma-rays. The gamma-ray emission can be naturally explained by the decay of neutral pions produced via the collision between high energy protons accelerated by the shock of Kesteven 41 and the adjacent MC. The electron energy budget would be too high for the SNR if the gamma-rays were produced via inverse Compton (IC) scattering off the Cosmic Microwave Background (CMB) photons.

G326.3-1.8 (also known as MSH 15-56) has been detected in radio as a middle-aged composite supernova remnant (SNR) consisting of a SNR shell and a pulsar wind nebula (PWN) which has been crushed by the reverse shock. With the recent Fermi-LAT data release Pass 8 providing increased acceptance and angular resolution, we investigate the morphology of this SNR to disentangle the PWN from the SNR contributions and understand the nature of the γ-ray emission. We thus perform a morphological and spectral analysis from 300 MeV to 300 GeV which highlights the contributions from these two components. The simplest interpretation is hadronic emission from the SNR and harder leptonic emission from the PWN.

Recent direct measurements of cosmic-ray (CR) light nuclei (protons, helium, and lithium) by AMS-02 have shown that the flux of each element has an unexpected hard component above ~300~GeV, and that the spectral indices of those components are almost the same (~2.5). This implies that there should be primary sources that produces CR lithium nuclei, which have been believed to be produced via spallation of heavier nuclei in the ISM (secondary origin). We propose the nearby Type Ia supernova following a nova eruption from a white dwarf as the origin of CR Li.

The Cosmic Ray (CR) physics has entered a new era driven by high precision measurements coming from direct detection (especially AMS-02 and PAMELA) and also from gamma-ray observations (Fermi-LAT). In this review we focus our attention on how such data impact the understanding of the supernova remnant paradigm for the origin of CRs. In particular we discuss advancement in the field concerning the three main stages of the CR life: the acceleration process, the escape from the sources and the propagation throughout the Galaxy. We show how the new data reveal a phenomenology richest than previously thought that could even challenge the current understanding of CR origin.

Supernova remnants (SNRs) are powerful particle accelerators. As a supernova (SN) blast wave propagates through the circumstellar medium (CSM), electrons and protons scatter across the shock and gain energy by entrapment in the magnetic field. The accelerated particles generate further magnetic field fluctuations and local amplification, leading to cosmic ray production. The wealth of data from Supernova 1987A is providing a template of the SN-CSM interaction, and an important guide to the radio detection and identification of core-collapse SNe based on their spectral properties. Thirty years after the explosion, radio observations of SNR 1987A span from 70 MHz to 700 GHz. We review extensive observing campaigns with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA), and follow-ups with other radio telescopes. Observations across the radio spectrum indicate rapid changes in the remnant morphology, while current ATCA and ALMA observations show that the SNR has entered a new evolutionary phase.

A large fraction of core-collapse supernovae are thought to result in the birth of a rotation-powered pulsar, which is later observable as a radio pulsar up to great ages. The birth properties of these pulsars, and in particular the distribution of their initial rotation periods, are however difficult to infer from studies of the radio pulsar population in our Galaxy. Yet the distributions of their birth properties is an important assumption for scenarios in which ultra-high-energy cosmic rays (UHECRs) originate in very young, extragalactic pulsars with short birth periods and/or high magnetic fields.

Using a model of the very young pulsar wind nebula’s dynamical and spectral evolution, with pulsar wind and accelerated particle parameters assumed similar to those inferred from modeling young pulsar wind nebulae (PWNe) in our Galaxy, we show that X-ray observations of supernovae, a few years to decades after the explosion, constitute a favored window to obtain meaningful constraints on the initial spin-down luminosity of the newly-formed pulsar. We examine the expected emerging PWN spectral component, taking into account the X-ray opacity of the expanding supernova ejecta, and find that it is typically best detectable in < 10 keV X-rays some years after the explosion. We use this framework to assess available X-ray observations and flux upper limits on supernovae, building on the work of Perna et al. (2008). We note that a resulting limit on spin-down luminosity corresponds univocally to a limit on the maximum magnetospheric acceleration potential, irrespective of the specific combination of magnetic field and rotation period that achieves it. We use available X-ray observations of supernovae to place constraints on the birth spin-down luminosity and period distribution of classical pulsars. We also examine the case of magnetars, born with much higher magnetic fields, and show that their much shorter initial spin-down time implies that any plausible signature of young magnetar wind nebulae can only be observed in harder X-ray or gamma-rays.

Supernova remnants (SNRs) are one of the most important sites where particles are accelerated with high efficiency and in a wide range of energies, becoming an important component of cosmic rays. A good test for this hypothesis will be possible using the data collected by next-generation radio and gamma-ray observatories, like the Square Kilometre Array (SKA) and the Cherenkov Telescope Array (CTA). Radio emission is fundamental to explore the SNR environment and to shed light on the physical processes involved in particle acceleration, providing direct links to high-energy physics. Two cases of SNRs recently studied in radio are presented, showing the importance of high-resolution radio images. An overview of SKA and its precursors is given with our ongoing preparation work. In particular, we present the EMU survey and the pathfinder project SCORPIO. Finally a direct view of the tight connection between SKA and CTA future studies of SNRs is provided.