Preprint typeset using LATEX style AASTeX6 v. 1.0
SECONDARY RADIO AND X-RAY EMISSIONS FROM GALAXY MERGERS
Chengchao Yuan1, Kohta Murase1,2, and Peter Mesz´ aros´ 1
1Department of Physics; Department of Astronomy and Astrophysics; Center for Particle and Gravitational Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 2Center for Gravitational Physics, Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan
ABSTRACT Shocks arising in galaxy mergers could accelerate cosmic-ray (CR) ions to TeV-PeV energies. While propagating in the intergalactic medium, these CRs can produce high-energy neutrinos, electron- positron pairs and gamma-rays. In the presence of intergalactic magnetic fields, the secondary pairs will radiate observable emissions through synchrotron radiation and inverse Compton scattering. In this paper, we demonstrate that these emissions can explain the radio and X-ray fluxes of merging galaxies such as NGC 660 and NGC 3256. Using our model in combination with the observations, we can constrain the gas mass, shock velocity, magnetic field and the CR spectral index s of these systems. For NGC 660 a single-zone model with a spectral index 2.1 ∼< s ∼< 2.2 is able to reproduce simultaneously the radio and X-ray observations, while a simple one-zone scenario with s ∼ 2 can describe the radio and a large fraction of X-ray observations of NGC 3256. Our work provides a useful approach for studying the dynamics and physical parameters of galaxy mergers, which can play an important part in future multi-messenger studies of similar and related extragalactic sources. Keywords: cosmic rays — galaxies: interactions — galaxies: individual (NGC 660, NGC 3256) — radio continuum: galaxies — X-rays: galaxies
1. INTRODUCTION electron-positron pairs produced by the CR interactions Star-forming galaxies including starbursts have been in such systems. Here the EIC is caused by scatter- considered as possible reservoirs of cosmic rays (CRs) ings with the cosmic microwave background (CMB), and sources of associated neutrinos and gamma rays infrared/optical starlight (SL) and extragalactic back- (e.g., Loeb & Waxman 2006; Thompson et al. 2007; ground light (EBL). In addition, since the radiation Murase et al. 2013), in which the CRs can be supplied spectrum of the merging galaxies is determined by the by not only supernovae but also hypernovae, superbub- dynamics of the galaxy interactions and the resulting bles and active galactic nuclei (Senno et al. 2015; Xiao physical conditions, this enables us to provide con- et al. 2016; Tamborra et al. 2014; Wang & Loeb 2016; straints on the magnetic field B, shock velocity vs, gas Lamastra et al. 2017; Liu et al. 2018). Interacting galax- mass Mg, etc. Different from Lisenfeld & Voelk(2010) ies, which may be accompanied by starburst activities, where shock-accelerated electrons are employed to de- have also been considered as additional accelerators of scribe the radio emissions of two colliding galaxies, UGC CRs (Kashiyama & M´esz´aros 2014; Yuan et al. 2018). 12914/5 and UGC 813/6, we present an alternative Under the conditions typical of galaxy merger systems model based on the secondary emission from inelastic pp collisions to reproduce simultaneously the radio and arXiv:1810.04155v4 [astro-ph.HE] 24 Jun 2019 synchrotron emission can extend from the radio band to the X-ray regime, while the inverse Compton scatter- X-ray observations of NGC 660 and NGC 3256. In gen- ing may be important in the ultraviolet (UV) and up to eral, secondary electrons are more natural to explain the beyond the X-ray band. electromagnetic emissions in merging galaxies. For the In this work we formulate a model which is ca- observed CRs, the electron acceleration efficiency, the pable of reproducing the radio and X-ray obser- fraction of plasma energy deposited to electrons, is at vations of specific systems using synchrotron and least two orders lower than the proton acceleration ef- −4 −2 synchrotron self-Compton (SSC) or external inverse ficiency, e.g. Ke/p = e/p ∼ 10 − 10 (Jones 2011; Compton (EIC) emissions from high-energy secondary Morlino & Caprioli 2012). This value is also consistent with the observations of Galactic supernova remnants. Furthermore, the recent particle-in-cell simulation shows −3 [email protected] a similar value, Ke/p ' 10 (Katz & Waxman 2008; 2 Yuan, Murase and Mesz´ aros´
Caprioli & Spitkovsky 2014; Park et al. 2015a). The ra- tio of the primary electrons (from shock accelerations) and the secondary electrons and positrons is approxi- mately γ, ν Ee,primary 6e −1 ' ∼< 10 . e± Ee,sec min[1, fpp,g]p Radio/X − ray R p, π0, π± 2 g where fpp,g is the effective pp optical depth in the emit- ting region. Therefore in our model with the typical electron/proton acceleration efficiencies, emission from primary electrons directly accelerated in shocks is sub- Diffusive Shock dominant compared to that from secondary electrons and positrons from pp collisions and pion decays. This is consistent with Murase et al.(2018) where they sug- gest that the secondary emissions overwhelm the pri- mary component in nearly proton calorimetric sources. Figure 1. Schematic figure showing the merger of two galax- It is possible that primary electrons can provide a non- ies. The shock was simplified as a straight line across the negligible contribution if K > 0.1, considering that dense core region. It is also in the core region where inter- e/p ∼ actions occur and neutrinos as well as electromagnetic radi- Ke/p is poorly constrained theoretically and observa- ation are produced. tionally for this system. In the following text, we focus on the primary electron/positron scenario and omit the The pions produced in the pp collisions between shock- primary electron contribution. accelerated CR ions and the galaxy gas generate, besides As a well-studied interacting system, NGC 660 is a high-energy neutrinos and γ-rays, also copious quanti- galaxy formed by the collision of two galaxies (van Driel ties of high-energy electron-positron pairs. These high- et al. 1995), which has been observed in both radio (e.g., energy leptons may produce observable synchrotron Douglas et al. 1996; Large et al. 1981; Condon et al. emissions while propagating inside the galactic mag- 2002, 1998; Dressel & Condon 1978; Bennett et al. 1986; netic fields. Here, considering the conservation of lep- Becker et al. 1991; Gregory & Condon 1991; Sramek ton numbers and muon decays, we approximate the to- 1975), microwave, infrared, UV and X-ray (e.g., Frater- tal electron-positron injection spectrum to be the same nali et al. 2004; Liu 2010; Brightman & Nandra 2011; with the neutrino production spectrum. Following the White et al. 2000) bands. Also, the magnetic field in the procedure in Yuan et al.(2018), the electron injection core region of NGC 660 is constrained in the range of spectrum can be written as 16 ± 5 µG through polarization studies (Drzazga et al. 2011). In this paper, we take NGC 660 as an example 2 1 2 dNν 1 −1 2 ε Ne(ε) = ε = pC Mgv and use our model to reproduce the radio, UV and X-ray 3 dε 12 s (1) fluxes. We also apply our model to constrain the shock × min [1, f ] , pp,g εp'20ε velocity and gas mass of the core region of NGC 660 by using the magnetic field 16 ± 5 µG as a precondition. where p is the CR ion acceleration efficiency (nor- To show that our model’s applicability can be extended mally fixed as 0.1), C = ln(εp,max/εp,min) is the normal- −2 to other similar systems, we also consider another well- ization coefficient for a ε spectrum, Mg is the gas mass studied galaxy formed through a merger, NGC 3256, as of the merging region, vs is the shock/collision veloc- a supplementary template. ity and fpp,g = κppcngσ(εp) min[tesc, tdyn] is pp optical This paper is organized as follows. In §2, we formu- depth inside the galaxy. In this expression, κpp = 0.5 is late the secondary electron-positron spectrum and calcu- the proton inelasticity, c is the speed of light, ng is the late resulting electromagnetic emissions, including syn- gas density, tesc is the escape time of CRs, tdyn ' Rg/vs chrotron radiation and SSC/EIC components. In §3 , we is the dynamic time of the merger and σ(εp) is the pp apply the formalism in §2 to the core regions of NGC cross section given by Kafexhiu et al.(2014). As galaxies 660 and NGC 3256. A summary and discussion, in- merge, strong shocks occur with a complicated morphol- cluding comparison with previous work in the context ogy over a galaxy scale, while merging cores of the two of starburst galaxies, is given in §4. galaxies lead to a dense core region. Particles are accel- erated by the shocks, and then will be distributed in a galaxy scale. The CRs diffusing in the core region will 2. SECONDARY ELECTRON SPECTRUM AND make neutrinos and gamma rays efficiently. In this work, ELECTROMAGNETIC EMISSIONS as a simplified approximation without covering the de- Secondary Radio and X-ray Emissions from Galaxy Mergers 3
2 dNγ 2 2 dNν tails of the shock structure, we assume that shocks are spectrum are correlated by ε = ε |ε =2ε . γ dεγ 3 ν dεν γ ν CR accelerators, which inject high energy CRs to the From energy conservation, we may approximately relate core region of the merging systems and initiate subse- the electron-positron injection rate to the gamma-ray quent interactions. Figure1 shows the schematic. After production rate, and the former spectrum can be writ- leaving the accelerator, the particles can propagate dif- ten as fusively or get advected away through galactic winds, 2 γγ 2 dNγ 1 2 dNν cut therefore the net escape rate is the sum of diffusion ε Ne (εe) = 2ε |εγ =2εe = ε |εν =εe , εe > εγγ /2. dεγ 3 dεν −1 −1 −1 rate and advection rate, e.g. tesc ≈ tdiff + tad . Al- (3) though the maximum CR energy εp,max and effective The total electron-positron injection spectrum is pp optical depth fpp,g depend on the geometry of the therefore the summation of Equations1 and3, or equiv- colliding galaxies, for simplicity and consistency, we as- alently we can introduce a modification factor χ(ε) = sume that the neutrinos are produced inside the core re- cut 1 + exp(−εγγ /2ε) to Equation1. gion of the interacting system and calculate the electro- With these preparatory work, we can now derive the magnetic radiation therein. This hypothesis is in good secondary electron-positron distributions and calculate agreement with the radio maps of NGC 660 and NGC the synchrotron and inverse Compton emissions. Con- 3256. Hence, to fully depict the physical condition of the sidering the dynamic time tdyn = Rg/vs, we have the core region, we introduce several quantities, the radius rate of lepton production Rg, the average magnetic field B as well as the previ- t−t Ne(ε)χ(ε) − dyn ously defined gas mass Mg and shock velocity vs. Us- Q(ε, t) = × min{1, e tesc }, (4) ing these parameters, we can write down the maximum tdyn CR energy, gas density and diffusion time explicitly as where the exponential factor describes the escape of CRs 3 vs 4 3 εp,max = 20 eBsRg c (Drury 1983), ng = Mg/( 3 πmpRg) after the dynamical time scale and is obtained through 2 and tdiff = Rg/(6Dg), respectively. Here, mp is the pro- the equation ∂N/∂t = −N/tesc. To get the electron ton mass, Dg is the diffusion coefficient and Bs is the distribution inside the galaxy, we solve the transport post-shock magnetic field which can be parameterized equation of a simplified leaky-box model 2 1 2 as a fraction of the ram pressure Bs /8π = 2 Bngmpvs ∂Ne Ne ∂ (Kashiyama & M´esz´aros 2014). As for the diffusion co- = Q(ε, t) − + [b(ε)Ne(ε, t)] (5) ∂t tesc ∂ε efficient Dg, we use a combined large and small angle diffusion equation as in Senno et al.(2015); Casse et al. where b(ε) is the electron energy loss rate due to (2001) and Yuan et al.(2018) and then it can be written synchrotron radiation, SSC/EIC and advection (bad ' explicitly as ε/tad). In our calculations, we assume Q and the dif- fusion coefficient Dg do not depend on the positions in R 2 D −1 the merging system. t '4.28 Myr g 0 diff 29 2 −1 In the synchrotron limit γ 1, the synchrotron ra- 3 kpc 10 cm s (2) e −1 diation power in the frequency range ω to ω + dω by h 1/2 2i × (ε/εc,g) + (ε/εc,g) one electron with Lorentz factor γe can be written in the well-known formula where D is defined by D ' cl /20, l ' 0.1R is √ 0 0 c c g 3e3B sin θ the coherence length of the magnetic field fluctuations p Psyn(ω, γe)dω = 2 F (X)dω (6) 2πmec and εc,g ≈ eBlc is the characteristic energy. As for the advection, the typical values of wind velocity in where θp is the angle between the electron velocity star-forming galaxies and star burst galaxies range from and the magnetic field, which is assumed to be π/2 in 500 km s−1 (Crocker 2012; Keeney et al. 2006) to 1500 our case, −1 ω 3 eB km s (Strickland & Heckman 2009). Here, we use a X = , ω = γ2 . −1 c e moderate value vw ≈ 1000 km s for interacting galax- ωc 2 mec ies since these galaxies may enter star-forming/starburst The function F (X) is given by phase. In this case we have the advection time t ' ad ∞ −1 Z 6 vw Rg Rg/vw ≈ 2.94 × 10 yr 1000 km s−1 3 kpc . F (X) = X K5/3(ξ)dξ. Inside the galaxy, the electron-positron injection spec- X trum can be modified due to additional injections via Then, it is straightforward to write down the inte- two-photon annihilation, γγ → e−e+, since the core re- grated radiation power gion can be opaque to high-energy gamma-ray photons Z 2 cut bsyn(ε) = Psyn(ω, ε/mec )dω. above a certain threshold energy εγγ . In the pion de- cay scenario, the gamma-ray spectrum and the neutrino 4 Yuan, Murase and Mesz´ aros´
ing. Now with the preparations on synchrotron radi- ation, we are able to solve the kinetic equation. One μ 1 1 μ special solution to the differential equation is the steady 1 μ state solution (∂Ne/∂t = 0), 1 1 1 1 −1