The variation in the intensity of cosmic rays at small angular scales is attributed to the interstellar turbulence in the vicinity of the Solar system. We show that {a turbulent origin of the small-scale structures implies that} the morphology of the observed cosmic-ray intensity skymap varies with our location in the interstellar turbulence. The gyroradius of cosmic rays is shown to be the length scale associated with an observable change in the skymap over a radian angular scale. The extent to which the intensity \mpo{at a certain} angular scale varies is proportional to the change in our location with a maximum change of about the amplitude of intensity variation at that scale in the existing skymap.} We suggest that for TeV cosmic rays a measurable variation could occur over a time scale of a decade due to the Earth's motion through the interstellar medium, if interstellar turbulence persists down to the gyroradius, \mpo{about $00\ \mu\mathrm{pc}$ for TeV-ish cosmic rays}. Observational evidence of the variability, or an absence of it, could provide a useful insight into the physical origin of the small-scale anisotropy.

Motivated by the ~7 % dipole anisotropy in the distribution of ultra-high energy cosmic-rays (UHECRs) above 8 EeV, we explore the anisotropy induced by the large scale structure, using constrained simulations of the local Universe and taking into account the effect of magnetic fields. The value of the intergalactic magnetic field (IGMF) is critical as it determines the UHECR cosmic horizon. We calculate the UHECR sky maps for different values of the IGMF variance and show the effect of the UHECR horizon on the observed anisotropy. The footprint of the local ($\lesssim350$ Mpc) Universe on the UHECR background, a small angular scale enhancement in the Northern Hemisphere, is seen. At 11.5 EeV (the median value of the energy bin at which the dipole has been reported), the LSS-induced dipole amplitude is $A_1\sim10\%$, for IGMF in the range [0.3-3] nG for protons, helium and nitrogen, compatible with the rms value derived from the cosmic power spectrum. However at these energies the UHECRs are also influenced by the Galactic Magnetic Field (GMF) and we discuss its effect on the LSS-induced anisotropy.

We explore the possibility that the recently detected dipole anisotropy in the arrival directions of~$>8$~EeV ultra-high energy cosmic-rays (UHECRs) arises due to the large-scale structure (LSS). We assume that the cosmic ray sources follow the matter distribution and calculate the flux-weighted UHECRs' RMS dipole amplitude taking into account the diffusive transport in the intergalactic magnetic field (IGMF). We find that the flux-weighted RMS dipole amplitude is $\sim8$% before entering the Galaxy. The amplitude in the [4-8] EeV is only slightly lower $\sim 5$%. The required IGMF is of the order of {5-30 nG}, and the UHECR sources must be relatively nearby, within $\sim$300 Mpc. The absence of statistically significant signal in the lower energy bin can be explained if the same nuclei specie dominates the composition in both energy bins and diffusion in the Galactic magnetic field (GMF) reduces the dipole of these lower rigidity particles. Photodisintegration of higher energy UHECRs could also reduce somewhat the lower energy dipole.

The Telescope Array (TA) shows a 20$^{\circ}$ hotspot as well as an excess of UHECRs above 50~EeV when compared with the Auger spectrum. We consider the possibility that both the TA excess and hotspot are due to a dominant source in the Northern sky. We carry out detailed simulations of UHECR propagation in both the intergalactic medium and the Galaxy, using different values for the intergalactic magnetic field. We consider two general classes of sources: transients and steady, adopting a mixed UHECR composition that is consistent with the one found by Auger. The spatial location of the sources is draw randomly. We generate Auger-like and TA-like data sets from which we determine the spectrum, the sky maps and the level of anisotropy. We find that, while steady sources are favored over transients, it is unlikely to account for all the currently available observational data. While we reproduce fairly well the Auger spectrum for the vast majority of the simulated data sets, most of the simulated data sets with a spectrum compatible with that of TA (at most a few percent depending on density model tested) show a much stronger anisotropy than the one observed. We find that the rare cases in which both the spectrum and the anisotropy are consistent require a steady source within $\sim 10$ Mpc, to account for the flux excess, and a strong extragalactic magnetic field $\sim 10$ nG, to reduce the excessive anisotropy.

GeV-TeV gamma-ray and PeV-EeV neutrino backgrounds provide a unique window on the nature of the ultra-high-energy cosmic-rays (UHECRs). We discuss the implications of the recent Fermi-LAT data regarding the extragalactic gamma-ray background (EGB) and related estimates of the contribution of point sources as well as IceCube neutrino data on the origin of the UHECRs. We calculate the diffuse flux of cosmogenic $\gamma$-rays and neutrinos produced during the UHECRs propagation and derive constraints on the possible cosmological evolution of UHECR sources. In particular, we show that the mixed-composition scenario which is in agreement with both (i) Auger measurements of the energy spectrum and composition up to the highest energies and (ii) the ankle-like feature in the light component detected by KASCADE-Grande, is compatible with both the Fermi-LAT measurements and with current IceCube limits.

GeV-TeV gamma-rays and PeV-EeV neutrino backgrounds provide a unique window on the nature of the ultra-high-energy cosmic-rays (UHECRs). We discuss the implications of the recent Fermi-LAT data regarding the extragalactic gamma-ray background (EGB) and related estimates of the contribution of point sources as well as IceCube neutrino data on the origin of the UHECRs. We calculate the diffuse flux of cosmogenic $\gamma$-rays and neutrinos produced by the UHECRs and derive constraints on the possible cosmological evolution of UHECR sources. In particular, we show that the mixed-composition scenario considered in \citet{Globus2015b}, which is in agreement with both (i) Auger measurements of the energy spectrum and composition up to the highest energies and (ii) the ankle-like feature in the light component detected by KASCADE-Grande, is compatible with both the Fermi-LAT measurements and with current IceCube limits. We also discuss the possibility for future experiments to detect associated cosmogenic neutrinos and further constrain the UHECR models, including possible subdominant UHECR proton sources.

The Telescope Array (TA) shows a 20$^{\circ}$ hotspot as well as an excess of UHECRs above 50~EeV when compared with the Auger spectrum. We consider the possibility that both the TA excess and hotspot are due to a dominant source in the Northern sky. We carry out detailed simulations of UHECR propagation in both the intergalactic medium and the Galaxy, using different values for the intergalactic magnetic field. We consider two general classes of sources: transients and steady, adopting a mixed UHECR composition that is consistent with the one found by Auger. The spatial location of the sources is draw randomly. We generate Auger-like and TA-like data sets from which we determine the spectrum, the sky maps and the level of anisotropy. We find that, while steady sources are favored over transients, it is unlikely to account for all the currently available observational data. While we reproduce fairly well the Auger spectrum for the vast majority of the simulated data sets, most of the simulated data sets with a spectrum compatible with that of TA (at most a few percent depending on density model tested) show a much stronger anisotropy than the one observed. We find that the rare cases in which both the spectrum and the anisotropy are consistent require a steady source within $\sim 10$ Mpc, to account for the flux excess, and a strong extragalactic magnetic field $\sim 10$ nG, to reduce the excessive anisotropy.

Stationary features are occasionally observed in AGN jets. A notable example is the HST-1 knot in M87. Such features are commonly interpreted as re-confinement shocks in hydrodynamic jets or focusing nozzles in Poynting jets. In this paper we compute the structure and Lorentz factor of a highly magnetized jet confined by external pressure having a profile that flattens abruptly at some radius. We find the development of strong oscillations upon transition from the steeper to the flatter pressure profile medium. Analytic formula is derived for the location of the nodes of these oscillations. We apply the model to the M87 jet and show that if the jet remains magnetically dominated up to sub-kpc scales, then focusing is expected. The location of the HST-1 knot can be reconciled with recent measurements of the pressure profile around the Bondi radius if the jet luminosity satisfies $L_j\simeq10^{43}$ erg/s. However, we find that magnetic domination at the collimation break implies a Lorentz factor in excess of $10^2$, atypical to FRI sources. A much lower value of the asymptotic Lorentz factor would require substantial loading close to the black hole. In that case HST-1 may be associated with a collimation nozzle of a hydrodynamic flow.

The collimation of a Poynting-flux dominated jet by a wind emanating from the surface of an accretion flow is computed using a semi-analytic model. The injection of the disk wind is treated as a boundary condition in the equatorial plane, and its evolution is followed by invoking a prescribed geometry of streamlines. Solutions are obtained for a wide range of disk wind parameters. It is found that jet collimation generally occurs when the total wind power exceeds about ten percents of the jet power. For moderate wind powers we find gradual collimation. For strong winds we find rapid collimation followed by focusing of the jet, after which it remains narrow over many Alfv\'en crossing times before becoming conical. We estimate that in the later case the jet's magnetic field may be dissipated by the current-driven kink instability over a distance of a few hundreds gravitational radii. We apply the model to M87 and show that the observed parabolic shape of the radio jet within the Bondi radius can be reproduced provided that the wind injection zone extends to several hundreds gravitational radii, and that its total power is about one third of the jet power. The radio spectrum can be produced by synchrotron radiation of relativistically hot, thermal electrons in the sheath flow surrounding the inner jet.

It is suggested that essentially all of the UHECRs we detect, including those at the highest energy, originate in our Galaxy. It is shown that even if the density of sources decreases with Galactic radius, then the anisotropy and composition can be understood. Inward anisotropy, as recently reported by the Auger collaboration can be understood as drift along the current sheet of UHECRs originating outside the solar circle, as predicted in Kumar and Eichler (2014), while those originating within the solar circle exit the Galaxy at high latitudes.

The dynamics of a dissipative Poynting dominated flow subject to a radiation drag due to Compton scattering of ambient photons by relativistic electrons accelerated in reconnecting current sheets is studied. It is found that the efficiency at which magnetic energy is converted to radiation is limited to a maximum value of $\epsilon_c=3l_{dis}\,\sigma_0/4(\sigma_0+1)$, where $\sigma_0$ is the initial magnetization of the flow and $l_{dis}\le1$ the fraction of initial Poynting flux that can dissipate. The asymptotic Lorentz factor satisfies $\Gamma_\infty\ge\Gamma_0(1+l_{dis}\,\sigma_0/4)$, where $\Gamma_0$ is the initial Lorentz factor. This limit is approached in cases where the cooling time is shorter than the local dissipation time. A somewhat smaller radiative efficiency is expected if radiative losses are dominated by synchrotron and SSC emissions. It is suggested that under certain conditions magnetic field dissipation may occur in two distinct phases: On small scales, asymmetric magnetic fields that are advected into the polar region and dragged out by the outflow dissipate to a more stable configuration. The dissipated energy is released predominantly as gamma rays. On much larger scales, the outflow encounters a flat density profile medium and re-collimates. This leads to further dissipation and wobbling of the jet head by the kink instability, as found recently in 3D MHD simulations. Within the framework of a model proposed recently to explain the dichotomy of radio loud AGN, this scenario can account for the unification of gamma-ray blazars with FRI and FRII radio sources.

We present a complete phenomenological model accounting for the evolution of the cosmic-ray spectrum and composition with energy, based on the available data over the entire spectrum. We show that there is no need to postulate any additional component, other than one single Galactic component depending on rigidity alone, and one extragalactic component, whose characteristics are similar to those derived from a study of particle acceleration at mildly relativistic shocks in a GRB environment (Globus et al., 2015). In particular, we show that the resulting cosmic ray spectrum and composition satisfy the various constraints derived from the current data in the Galactic/extragalactic transition region, notably from the measurements of KASCADE Grande and Auger. Finally, we derive some generic features that a working phenomenological scenario may exhibit to give a global account of the cosmic ray data with a minimum number of free parameters.

We study the acceleration of CR protons and nuclei at GRB internal shocks. Physical quantities and their time evolution are estimated using the internal shock modeling implemented by Daigne & Mochkovitch 1998. We consider different hypotheses about the way the energy dissipated at internal shocks is shared between accelerated CR, e- and B field. We model CR acceleration at mildly relativistic shocks, including all the significant energy loss processes. We calculate CR and neutrino release from single GRBs, assuming that nuclei heavier than protons are present in the relativistic wind. Protons can only reach maximum energies of ~ 10^19.5 eV, while intermediate and heavy nuclei are able to reach values of ~ 10^20 eV and above. The spectra of nuclei escaping from the acceleration site are found to be very hard while the combined spectrum of protons and neutrons is much softer. We calculate the diffuse UHECR flux expected on Earth using the GRB luminosity function from Wanderman & Piran 2010. Only the models assuming that the prompt emission represent a very small fraction of the energy dissipated at internal shocks, and that most of this dissipated energy is communicated to accelerated CR, are able to reproduce the magnitude of the UHECR flux observed. For these models, the observed shape of the UHECR spectrum can be well reproduced and the evolution of the composition is compatible with the trend suggested by Auger. We discuss implications of the softer proton component for the GCR to EGCR transition in the light of the recent composition analyses (KASCADE-Grande experiment). The associated secondary particle diffuse fluxes do not upset any current observational limit. Diffuse neutrino flux from GRB sources should however be detected with the lifetime of neutrino observatories.

We construct a semi-analytic model for MHD flows in Kerr geometry, that incorporates energy loading via neutrino annihilation on magnetic field lines threading the horizon. We compute the double-flow structure for a wide range of energy injection rates, and identify the different operation regimes. At low injection rates the outflow is powered by the spinning black hole via the Blandford-Znajek mechanism, whereas at high injection rates it is driven by the pressure of the plasma deposited on magnetic field lines. In the intermediate regime both processes contribute to the outflow formation. The parameter that quantifies the load is the ratio of the net power injected below the stagnation radius and the maximum power that can be extracted magnetically from the black hole.

We present a model for the spine of relativistic MHD outflows in the Kerr geometry. Meridional self-similarity is invoked to derive semi-analytical solutions close to the polar axis. The study of the energy conservation along a particular field line gives a simple criterion for the collimation of jets. Such parameter have already been derived in the classical case by Sauty et al. 1999 and also extended to the Schwarzschild metric by Meliani et al. 2006. We generalize the same study to the Kerr metric. We show that the rotation of the black hole increases the magnetic self-confinement.

We discuss the role plasma injection plays in the formation of outflows in Kerr spacetime. Using a model for the double flow established in the polar region of a rotating black hole, we study the interplay between the different processes that can power the outflow. In particular, we find two types of flows with distinct properties that depend on the rate at which energy is deposited in the magnetosphere. We discuss the implications of this result for gamma ray bursts outflows.

The effect of mass and energy loading on the efficiency at which energy can be extracted magnetically from a Kerr black hole is explored, using a semi-analytic, ideal MHD model that incorporates plasma injection on magnetic field lines. We find a critical load below which the specific energy of the plasma inflowing into the black hole is negative, and above which it is positive, and identify two types of flows with distinct properties; at sub-critical loads a magnetic outflow is launched from the ergosphere, owing to extraction of the black hole spin energy, as originally proposed by Blandford and Znajek. At super-critical loads the structure of the flow depends on the details of the injection process. In cases where the injected plasma is relativistically hot, a pressure driven, double trans-magnetosonic flow is launched from a stagnation point located outside the ergosphere, between the inner and outer light cylinders. Some fraction of the energy deposited in the magnetosphere is then absorbed by the black hole and the rest emerges at infinity in the form of a relativistic outflow. When the injected plasma is cold an outflow may not form at all. We discuss the implications of our results to gamma ray bursts and active galactic nuclei.

The structure of a hydrodynamic, double transonic flow driven by neutrino annihilation in the polar region of a Schwarzschild black hole is computed for different energy deposition profiles. The requirement that both, the inflow into the black hole and the outflow to infinity pass smoothly through their sonic points fixes the stagnation radius and stagnation pressure. The asymptotic power of the outflow is shown to be the integral of the energy deposition rate above the stagnation radius. The outflow production efficiency depends on the energy deposition profile, and is generally higher for shallower profiles. Using recent calculations of the neutrino annihilation rate, we estimate that over 50 percents of the total energy deposited above the horizon can emerge in the form of a relativistic outflow at infinity. The continuous creation of plasma during the expansion of the outflow leads to generation of a large specific entropy. This has important implications for the prompt photospheric emission in GRBs.

Context: A large sample of T Tauri stars exhibits optical jets, approximately half of which rotate slowly, only at ten per cent of their breakup velocity. The disk-locking mechanism has been shown to be inefficient to explain this observational fact. Aims: We show that low mass accreting T Tauri stars may have a strong stellar jet component that can effectively brake the star to the observed rotation speed. Methods: By means of a nonlinear separation of the variables in the full set of the MHD equations we construct semi- analytical solutions describing the dynamics and topology of the stellar component of the jet that emerges from the corona of the star. Results: We analyze two typical solutions with the same mass loss rate but different magnetic lever arms and jet radii. The first solution with a long lever arm and a wide jet radius effectively brakes the star and can be applied to the visible jets of T Tauri stars, such as RY Tau. The second solution with a shorter lever arm and a very narrow jet radius may explain why similar stars, either Weak line T Tauri Stars (WTTS) or Classical T Tauri Stars (CTTS) do not all have visible jets. For instance, RY Tau itself seems to have different phases that probably depend on the activity of the star. Conclusions: First, stellar jets seem to be able to brake pre-main sequence stars with a low mass accreting rate. Second, jets may be visible only part time owing to changes in their boundary conditions. We also suggest a possible scenario for explaining the dichotomy between CTTS and WTTS, which rotate faster and do not have visible jets.

We extend previous studies of mixed-composition extragalactic cosmic-ray source models, by investigating the influence of a non-negligible extragalactic magnetic field on the propagated cosmic-ray spectrum and composition. We study the transport of charged particles in turbulent fields and the transition from a ballistic to a diffusive propagation regime. We introduce a method allowing a fast integration of the particle trajectories, which allows us to calculate extragalactic cosmic-ray spectra in the general case, without using either the diffusive or the rectilinear approximation. We find that the main features of the mixed-composition models -- regarding the interpretation of the ankle and the non-monotonous evolution of the average cosmic-ray mass -- remain essentially unchanged as long as the magnetic field intensity does not exceed a few nG.