An interacting molecular cloud scenario for production of gamma-rays and neutrinos from MAGIC J1835-069, and MAGIC J1837-073

Recently the MAGIC telescope observed three TeV gamma-ray extended sources in the galactic plane in the neighborhood of radio SNR G24.7+0.6. Among them, the PWN HESS J1837-069 was detected earlier by the HESS observatory during its first galactic plane survey. The other two sources, MAGIC J1835-069 and MAGIC J1837-073 are detected for the first time at such high energies. Here we shall show that the observed gamma-rays from the SNR G24.7+0.6 and the HESS J1837-069 can be explained in terms of hadronic interactions of the PWN/SNR accelerated cosmic rays with the ambient matter. We shall further demonstrate that the observed gamma-rays from the MAGIC J1837$-$073 can be interpreted through hadronic interactions of runaway cosmic-rays from PWN HESS J1837-069 with the molecular cloud at the location of MAGIC J1837-073. No such association has been found between MAGIC J1835$-$069 and SNR G24.7+0.6 or PWN HESS J1837$-$069. We have examined the maximum energy attainable by cosmic-ray particles in the SNR G24.7+0.6/ PWN HESS J1837-069 and the possibility of their detection with future gamma-ray telescopes. The study of TeV neutrino emissions from the stated sources suggests that the HESS J1837$-$069 should be detected by IceCube Gen-2 neutrino telescope in a few years of observation.


I. INTRODUCTION
While studying a star-forming region in the inner part of the galactic plane centred around the supernova remnant (SNR) G24.7+0.6 (3FHL J1834.1 − 0706 / FGES J1834.1 − 0706), which a middle-aged composite SNR, the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) Telescopes recently discovered a new extended gamma-ray source MAGIC J1835−069 at energies above 150 GeV to till 5 TeV in the neighborhood of the supernova remnant (SNR) G24.7+0.6 [1]. Besides, the MAGIC telescope also detected high energy gamma rays from two other sources in the region namely HESS J1837−069, and MAGIC J1837−073. In fact, the HESS J1837−069 is the brightest object in the region which was discovered by the High Energy Stereoscopic System (HESS) observatory during its first galactic plane scan [2]. The MAGIC J1837−073 is adjacent to HESS J1837−069 and coincides with the HESS hotspot in the south of HESS J1837−069.
The MAGIC observations suggest that the center position of the extended object MAGIC J1835−069 is located 0.34 • away from the center of the SNR G24.7+0. 6. The spectral studies of the stated two sources suggest that the statistical significance for a common origin of the MAGIC J1835−069 and the SNR G24.7+0.6 is less than 1.5 . The gamma-ray spectral measurements done by Fermi-LAT and MAGIC observatories suggest that the two sources are likely to be correlated [1]. The observation of 13 CO J = 1 − 0 line at 110 GHz favors for molecular cloud origin of observed gamma-rays from MAGIC J1835−069 [3,4].
The gamma-ray extended source HESS J1837−069 was initially classified as an unidentified object [2]. It coincides with the x-ray source AX J1838.0-0655. A young and energetic 70.5 ms rotation powered pulsar PSR J1838-0655 in the region * pbanik74@yahoo.com † aru_bhadra@yahoo.com was discovered by the Rossi X-Ray Timing Explorer (RXTE) [5]. The high resolution of Chandra observation revealed that AX J1838.0-0655 contains a bright point source surrounded by a centrally peaked nebula [5]. Following these discoveries, the source HESS J1837−069 is interpreted as a pulsar wind nebula [6,7]. The size of the AX J1838.0-0655 in x-rays is smaller than the HESS J1837−069 observed in gamma rays which are generally interpreted as the diffusion of high energy particles in the surrounding medium. The origin of GeV-TeV gamma rays from MAGIC J1837−073 is so far not discussed in the literature.
Under the circumstances, the main objectives of the present work are as follows: First, we would like to examine whether the spectral behavior of observed gamma-rays from the HESS J1837−069 and SNR G24.7+0.6 can be explained by hadronic interactions of the PWN/ SNR accelerated cosmic rays with the ambient (proton) matter. Secondly, with the same cosmic ray production spectrum which is required to explain the observed gamma-rays directly from HESS J1837−069 and SNR G24.7+0.6, we would like to explore whether one can consistently explain the observed EM SED from the sources MAGIC J1837−073 and MAGIC J1835−069 with a scenario of cosmic-ray-illuminated cloud interacting with the PWN HESS J1837−069/SNR G24.7+0.6. We would like to determine the corresponding neutrino fluxes from all the stated sources and subsequently study their detection likelihood by the present/upcoming neutrino observatories.
The plan of the paper is the following: In the next section, we shall describe the methodology for evaluating the gammaray and neutrino fluxes from the PWN/SNR and the molecular cloud interacting with the PWN/SNR. The gas distributions and subsequently masses of the regions under the study assuming MAGIC J1835−069, and MAGIC J1837−073 located in the vicinity of PWN HESS J1837−069/SNR G24.7+0.6 are estimated in the same section. In sect. 3, the numerical estimated fluxes of hadronically produced gamma-rays and neutrinos from the PWN/SNR and surrounding molecular clouds over the GeV to TeV energy range along with the observations. The findings of the present work will be discussed in Sect. 4 and we shall conclude finally in the same section.

II. THE GAMMA RAYS AND NEUTRINOS PRODUCTION METHODOLOGY
The production spectrum of cosmic rays at PWN/SNR can be represented by a power law and is given by where is the spectral index, is the energy of the cosmic rays. The normalization constant is related with the supernova explosion energy via the relation where denotes the efficiency of supernova explosion energy into cosmic rays, , and , represents the minimum and maximum energies of accelerated cosmic ray protons.
The high energy gamma-ray flux from the PWN/SNR can be explained by the interaction of the shock accelerated cosmic rays with the ambient matter (protons) of density . Such interaction will produce gamma-rays and neutrinos after subsequent decay of neutral and charged pions respectively along with the other secondary particles. Here we follow Banik & Bhadra (2017) [8] to estimate the differential flux of high energy gamma rays reaching the Earth from the SNR G24.7+0.6. The corresponding flux of high energy neutrinos reaching the Earth from the direction of the SNR can be written as where ( ) represents the resulting neutrino emissivity which can be obtained by following [9,10] and denotes the distance between the PWN/SNR and the Earth.
After escaped away from the PWN/SNR, the accelerated cosmic ray protons diffuse in the interstellar medium with a diffusion coefficient given by [11] = 0 ( 10 ) where 0 represents the diffusion coefficient of accelerated protons at 10 GeV energy (∼ 10 28 cm 2 s −1 in our Galactic medium) and indicates a constant having value between 0.3 to 0.7 [12,13]. However, a slow diffusion has been inferred in the dense gaseous medium of molecular clouds [14,15]. The measurement of Boron to Carbon Flux Ratio in Cosmic Rays over the rigidity extent of 1.9 GV to 2.6 TV by the Alpha Magnetic Spectrometer (AMS-02) on the International Space Station found that is around 0.33 [16]. Becker et al. (2016) [17] have recently demonstrated that the cosmic ray budget can be likened well for a diffusion coefficient that is close to ∝ 0.3 by studying 21 SNRs, those are well-studied from radio wavelengths up to gamma-ray energies. On the other hand, if cosmic rays originated at SNRs of the galaxy, a more powerful diffusion with = 0.5 cannot generate the observed cosmic ray energy spectrum especially at the highenergy (TeV) component of the spectrum [17]. The probability density of finding a cosmic ray particle at a given radius from a PWN/SNR for an impulsive injection spectrum, after including the energy loss due to hadronic −interaction of cosmic rays in PWN/SNR, can be written as [14,18,19] ( , , ) = 1 where denotes the age of the supernova explosion/PWN, = stands for represents the cosmic ray energy loss time-scale due to hadromic interaction of cosmic rays with ambient matter of the PWN/SNR, ∼ 0.45 denotes the inelasticity of the interaction, and [9] being the interaction cross-section [19]. Here, is the radius of the sphere up to which the cosmic ray particles have time to travel after their injection [19]. Therefore, the intensity of cosmic rays at a distance from the source can be written as [10] where it is assumed that cosmic rays are diffusing from a point source.
To explain the MAGIC collaboration observed high energy gamma-rays from the source MAGIC J1835−069 and MAGIC J1837−073, we consider a dense molecular cloud at the center of this source. The cloud is continuously illuminated by runaway relativistic cosmic ray particles accelerated at the nearby SNR G24.7+0.6 or the PWN HESS J1837−069. We follow Banik & Bhadra (2017) [8] and Banik & Bhadra (2018) [10] to find out the differential flux of high energy gamma rays and neutrinos reaching the Earth from the source MAGIC J1835−069 and MAGIC J1837−073 after their production in the interaction of runaway relativistic cosmic ray particles with molecular clouds.
Since our objective is to examine whether the observed gamma rays from MAGIC J1835−069 and MAGIC J1837−073 are correlated with the two known high energy gamma rays SNR/PWN of nearby regions (SNR G24.7+0.6 and HESS J1837−069), we exclude the leptonic origin scenario of gamma rays from SNR G24.7+0.6 and HESS J1837−069. The decay of charged pions produced in pp interactions at the molecular cloud region will finally lead to electrons and positrons along with the neutrinos. The so produced electrons will give synchrotron photons while traveling through the magnetic field of the molecular cloud. However, the density of soft photons in the cloud region (synchrotron photons) will be low because of the weak magnetic field of the clouds. As a result, the inverse Compton process is not significant in the present scenario. The interstellar molecular clouds are a complex mixture of gas at certain velocities. The major constituent of molecular clouds is hydrogen which may remain in three different phases, the neutral atomic hydrogen (HI), the ionized hydrogen (HII), and the molecular hydrogen (H 2 ). The cold H 2 is not directly observable in emission. Carbon monoxide (CO) is widely used as tracer for H 2 , the intensity of rotational transition of CO relates with the abundance of H 2 . The column density of atomic hydrogen can be estimated straightway from the 21 cm radio frequency transition between the two hyperfine levels of the ground state of HI. The regions of HII is characterized by thermal free-free radiation and recombination lines.
To estimate the column density for H 2 , we take data for 12 CO = 1 − 0 emission from the carbon monoxide (CO) composite survey [20]. We evaluate the column density of molecular hydrogen by assuming a linear relationship 2 = × , where represent the velocity-integrated brightness temperature of CO 2.6-mm line [21,22]. Here, the conversion factor, is assumed to be 3 × 10 20 cm −2 K −1 km −1 s which is consistent with the suggested conversion factor of ≈ (2 − 4) × 10 20 cm −2 K −1 km −1 s from different observations for the molecular clouds very near to galactic plane [20,23]. We use the circular rotation model for the milky-way Galaxy as given in Ranasinghe & Leahy (2018) [24] and the rotational curve data for gas velocities at different Galactocentric distance [25] to correlate measured radial velocity of gases with kinematic distance. We obtained a radial velocity of the molecular hydrogen gas about 88 km s −1 and 31.2 km s −1 which corresponds to 6.6 kpc and 3.5 kpc distance from us respectively. Therefore we choose a radial velocity range of = 78 − 98 km s −1 (or = 21 − 41 km s −1 ) to be integrated to determine the column density of molecular hydrogen at a distance of about 6.6 Kpc (or 3.5 Kpc) and obtained column density map of gases are shown in the Fig. 1.
We obtain the column density of HII from the Planck observed free-free emission at 30 GeV frequency [26]. First, we use the conversion factor at 30 GHz as given in table 3 of Planck Collaboration X 2016 [26] to convert the measured brightness temperature into free−free intensity ( ). Then we follow the equation (5) in Sodroski et al. (1997) [27] to transform the free−free intensity into column density. To calculate the H II column density, we adopted an electron temperature of = 8000 K and the effective density of electrons of = 10 cm −3 [27].
The column density of HI in the aforementioned regions is estimated under the consideration of the optically thin limit [21,28] from the data cube of Galactic HI 4 survey (HI4PI) [29]. To evaluate gas column density distributions of HI, we integrated the gases over the radial velocities = 78 − 98 km s −1 for 6.6 kpc and = 21 − 41 km s −1 for 3.5 kpc.
The total mass within a molecular cloud/source has been estimated from the expression [7] represents the mass of the hydrogen atom, is the column density of hydrogen in ith phase ( ≡ 2, , ), denotes to the angular area, and indicates the distance of the molecular cloud/source.
The total mass of hydrogen gas content of the sources SNR G24.7+0.6 (at 3.5 kpc), HESS J1837−069 (at 6.6 kpc) and two possible associated sources, i.e MAGIC J1835−069 and MAGIC J1837−073 either at 3.5 kpc or 6.6 kpc are shown in Table I and Table II respectively.

III. RESULTS AND DISCUSSION
Before the report by the MAGIC collaboration, Katsuta et al. (2017) [21] carried out a study of the complex region around the SNR G24.7+0.6 with the Fermi-LAT telescope and found GeV gamma-ray emission from two elliptical extended regions   [1]. In addition, through X-ray observations of the region with XMM−Newton [21] found a young massive OB association/ cluster, G25.18+0.26 which is proposed to be associated with the extended GeV emissions (G25A and G25B). According to the MAGIC collaboration the emission detected at VHE with MAGIC is unlikely to be correlated with the OB association/ cluster G25.18+0.26 detected in X-rays as MAGIC only detects emission from the G25A1 component [1]. As mentioned above we have estimated high energy gamma rays to explain the observed gamma-ray spectrum from the sources, SNR G24.7+0.6 and HESS J1837−069 in the hadronic interaction of cosmic rays with ambient matter. We have also investigated below the observed gamma-ray spectrum of the sources, MAGIC J1835−069 and MAGIC J1837−073 in terms of the hadronic interaction of cosmic rays with giant molecular clouds at the centre of these sources which are illuminated by the accelerated cosmic rays either from SNR G24.7+0.6 or HESS J1837−069. The estimation of gamma ray flux along with observed data and consequently corresponding predicted high energy neutrino fluxes from the said sources are discussed bellow.
A. SNR G24.7+0.6 The SNR G24.7+0.6 is a radio and gamma-ray supernova remnant advancing in a dense medium [1,[30][31][32][33][34]. A recent study based on analysis of the 1420 MHz radio continuum data, H I-line absorption spectra, and H I channel maps which is obtained from the VLA (Very Large Array) Galactic Plane Survey (VGPS) [35], suggests a new distance about ∼ 3.5 ± 0.2 kpc of the SNR from us [24]. Using the new distance ( ) measurement of the SNR [24], we have estimated the age of the SNR to be = 2.09 × 10 4 year from 2 − relation as given in [36] where (1 ) = 20 is the observed flux density from the SNR at 1 GHz frequency [30]. Using this FIG. 2. The estimated differential energy spectrum of gamma-rays reaching the Earth from the SNR G24.7+0.6. The red dotted line represents the gamma-ray flux produced due to interaction of SNR accelerated cosmic rays with the ambient matter of the SNR. The gray dotted line denotes the background gamma ray flux from the SNR due to interaction the cosmic rays galactic background with SNR material. The black continuous line represents the estimated overall differential EM SED from the SNR G24.7+0.6. The blue dash-double-dotted and brown dash-single-dotted lines represent the detection sensitivity of the CTA detector for 1000 hours and the MAGIC detector for 31 hour, respectively. distance of the SNR, the size of the SNR has been estimated to be (30.5 ± 1.7) × (15.3 ± 0.9) pc [24] and therefore mean radius of the SNR will be ∼ 11.45 pc. The mean radius in the 1420 MHz radio continuum image can be considered as the shock wave radius . Here we apply a basic Sedov model [37] and the shock wave radius in the Sedov stage can be written as where is the ambient medium (proton) density in cm −3 units, is the supernova explosion energy in units of 0.75 × 10 51 erg and 4 = /10 4 is the age of the SNR in units of 10 4 year. For an age of the SNR about ∼ 2.09 × 10 4 year and shock wave radius = 11.45 pc, we can adopt a typical supernova explosion energy of = 1.8 × 10 51 erg and an ambient matter density of the SNR of = 19 cm −3 in order to explain observed gamma-rays from the SNR G24.7+0.6.
Here we have explored whether the observed high energy gamma-ray spectrum from the SNR G24.7+0.6 can be interpreted through hadronic interactions of cosmic rays (mainly protons) accelerated by the SNR with the ambient (proton) matter. The primary cosmic-ray production spectrum has been obtained assuming a typical supernova explosion energy of we have considered an ambient matter density of the SNR of = 19 cm −3 as mentioned above in order to match the observed gamma-ray spectrum. We have found that a ∼ 10% efficiency of conversion of supernova explosion energy to cosmic ray (proton) energy can explain well the observed experimental data. The estimated differential gamma-ray flux reaching the Earth from the SNR is displayed in the Fig. 2 along with the observations. The MAGIC detector sensitivity for 31 hours of observation with the medium zenith angle (30 • − 45 • ) has been obtained from the MAGIC detector sensitivity for 50 hours as given in [38] by assuming 1/observational-time dependency of the sensitivity for short observation times [38]. The non-detection of gamma-rays from the SNR G24.7+0.6 by MAGIC observatory restricts the maximum energy , attainable by a cosmic ray proton in the SNR shock to 1 TeV. Such a choice of maximum energy of cosmic rays leads our estimation of gamma-ray flux compatible with the Fermi-LAT observation and non-detection by the MAGIC detector (considering sensitivity of 31 h for detection of a point source with 5 significance) as shown in the figure. Future large-area telescopes like CTA [39] and LHAASO [40], which have unprecedented sensitivity up to 100 TeV energies, should be able to identify the maximum attainable energy up to which cosmic rays can accelerate in the supernova remnant. We have also included the effect of background gamma-ray flux produced due to the interaction of diffuse cosmic rays in the galaxy with the matter spread over by the progenitor star after a supernova explosion. Here we use our estimated mass of gas content of the SNR of ∼ 5.3 × 10 4 M ⊙ to estimate the background gamma-ray flux. The Model fitting parameters for observed gamma-rays from SNR G24.7+0.6 are shown in Table III.
Because of the low value of maximum energy attainable by cosmic rays in the source, the SNR G24.7+0.6 should not generate neutrinos above 1 TeV energies.

MAGIC J1835−069
The gamma-ray source MAGIC J1835−069 has been detected with a peak significance of 11 and resolved from SNR G24.7+0.6 at a level of 13 for the first time [1]. Its gammaray spectrum is well-represented by a power-law function with a spectral photon index of 2.74 ± 0.08 [1]. The source positions between two extended sources detected above 10 GeV by Fermi-LAT, one of which is FGES J1836.5−0652 and the other is FGES J1834.1−0706, which are associated with HESS J1837−069 and SNR G24.7+0.6 respectively. The centre of MAGIC J1835−069 is located at = 0.34 • away with respect to the center of the radio SNR G24.7+0.6. MAGIC collaboration suggested that the detected gamma-ray emission can be explained as the result of proton−proton interaction between the cosmic rays accelerated by SNRs and the CO-rich surrounding [1].
To explain the observed high gamma-ray spectrum from the unknown source MAGIC J1835−069, we have considered MAGIC J1835−069 as a giant molecular cloud that is illuminated by the runaway cosmic rays accelerated in the SNR G24.7+0.6. Here we have taken the distance of MAGIC J1835−069 from the SNR G24.7+0.6 as = 30 pc following the MAGIC observation. Here we take the mass of the molecular cloud MAGIC J1835−069 4.6 × 10 4 M ⊙ which is estimated considering the source distance at 3.5 kpc to explain gamma-ray observations. The high energy cosmic rays accelerated in a DSA mechanism at the shock front of the SNR G24.7+0.6 may escape and finally interact hadronically with the matter (dominantly protons) of molecular cloud MAGIC J1835−069 after a diffusive propagation over a distance of . Consequently, such hadronic interaction will produce high energy gamma-rays and neutrinos from the molecular cloud after subsequent decays of neutral and charged pions. The estimated differential gammaray flux reaching the Earth from MAGIC J1835−069 is shown in the left-hand panel of Fig. 3 along with the observations. Here we have adopted a diffusion constant of 0 = 10 27 cm 2 s −1 with = 0.33 for cosmic rays. Our findings suggest that the observed GeV-TeV gamma-rays from the location of MAGIC J1835−069 can not be interpreted by the presented model for any reasonable fitting parameters. The present analysis thus does not support the hypothesis of possible association of the unknown source MAGIC J1835−069 with the SNR G24.7+0.6. The possible correlation of MAGIC J1835−069 with HESS J1837−069 is examined below.

MAGIC J1837−073
The detailed production model for GeV-TeV gamma-rays from MAGIC J1837-073 (3FGL J1837.6-0717 [41]) is not available in the literature yet. An analysis similar to the MAGIC J1835−069 as reported above suggests that this unknown TeV source is not associated with the SNR G24.7+0.6. The estimated differential gamma-ray flux reaching the Earth from MAGIC J1837−073 is shown in the right-hand panel of Fig. 3 along with the observations where we adopt the distance of MAGIC J1835−073 from the SNR G24.7+0.6 to be = 65 pc and a diffusion constant of 0 = 3 × 10 27 cm 2 s −1 with = 0.33 for cosmic rays. Here we take the mass of the molecular cloud as ∼ 2.2 × 10 4 M ⊙ which is estimated as stated above considering cloud distance of 3.5 kpc.

B. HESS J1837−069
The H.E.S.S. Galactic Plane Survey (HGPS; [42]) in this region above 500 GeV shows a large and bright source, dubbed HESS J1837−069 [2,43]. The source HESS J1837−069 is found to be associated with a spin-down pulsar PSR J1838−065 of characteristic age ≈ 2.3 × 10 4 year, surface dipole magnetic field ≈ 1.9 × 10 12 G and pulsation period ≈ 70.5 ms [5]. Using the relation given in [44,45], we find that the initial period of the pulsar was about few millisecond only and therefore initial rotation energy of the pulsar was of the order of 10 51 ergs which satisfies the required energy budget to interpret the observed radiations in hadronic interactions of so accelerated cosmic rays. A fraction of rotational energy released during slowing down of pulsar can be assumed to accelerate cosmic ray protons extracted from the surface of the neutron star [44][45][46]. The accelerated primary cosmic rays may interact with the ambient medium and produce high energy gamma-rays and neutrinos. The primary cosmic-ray production spectrum has inferred by assuming released rotational energy of = 3 × 10 51 erg with a spectral index of = −1.9 using Eq. 2. Here, we consider an ambient matter density of the nebula of = 8.7 cm −3 to meet the observed gamma-ray spectrum. We notice that the efficiency of conversion of rotational energy to cosmic ray energy of = 10% for primary cosmic ray protons can explain well the observed data.
The expected differential gamma-ray flux reaching the Earth from the PWN has presented in Fig. 4 along with the observations. We find that this hadronic model can also explain the observed SED from the source considering the maximum energy attainable by a cosmic ray proton in the PWN to be , = 70 TeV. As shown in Fig. 4, future telescopes like e−ASTROGAM [47], CTA [39] and LHAASO [40], which have higher sensitivity than the present generation gamma ray telescopes, should able to detect the EM SED of the source in a wide energy band and hence can help to identify the true physical model that leads to its EM emission. We have also considered the effect of background gamma-ray flux produced due to the interaction of galactic diffuse cosmic rays with the matter of the PWN. The mass of the gas content for the source is found out to be ∼ 2.8 × 10 5 M ⊙ which is adopted to explain the observed gamma-ray spectrum. Along with the EM SED matching by modeling, we have also reproduced the observed integral gamma-ray flux above 200 GeV from the source as reported in [2] in order to precise determination of cosmic ray production spectral index and also the maximum attainable energy . We have also determined the neutrino flux originated due to the decay of charged pions produced in the hadronic interaction of cosmic rays accelerated by HESS J1837−069 with the ambient matter (proton). The expected all flavor neutrino flux reaching at the Earth from the source is shown in the right-hand panel of Fig. 4. We have estimated the expected muon-neutrino events in the IceCube detector from the PWN  [48] and it is found out to be ∼ 1 event in one year of observations above 1 TeV energies. We found that the estimated all flavor neutrino flux around 1 TeV energies is well within the reach of the sensitivity of IceCube-Gen2, a future extension of IceCube detector [49].

MAGIC J1835−069
The projected distance of the gamma-ray source MAGIC J1835−069 from the pulsar associated with HESS J1837−069 (for a distance of 6.6 kpc) is reported to be more than ∼ 65 pc. As discussed above, it is unlikely that MAGIC J1835−069 is correlated with SNR G24.7+0.6. Here we have explored the possibility of molecular cloud origin of the observed gamma rays from MAGIC J1835−069, associated with the source HESS J1837−069. The energy density of runaway accelerated cosmic rays at the location of MAGIC J1835−069 coming from HESS J1837−069 can be obtained from Eq. 5 consistently where we have adopted the age of HESS J1837−069 to be = 2.3 ×10 4 years [5], a diffusion constant of 0 = 2 ×10 27 cm 2 s −1 with = 0.33 for cosmic rays. The distance of the molecular cloud from HESS J1837−069 is taken as 65 pc following the MAGIC observation and the mass of the cloud is considered as = 1.9 × 10 5 M ⊙ (estimated for a distance of 6.6 kpc). The estimated differential gamma-ray flux reaching the Earth from MAGIC J1835−069 produced in hadronic interaction of these runaway cosmic rays with molecular cloud materials (protons) is shown in the Fig. 5 along with the observations. Our findings suggest that MAGIC J1835−069 is also unlikely to be associated with HESS J1837−069.

MAGIC J1837−073
The energy density of escaped accelerated cosmic rays at the position of MAGIC J1837−073 coming from HESS J1837−069 can be obtained from Eq. 5 where we consider a diffusion constant of 0 = 2 × 10 27 cm 2 s −1 with = 0.33 for cosmic rays. The distance of the molecular cloud from HESS J1837−069 and mass of the cloud are adopted to be 38 pc and = 2.3 × 10 5 M ⊙ respectively. The estimated differential gamma-ray flux reaching the Earth from MAGIC J1837−073 created in hadronic interaction of these runaway cosmic rays with molecular cloud materials (protons) is shown in the left panel of Fig. 6 along with the Fermi-LAT and MAGIC observations. Our findings suggest that MAGIC J1835−073 is likely to be associated with HESS J1837−069. We have also estimated the neutrino flux originated due to the decay of charged pions created in the hadronic interaction of cosmic rays accelerated by HESS J1837−069 with the molecular cloud. The expected all flavor neutrino flux reaching at the Earth from the source is shown in the right-hand panel of Fig. 6. We have estimated the corresponding expected muon-neutrino events in the IceCube detector from the SNR following [48] and it is found out to be = 0.21 events in 1 year of observations above 1 TeV energies. We found that the estimated all flavor neutrino flux around 1 TeV energies may be detected by IceCube-Gen2 [49].

IV. CONCLUSION
On the basis of our present analysis we conclude the followings: i) The Fermi-LAT has detected the SNR G24.7+0.6 in the GeV energy range but the MAGIC telescope did not detect any TeV gamma rays from the source, instead they found a new TeV source, the MAGIC J1835−069, in the neighborhood of the SNR G24.7+0.6 resolving it from SNR G24.7+0.6 at a level of 13 . The observed gamma ray emission from the SNR G24.7+0.6 can be interpreted in terms of hadronic interaction of the SNR accelerated cosmic rays with the ambient matter but the non-detection of TeV gamma rays from the source restricts the maximum energy attainable by a cosmic ray particle in the source to just 1.3 TeV. Consequently, there will be no TeV neutrinos from the source.
ii) The present analysis suggests that the GeV-TeV emission of HESS J1837−069 can well be described by hadronic interaction of the PWN accelerated cosmic rays with the ambient matter. The maximum energy attainable by a cosmic ray proton in the source turns out to be 70 TeV which can be verified in the future by the upcoming CTA and LHAASO detectors. The estimated all flavor neutrino flux above 1 TeV energy is found to be well within the reach of the sensitivity of IceCube-Gen2, a future extension of IceCube detector.
iii) The newly detected source MAGIC J1835−069 is unlikely to be a cosmic ray illuminated molecular cloud. The present analysis shows that the runway cosmic rays either from the SNR G24.7+0.6 or the PWN HESS J1837−069 fall short to reproduce the observed TeV emission from MAGIC J1835−069 by interacting with molecular cloud at the location of MAGIC J1835−069. Instead of considering the SNR G24.7+0.6 and MAGIC J1835−069 as a resolved separate sources recently Sun et al. [7] have considered the gamma ray emission from the two sources together (over an extended region) and described the observed GeV-TeV gamma rays from the region containing SNR G24.7+0.6 and MAGIC J1835−069 by hadronic interaction of cosmic rays with ambient matter (cloud). Since the MAGIC telescope resolved MAGIC J1835−069 from the SNR G24.7+0.6 with high significance we have here treated the stated sources separately. It may happen that MAGIC J1835−069 is the major part of the ejecta of supernova explosion that leads to G24.7+0. 6. iv) The newly detected TeV gamma ray source MAGIC J1837−073 appears to be associated with the PWN HESS J1837−069.
The gamma ray emission from MAGIC J1837−073 can well be interpreted by hadronic interactions of runway cosmic rays from HESS J1837−069 with the molecular cloud at the location of MAGIC J1837−073. The flux of TeV neutrinos from the source is, however, too small to be detected by the present or any upcoming neutrino telescope.