High-energy Gamma-ray Emission Research Articles

The Galactic Population of Pulsar Wind Nebulae and the Contribution of Its Unresolved Component to the Diffuse High-Energy Gamma-ray Emission

In this work, we provide a phenomenological description of the population of galactic TeV pulsar wind nebulae (PWNe) based on suitable assumptions for their space and luminosity distribution. We constrain the general features of this population by assuming that it accounts for the majority of bright sources observed by H.E.S.S. Namely, we determine the maximal luminosity and fading time of PWNe (or, equivalently, the initial period and magnetic field of the pulsar powering the observed emission) by performing a statistical analysis of bright sources in the H.E.S.S. galactic plane survey. This allows us to estimate the total luminosity and flux produced by galactic TeV PWNe. We also evaluate the cumulative emission from PWNe that cannot be resolved by H.E.S.S., showing that this contribution can be as large as ∼40% of the total flux from resolved sources. We argue that also in the GeV domain, a relevant fraction of this population cannot be resolved by Fermi-LAT, providing a non-negligible contribution to the large-scale diffuse emission in the inner galaxy. This additional component could naturally account for a large part of the spectral index variation observed by Fermi-LAT, weakening the evidence of cosmic ray spectral hardening in the inner galaxy. Finally, the same result is obtained for PeV energy, for which the sum of the diffuse component, due to unresolved PWNe, and the truly diffuse emission well saturates the recent Tibet AS-γ data, without the need to introduce a progressive hardening of the cosmic-ray spectrum toward the galactic centre.

The Upcoming GAMMA-400 Experiment

The upcoming GAMMA-400 experiment will be implemented aboard the Russian astrophysical space observatory, which will be operating in a highly elliptical orbit over a period of 7 years to provide new data on gamma-ray emissions and cosmic-ray electron + positron fluxes, mainly from the galactic plane, the Galactic Center, and the Sun. The main observation mode will be a continuous point-source mode, with a duration of up to ~100 days. The GAMMA-400 gamma-ray telescope will study high-energy gamma-ray emissions of up to several TeV and cosmic-ray electrons + positrons up to 20 TeV. The GAMMA-400 telescope will have a high angular resolution, high energy and time resolutions, and a very good separation efficiency for separating gamma rays from the cosmic-ray background and the electrons + positrons from protons. A distinctive feature of the GAMMA-400 gamma-ray telescope is its wonderful angular resolution for energies of >30 GeV (0.01° for Eγ = 100 GeV), which exceeds the resolutions of space-based and ground-based gamma-ray telescopes by a factor of 5–10. GAMMA-400 studies can reveal gamma-ray emissions from dark matter particles’ annihilation or decay, identify many unassociated, discrete sources, explore the extended sources’ structures, and improve the cosmic-ray electron + positron spectra data for energies of >30 GeV.

Statistical Relationship between Long-duration High-energy Gamma-Ray Emission and Solar Energetic Particles

Large solar eruptions are often associated with long-duration γ-ray emission extending well above 100 MeV. While this phenomenon is known to be caused by high-energy ions interacting with the solar atmosphere, the underlying dominant acceleration process remains under debate. Potential mechanisms include continuous acceleration of particles trapped within large coronal loops or acceleration at coronal mass ejection (CME)-driven shocks, with subsequent back-propagation toward the Sun. As a test of the latter scenario, previous studies have explored the relationship between the inferred particle population producing the high-energy γ-rays and the population of solar energetic particles (SEPs) measured in situ. However, given the significant limitations on available observations, these estimates unavoidably rely on a number of assumptions. In an effort to better constrain theories of the γ-ray emission origin, we reexamine the calculation uncertainties and how they influence the comparison of these two proton populations. We show that, even accounting for conservative assumptions related to the γ-ray flare, SEP event, and interplanetary scattering modeling, their statistical relationship is only poorly/moderately significant. However, though the level of correlation is of interest, it does not provide conclusive evidence for or against a causal connection. The main result of this investigation is that the fraction of the shock-accelerated protons required to account for the γ-ray observations is >20%–40% for six of the 14 eruptions analyzed. Such high values argue against current CME-shock origin models, predicting a <2% back-precipitation; hence, the computed number of high-energy SEPs appears to be greatly insufficient to sustain the measured γ-ray emission.

Gamma-ray emission from spectrally resolved cosmic rays in galaxies

ABSTRACT Cosmic rays (CRs) are ubiquitous in the interstellar medium (ISM) of nearby galaxies, but many of their properties are not well-constrained. Gamma-ray observations provide a powerful tool in this respect, allowing us to constrain both the interaction of CR protons with the ISM and their transport properties. To help better understand the link between observational signatures and CR physics, we use a series of magneto-hydrodynamical arepo simulations of isolated galaxies performed using spectrally resolved CR transport in every computational cell, with subsequent gamma-ray emission calculated using the Crayon + (Cosmic RAY emissiON) code. In each of our simulated haloes, modelling the energy-dependent spatial diffusion of CRs leads to a more extended distribution of high-energy (~100 GeV) gamma rays compared to that predicted by a ‘grey’ steady-state model, which is especially visible in the corresponding emission maps and radial profiles. Despite this, the total gamma-ray spectra can often be well approximated by the steady-state model, although recovering the same spectral index typically requires a minor variation of the energy dependence of the diffusion coefficient. Our simulations reproduce the observed shapes of the gamma-ray spectra of nearby star-forming galaxies and also match recent observations of the far infrared–gamma-ray relation. We find, however, that the spectrally resolved model yields marginally smaller luminosities for lower star formation rates compared to grey simulations of CRs. Our work highlights the importance of modelling spectrally resolved CR transport for an accurate prediction of spatially resolved high-energy gamma-ray emission, as will be probed by the upcoming Cherenkov Telescope Array observatory.

VERITAS Discovery of Very High Energy Gamma-Ray Emission from S3 1227+25 and Multiwavelength Observations

We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on 2015 May 15 with the Fermi-Large Area Telescope (LAT). A combined 5 hr VERITAS exposure on May 16 and 18 resulted in a strong 13σ detection with a differential photon spectral index, Γ = 3.8 ± 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target-of-opportunity observations with Swift, optical photometry, polarimetry, and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of τ obs = 6.2 ± 0.9 hr, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cutoff. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model.

High-energy gamma-ray emission powered by a young protostar: the case of S255 NIRS 3

ABSTRACT Evidence of efficient acceleration of cosmic rays in massive young stellar objects has been recently reported. Among these massive protostars, S255 NIRS 3, for which extreme flaring events associated with radio jets have been detected, is one of the best objects to test this hypothesis. We search for gamma-ray emission associated with this object in Fermi-LAT data and inspect the gas content in different molecular lines using the MWISP survey. A GeV source dubbed 4FGL J0613.1+1749c lies on top of the MYSO region, where two filamentary ∼10 pc CO structures extend along the same direction of the subparsec radio jets. We investigate the spectrum, morphology, and light curve of the gamma-ray source and compare it with the theoretical emission expected from hadronic and leptonic populations accelerated in the radio jets. We argue that the gamma-ray source could be powered by particles accelerated in the S255 NIRS 3 jets, radiating via bremsstrahlung or proton–proton interaction, and with a synchrotron component shinning in radio from primary or secondary electrons in the case of a leptonic or hadronic population.

An updated view and perspectives on high-energy gamma-ray emissionfrom SGR J1935+2154 and its environment

SGR J1935+2154 was discovered in 2016 and is currently one of the most burst-active Soft Gamma-ray Repeaters (SGR), having emitted many X-ray bursts in recent years. In one of our previous articles, we investigated the contribution to high-energy and very high-energy gamma-ray emission (VHE, E > 100 GeV) due to cosmic-ray acceleration of SNR G57.2+0.8 hosting SGR J1935+2154 using the GALPROP propagation code. However, follow-up observations of SGR 1935+2154 were made for 2 hours on April 28, 2020, using the High Energy Stereoscopic System (H.E.S.S.). The observations coincide with X-ray bursts detected by INTEGRAL and Fermi/Gamma-ray Burst Monitor (GBM). These are the first high-energy gamma-ray observations of an SGR in a flaring state, and upper limits on sustained and transient emission have been derived. Now that new H.E.S.S. observations have been made, it is interesting to update our model with respect to these new upper limits. We extend our previous results to a more general situation using the new version of GALPROP. We obtain a hadronic model that confirms the results discussed by H.E.S.S. . This leads to an optimistic prospect that cosmic ray gamma rays from SGR J1935+2154 can contribute to the overall gamma energy density distribution and in particular to the diffusion gamma rays from the Galactic center.

Constraints on the Very High Energy Gamma-Ray Emission from Short GRBs with HAWC

Many gamma-ray bursts (GRBs) have been observed from radio wavelengths, and a few at very high energies (VHEs, >100 GeV). The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is well suited to study transient phenomena at VHEs owing to its large field of view and duty cycle. These features allow for searches of VHE emission and can probe different model assumptions of duration and spectra. In this paper, we use data collected by HAWC between 2014 December and 2020 May to search for emission in the energy range from 80 to 800 GeV coming from a sample of 47 short GRBs that triggered the Fermi, Swift, and Konus satellites during this period. This analysis is optimized to search for delayed and extended VHE emission within the first 20 s of each burst. We find no evidence of VHE emission, either simultaneous or delayed, with respect to the prompt emission. Upper limits (90% confidence level) derived on the GRB fluence are used to constrain the synchrotron self-Compton forward-shock model. Constraints for the interstellar density as low as 10−2 cm−3 are obtained when assuming z = 0.3 for bursts with the highest keV fluences such as GRB 170206A and GRB 181222841. Such a low density makes observing VHE emission mainly from the fast-cooling regime challenging.

Search for Gamma-ray Emission from Accretion Flares of Tidal Disruption Events Possibly Associated with the IceCube Neutrinos

Outflows or disk-coronas generated in tidal disruption events (TDEs) of supermassive black holes have been suggested as possible sites of high-energy neutrinos. Three TDEs (AT2019dsg, AT2019fdr and AT2019aalc) have been claimed to be associated with high-energy astrophysical neutrinos in multi-messenger follow-ups. No GeV photons have been detected accompanying the neutrino for the three sources. In this work, we searched for the high-energy gamma-ray emission from a larger sample of TDE candidates observed by the Zwicky Transient Facility (ZTF). No significant GeV emission was observed, and the upper limits of the gamma-ray emission flux are reported. We then performed a stacking analysis for the sample sources and found that the collective gamma-ray emission of this class of sources was also not bright enough to be detected by the Fermi Large Area Telescope (Fermi-LAT). The nondetection of the high-energy gamma-ray emission from the sample TDEs could be due to the fact that the high-energy gamma rays are absorbed by soft photons in the source. Using a model-based hypothesis, the upper limit on the emission radius of the neutrino production is obtained for these TDEs: R&lt;1016 cm for typical TDE parameter values.

VLBI observations of GRB 201015A, a relatively faint GRB with a hint of very high-energy gamma-ray emission

Context. A total of four long-duration gamma-ray bursts (GRBs) have been confirmed at very high-energy (≥100GeV) with high significance, and any possible peculiarities of these bursts will become clearer as the number of detected events increases. Multi-wavelength follow-up campaigns are required to extract information on the physical conditions within the jets that lead to the very high-energy counterpart, hence they are crucial to reveal the properties of this class of bursts. Aims. GRB 201015A is a long-duration GRB detected using the MAGIC telescopes from ~40 s after the burst. If confirmed, this would be the fifth and least luminous GRB ever detected at these energies. The goal of this work is to constrain the global and microphysical parameters of its afterglow phase, and to discuss the main properties of this burst in a broader context. Methods. Since the radio band, together with frequent optical and X-ray observations, proved to be a fundamental tool for overcoming the degeneracy in the afterglow modelling, we performed a radio follow-up of GRB 201015A over 12 different epochs, from 1.4 days (2020 October 17) to 117 days (2021 February 9) post-burst, with the Karl G. Jansky Very Large Array, e-MERLIN, and the European VLBI Network. We include optical and X-ray observations, performed respectively with the Multiple Mirror Telescope and the Chandra X-ray Observatory, together with publicly available data, in order to build multi-wavelength light curves and to compare them with the standard fireball model. Results. We detected a point-like transient, consistent with the position of GRB 201015A until 23 and 47 days post-burst at 1.5 and 5 GHz, respectively. No emission was detected in subsequent radio observations. The source was also detected in optical (1.4 and 2.2 days post-burst) and in X-ray (8.4 and 13.6 days post-burst) observations. Conclusions. The multi-wavelength afterglow light curves can be explained with the standard model for a GRB seen on-axis, which expands and decelerates into a medium with a homogeneous density. A circumburst medium with a wind-like profile is disfavoured. Notwithstanding the high resolution provided by the VLBI, we could not pinpoint any expansion or centroid displacement of the outflow. If the GRB is seen at the viewing angle θ that maximises the apparent velocity βapp (i.e. θ ~ βapp-1), we estimate that the Lorentz factor for the possible proper motion is Гα ≤ 40 in right ascension and Гδ ≤ 61 in declination. On the other hand, if the GRB is seen on-axis, the size of the afterglow is ≤5pc and ≤16pc at 25 and 47 days. Finally, the early peak in the optical light curve suggests the presence of a reverse shock component before 0.01 days from the burst.

Evidence of TeV halos around millisecond pulsars

Using data from the HAWC gamma-ray Telescope, we have studied a sample of 37 millisecond pulsars (MSPs), selected for their spindown power and proximity. From among these MSP, we have identified four which favor the presence of very high-energy gamma-ray emission at a level of $(2\Delta \ln \mathcal{L})^{1/2} \ge 2.5$. Adopting a correlation between the spindown power and gamma-ray luminosity of each pulsar, we performed a stacked likelihood analysis of these 37 MSPs, finding that the data supports the conclusion that these sources emit very high-energy gamma-rays at a level of $(2\Delta \ln \mathcal{L})^{1/2} = 4.24$. Among sets of randomly selected sky locations within HAWC's field-of-view, less than 1\% of such realizations yielded such high statistical significance. Our analysis suggests that MSPs produce very high-energy gamma-ray emission with a similar efficiency to that observed from the Geminga TeV-halo, $\eta_{\rm MSP} = (0.39-1.08) \times \eta_{\rm Geminga}$. This conclusion poses a significant challenge for pulsar interpretations of the Galactic Center gamma-ray excess, as it suggests that any population of MSPs potentially capable of producing the GeV excess would also produce TeV-scale emission in excess of that observed by HESS from this region. Future observations by CTA will be able to substantially clarify this situation.

Probing the High-energy Gamma-Ray Emission Mechanism in the Vela Pulsar via Phase-resolved Spectral and Energy-dependent Light-curve Modeling

Abstract Recent kinetic simulations sparked a debate regarding the emission mechanism responsible for pulsed GeV γ-ray emission from pulsars. Some models invoke curvature radiation, while other models assume synchrotron radiation in the current sheet. We interpret the curved spectrum of the Vela pulsar as seen by H.E.S.S. II (up to ∼100 GeV) and the Fermi Large Area Telescope to be the result of curvature radiation due to primary particles in the pulsar magnetosphere and current sheet. We present phase-resolved spectra and energy-dependent light curves using an extended slot gap and current-sheet model, invoking a step function for the accelerating electric field as motivated by kinetic simulations. We include a refined calculation of the curvature radius of particle trajectories in the lab frame, impacting the particle transport, predicted light curves, and spectra. Our model reproduces the decrease of the flux of the first peak relative to the second one, evolution of the bridge emission, near-constant phase positions of peaks, and narrowing of pulses with increasing energy. We can explain the first of these trends because we find that the curvature radii of the particle trajectories in regions where the second γ-ray light-curve peak originates are systematically larger than those associated with the first peak, implying that the spectral cutoff of the second peak is correspondingly larger. However, an unknown azimuthal dependence of the E field, as well as uncertainty in the precise spatial origin of the GeV emission, precludes a simplistic discrimination of emission mechanisms.

Red Dwarfs as Sources of Cosmic Rays

The Pamela, Fermi, AMS-02 spectrometer experimental data cannot be explained using the diffusive models of propagation of cosmic-rays accelerated at the supernova shocks and require the existence of nearby sources of cosmic rays at the distances less than one kpc. These sources could explain the growth of the ratio of galactic positrons to electrons with an energy increase, the complex dependence of the exponent of the proton and alpha spectra on the energy, the cosmic ray anomaly component origin. We consider active dwarf stars as possible sources of galactic cosmic rays in the energy range up to ~ 1014 eV. The generation of high-energy cosmic rays should be accompanied by high-energy gamma-ray emission, which may be detected. The TeV gamma-ray emission mostly of flaring type from active red dwarf stars V388 Cas, V547 Cas, V780 Tau, V962 Tau, V1589 Cyg, GJ 3684, GJ 1078 and GL 851.1 was detected recently with SHALON long-term observations. This result confirms that active dwarf stars are also the sources of high-energy galactic cosmic rays.

Searching signals of dark matter from unidentified Fermi-LAT objects with H.E.S.S

Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos, as shown from cosmological N-body simulations. These subhalos could shine in gamma-rays and be eventually detected as unidentified sources in gamma-ray surveys. From a thorough selection of unidentified Fermi-LAT Objects (UFOs), we observe four UFOs with H.E.S.S. and we search for very high-energy (VHE, E ≥ 100 GeV) gamma-ray emission. Considering dark matter masses above a few hundred GeV, the observed UFOs could be identified as dark matter subhalos, given their hard gamma-ray spectra in the few-ten-to-hundred GeV energy range. Since no significant very-high-energy gamma-ray emission is detected in any of the four UFOs dataset nor in the combined one, we derive constraints on the product of the velocity-weighted annihilation cross-section 〈συ〉 by the J-factor for the dark matter models. We derive 95% CL upper limits on 〈συ〉 J in W + W − and τ+τ− annihilation channels for the TeV dark matter particles. Considering thermal WIMPs, we derive constraints on the J-factors from the H.E.S.S. observations. Assuming model-dependent predictions from cosmological N-body simulations on the J-factor distribution for Milky Way-sized galaxies, the dark matter models with masses greater than 0.3 TeV for the UFO emissions can be ruled out at high confidence level.

Redshift determination of blazars for the Cherenkov Telescope Array

AbstractBlazars are the most numerous type of observed high-energy gamma-ray emitters. However, their emission mechanisms and population properties are still not well-understood. Crucial to this understanding are their cosmological redshifts, which are often not easy to obtain. This presents a great challenge to the next-generation ground-based observatory for very-high-energy gamma rays, the Cherenkov Telescope Array (CTA), which aims to detect a large number of distant blazars to study their intrinsic emission properties and to place tight constraints on the extragalactic background light density, amongst others. The successful investigation of these subjects needs a precise redshift determination. Motivated by these challenges, the CTA redshift task force initiated more than 3 years ago a spectroscopic observing program using some of the largest optical and infrared telescopes to measure the redshifts of a large fraction of blazars that are likely to be detected with CTA. In this proceedings, we give an overview of the CTA redshift task force, discuss some of the difficulties associated with measuring the redshifts of blazars and present our sample selection and observing strategies. We end the proceedings with reporting selected results from the program, the on-going collaborative efforts and our plans for the future.

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

ABSTRACT We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016–2020. Optical, X-ray, and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very high energy (VHE; ≳100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by The Owens Valley Radio Observatory indicates the multizone nature of the emission from this object. We constrain the VHE duty cycle of the source to be &amp;lt;16 2014-like flares per year (95 per cent confidence). For the first time for this source, a broad-band low-state spectral energy distribution is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 ± 0.93stat) × 1021 cm−2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modelled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad-line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet.

High-energy gamma-ray emission from SNR G57.2+0.8 hosting SGR J1935+2154

In recent years some efforts have been made to identify active sources capable of accelerating particles up to 1015 eV, known as PeVatrons. Measurements of TeV (1012 eV) gamma-rays from supernova remnants (SNRs) have shown that efficient particle acceleration can occur in SNR diffusive shocks. In this paper, we obtain the contribution to the high energy and very high-energy gamma-ray (VHE, E > 100 GeV) emission due to cosmic-ray acceleration from SNR G57.2+0.8 hosting the Soft Gamma Repeater (SGR) J1935+2154 with the use of the GALPROP code. To do so, we take into account the SNR + SGR association as a single source close to the Galactic center.We propose that the above setting can provide a more comprehensive scenario for the generation of GeV-TeV gamma-rays. We also discuss the contribution from the SNR G57.2+0.8 and SGR J1935+2154 region to the diffusive TeV energy gamma-ray emission from the Galactic center.

Exploring the possible gluon condensation signature in gamma-ray emission from pulsars

The very high energy gamma-ray emissions from pulsars are usually considered to be dominated by leptonic scenario since the hadronic flux is weak. We point out that the gluon condensation predicted by a nonlinear QCD evolution equation may greatly enhance the cross sections of proton-target interactions and give rise to a characteristic broken power law in the gamma-ray spectra. The result is used to explore the gluon condensation signature in the observed gamma-ray spectra from pulsars.

Nearby SN-associated GRB 190829A: Environment, Jet Structure, and VHE Gamma-Ray Afterglows

Abstract We present a self-consistent paradigm to interpret the striking features of nearby low-luminosity GRB 190829A. Its prompt gamma-ray lightcurve has two separated pulses. We propose that the interaction of the hard prompt gamma-ray photons ( E p = 624 − 303 + 2432 keV) of its initial pulse with the dusty medium (A V = 2.33) not only results in the second soft gamma-ray pulse (E p ∼ 12 keV) but also makes a pre-accelerated e ± -rich medium shell via the γγ annihilation. In this paradigm, we show that the observed radio, optical, and X-ray afterglow lightcurves are well-fitted with the forward shock model. Its jet is almost isotropic (θ j &gt; 1.0 rad) with a Lorentz factor of ∼35, and the electron density of the e ± -rich medium shell is ∼15 cm−3, which is about seven times higher than the electron density of its normal surrounding medium. The GRB ejecta catches up with and propagates into the e ± -rich medium shell at a region of R = (4.07–6.46) × 1016 cm. This results in a bright afterglow bump at ∼103 seconds after the GRB trigger. The predicted very high energy (VHE) gamma-ray emission from the synchrotron self-Compton process agrees with the H.E.S.S. observation. The derived broadband spectral energy distribution shows that GRB 190829A like nearby GRBs would be promising targets of the VHE gamma-ray telescopes, such as H.E.S.S., MAGIC, and CTA (Cerenkov Telescope Arrays).

Prospects of newly detecting nearby star-forming galaxies by the Cherenkov Telescope Array

ABSTRACT Prospects of the Cherenkov Telescope Array (CTA) for the study of very high energy gamma-ray emission from nearby star-forming galaxies are investigated. In the previous work, we constructed a model to calculate luminosity and energy spectrum of pion-decay gamma-ray emission produced by cosmic ray interaction with the interstellar medium (ISM), from four physical quantities of galaxies [star formation rate (SFR), gas mass, stellar mass, and effective radius]. The model is in good agreement with the observed GeV–TeV emission of several nearby galaxies. Applying this model to nearby galaxies that are not yet detected in TeV (mainly from the KINGFISH catalogue), their hadronic gamma-ray luminosities and spectra are predicted. We identify galaxies of the highest chance of detection by CTA, including NGC 5236, M33, NGC 6946, and IC 342. Concerning gamma-ray spectra, NGC 1482 is particularly interesting because our model predicts that this galaxy is close to the calorimetric limit and its gamma-ray spectral index in GeV–TeV is close to that of cosmic ray protons injected into ISM. Therefore, this galaxy may be detectable by CTA even though its GeV flux is below the Fermi Large Area Telescope sensitivity limit. In the TeV regime, most galaxies are not in the calorimetric limit, and the predicted TeV flux is lower than that assuming a simple relation between the TeV luminosity and SFR of M82 and NGC 253, typically by a factor of 15. This means that a more sophisticated model beyond the calorimetric limit assumption is necessary to study TeV emission from star-forming galaxies.