Since the detection of non-thermal radio emission from the bow shock of the massive runaway star BD +43$^{\circ}$3654 simple models have predicted high-energy emission, at X and gamma-rays, from these Galactic sources. Observational searches for this emission so far give no conclusive evidence but a few candidates at gamma rays. In this work we aim at developing a more sophisticated model for the non-thermal emission from massive runaway star bow shocks. The main goal is to establish whether these systems are efficient non-thermal emitters, even if they are not strong enough to be yet detected. For modeling the collision between the stellar wind and the interstellar medium we use 2D hydrodynamic simulations. We then adopt the flow profile of the wind and the ambient medium obtained with the simulation as the plasma state for solving the transport of energetic particles injected in the system, and the non-thermal emission they produce. For this purpose we solve a 3D (2 spatial + energy) advection-diffusion equation in the test-particle approximation. We find that a massive runaway star with a powerful wind converts 0.16-0.4% of the power injected in electrons into non-thermal emission, mostly produced by inverse Compton scattering of dust-emitted photons by relativistic electrons, and secondly by synchrotron radiation. This represents a fraction of $\sim$ $10^{-5}-10^{-4}$ of the wind kinetic power. Given the better sensibility of current instruments at radio wavelengths theses systems are more prone to be detected at radio through the synchrotron emission they produce rather than at gamma energies.

The bowshocks of runaway stars had been theoretically proposed as gamma-ray sources. However, this hypothesis has not been confirmed by observations up to date. In this paper, we present two runaway stars (Lambda Cephei and LS 2355) whose bowshocks are coincident with the unidentified Fermi gamma-ray sources 3FLG J2210.1+5925 and 3FGL J1128.7-6232, respectively. After performing a cross-correlation between different catalogues at distinct wavelengths, we found that these bowshocks are the most peculiar objects in the Fermi position ellipses. Then, we computed the inverse Compton emission and fitted the Fermi data in order to test the viability of both runaway stars as potential counterparts of the two high-energy sources. We obtained very reasonable values for the fitted parameters of both stars. We also evaluated the possibility for the source 3FGL J1128.7-6232, which is positionally coincident with an HII region, to be the result of background cosmic-rays protons interacting with the matter of the cloud, as well as the probability of a pure chance association. We conclude that the gamma rays from these Fermi sources might be produced in the bowshocks of the considered runaway stars. In such a case, these would be the first sources of this class ever detected at gamma rays.

High-velocity clouds (HVCs) are interstellar clouds of atomic hydrogen that do not partake of the Galactic rotation and have velocities of a several hundred kilometers per second. A considerable number of these clouds are falling down towards the Galactic disk. HVCs form large and massive complexes, so their collisions with the disk must release a great amount of energy into the interstellar medium. The cloud-disk interaction produces two shocks, one propagates through the cloud and the other through the disk; the properties of these shocks depend mainly on the cloud velocity and the disk-cloud density ratio. In this work we study the conditions necessary for these shocks to accelerate particles by diffusive shock acceleration and the produced non-thermal radiation. We analyze particle acceleration in both the cloud and disk shocks. Solving a time-dependent 2-D transport equation for both relativistic electrons and protons we obtain particle distributions and non-thermal spectral energy distributions. In a shocked cloud significant synchrotron radio emission is produced along with soft gamma rays. In the case of acceleration in the shocked disk, the non-thermal radiation is stronger; the gamma rays, of leptonic origin, might be detectable with current instruments. A large number of protons are injected into the Galactic interstellar medium, and locally exceed the cosmic-ray background. We conclude that under adequate conditions the contribution from HVC-disk collisions to the galactic population of relativistic particles and the associated extended non-thermal radiation might be important.

Massive runaway stars produce bow shocks through the interaction of their winds with the interstellar medium, with the prospect for particle acceleration by the shocks. These objects are consequently candidates for non-thermal emission. Our aim is to investigate the X-ray emission from these sources. We observed with XMM-Newton a sample of 5 bow shock runaways, which constitutes a significant improvement of the sample of bow shock runaways studied in X-rays so far. A careful analysis of the data did not reveal any X-ray emission related to the bow shocks. However, X-ray emission from the stars is detected, in agreement with the expected thermal emission from stellar winds. On the basis of background measurements we derive conservative upper limits between 0.3 and 10 keV on the bow shocks emission. Using a simple radiation model, these limits together with radio upper limits allow us to constrain some of the main physical quantities involved in the non-thermal emission processes, such as the magnetic field strength and the amount of incident infrared photons. The reasons likely responsible for the non-detection of non-thermal radiation are discussed. Finally, using energy budget arguments, we investigate the detectability of inverse Compton X-rays in a more extended sample of catalogued runaway star bow shocks. From our analysis we conclude that a clear identification of non-thermal X-rays from massive runaway bow shocks requires one order of magnitude (or higher) sensitivity improvement with respect to present observatories.

Fast radio bursts, or FRBs, are transient sources of unknown origin. Recent radio and optical observations have provided strong evidence for an extragalactic origin of the phenomenon and the precise localization of the repeating FRB 121102. Observations using the Karl G. Jansky Very Large Array (VLA) and very-long-baseline interferometry (VLBI) have revealed the existence of a continuum non-thermal radio source consistent with the location of the bursts in a dwarf galaxy. All these new data rule out several models that were previously proposed, and impose stringent constraints to new models. We aim to model FRB 121102 in light of the new observational results in the active galactic nucleus (AGN) scenario. We propose a model for repeating FRBs in which a non-steady relativistic $e^\pm$-beam, accelerated by an impulsive magnetohydrodynamic (MHD)-driven mechanism, interacts with a cloud at the centre of a star-forming dwarf galaxy. The interaction generates regions of high electrostatic field called cavitons in the plasma cloud. Turbulence is also produced in the beam. These processes, plus particle isotropization, the interaction scale, and light retardation effects, provide the necessary ingredients for short-lived, bright coherent radiation bursts. The mechanism studied in this work explains the general properties of FRB 121102, and may also be applied to other repetitive FRBs. Coherent emission from electrons and positrons accelerated in cavitons provides a plausible explanation of FRBs.

Fast magnetic reconnection may occur in different astrophysical sources, producing flare-like emission and particle acceleration. Currently, this process is being studied as an efficient mechanism to accelerate particles via a first-order Fermi process. In this work we analyse the acceleration rate and the energy distribution of test particles injected in three-dimensional magnetohydrodynamical (MHD) domains with large-scale current sheets where reconnection is made fast by the presence of turbulence. We study the dependence of the particle acceleration time with the relevant parameters of the embedded turbulence, i.e., the Alfv\'en speed $V_{\rm A}$, the injection power $P_{\rm inj}$ and scale $k_{\rm inj}$ ($k_{\rm inj} = 1/l_{\rm inj}$). We find that the acceleration time follows a power-law dependence with the particle kinetic energy: $t_{acc} \propto E^{\alpha}$, with $0.2 < \alpha < 0.6$ for a vast range of values of $c / V_{\rm A} \sim 20 - 1000$. The acceleration time decreases with the Alfv\'en speed (and therefore with the reconnection velocity) as expected, having an approximate dependence $t_{acc} \propto (V_{\rm A} / c)^{-\kappa}$, with $\kappa \sim 2.1- 2.4$ for particles reaching kinetic energies between $1 - 100 \, m_p c^2$, respectively. Furthermore, we find that the acceleration time is only weakly dependent on the $P_{\rm inj}$ and $l_{\rm inj}$ parameters of the turbulence. The particle spectrum develops a high-energy tail which can be fitted by a hard power-law already in the early times of the acceleration, in consistency with the results of kinetic studies of particle acceleration by magnetic reconnection in collisionless plasmas.

Fast Magnetic Reconnection is currently regarded as an important process also beyond the solar system, specially in magnetically dominated regions of galactic and extragalactic sources like the surrounds of black holes and relativistic jets. In this lecture we discuss briefly the theory of fast magnetic reconnection, specially when driven by turbulence which is very frequent in Astrophysical flows, and its implications for relativistic particle acceleration. Then we discuss these processes in the context of the sources above, showing recent analytical and multidimensional numerical MHD studies that indicate that fast reconnection can be a powerful process to accelerate particles to relativistic velocities, produce the associated high energy non-thermal emission, and account for efficient conversion of magnetic into kinetic energy in these flows.

Galactic cosmic rays are believed to be mostly accelerated at supernova shocks. However, the interstellar magnetic field is too weak to efficiently accelerate galactic cosmic rays up to the highest energies, i.e. $10^{15}$ eV. A stronger magnetic field in the pre-shock region could provide the efficiency required. Bell's cosmic-ray nonresonant streaming instability has been claimed to be responsible for the amplification of precursor magnetic fields. However, an alternative mechanism has been proposed in which the cosmic-ray pressure gradient forms the shock precursor and drives turbulence, amplifying the magnetic field via the small-scale dynamo. A key ingredient for the mechanism to operate are the inhomogeneities present in the interstellar medium (ISM). These inhomogeneities are the consequence of turbulence. In this work we explore the magnetic field amplification in different ISM conditions through 3D MHD numerical simulations.

Fast radio bursts are mysterious transient sources likely located at cosmological distances. The derived brightness temperatures exceed by many orders of magnitude the self-absorption limit of incoherent synchrotron radiation, implying the operation of a coherent emission process. We propose a radiation mechanism for fast radio bursts where the emission arises from collisionless Bremsstrahlung in strong plasma turbulence excited by relativistic electron beams. We discuss possible astrophysical scenarios in which this process might operate. The emitting region is a turbulent plasma hit by a relativistic jet, where Langmuir plasma waves produce a concentration of intense electrostatic soliton-like regions (cavitons). The resulting radiation is coherent and, under some physical conditions, can be polarised and have a power-law distribution in energy. We obtain radio luminosities in agreement with the inferred values for fast radio bursts. The timescale of the radio flare in some cases can be extremely fast, of the order of $10^{-3}$ s. The mechanism we present here can explain the main features of fast radio bursts and is plausible in different astrophysical sources, such as gamma-ray bursts and some Active Galactic Nuclei.

Recent studies have indicated that cosmic ray acceleration by a first-order Fermi process in magnetic reconnection current sheets can be efficient enough in the surrounds of compact sources. In this work, we discuss this acceleration mechanism operating in the core region of galactic black hole binaries (or microquasars) and show the conditions under which this can be more efficient than shock acceleration. In addition, we compare the corresponding acceleration rate with the relevant radiative loss rates obtaining the possible energy cut-off of the accelerated particles and also compute the expected spectral energy distribution (SED) for two sources of this class, namely Cygnus X-1 and Cygnus X-3, considering both leptonic and hadronic processes. The derived SEDs are comparable to the observed ones in the low and high energy ranges. Our results suggest that hadronic non-thermal emission due to photo-meson production may produce the very high energy gamma-rays in these microquasars.

Giant molecular clouds (GMCs) are a new population of gamma-ray sources, being the target of cosmic rays (CRs) -- locally accelerated or not --. These clouds host very young stellar clusters where massive star formation takes place. Eventually, some of the stars are ejected from the clusters, becoming runaway stars. These stars move supersonically through the cloud and develop bowshocks where particles can be accelerated up to relativistic energies. As a result, the bowshocks present non-thermal emission, and inject relativistic protons in the cloud. These protons diffuse in the GMC interacting with the matter. We present a model for the non-thermal radiation generated by protons and secondary pairs accelerated in the bowshocks of massive runaways stars within young GMCs. We solve the transport equation for primary protons and secondary pairs as the stars move through the cloud. We present non-thermal emissivity maps in radio and in gamma rays as a function of time. We obtain X-ray luminosities of the order $\sim$ ${ 10^{32}}$ erg~s$^{-1}$ and gamma-ray luminosities $\sim$ $10^{34}$ erg~s$^{-1}$. We conclude that, under some assumptions, relativistic protons from massive runaway stars interacting with matter in GMCs give rise to extended non-thermal sources.

After the so-called cosmic recombination, the expanding universe entered into a period of darkness since most of the matter was in a neutral state. About a billion years later, however, the intergalactic space was once again ionized. The process, known as the cosmic reionization, required the operation of mechanisms that are not well understood. Among other ionizing sources, Population III stars, mini-quasars, and X-ray emitting microquasars have been invoked. In this article we propose that primordial cosmic rays, accelerated at the termination points of the jets of the first microquasars, may have contributed to the reionization of the intergalactic space as well. For this we quantify the ionization power of cosmic rays (electrons and protons) in the primordial intergalactic medium using extensive particle cascade simulations. We establish that, depending on the fraction of electrons to protons accelerated in the microquasar jets, cosmic rays should have contributed to the reionization of the primordial intergalactic medium as much as X-rays from microquasar accretion disks. If the primordial magnetic field was of the order of $10^{-17}$ G, as some models suggest, cosmic rays had an important role in ionizing the neutral material far beyond the birth places of the first stars.

Runaway stars are ejected from their formation sites well within molecular cores in giant dark clouds. Eventually, these stars can travel through the molecular clouds, which are highly inhomogeneous. The powerful winds of massive runaway stars interact with the medium forming bowshocks. Recent observations and theoretical modelling suggest that these bowshocks emit non-thermal radiation. As the massive stars move through the inhomogeneous ambient gas the physical properties of the bowshocks are modified, producing changes in the non-thermal emission. We aim to compute the non-thermal radiation produced in the bowshocks of runaway massive stars when travelling through a molecular cloud. We calculate the non-thermal emission and absorption for two types of massive runaway stars, an O9I and an O4I, as they move through a density gradient. We present the spectral energy distributions for the runaway stars modelled. Additionally, we obtain light curves at different energy ranges. We find significant variations in the emission over timescales of $\sim$ 1 yr. We conclude that bowshocks of massive runaway stars, under some assumptions, might be variable gamma-ray sources, with variability timescales that depend on the medium density profile. These objects might constitute a population of galactic gamma-ray sources turning on and off within years.

HD 195592 is an O-type super-giant star, known as a well-established runaway. Recently, a Fermi gamma-ray source (2FGL J2030.7+4417) with a position compatible with that of HD 195592 has been reported. Our goal is to explore the scenario where HD 195592 is the counterpart of the Fermi gamma-ray source modeling the non-thermal emission produced in the bowshock of the runaway star. We calculate the spectral energy distribution of the radiation produced in the bowshock of HD 195592 and we compare it with Fermi observations of 2FGL J2030.7+4417. We present relativistic particle losses and the resulting radiation of the bowshock of HD 195592 and show that the latter is compatible with the detected gamma-ray emission. We conclude that the Fermi source 2FGL J2030.7+4417 might be produced, under some energetic assumptions, by inverse Compton up-scattering of photons from the heated dust in the bowshock of the runaway star. HD 195592 might therefore be the very first object detected belonging to the category of gamma-ray emitting runaway massive stars, whose existence has been recently predicted.

Runaway stars produce shocks when passing through interstellar medium at supersonic velocities. Bow shocks have been detected in the mid-infrared for several high-mass runaway stars and in radio waves for one star. Theoretical models predict the production of high-energy photons by non-thermal radiative processes in a number sufficiently large to be detected in X-rays. To date, no stellar bow shock has been detected at such energies. We present the first detection of X-ray emission from a bow shock produced by a runaway star. The star is AE Aur, which was likely expelled from its birthplace by the encounter of two massive binary systems and now is passing through the dense nebula IC 405. The X-ray emission from the bow shock is detected at 30" to the northeast of the star, coinciding with an enhancement in the density of the nebula. From the analysis of the observed X-ray spectrum of the source and our theoretical emission model, we confirm that the X-ray emission is produced mainly by inverse Compton upscattering of infrared photons from dust in the shock front.

Runaway massive stars are O- and B-type stars with high spatial velocities with respect to the interstellar medium. These stars can produce bowshocks in the surrounding gas. Bowshocks develop as arc-shaped structures, with bows pointing to the same direction as the stellar velocity, while the star moves supersonically through the interstellar gas. The piled-up shocked matter emits thermal radiation and a population of locally accelerated relativistic particles is expected to produce non-thermal emission over a wide range of energies. We aim to model the non-thermal radiation produced in these sources. Under some assumptions, we computed the non-thermal emission produced by the relativistic particles and the thermal radiation caused by free-free interactions, for O4I and O9I stars. We applied our model to Zeta Oph (HD 149757), an intensively studied massive star seen from the northern hemisphere. This star has spectral type O9.5V and is a well-known runaway. Spectral energy distributions of massive runaways are predicted for the whole electromagnetic spectrum. We conclude that the non-thermal radiation might be detectable at various energy bands for relatively nearby runaway stars, especially at high-energy gamma rays. Inverse Compton scattering with photons from the heated dust gives the most important contribution to the high-energy spectrum. This emission approaches Fermi sensitivities in the case of Zeta Oph.

T Tauri stars are young, low mass, pre-main sequence stars surrounded by an accretion disk. These objects present strong magnetic activity and powerful magnetic reconnection events. Strong shocks are likely associated with fast reconnection in the stellar magnetosphere. Such shocks can accelerate particles up to relativistic energies. We aim at developing a simple model to calculate the radiation produced by non-thermal relativistic particles in the environment of T Tauri stars. We want to establish whether this emission is detectable at high energies with the available or forthcoming gamma-ray telescopes. We assume that particles (protons and electrons) pre-accelerated in reconnection events are accelerated at shocks through Fermi mechanism and we study the high-energy emission produced by the dominant radiative processes. We calculate the spectral energy distribution of T Tauri stars up to high-energies and we compare the integrated flux obtained with that from a specific Fermi source, 1FGL J1625.8-2429c, that we tentatively associate with this kind of young stellar objects (YSOs). We suggest that under reasonable general conditions nearby T Tauri stars might be detected at high energies and be responsible for some unidentified Fermi sources on the Galactic plane.

T Tauri stars are young, low mass, pre-main sequence stars surrounded by an accretion disk. These objects present strong magnetic activity and powerful magnetic reconnection events. Strong shocks are likely produced in the stellar magnetosphere by such events and charged particles accelerated up to relativistic energies. We present a simple model for the non-thermal radiation generated by high-energy particles in T Tauri stars. We discuss whether this emission is detectable at high energies with the available gamma-ray telescopes.

The high-mass microquasar Cygnus X-1, the best-established candidate for a stellar-mass black hole in the Galaxy, has been detected in a flaring state at very high energies (VHE), E > 200 GeV, by the Atmospheric Cherenkov Telescope MAGIC. The flare occurred at orbital phase 0.91, where phase 1 is the configuration with the black hole behind the companion high-mass star, when the absorption of gamma-ray photons by photon-photon annihilation with the stellar field is expected to be highest. We aim to set up a model for the high-energy emission and absorption in Cyg X-1 that can explain the nature of the observed gamma-ray flare. We study the gamma-ray opacity due to pair creation along the whole orbit, and for different locations of the emitter. Then we consider a possible mechanism for the production of the VHE emission. We present detailed calculations of the gamma-ray opacity and infer from these calculations the distance from the black hole where the emitting region was located. We suggest that the flare was the result of a jet-clump interaction where the decay products of inelastic proton-proton collisions dominate the VHE outcome. We are able to reproduce the spectrum of Cyg X-1 during the observed flare under reasonable assumptions. The flare may be the first event of jet-cloud interaction ever detected at such high energies.