Radiation from the impact of broad-line region clouds onto AGN accretion disks
Ana Laura Müller Gustavo E. Romero
astro-ph.HE
 


   

   

   

   Instituto Argentino de Radioastronomía (CCT-La Plata, CONICET; CICPBA), 
               C.C. No. 5, 1894, Villa Elisa, Argentina

              almuller@iar.unlp.edu.ar
         Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, 
             Paseo del Bosque, 1900, La Plata, Argentina

             romero@iar.unlp.edu.ar
         Institute for Nuclear Physics (IKP), Karlsruhe Institute of Technology (KIT), Germany
         Instituto de Tecnologías en Detección y Astropartículas (CNEA, CONICET, UNSAM), Buenos Aires, Argentina

             
         
   
   Radiation from BLR clouds onto AGN accretion disks
   A. L. Müller and G. E. Romero
   
         

   





 
  
  
   Active galactic nuclei are supermassive black holes surrounded by an accretion disk, two populations of clouds, bipolar jets, and a dusty torus. The clouds move in Keplerian orbits at high velocities. In particular, the broad-line region (BLR) clouds have velocities ranging from 1000 to 10000 km s^-1. Given the extreme proximity of these clouds to the supermassive black hole, frequent collisions with the accretion disk should occur.
  
   The impact of BLR clouds onto the accretion disk can produce strong shock waves where particles might be accelerated. The goal of this work is to investigate the production of relativistic particles, and the associated non-thermal radiation in these events. In particular, we apply the model we develop to the Seyfert galaxy NGC 1068.
  
   We analyze the efficiency of diffusive shock acceleration in the shock of colliding clouds of the BLR with the accretion disk. We calculate the spectral energy distribution of photons generated by the relativistic particles and estimate the number of simultaneous impacts needed to explain the gamma radiation observed by Fermi in Seyfert galaxies.
  
   We find that is possible to understand the measured gamma emission in terms of the interaction of clouds with the disk if the hard X-ray emission of the source is at least obscured between 20% and 40%. The total number of clouds contained in the BLR region might be between 3×10^8 and 6×10^8, which are values in good agreement with the observational evidence. The maximum energy achieved by the protons (∼ PeV) in this context allows the production of neutrinos in the observing range of IceCube.
  
   

   

Radiation from the impact of broad-line region clouds onto AGN accretion disks
    A. L. Müller
          1,3,4
          G. E. Romero1,2 
          
          
    December 30, 2023
=================================================================================






§ INTRODUCTION


Active galactic nuclei (AGNs) are formed by accreting supermassive black holes (BHs). Their characteristic emission is produced by a very compact region and covers a wide range of frequencies <cit.>. From an observational point of view, objects defined as AGNs are actually very diverse. A first distinction can be made between radio-quiet and radio-loud AGNs. The latter are bright radio objects whose radiation in that band is several orders of magnitude larger than the typical emission of the radio-quiet nuclei. These two main groups can be additionally divided taking into account a variety of characteristics (e.g., the alignment of the jet with the line of sight, the intensity of the lines in the spectra, their luminosity; see <cit.>). The heterogeneity of AGNs can be understood in terms of a unified model adjusting the orientation to the observer and the values of parameters related to the central BH <cit.>. In unified pictures, an AGN is essentially a supermassive BH surrounded by a subparsec accretion disk and a dusty torus. Inside the torus two populations of clouds move in Keplerian orbits: the broad-line region (BLR) and the narrow-line (NLR) region clouds (see ). In the case of radio-loud AGNs, the system also includes  a relativistic jet emitting  synchrotron radiation.

The ultraviolet (UV) and optical spectra of some subclasses of AGNs have prominent broad emission lines (e.g., Seyfert 1). The gas producing these lines should be contained in a central region close to the BH. The structure of this zone is modeled as a group of clouds orbiting in random directions, but with velocities in the range of  to  <cit.>. The electron number density of the BLR clouds ranges typically from  to  and the gas is completely photoionized by the disk radiation. The BLR reprocesses around 10% of the disk luminosity and re-emits lines with a mean energy of  and a typical photon density of  <cit.>. 

The nucleus of Seyfert 2 galaxies is typically obscured by the dusty torus. Therefore, the BLR appears partially hidden but still detectable in the spectropolarimetric data <cit.>. Only in the low-luminosity Seyfert 2 AGNs the existence of the BLR has not been confirmed by observational data <cit.>.

Since in the standard AGN model the BLR clouds co-exist with the accretion disk, and given the strong evidence of infall motion <cit.>, direct collisions between these clouds and the disk should occur. Similarly, the interaction of stars and BHs with the accretion disk has been analyzed before by many authors, but with emphasis on the AGN fuelling consequences, the thermal emission, and the gravitational waves produced in the impacts <cit.>. 

In this work, we study the possibility of accelerating particles by first-order Fermi mechanism in the shock waves produced by the impacts of BLR clouds with the accretion disk (Section <ref>). In Sections <ref>, <ref>, and <ref>, we present estimates of the cosmic ray acceleration inside the shocked cloud and model the non-thermal emission. Finally, we apply our model to the Seyfert galaxy NGC 1068 (Section <ref>) and discuss the contribution of these impacts to the high-energy radiation in Section <ref>.

 




§ BASIC MODEL


  We assume a standard AGN with a central Schwarzschild BH of mass  surrounded by a Shakura-Sunyaev accretion disk <cit.>. The disk extends from the last stable orbit[R_ ISCO=6 R_ g∼ 3 × 10^-5 pc] to hundreds of thousands of gravitational radii <cit.>. We adopt standard values for the accretion efficiency and viscosity parameters, namely η_ accre=0.1 and α_ accre∼ 0.1, respectively <cit.>. The bolometric luminosity is ∼ 7×10^45 erg s^-1, i.e., λ_ Edd=L_ bol/L_ Edd=0.7. We calculate the characteristic values of the different parameters of the disk at each radius using the expressions provided by <cit.>. The spectrum of the accretion disk is obtained integrating the Planck function over the surface area. Each disk ring has a characteristic temperature. From the resulting expression, the total luminosity L_ disk is calculated integrating the spectrum over the whole energy range of the emission.
 
  The clouds in the BLR are considered to be spherical and homogeneous with a radius  <cit.>. The parameters adopted for an average cloud in our model are shown in Table <ref>. The evidence obtained from many observational studies indicates that the clouds existing in the BLR move in Keplerian orbits with velocities between 10^3 and  <cit.>. We adopt a scenario where the cloud velocity is . This speed corresponds to a circular Keplerian orbit of radius , which is the distance from the galactic center to the place where the impact of the cloud on the equatorial disk occurs. The relevant physical properties of the disk at that radius are shown in Table <ref>.
  
  
    
  
  
  The size of the BLR ranges typically from 0.01 to  <cit.>. One way to estimate this quantity is through reverberation studies <cit.>. The values obtained can differ by about an order of magnitude using different emission lines <cit.>. Therefore, it is necessary to account for a wide range of radii (e.g., from 10^-3 to 0.1 pc) to reproduce the line pattern attributed to a BLR <cit.>. We  assume that the BLR is a thin shell, whose internal radius is ∼ r, whereas the external one is given by R_ BLR=√(0.1 L_ disk/(π U_ BLR c)) <cit.>. 
  
     

  The cloud moves supersonically. The collision of the cloud with the accretion disk produces two shock waves: a forward shock propagating through the disk and a reverse shock propagating through the cloud. The velocities of the shocks are calculated with the expressions presented in <cit.>, whereas the values of the physical parameters in the shocked regions are obtained using the equations for strong adiabatic shocks deduced from the Rankine-Hugoniot relations <cit.>. Similar collisions between high-velocity clouds and galactic disks have been studied by several authors <cit.>. 
  
  Adiabatic shocks can be defined demanding that their cooling length R_Λ is greater than the length of the traversed medium (i.e., the cloud radius and the width of the disk). We calculate the cooling length using the following expression <cit.>:
  
  
    R_Λ=5.06×10^-29(U/km s^-1)^3 A/(n/ 
    cm^-3) (Λ(T)/erg cm^3 s^-1) pc

  with
   
    T=22 A ( U/km s^-1)^2 K,

  
  where U is the shock velocity with respect to the undisturbed medium of density n and A is a parameter which depends on the conditions of the unshocked gas; its value is 0.6 if the medium is ionized or 1.6 if it is neutral. In addition, the function Λ(T) is the cooling rate <cit.>. 

  The gas in the acceleration region should not be magnetically dominated, otherwise the medium becomes mechanically incompressible and a strong shock cannot exist <cit.>. Consequently, the magnetic energy of the medium (u_ B) must be in subequipartition with the kinetic energy of the shocked gas (u_ g), and the magnetization parameter β=u_ B/u_ g becomes less than 1. Taking this into account, we assume a modest value of 0.1 for β in order to grant effective shock generation and derive the magnetic field in the cloud from

 
  
    u_ B=β u_ g


  
    B^2/8π=0.1 9/8n_ c m_ p U^2.


  Table <ref> shows the values of the physical parameters for the shocked media.
  
    
  
  Since the shock moving through the disk turns out to be radiative, it is not efficient enough to accelerate particles. Therefore, we study the cosmic ray production only in the shock that propagates through the cloud. The collision ends after , when the shock finally reaches the total length of the cloud. After this time, hydrodynamic instabilities may become important and destroy the cloud <cit.>. In the case of magnetized clouds, it is possible for them to survive up to ∼ 4 t_ coll or even longer <cit.>.  

     


§ PARTICLE ACCELERATION AND ENERGY LOSSES


   First-order Fermi mechanisms can operate in scenarios with strong adiabatic shock waves <cit.>. The acceleration rate for a particle of energy E and charge q in a region with a magnetic field B and where diffusive shock acceleration (DSA) takes place is 
   
   
    t_ acc=E( dE/dt)^-1=η^-1E/q c B.


   Here  is the efficiency of the process. Under conditions of the  first-order Fermi mechanism <cit.>
   
   
    η^-1∼ 20 D/r_ g c( c/V_ s)^2,


   where D is the diffusion coefficient of the medium and  is the gyroradius of the particle. We assume that the diffusion proceeds in the Bohm regime, which means that .
   
   
    Given that the acceleration can be suppressed in very high-density media, it is necessary to evaluate the importance of the Coulomb and ionization losses suffered by the particles <cit.>. In order to evaluate this, we calculate the corresponding cooling times using the expressions provided by <cit.>

      
    t^e_ ion=1.3×10^8 (n_ e/ cm^-3)^-1 (E/ eV) [ ln(E/ eV/n_e/ cm^-3)+61.15 ]^-1 s

     
      
    t^p_ ion=    3.2×10^6 Z^-2 (n_ e/ cm^-3)^-1 (E/ eV)
       ×(β^2/2.34×10^-5 x_m^3+β^3)^-1 Θ(β-7.4×10^-4 x_m) s,

      where Θ is the Heaviside function,  β=√(1-γ^-2) (with γ the Lorentz factor of the particle), and x_m=(T_e/2×10^6  K)^1/2.
    

    The relativistic particles injected lose energy due to the interaction with the matter, photon, and magnetic fields of the cloud. We consider the synchrotron losses (sync) produced by the interaction of the electrons with the magnetic field and the relativistic Bremsstrahlung losses (BS) produced by the interaction of the same particles with the ionized hot matter of the cloud. We also calculate the inverse Compton (IC) upscattering of the photons from the BLR, the accretion disk, and the synchrotron radiation (SSC). The local emission from the disk is approximately a blackbody, whose temperature is T_ disk=1970.7 K. On the other hand, the BLR radiation is a monochromatic photon field with  and . For the protons, the most relevant radiative process is the proton-proton inelastic collisions (pp). We calculate the cooling timescales associated with these processes, using the expressions presented by <cit.> (see Eqs. 5-12). 
     
      We also take into account the fact that the particles can escape from the region of acceleration because of diffusion. The cooling rate for this process is <cit.>
      
      
    t_ diff=X^2/D,

      
      where D is the diffusion coefficient of the medium (i.e., the Bohm diffusion coefficient in our model) and X is the characteristic size of the acceleration region. We assume X=R_ c.
  
      Another non-radiative process that we include is the adiabatic loss, i.e., the energy loss due to the work done by the particles expanding the shocked cloud matter. The cooling timescale is given by <cit.>
      
    t_ adi=4.9 × 10^-8(U/10^3 km s^-1)^-1 (R_ c/ cm) s.

      
      The maximum energy that the particles can reach before they escape from the acceleration region is constrained by the Hillas criterion  <cit.>. The maximum value according to this criterion is .
      
      In Fig. <ref> and Fig. <ref> we show the cooling timescales together with the diffusion and acceleration rates for the electrons and protons. The maximum energy for each kind of particle can be inferred looking at the point where the acceleration rate is equal to the total cooling and/or escape rate. 
      
      Synchrotron dominates the energy losses for the electrons over the whole energy range, whereas the IC losses become negligible (see Fig. <ref>). This is expected because the magnetic energy density in the cloud  is much higher than the blackbody radiation energy density  of the disk and the photon density of the BLR . For protons, the pp dominates the energy losses (see Fig. <ref>). Consequently, the maximum energies are  and  for electrons and protons, respectively. The Hillas criterion is satisfied by particles of such energies, and thus   are the maximum energies that the particles can reach in our scenario. 
      
      Another important result shown by these plots is that, after the end of the collision, the produced cosmic rays  cool down locally and do not propagate. Electrons  lose all their energy almost immediately (), whereas the protons will lose it after ∼ 3×10^4 s. This timescale is comparable to t_ coll, thus the accelerated hadrons, and the secondary particles created by them, will emit longer than the primary leptons. 
      
       



§ PARTICLE DISTRIBUTIONS


We solve the transport equation for relativistic particles,


    ∂ N_e,p(E, t)/∂ t+∂[b(E) N_e,p(E,t)]/∂ E+ N_e,p(E, t)/t_esc=Q_e,p(E) 
    ,


to find the particle distributions <cit.>. In this expression, Q_e,p(E) represents the injection term,  the sum of all the energy losses, and t_esc the escape time of the particles (i.e., the diffusion). For the injection function, we assume a power law with an exponential cutoff Q_e,p(E)=K_e,p E^-2 exp(-E/E^e/p_ max), which corresponds to the injection produced by a DSA in adiabatic strong shocks. Given that the particle distribution of electrons reaches the steady state at once and the proton distribution does the same in only ∼ 10^4 s, we study the steady state solution of the transport equation. The time interval where this solution is valid is t_ ss=t_ coll-10^4 s=2.4×10^4 s.

  The power injected per impact can be calculated as  <cit.>. The kinetic energy obtained using the set of parameters of our model is . We assume that 10% of the energy available is used to accelerate particles up to relativistic energies. Therefore, the power available to accelerate electrons and protons in the cloud is . How this luminosity is divided between the electrons (L_e) and protons (L_p) is uncertain. We consider two situations: energy equally distributed between the two particle types  and 100 times the energy injected in electrons to protons .
      


§ SPECTRAL ENERGY DISTRIBUTIONS


Using the particle energy distributions obtained in the previous section (Section <ref>), we calculate the spectral energy distribution (SED) taking into account all the radiative processes mentioned and correcting the result by absorption. To this end, we suppose that the emission region is a spherical cap with height X, so its volume is V_ c=π X^2 (3 R_ c-X)/3.



 §.§ Radiative processes


In the case of the synchrotron emission, we use the expressions provided by <cit.>. Then the synchrotron luminosity emitted by a distribution of electrons N_ e(E_ e) can be calculated as


    L_ sync(E_γ)=E_γ V_ c κ_ SSA(E_γ)∫_E^e_ min^E^e_ max N_ e(E) P_ sync(E,E_γ) dE,


with


    P_ sync(E_ e, E_γ)=√(3) e^3 B/h m_ e c^2 E_γ/E_ c∫_E_γ/E_ c^∞ K_5/3(ζ) dζ


and


    E_ c=3/4 πe h B/m_e c(E/m_e c^2)^2.


Here, K_5/3(ζ) is a modified Bessel function. Defining , we use that φ ∫_φ^∞K_5/3(ζ) dζ≈1.85 φ^1/3e^-φ if 0.1≤φ≤ 10 <cit.>. The coefficient κ_ SSA(E_γ) is the correction due to synchrotron self-absorption (SSA):


    κ_ SSA(E_γ)=1-e^-τ_ SSA(E_γ)/τ_ SSA(E_γ).


The expression for the optical depth τ_ SSA can be found in <cit.>. 

We calculate the IC emission and Bremsstrahlung[There is a typo in Eq. 31. The 4π in the denominator should be removed.] using the expressions presented by <cit.> (Eqs. 28-33). To estimate the gamma luminosity generated by pp inelastic collisions, we follow the procedure given by <cit.> (see Section IV and V). Following this approach, the emissivity produced by protons with E_p < 100 GeV is obtained using the δ-functional approximation <cit.>, whereas for E_p > 100 GeV corrections accounting for the production of charge particles are introduced. Finally, we also calculate and include the thermal contribution from the accretion disk. 




 §.§ Absorption and secondary particles


The interaction of the gamma photons generated by pp collisions with the UV photons from the BLR, and with the optical photons coming from the accretion disk, inject secondary electron-positron pairs.

The optical depth for gamma rays propagating in this scenario can be calculated with the expression for the total cross section provided by <cit.>, being the threshold condition for pair production . Since we assume ϵ=10eV for the BLR photons, gamma rays with E_γ>30GeV satisfy this condition. In the case of the absorption by accretion disk photons, the threshold is exceeded by gamma photons with E_γ>1.5TeV.

The injection of secondary particles Q^ sec_ e(E_ e) (in units of erg^-1 s^-1 cm^-3) produced in photon-photon interactions, if , can be approximated as <cit.>

 
    Q^ sec_ e(E_ e)=   3/32c σ_ T/m_ e c^2_γ_ e^∞_ϵ_γ/4 γ_ e (ϵ_γ-γ_ e)^∞ dϵ_γ dω n_γ(ϵ_γ)/ϵ_γ^3 n_ ph(ω)/ω^2
       ×{4 ϵ_γ^2/γ_ e (ϵ_γ-γ_ e)ln[ 4 γ_ e ω (ϵ_γ-γ_ e)/ϵ_γ]-8 ϵ_γ ω. 
       . + 2 (2 ϵ_γ ω-1) ϵ_γ^2/γ_ e (ϵ_γ-γ_ e)-( 1-1/ϵ_γ ω) ϵ_γ^4/γ_ e^2 (ϵ_γ-γ_ e)^2},
 

 where γ_ e=E_ e/(m_ e c^2), ϵ_γ=E_γ/(m_ e c^2), and ω=ϵ/(m_ e c^2). These particles interact and emit by the same processes as the primary electrons. Considering that the synchrotron radiation dominates the cooling of the electrons, we only calculate this emission for the secondaries.



 §.§ Results


Figures <ref> and <ref> show the SEDs obtained for L_p/L_e=1 and L_p/L_e=100, respectively. We find that the luminosity at the lowest frequencies (radio) is particularly sensitive to this ratio. The radiation from primary leptons dominates in this part of the spectrum only if the power that goes to the protons is significantly less than 100 L_e. The radio luminosities in the cases of L_p/L_e=1 and L_p/L_e=100 differ by about a factor 10.

The optical region of the spectrum is dominated by the thermal radiation from the accretion disk, whereas the high-energy part is non-thermal emission produced as a consequence of the acceleration of hadrons. Most of the gamma emission generated in pp collisions is absorbed and converted to secondary particles. The synchrotron radiation of these secondaries prevails in the energy range from  to , having a maximum of  at around .

The hard X-rays and gamma luminosity produced by the collision of one BLR cloud are several orders of magnitude below the values typically detected in AGNs by Swift, INTEGRAL, and Fermi. Therefore, a single  event is not expected to be observed as a flare. The very high-energy gamma-ray tail of a single impact might be detected in the future in nearby sources by the forthcoming Cherenkov Telescope Array (CTA). However, since these photons can be easily absorbed if they travel through a dense visible or IR photon field (e.g., from a stellar association or the emission of the dusty torus), they might be strongly attenuated. Nevertheless, we note that the slope of the SED agrees very well with the observational data of a few galaxies like NGC 1068, NGC 4945, and Circinus <cit.>. Given that the total number of BLR clouds may be around 10^8 or more, it seems more realistic to think about multiple simultaneous collisions, in which case the observed luminosity will be the sum of the individual events. For this reason, in the next section we apply our model to NGC 1068 and discuss the possibility of simultaneous impacts.

       




§ APPLICATION TO NGC 1068


NGC 1068 is a spiral edge-on galaxy in the constellation Cetus whose distance to Earth is D∼14.4 Mpc <cit.>. This object is classified as a Seyfert 2 galaxy and inspired the AGN unified model <cit.>. Its bolometric luminosity is estimated to be  <cit.>. Although it is considered a star-forming galaxy, its emission can not be completely explained using a starburst model only <cit.>.

We apply our model to NGC 1068, using the parameters provided by <cit.> (see Table <ref>). For the BLR cloud, we assume the same parameters adopted previously. We suppose that the total luminosity will be the radiation produced by a single event multiplied by a number N_ events, which is the number of simultaneous events. We fix N_ events requiring to match the total gamma emission observed by Fermi in the range from 100 MeV to 100 GeV, which is  <cit.>. We discuss the validity of this assumption in Section <ref>.

We compare the multiple-event SED with the radio observations taken with VLBA <cit.> and ALMA <cit.>, and the gamma-ray spectra produced with the last Fermi catalog (8 yr) <cit.> (see Table <ref>). The radio observations are considered upper limits because the fluxes reported correspond to regions whose sizes are far larger than the region we are modeling. Furthermore, the data at 256 GHz is the integrated luminosity in a region of 9.1 pc <cit.>. On the other hand, the spatial resolution of the observation at 694 GHz is 4 pc and consequently the thermal emission from the dusty torus is also included <cit.>. 

<cit.>, based on NuSTAR observations, suggested that even the hard X-ray luminosity of NGC 1068 is obscured because of its Compton-thick nature. This scenario was recently reviewed and confirmed by <cit.>. This implies that the measured X-ray emission is not intrinsic, but transmitted by reflections. In this situation, the intrinsic radiation in the source is higher than observed.


 

 

Considering a magnetization of 10% (Eq. <ref>), we find in the case of L_p/L_e=100 that the required number of simultaneous N_ events to match the gamma luminosity observed by Fermi is ∼ 2.8×10^3. The luminosity in the range of the Swift data is , which is more than twice the emission measured, implying that the source should be ∼ 60% obscured if the contribution of other sources in that band is negligible. The radio flux at 256 GHz is overestimated by about 12% (see Fig. <ref>). In consequence, we calculate the SED for  and  to see whether higher magnetic fields improve the results. The corresponding magnetization ratios, maximum energies for the particles, and the luminosity in some bands are shown in Table <ref> for the two scenarios. The corresponding SEDs are presenting in Fig. <ref> and Fig. <ref>.

 

 

 

 





 We see  in all the cases that the VLBA limit is not exceeded because the radiation is strongly attenuated by SSA (see Fig. <ref>). With , the total number of events required to reach the total gamma emission measured by Fermi is 1.9×10^3, whereas with  G the presence of 1.5×10^3 simultaneous impacts is enough (see Table <ref> and Fig. <ref>). The hard X-ray total luminosity predicted is  for , and  for . Therefore, the obscuration of the source should be at least between 20% and 40%. These values can be increased by the existence of other sources emitting hard X-rays (e.g., a corona). Coronae have characteristic temperatures of ∼10^9 K and emit X-rays by Comptonization of photons from the accretion disk <cit.>. The expected luminosity of coronae can be similar or even up to a few of orders of magnitude higher than  produced in our scenario, in which case the obscuration percentage raises. The detection of the Fe K-alpha line in NGC 1068 suggests the presence of a corona in this source, but the evidence is still  not conclusive <cit.>. With these magnetic field values, the upper limits imposed by the observations with ALMA are not violated (see Fig. <ref> and Table <ref>).





§ DISCUSSION


In Section <ref>, we assume between 1.5×10^3 and 2.8×10^3 simultaneous events in order to achieve the emission observed by Fermi. Many authors found that the number of clouds in the BLR should be ∼ 10^8 or even larger <cit.>. <cit.> found that the total number could go up to 10^18 depending on the value of the filling factor and the optical depth. 

Under the assumption that the clouds are uniformly distributed, their number per unit of volume (n_ clouds) will be just the total number of clouds (N_ clouds) divided by the volume of the BLR (V_ BLR). Calculating the characteristic radius as mentioned previously, for NGC 1068 we find . Since the BLR can be thought of as a thin shell extending from r to R_ BLR (see Table <ref>), the resulting volume is . On the other hand, the number of impacts per unit of time is Ṅ_ clouds=n_ clouds v_ c π r^2. Requiring that N_ events=Ṅ_ clouds t_ ss=1.5×10^3, Ṅ_ clouds becomes 0.06, which means N_ clouds should be ∼ 3×10^8. This value agrees very well with the number of clouds estimated from the observations. The characteristic luminosity fluctuation would be √(N_ events)/N_ events∼ 2.6% in t_ ss/N_ events∼ 16 s, assuming Poisson statistics <cit.>. Increasing the number of events  to 2.8×10^3, the total number of clouds in the BLR should be N_ clouds=6×10^8. The variability expected in this case is √(N_ events)/N_ events∼ 1.8% in t_ ss/N_ events∼ 9 s. Therefore, the emission produced by these events for any of the magnetic field values considered will be detected as continuous and our previous analysis becomes valid.

Long-term variations in the X-ray luminosity of radio-quiet AGNs, usually understood as changes in the size and properties of the corona, are not predicted by this model <cit.>. Nevertheless, the existence of a corona is not incompatible with the model here presented. Fluctuations in the observed X-ray emission from the impacts could be produced by variations in the absorbers in the line of sight. Strong modifications of the local environment, for instance due to a change in the accretion regimen, could also result in alterations of the X-ray luminosity, but the gamma emission should also be  affected.

Considering that in pp collisions the neutrinos produced by charged pions take ∼ 5% of the energy of the relativistic proton <cit.>, it is possible to have neutrinos with energies in the detection range of IceCube. Therefore, and given the maximum energies achieved by the particles, the processes presented in this work might contribute to the spectrum reported by <cit.>.





§ SUMMARY AND CONCLUSIONS


In this paper we investigated the acceleration of particles and the non-thermal emission produced by the collision of broad-line region clouds and the accretion disk in active galactic nuclei. We proposed as the acceleration mechanism the diffusive shock acceleration and calculated the maximum energies that the particles can reach. We found that, depending on the strength of the magnetic field, electrons can be accelerated up to 36 GeV, whereas the proton maximum energy rises  to ∼ 4 PeV. The energy losses for electrons are dominated by synchrotron, whereas pp interactions dominate the cooling for the protons. The accelerated particles cool down locally and do not escape from the source.

We found that the emission of a single event cannot be detected as a flare, whereas the luminosity of multiple simultaneous events can explain the gamma radiation of NGC 1068 if its nucleus is at least obscured between 20% and 40% at hard X-ray frequencies. The high-energy gamma photons produced in pp inelastic collisions are absorbed in the BLR radiation field, injecting secondary electrons. These secondaries emit synchrotron radiation in the detection range of Fermi.

The number of simultaneous events needed to account for the gamma rays observed varies between 1.5×10^3 and 2.8×10^3, depending on the magnetic field assumed. These numbers are feasible if the total number of BLR clouds is between 3×10^8 and 6×10^8. The variability of luminosity in time generated because of the superposition of sources is too small to be detected. Given the maximum energies achieved by the protons, neutrinos with energies in the detection range of IceCube might be created in the collision of BLR clouds with accretion disks.

All in all, the model presented is an attractive alternative to explain the high-energy emission in active systems deprived of powerful jets. Further observations with the next generation of X-ray and gamma satellites <cit.> might contribute to validating and distinguishing our model from other possible scenarios <cit.>.




 We would like to thank the anonymous reviewer for suggestions and comments. ALM thanks to S. del Palacio for fruitful discussions, and to A. Streich for his technical support. This work was supported by the Argentine agencies CONICET (PIP 2014-00338), ANPCyT (PICT 2017-2865) and the Spanish Ministerio de Economía y Competitividad (MINECO/FEDER, UE) under grant  AYA2016-76012-C3-1-P.
















aa 
 


