Source author record

Stefano Gabici

Stefano Gabici appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2025arXiv

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Black holes (BH), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to tera-electronvolt (TeV) gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectrum extending to the ultrahigh-energy (UHE; photon energy $E>100$ TeV) band and one microquasar (Cygnus X-1) of spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory (LHAASO). Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles out of 1 peta-electronvolt (PeV), suggesting microquasars to be important contributors to Galactic cosmic rays especially around the `knee' region.

preprint2023arXiv

Star clusters as cosmic ray accelerators

Massive stars blow powerful winds and eventually explode as supernovae. By doing so, they inject energy and momentum in the circumstellar medium, which is pushed away from the star and piles up to form a dense and expanding shell of gas. The effect is larger when many massive stars are grouped together in bound clusters or associations. Large cavities form around clusters as a result of the stellar feedback on the ambient medium. They are called superbubbles and are characterised by the presence of turbulent and supersonic gas motions. This makes star clusters ideal environments for particle acceleration, and potential contributors to the observed Galactic cosmic ray intensity. The acceleration of particles at star clusters and in their surroundings may provide a major contribution to the observed CR flux. Moreover, it may explain the fine structures observed in the chemical composition of these particles, and possibly provide a solution to the puzzle of the origin of cosmic rays of energies in the PeV range and beyond.

preprint2022arXiv

Cosmic ray production in superbubbles

We compute the production of cosmic rays in the dynamical superbubble produced by a cluster of massive stars. Stellar winds, supernova remnants and turbulence are found to accelerate particles so efficiently that the nonlinear feedback of the particles must be taken into account in order to ensure that the energy balance is not violated. High energy particles do not scatter efficiently on the turbulence and escape quickly after each supernova explosion, which makes both their intensity inside the bubble and injection in the interstellar medium intermittent. On the other hand, the stochastic acceleration of low energy particles hardens the spectra at GeV energies. Because cosmic rays damp the turbulence cascade, this hardening is less pronounced when nonlinearities are taken into account. Nevertheless, spectra with hard components extending up to 1 to 10 GeV and normalised to an energy density of 1 to 100 eV cm$^{-3}$ are found to be typical signatures of cosmic rays produced in superbubbles. Efficient shock reacceleration within compact clusters is further shown to produce hard, slightly concave spectra, while the presence of a magnetised shell is shown to enhance the confinement of cosmic rays in the bubble and therefore the collective plasma effects acting on them. We eventually estimate the overall contribution of superbubbles to the galactic cosmic ray content and show typical gamma-ray spectra expected from hadronic interactions in superbuble shells. In both cases, a qualitative agreement with observations is obtained.

preprint2021arXiv

The theory of cosmic-ray scattering on pre-existing MHD modes meets data

We present a comprehensive study about the phenomenological implications of the theory describing Galactic cosmic-ray scattering onto magnetosonic and Alfvénic fluctuations in the $\mathrm{GeV} - \mathrm{PeV}$ domain. We compute a set of diffusion coefficients from first principles, for different values of the Alfvénic Mach number and other relevant parameters associated to both the Galactic halo and the extended disk, taking into account the different damping mechanisms of turbulent fluctuations acting in these environments. We confirm that the scattering rate associated to Alfvénic turbulence is highly suppressed if the anisotropy of the cascade is taken into account. On the other hand, we highlight that magnetosonic modes play a dominant role in Galactic confinement of cosmic rays up to $\mathrm{PeV}$ energies. We implement the diffusion coefficients in the numerical framework of the {\tt DRAGON} code, and simulate the equilibrium spectrum of different primary and secondary cosmic-ray species. We show that, for reasonable choices of the parameters under consideration, all primary and secondary fluxes at high energy (above a rigidity of $\simeq 200 \, \mathrm{GV}$) are correctly reproduced within our framework, in both normalization and slope.

preprint2020arXiv

Particle acceleration at colliding shock waves

We model the diffusive shock acceleration of particles in a system of two colliding shock waves and present a method to solve the time-dependent problem analytically in the test-particle approximation and high energy limit. In particular, we show that in this limit the problem can be analysed with the help of a self-similar solution. While a number of recent works predict hard ($E^{-1}$) spectra for the accelerated particles in the stationary limit, or the appearance of spectral breaks, we found instead that the spectrum of accelerated particles in a time-dependent collision follows quite closely the canonical $E^{-2}$ prediction of diffusive shock acceleration at a single shock, except at the highest energy, where a hardening appears, originating a bumpy feature just before the exponential cutoff. We also investigated the effect of the reacceleration of pre-existing cosmic rays by a system of two shocks, and found that under certain conditions spectral features can appear in the cutoff region. Finally, the mathematical methods presented here are very general and could be easily applied to a variety of astrophysical situations, including for instance standing shocks in accretion flows, diverging shocks, backward collisions of a slow shock by a faster shock, and wind-wind or shock-wind collisions.

preprint2019arXiv

Runaway O-star Bow Shocks as Particle Accelerators? The Case of AE Aur revisited

We present results of our Chandra/ACIS observations of the field centered on the fast, runaway O star AE Aur and its bow shock. Previous XMM-Newton observations revealed an X-ray "blob" near the IR arc tracing the bow shock, possibly a nonthermal source consistent with models of Inverse Compton scattering of dust IR photons by electrons accelerated at the shock. The new, subarcsecond resolution Chandra data, while confirming the presence of the XMM-Newton source, clearly indicate that the latter is neither extended nor coincident with the IR arc and strongly suggest it is a background AGN. Motivated by results published for the bow shock of BD+43 3654, we extended our study to the radio domain, by analyzing archival EVLA data. We find no radio emission from the AE Aur bow shock either. The corresponding upper limits for the absorbed (unabsorbed) X-ray flux of 5.9(7.8)x10^-15 erg/cm^2/s (3 sigma) and, in the radio range, of 2 mJy (1.4 GHz), and 0.4 mJy (5.0 GHz), are used to put constraints on model predictions for particle acceleration within the bow shock. In the "classical" framework of Diffusive Shock Acceleration, we find that the predicted X-ray and radio emission by the bow shock is at least two orders of magnitude below the current upper limits, consistent with the systematic non-detections of up to 60 stellar bow shocks. The only exception so far remains that of BD+43 3654, probably the result of its very large mass-loss rate among runaway O stars.

preprint2016arXiv

Non-linear diffusion of cosmic rays escaping from supernova remnants - I. The effect of neutrals

Supernova remnants are believed to be the main sources of galactic Cosmic Rays (CR). Within this framework, particles are accelerated at supernova remnant shocks and then released in the interstellar medium. The mechanism through which CRs are released and the way in which they propagate still remain open issues. The main difficulty is the high non-linearity of the problem: CRs themselves excite the magnetic turbulence that confines them close to their sources. We solve numerically the coupled differential equations describing the evolution in space and time of the escaping particles and of the waves generated through the CR streaming instability. The warm ionized and warm neutral phases of the interstellar medium are considered. These phases occupy the largest fraction of the disc volume, where most supernovae explode, and are characterised by the significant presence of neutral particles. The friction between those neutrals and ions results in a very effective wave damping mechanism. It is found that streaming instability affects the propagation of CRs even in the presence of ion-neutral friction. The diffusion coefficient can be suppressed by more than a factor of $\sim 2$ over a region of few tens of pc around the remnant. The suppression increases for smaller distances. The propagation of $\approx 10$ GeV particles is affected for several tens of kiloyears after escape, while $\approx 1$ TeV particles are affected for few kiloyears. This might have a great impact on the interpretation of gamma-ray observations of molecular clouds located in the vicinity of supernova remnants.

preprint2015arXiv

CRIME - cosmic ray interactions in molecular environments

Molecular clouds act as targets for cosmic rays (CR), revealing their presence through either gamma-ray emission due to proton-proton interactions, and/or through the ionization level in the cloud, produced by the CR flux. The ionization rate is a unique tool, to some extent complementary to the gamma-ray emission, in that it allows to constrain the CR spectrum especially for energies below the pion production rate ($\approx 280$ MeV). Here we study the effect of ionization on $H_2$ clouds due to both CR protons and electrons, using the fully relativistic ionization cross sections, which is important to correctly account for the contribution due to relativistic CRs. The contribution to ionization due to secondary electrons is also included self-consistently. The whole calculation has been implemented into a numerical code which is publicly accessible through a web-interface. The code also include the calculation of gamma-ray emission once the CR spectrum

preprint2015arXiv

Gamma-ray observations of blazars and the intergalactic magnetic field spectrum

Very-high energy observations of blazars can be used to constrain the strength of the intergalactic magnetic field. A simplifying assumption which is often made is that of a magnetic field of constant strength composed by randomly oriented and identical cells. In this paper, we demonstrate that a more realistic description of the structure of the intergalactic magnetic field is indeed needed. If such a description is adopted, the observational bounds on the field strength are significantly affected in the limit of short field correlation lengths: in particular, they acquire a dependence on the magnetic field power spectrum. In the case of intergalactic magnetic fields which are generated causally, for which the magnetic field large scale spectral index is $n_B\geq 2$ and even, the observational lower bound becomes more constraining by about a factor 3. If instead $-3<n_B<-2$, the lower bound is significantly relaxed. Such magnetic fields with very red spectra can in principle be produced during inflation, but remain up to now speculative.

preprint2015arXiv

High-energy gamma-ray and neutrino backgrounds from clusters of galaxies and radio constraints

Cosmic-ray protons accumulate for cosmological times in clusters of galaxies as their typical radiative and diffusive escape times are longer than the Hubble time. Their hadronic interactions with protons of the intra-cluster medium generate secondary electrons, gamma-rays and neutrinos. We here estimate the contribution from clusters to the diffuse gamma-ray and neutrino backgrounds. We model the cluster population by means of their mass function, using a phenomenological luminosity-mass relation applied to all clusters, as well as a detailed semi-analytical model. Additionally, we consider observationally-motivated values for the cluster magnetic field. This is a crucial parameter since the observed radio counts due to synchrotron emission by secondary electrons need to be respected. For a choice of parameters respecting all current constraints, and assuming a spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the extragalactic gamma-ray background observed by Fermi and to the IceCube flux. They account for less than 1% for spectral indices <=-2. The neutrino flux observed by IceCube can be reproduced without violating constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with energy. In the case of proton-photon interactions, we find that very likely protons do not reach sufficiently high energies to produce neutrinos in clusters. We argue that our results are optimistic due to our assumptions, and that clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds. Finally, we find that the cluster contribution to the angular fluctuations in the gamma-ray background is subdominant, less than 10%. [abridged]

preprint2015arXiv

Neutrinos from Clusters of Galaxies and Radio Constraints

Cosmic-ray (CR) protons can accumulate for cosmological times in clusters of galaxies. Their hadronic interactions with protons of the intra-cluster medium (ICM) generate secondary electrons, gamma-rays and high-energy neutrinos. In light of the high-energy neutrino events recently discovered by the IceCube observatory, we estimate the contribution from galaxy clusters to the diffuse gamma-ray and neutrino backgrounds. For the first time, we consistently take into account the synchrotron emission generated by secondary electrons and require the clusters radio counts to be respected. For a choice of parameters respecting current constraints from radio to gamma-rays, and assuming a proton spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the IceCube flux, and much less to the total extragalactic gamma-ray background observed by Fermi. They account for less than 1% for spectral indexes <-2. The high-energy neutrino flux observed by IceCube can be reproduced without violating radio constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the high-energy IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with the particle energy. We stress that our results are valid for all kind of sources injecting CR protons into the ICM, and that, while IceCube can test the most optimistic scenarios for spectral indexes >=-2.2 by stacking few nearby massive objects, clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds in any realistic scenario.

preprint2015arXiv

Non-linear Cosmic Ray propagation close to the acceleration site

Recent advances on gamma-ray observations from SuperNova Remnants and Molecular Clouds offer the possibility to study in detail the properties of the propagation of escaping Cosmic Rays (CR). However, a complete theory for CR transport outside the acceleration site has not been developed yet. Two physical processes are thought to be relevant to regulate the transport: the growth of waves caused by streaming instability, and possible wave damping mechanisms that reduce the growth of the turbulence. Only a few attempts have been made so far to incorporate these mechanisms in the theory of CR diffusion. In this work we present recent advances in this subject. In particular, we show results obtained by solving the coupled equations for the diffusion of CRs and the evolution of Alfven waves. We discuss the importance of streaming instabilities and wave damping in different ISM phases.

preprint2015arXiv

Observations of supernova remnants and molecular clouds from the mm to the gamma-ray domain: bridging low and high energy cosmic rays

New evidence that cosmic rays (hadronic component) are accelerated by supernova remnant shocks all the way from low energies to high energies, has come from recent works combining gamma-ray observations in the sub-GeV to TeV domain on the one hand, and in the submm-mm domain on the other hand. These observations concern the specific cases of supernova remnants interacting with molecular cloud complexes, that have long been suspected to be ideal laboratories to study in situ cosmic ray acceleration and diffusion. Indeed, enhanced gamma-ray emission from neutral pion decay, as well as enhanced ionization (both by at least one order of magnitude with respect to average galactic values) have been observed in several regions of massive star formation housing supernova remnants interacting with molecular cloud complexes. This paper summarizes the main physical and chemical processes at work, as well as recent observations, that further support the paradigm of cosmic ray acceleration by supernova remnants all the way from the MeV domain up to several tens of TeV, although much work remains to be done to understand cosmic ray penetration and diffusion inside and around molecular clouds, and reveal the actual spectrum of the accelerated cosmic rays.

preprint2015arXiv

On the electron-ion temperature ratio established by collisionless shocks

Astrophysical shocks are often collisionless shocks. An open question about collisionless shocks is whether electrons and ions each establish their own post-shock temperature, or whether they quickly equilibrate in the shock region. Here we provide simple relations for the minimal amount of equilibration to expect. The basic assumption is that the enthalpy-flux of the electrons is conserved separately, but that all particle species should undergo the same density jump across the the shock. This assumption results in an analytic treatment of electron-ion equilibration that agrees with observations of collisionless shocks: at low Mach numbers ($<2$) the electrons and ions are close to equilibration, whereas for Mach numbers above $M \sim 60$ the electron-ion temperature ratio scales with the particle masses $T_e/T_i = m_e/m_i$. In between these two extremes the electron-ion temperature ratio scales as $T_e/T_i \propto 1/M_s^2$. This relation also hold if adiabatic compression of the electrons is taken into account. For magnetised plasmas the compression is governed by the magnetosonic Mach number, whereas the electron-ion temperatures are governed by the sonic Mach number. The derived equations are in agreement with observational data at low Mach numbers, but for supernova remnants the relation requires that the inferred Mach numbers for the observations are over- estimated, perhaps as a result of upstream heating in the cosmic-ray precursor. In addition to predicting a minimal electron/ion temperature ratio, we also heuristically incorporate ion-electron heat exchange at the shock, quantified with a dimensionless parameter $ξ$. Comparing the model to existing observations in the solar system and supernova remnants suggests that the data are best described by $ξ \sim 5$ percent. (Abridged abstract.)

preprint2014arXiv

A Galactic Halo Origin of the Neutrinos Detected by IceCube

Recent IceCube results suggest that the first detection of very high energy astrophysical neutrinos have been accomplished. We consider these results at face value in a Galactic origin context. Emission scenarios from both the Fermi bubble and broader halo region are considered. We motivate that such an intensity of diffuse neutrino emission could be Galactic in origin if it is produced from an outflow into the halo region. This scenario requires cosmic ray transport within the outflow environment to be different to that inferred locally within the disk and that activity in the central part of the Galaxy accelerates cosmic rays to trans-"knee" energies before they escape into an outflow. The presence of a large reservoir of gas in a very extended halo around the Galaxy, recently inferred from X-ray observations, implies that relatively modest acceleration power of $10^{39}$ erg s$^{-1}$ in PeV energy cosmic rays may be sufficient to explain the observed neutrino flux. Such a luminosity is compatible with that required to explain the observed intensity of cosmic rays around the "knee".

preprint2014arXiv

Cosmic ray induced ionisation of a molecular cloud shocked by the W28 supernova remnant

Cosmic rays are an essential ingredient in the evolution of the interstellar medium, as they dominate the ionisation of the dense molecular gas, where stars and planets form. However, since they are efficiently scattered by the galactic magnetic fields, many questions remain open, such as where exactly they are accelerated, what is their original energy spectrum, and how they propagate into molecular clouds. In this work we present new observations and discuss in detail a method that allows us to measure the cosmic ray ionisation rate towards the molecular clouds close to the W28 supernova remnant. To perform these measurements, we use CO, HCO$^+$, and DCO$^+$ millimetre line observations and compare them with the predictions of radiative transfer and chemical models away from thermodynamical equilibrium. The CO observations allow us to constrain the density, temperature, and column density towards each observed position, while the DCO$^+$/HCO$^+$ abundance ratios provide us with constraints on the electron fraction and, consequently, on the cosmic ray ionisation rate. Towards positions located close to the supernova remnant, we find cosmic ray ionisation rates much larger ($\gtrsim100$) than those in standard galactic clouds. Conversely, towards one position situated at a larger distance, we derive a standard cosmic ray ionisation rate. Overall, these observations support the hypothesis that the $γ$ rays observed in the region have a hadronic origin. In addition, based on CR diffusion estimates, we find that the ionisation of the gas is likely due to $0.1 -1$ GeV cosmic rays. Finally, these observations are also in agreement with the global picture of cosmic ray diffusion, in which the low-energy tail of the cosmic ray population diffuses at smaller distances than the high-energy counterpart.

preprint2014arXiv

Hadronic gamma-rays from RX J1713.7-3946?

RXJ1713.7-3946 is a key object to check the supernova remnant paradigm of the origin of Galactic cosmic rays. While the origin of its gamma-ray emission (hadronic versus leptonic) is still debated, the hard spectrum at GeV energies reported by the Fermi collaboration is generally interpreted as a strong argument in favor of a leptonic scenario. On the contrary, we show that hadronic interactions can naturally explain the gamma-ray spectrum if gas clumps are present in the supernova remnant shell. The absence of thermal X-rays from the remnant fits well within this scenario.

preprint2014arXiv

Propagation of cosmic rays into diffuse clouds

We study the capability of low-energy cosmic rays (CR) to penetrate into diffuse clouds when they move from the hot ionized plasma to a cool cloud embedded in that plasma. The spectrum of CR inside a cloud can be remarkably different from the the one present in the hot interstellar medium because when CRs pass through a dense cloud of matter, they suffer energy losses due to ionization and nuclear interactions. Hence there is a net flux of CRs towards the cloud that can excite Alfvén waves. In turn, self-excited Alfvén waves enhances the diffusion of CRs near the edge of the cloud, forcing CRs to spend more time in this layer and increasing the amount of energy losses. The final effect is that the flux of CR entering into the cloud is strongly suppressed below an energy threshold whose value depends on ambient parameters. For the first time we use the full kinetic theory to describe this problem, coupling CRs and Alfvén waves through the streaming instability, and including the damping of the waves due to ion-neutral friction and the CR energy losses due to ionization and pion production. Differently from previous approaches, this method allow us to describe the full CR spectrum and the magnetic turbulence spectrum in the transition region between hot interstellar medium and the cloud. For the typical size and density of diffuse clouds in the Galaxy, we find that the flux of particles entering the cloud is strongly reduced below E~10-100 MeV, while for larger energies the CR spectrum remain unaltered.

preprint2013arXiv

Interaction of escaping cosmic rays with molecular clouds

The study of the gamma-ray radiation produced by cosmic rays that escape their accelerators is of paramount importance for (at least) two reasons: first, the detection of those gamma-ray photons can serve to identify the sources of cosmic rays and, second, the characteristics of that radiation give us constraints on the way in which cosmic rays propagate in the interstellar medium. This paper reviews the present status of the field.

preprint2012arXiv

Anisotropic CR diffusion and gamma-ray production close to supernova remnants, with an application to W28

Cosmic rays that escape their acceleration site interact with the ambient medium and produce gamma rays as the result of inelastic proton-proton collisions. The detection of such diffuse emission may reveal the presence of an accelerator of cosmic rays, and also constrain the cosmic ray diffusion coefficient in its vicinity. Preliminary results in this direction have been obtained in the last years from studies of the gamma-ray emission from molecular clouds located in the vicinity of supernova remnants, which are the prime candidate for cosmic ray production. Hints have been found for a significant suppression of the diffusion coefficient with respect to the average one in the Galaxy. However, most of these studies rely on the assumption of isotropic diffusion, which may not be very well justified. Here, we extend this study to the case in which cosmic rays that escape an accelerator diffuse preferentially along the magnetic field lines. As a first approximation, we further assume that particles are strongly magnetized and that their transport perpendicular to the magnetic field is mainly due to the wandering of the field lines. The resulting spatial distribution of runaway cosmic rays around the accelerator is, in this case, strongly anisotropic. An application of the model to the case of the supernova remnant W28 demonstrates how the estimates of the diffusion coefficient from gamma-ray observations strongly depend on the assumptions made on the isotropy (or anisotropy) of diffusion. For higher levels of anisotropy of the diffusion, larger values of the diffusion coefficient are found to provide a good fit to data. Thus, detailed models for the propagation of cosmic rays are needed in order to interpret in a correct way the gamma-ray observations.

preprint2012arXiv

Cosmic rays and molecular clouds

This paper deals with the cosmic-ray penetration into molecular clouds and with the related gamma--ray emission. High energy cosmic rays interact with the dense gas and produce neutral pions which in turn decay into two gamma rays. This makes molecular clouds potential sources of gamma rays, especially if they are located in the vicinity of a powerful accelerator that injects cosmic rays in the interstellar medium. The amplitude and duration in time of the cosmic--ray overdensity around a given source depend on how quickly cosmic rays diffuse in the turbulent galactic magnetic field. For these reasons, gamma-ray observations of molecular clouds can be used both to locate the sources of cosmic rays and to constrain the properties of cosmic-ray diffusion in the Galaxy.

preprint2011arXiv

Cosmic ray escape from supernova remnants

Galactic cosmic rays are believed to be accelerated at supernova remnants via diffusive shock acceleration. Though this mechanism gives fairly robust predictions for the spectrum of particles accelerated at the shock, the spectrum of the cosmic rays which are eventually injected in the interstellar medium is more uncertain and depends on the details of the process of particle escape from the shock. Knowing the spectral shape of these escaping particles is of crucial importance in order to assess the validity of the supernova remnant paradigm for cosmic ray origin. Moreover, after escaping from a supernova remnant, cosmic rays interact with the surrounding ambient gas and produce gamma rays in the vicinity of the remnant itself. The detection of this radiation can be used as an indirect proof of the fact that the supernova remnant was indeed accelerating cosmic rays in the past.

preprint2011arXiv

Molecular Clouds as Cosmic Ray Laboratories

We will here discuss how the gamma-ray emission from molecular clouds can be used to probe the cosmic ray flux in distant regions of the Galaxy and to constrain the highly unknown cosmic ray diffusion coefficient. In particular we will discuss the GeV to TeV emission from runaway cosmic rays penetrating molecular clouds close to young and old supernova remnants and in molecular clouds illuminated by the background cosmic ray flux.

preprint2009arXiv

Broad-band nonthermal emission from molecular clouds illuminated by cosmic rays from nearby supernova remnants

Molecular clouds are expected to emit non-thermal radiation due to cosmic ray interactions in the dense magnetized gas. Such emission is amplified if a cloud is located close to an accelerator of cosmic rays and if energetic particles can leave the accelerator site and diffusively reach the cloud. We consider here the situation in which a molecular cloud is located in the proximity of a supernova remnant which is efficiently accelerating cosmic rays and gradually releasing them in the interstellar medium. We calculate the multiwavelength spectrum from radio to gamma rays which is emerging from the cloud as the result of cosmic ray interactions. The total energy output is dominated by the gamma ray emission, which can exceed the emission in other bands by an order of magnitude or more. This suggests that some of the unidentified TeV sources detected so far, with no obvious or very weak counterparts in other wavelengths, might be in fact associated with clouds illuminated by cosmic rays coming from a nearby source. Moreover, under certain conditions, the gamma ray spectrum exhibit a concave shape, being steep at low energies and hard at high energies. This fact might have important implications for the studies of the spectral compatibility of GeV and TeV gamma ray sources.

preprint2007arXiv

Particle Acceleration in Supernova Remnants and the Production of Thermal and Nonthermal Radiation

If highly efficient, cosmic ray production can have a significant effect on the X-ray emission from SNRs as well as their dynamical evolution. Using hydrodynamical simulations including diffusive shock acceleration, we produce spectra for both the thermal and nonthermal forward shock emission. For a given ambient density and explosion energy, we find that the position of the forward shock at a given age is a strong function of the acceleration efficiency, providing a signature of cosmic-ray production. Using an approximate treatment for the ionization state of the plasma, we investigate the effects of slow vs. rapid heating of the postshock electrons on the ratio of thermal to nonthermal X-ray emission at the forward shock. We also investigate the effects of magnetic field strength on the observed spectrum for efficient cosmic-ray acceleration. The primary effect of a large field is a considerable flattening of the nonthermal spectrum in the soft X-ray band. Spectral index measurements from X-ray observations may thus be indicators of the postshock magnetic field strength. The predicted gamma-ray flux from inverse-Compton (IC) scattering and neutral pion decay is strongly affected by the ambient conditions and, for the particular parameters used in our examples, the IC emission at E ~ 1 TeV exceeds that from pion decay, although at both lower and higher energies this trend is reversed for cases of high ambient density. More importantly, high magnetic fields produce a steepening of the electron spectrum over a wide energy range which may make it more difficult to differentiate between IC and pion-decay emission solely by spectral shape.

preprint2004arXiv

Alfvenic reacceleration of relativistic particles in galaxy clusters: MHD waves, leptons and hadrons

There is a growing evidence that extended radio halos are most likely generated by electrons reaccelerated via some kind of turbulence generated in the cluster volume during major mergers. It is well known that Alfvén waves channel most of their energy flux in the acceleration of relativistic particles. Much work has been done recently to study this phenomenon and its consequences for the explanation of the observed non-thermal phenomena in clusters of galaxies. We investigate here the problem of particle-wave interactions in the most general situation in which relativistic electrons, thermal protons and relativistic protons exist within the cluster volume. The interaction of all these components with the waves, as well as the turbulent cascading and damping processes of Alfvén waves, are treated in a fully time-dependent way. This allows us to calculate the spectra of electrons, protons and waves at any fixed time. The {\it Lighthill} mechanism is invoked to couple the fluid turbulence, supposedly injected during cluster mergers, to MHD turbulence. We find that present observations of non-thermal radiation from clusters of galaxies are well described within this approach, provided the fraction of relativistic hadrons in the intracluster medium (ICM) is smaller than 5-10 %.