Galactic Emission Research Articles

Maser Activity of Organic Molecules toward Sgr B2(N)

At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH2NH, HNCNH, and HCOOCH3. However, the lack of spatial distribution information on these new maser species has prevented us from assessing the excitation conditions of the maser emission as well as their pumping mechanisms. Here, we present a mapping study toward Sgr B2 north (N) to locate the region where the complex maser emission originates. We report the first detection of the Class I methanol (CH3OH) maser at 84 GHz and the first interferometric map of the methanimine (CH2NH) maser at 5.29 GHz toward this region. In addition, we present a tool for modeling and fitting the unsaturated molecular maser signals with non-LTE radiative transfer models and Bayesian analysis using the Markov Chain Monte Carlo approach. These enable us to quantitatively assess the observed spectral profiles. The results suggest a two-chain-clump model for explaining the intense CH3OH Class I maser emission toward a region with low continuum background radiation. By comparing the spatial origin and extent of maser emission from several molecular species, we find that the 5.29 GHz CH2NH maser has a close spatial relationship with the 84 GHz CH3OH Class I masers. This relationship serves as observational evidence to suggest a similar collisional pumping mechanism for these maser transitions.

A new analytical model of the cosmic-ray energy flux for Galactic diffuse radio emission

Low-frequency radio observations of diffuse synchrotron radiation offer a unique vantage point from which to investigate the intricate relationship between gas and magnetic fields in the formation of structures within the Galaxy, spanning from the diffuse interstellar medium (ISM) to star-forming regions. Achieving this pivotal objective hinges on a comprehensive understanding of cosmic-ray properties; these dictate the effective energy distribution of relativistic electrons, which are primarily responsible for the observable synchrotron radiation. Notably, cosmic-ray electrons (CRe) with energies of between 100 MeV and 10 GeV play a crucial role in determining the majority of the sky brightness below the GHz range. However, their energy flux (je) remains elusive because of solar modulation. We propose a way to derive observational constraints on this energy gap of interstellar CRe through the brightness temperature spectral index of low-frequency radio emission, here denoted βobs. We introduce a new parametric analytical model that fits available data for je in accordance with the βobs values measured in the literature between 50 MHz and 1 GHz for diffuse emission in the Milky Way. Our model accounts for multiple observations considering magnetic-field strengths consistent with existing measurements below 10 μG. We present a first all-sky map of the average component of the magnetic field perpendicular to the line of sight and validate our methodology against state-of-the art numerical simulations of the diffuse ISM. This research makes headway in modeling Galactic diffuse emission with a practical, parametric form. It provides essential insights that will help preparations for the imminent arrival of the Square Kilometre Array.

Constraints on the Origins of the Galactic Neutrino Flux Detected by IceCube

Galactic and extragalactic objects in the Universe are sources of high-energy neutrinos that may contribute to the astrophysical neutrino signal seen by IceCube. Recently, a study done using cascade-like events seen by IceCube reported neutrino emission from the Galactic plane with >4σ significance. In this work, we put a lower limit on the number of Galactic sources required to explain this emission. To achieve this, we use a simulation package created to simulate point sources in the Galaxy along with the neutrino and gamma-ray flux emissions originating from them. Along with using past IceCube discovery potential curves, we also account for Eddington bias effects due to Poisson fluctuations in the number of detected neutrino events. We present a toy Monte Carlo simulation to show that there must be at least eight sources, each with luminosity less than 1.6 × 1035 erg s−1, responsible for the Galactic neutrino emission. Our results constrain the number of individual point-like emission regions, which apply both to discrete astrophysical sources and to individual points of diffuse emission.

Constraining the Inner Galactic DM Density Profile with H.E.S.S.

In this short review, corresponding to a talk given at the conference “Cosmology 2023 in Miramare”, we combine an analysis of five regions observed by H.E.S.S. in the Galactic Center, intending to constrain the Dark Matter (DM) density profile in a WIMP annihilation scenario. For the analysis, we include the state-of-the-art Galactic diffuse emission Gamma-optimized model computed with DRAGON and a wide range of DM density profiles from cored to cuspy profiles, including different kinds of DM spikes. Our results are able to constrain generalized NFW profiles with an inner slope γ≳1.3. When considering DM spikes, the adiabatic spike is completely ruled out. However, smoother spikes given by the interactions with the bulge stars are compatible if γ≲0.8, with an internal slope of γsp-stars=1.5.

Robust inference of the Galactic Centre gamma-ray excess spatial properties

ABSTRACT The gamma-ray Fermi-LAT Galactic Centre excess (GCE) has puzzled scientists for over 15 yr. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and point sources are masked. Our main aim is to compare a spherically symmetric morphology – potentially arising from the annihilation of dark matter (DM) particles – with a boxy morphology – expected if faint unresolved sources in the Galactic bulge dominate the excess emission. Recent claims favouring a DM-motivated template for the GCE are shown to rely on a specific Galactic bulge template, which performs worse than other templates for the Galactic bulge. We find that a non-parametric model of the Galactic bulge derived from the VISTA Variables in the Via Lactea survey results in a significantly better fit for the GCE than DM-motivated templates. This result is independent of whether a galprop-based model or a more non-parametric ring-based model is used to describe the diffuse Galactic emission. This conclusion remains true even when additional freedom is added in the background models, allowing for non-parametric modulation of the model components and substantially improving the fit quality. When adopted, optimized background models provide robust results in terms of preference for a boxy bulge morphology for the GCE, regardless of the mask applied to the Galactic plane.

On the Galactic radio signal from stimulated decay of axion dark matter

We study the full-sky distribution of the radio emission from the stimulated decay of axions which are assumed to compose the dark matter in the Galaxy. Besides the constant extragalactic and CMB components, the decays are stimulated by a Galactic radio emission with a spatial distribution that we empirically determine from observations. We compare the diffuse emission to the counterimages of the brightest supernovæ remnants, and take into account the effects of free-free absorption. We show that, if the dark matter halo is described by a cuspy NFW profile, the expected signal from the Galactic center is the strongest. Interestingly, the emission from the Galactic anti-center provides competitive constraints that do not depend on assumptions on the uncertain dark matter density in the inner region. Furthermore, the anti-center of the Galaxy is the brightest spot if the Galactic dark matter density follows a cored profile. The expected signal from stimulated decays of axions of mass ma ∼ 10-6 eV is within reach of the Square Kilometer Array for an axion-photon coupling gaγ ≳ (2-3) × 10-11 GeV-1.

Galactic Diffuse Emission from Radio to Ultra-high-energy γ-Rays in Light of Up-to-date Cosmic-Ray Measurements

Cosmic rays (CRs) travel throughout the Galaxy, leaving traces from radio to ultra-high-energy γ-rays due to interactions with the interstellar gas, radiation field, and magnetic field. Therefore, it is necessary to utilize multiwavelength investigations on the Galactic diffuse emission to shed light on the physics of CR production and propagation. In this work, we present a spatially dependent propagation scenario, taking account of a local source contribution, while making allowances for an additional CR component freshly accelerated near their sources. In this picture, after reproducing the particle measurements at the solar system, we calculated the intensity and compared the spectral energy distribution to observations from Fermi-LAT and LHAASO-KM2A in the γ-ray band, and from WMAP and Planck among other radio surveys at lower energies. Multiband data considered in conjunction, the former comparison exhibits sufficiently good consistency in favor of our model, while the latter calls for improvement in data subtraction and processing. From this standpoint, there remains potential for advanced observations at energies from milli-eVs to MeVs toward the Galactic plane, in order to evaluate our model further and more comprehensively in the future.

Galactic diffuse neutrino emission from sources beyond the discovery horizon

Galactic diffuse neutrino emission from sources beyond the discovery horizon

Modeling the spectral energy distribution of starburst galaxies

Context. Analyzing multiwavelength observations of galaxies from the far-ultraviolet to the millimeter domains provides a wealth of information on the physical properties of galaxies and their evolution across cosmic time. Existing or upcoming ground-based or space-borne facilities with enhanced sensitivities and resolutions open an unprecedented window on the galaxy evolution in the early Universe. However, the derivation of galaxy properties from nebular emission lines is not trivial because the interstellar medium in a galaxy may be patchy, and emission might originate both from starburst emission regions and from partially covered photon-dominated regions. Aims. We model both the nebular continuum emission and the line emission of the spectral energy distribution for galaxies exhibiting both a HII region-like emission and emission like that from a photon-dominated regions to account for the partial shielding of the starburst emission region by dense clouds. Methods. Nebular galactic emission was modeled from far-ultraviolet to millimeter ranges in a two-sector model with an HII region and a photon-dominated region. The partial overlap of the HII region by the photon-dominated region was accounted for by a covering factor. We generated grids of emission spectra using the Cloudy photoionization code for our two-sector model. Results. We compared our models with spectral lines from different samples of galaxies for which we mixed characteristic emission from starburst regions and denser regions. We show that the infrared line ratios can constrain the density, metallicity, photoionization parameter, and the covering factor. We also built infrared diagnostic diagrams based on different infrared line ratios in which the galaxy location contains information about its physical conditions. Conclusions. The two-sector model that couples starburst emission regions and photon-dominated regions can span the existing observations. We implement the resulting emission line libraries in the CIGALE galaxy spectral energy distribution code to help interpret spectrophotometric observations.

Origin and Composition of the Galactic Diffuse X-Ray Emission Spectra by Unresolved X-Ray Sources

Galactic diffuse X-ray emission (GDXE) can be spatially segmented into Galactic center X-ray emission (GCXE), Galactic ridge X-ray emission (GRXE), and Galactic bulge X-ray emission (GBXE). The X-ray spectra of GDXE are expressed by the assembly of compact X-ray sources, which are either white dwarfs (WDs) or X-ray active stars consisting of binaries with late-type stars. WDs have either a strong magnetic field or a weak magnetic field. WDs and X-ray active stars are collectively called compact X-ray stars. However, spectral fittings by the assembly of all compact X-ray stars for GCXE, GRXE, and GBXE are rejected, leaving significant excess near the energies of the Kα, Heα, and Lyα lines. These excesses are found in the collisional ionization equilibrium (CIE) plasma. Thus, the spectra of GRXE and GBXE are improved by adding CIE supernova remnants (SNRs). However, the GCXE spectrum is still unacceptable, with significant data excess due to radiative recombination emission (recombining plasma (RP)). The GCXE fit is then significantly improved by adding aged RP-SNRs. Aged RP-SNRs are made by a past big flare of Sgr A* emitting either hard X-rays or low-energy cosmic rays. The big flares may excite Fe and Ni atoms in cold diffuse gas (cold matter (CM)) and emit fluorescent X-ray lines. The CIE-SNRs, RP-SNRs, and CM are called diffuse X-ray sources. This paper presents the spectral fits by the assembly of all the compact and diffuse X-ray sources together with high-quality spectra and a combined fit among all the GDXE of GCXE, GRXE, and GBXE. This provides the first scenario to quantitatively and comprehensively predict the origin of the GDXE spectra.

The history of the Milky Way: The evolution of star formation, cosmic rays, metallicity, and stellar dynamics over cosmic time

Abstract We study the long-term evolution of the Milky Way (MW) over cosmic time by modeling the star formation, cosmic rays, metallicity, stellar dynamics, outflows, and inflows of the galactic system to obtain various insights into the galactic evolution. The mass accretion is modeled by the results of cosmological N-body simulations for the cold dark matter. We find that the star formation rate is about half the mass accretion rate of the disk, given the consistency between observed Galactic diffuse X-ray emissions (GDXEs) and possible conditions driving the Galactic wind.Our model simultaneously reproduces the quantities of star formation rate, cosmic rays, metals, and the rotation curve of the current MW. The most important predictions of the model are that there is an unidentified accretion flow with a possible number density of ∼10−2 cm−3 and that part of the GDXEs originates from a hot, diffuse plasma which is formed by consuming about $10\%$ of supernova explosion energy. The latter is the science case for future X-ray missions: XRISM, Athena, and so on. We also discuss further implications of our results for the planet formation and observations of external galaxies in terms of multi-messenger astronomy.

Galactic Gamma-Ray Diffuse Emission at TeV Energies with HAWC Data

Galactic gamma-ray diffuse emission (GDE) is emitted by cosmic rays (CRs), ultra-relativistic protons, and electrons, interacting with gas and electromagnetic radiation fields in the interstellar medium. Here we present the analysis of teraelectronvolt diffuse emission from a region of the Galactic plane over the range in longitude of l ∈ [43°, 73°], using data collected with the High Altitude Water Cherenkov (HAWC) detector. Spectral, longitudinal, and latitudinal distributions of the teraelectronvolt diffuse emission are shown. The radiation spectrum is compatible with the spectrum of the emission arising from a CR population with an index similar to that of the observed CRs. When comparing with the DRAGON base model, the HAWC GDE flux is higher by about a factor of 2. Unresolved sources such as pulsar wind nebulae and teraelectronvolt halos could explain the excess emission. Finally, deviations of the Galactic CR flux from the locally measured CR flux may additionally explain the difference between the predicted and measured diffuse fluxes.

Probing the Galactic Diffuse Continuum Emission with COSI

In 2016, the Compton Spectrometer and Imager (COSI) had a successful 46 day flight on board NASA’s Super Pressure Balloon platform. In this work, we report measurements of the Galactic diffuse continuum emission (GDCE) observed toward the inner Galaxy during the flight, which in the COSI energy band (0.2–5 MeV) is primarily generated from inverse Compton radiation. Within uncertainties, we find overall good agreement with previous measurements from INTEGRAL/SPI and COMPTEL. Based on these initial findings, we discuss the potential for further probing the GDCE with the 2016 COSI balloon data, as well as prospects for the upcoming satellite mission.

The Milky Way revealed to be a neutrino desert by the IceCube Galactic plane observation

The Galactic diffuse emission (GDE) is formed when cosmic rays leave the sources where they were accelerated, diffusively propagate in the Galactic magnetic field and interact with the interstellar medium and interstellar radiation field. GDE in γ-rays (GDE-γ) has been observed up to subpetaelectronvolt energies, although its origin may be explained by either cosmic-ray nuclei or electrons. Here we show that the γ-rays accompanying the high-energy neutrinos recently observed by the IceCube Observatory from the Galactic plane have a flux that is consistent with the GDE-γ observed by the Fermi-LAT and Tibet ASγ experiments around 1 TeV and 0.5 PeV, respectively. The consistency suggests that the diffuse γ-ray emission above ~1 TeV could be dominated by hadronuclear interactions, although a partial leptonic contribution cannot be excluded. Moreover, by comparing the fluxes of the Galactic and extragalactic diffuse emission backgrounds, we find that the neutrino luminosity of the Milky Way is one-to-two orders of magnitude lower than the average of distant galaxies. This finding implies that our Galaxy has not hosted the type of neutrino emitters that dominates the isotropic neutrino background at least in the past few tens of kiloyears.

The intensity of diffuse galactic emission reflected by meteor trails

ABSTRACT We calculate the reflection of diffuse galactic emission by meteor trails and investigate its potential relationship to meteor radio afterglow (MRA). The formula to calculate the reflection of diffuse galactic emission is derived from a simplified case, assuming that the signals are mirrored by the cylindrical overdense ionization trail of meteors. The overall observed reflection is simulated through a ray tracing algorithm together with the diffuse galactic emission modelled by the Global Sky Model sky model. We demonstrate that the spectrum of the reflected signal is broad-band and follows a power law with a negative spectral index of around −1.3. The intensity of the reflected signal varies with local sidereal time and the brightness of the meteor and can reach 2000 Jy. These results agree with some previous observations of MRAs. Therefore, we think that the reflection of galactic emission by meteor trails can be a possible mechanism causing MRAs, which is worthy of further research.

Uncertainties of the 30–408 MHz Galactic emission as a calibration source for radio detectors in astroparticle physics

Context. Arrays of radio antennas have proven to be successful in astroparticle physics with the observation of extensive air showers initiated by high-energy cosmic rays in the Earth’s atmosphere. Accurate determination of the energy scale of the primary particles’ energies requires an absolute calibration of the radio antennas for which, in recent years, the utilization of the Galactic emission as a reference source has emerged as a potential standard. Aims. To apply the “Galactic calibration” a proper estimation of the systematic uncertainties on the prediction of the Galactic emission from sky models is necessary, which we aim to quantify on a global level and for the specific cases of selected radio arrays. We further aim to determine the influence of additional natural radio sources on the Galactic calibration. Methods. We compared seven different sky models that predict the full-sky Galactic emission in the frequency range from 30 to 408 MHz. We made an inventory of the reference maps on which they rely and used the output of the models to determine their global level of agreement. We subsequently took typical sky exposures and the frequency bands of selected radio arrays into account and repeated the comparison for each of them. Finally, we studied and discuss the relative influence of the quiet Sun, the ionosphere, and Jupiter. Results. We find a systematic uncertainty of 14.3% on the predicted power from the Galactic emission, which scales to approximately half of that value as the uncertainty on the determination of the energy of cosmic particles. When looking at the selected radio arrays, the uncertainty on the predicted power varies between 11.7% and 21.5%. The influence of the quiet Sun turns out to be insignificant at the lowest frequencies but increases to a relative contribution of ~30% around 400 MHz.

Constraints on the magnetic field in the intercluster bridge A399–A401

Galaxy cluster mergers are natural consequences of structure formation in the Universe. Such events involve the dissipation of a large amount of energy (∼1063 erg) during the process. Part of this energy can be channelled in particle acceleration and magnetic field amplification, enhancing non-thermal emission of the intra- and intercluster environment. Recently, low-frequency observations led to the detection of a bridge of diffuse synchrotron emission connecting two merging galaxy clusters, Abell 399 and Abell 401. This result provides clear observational evidence of relativistic particles and magnetic fields in between clusters. In this work, we used LOw Frequency ARray (LOFAR) observations at 144 MHz to study the polarised emission in the A399–A401 bridge region for the first time. No polarised emission was detected from the bridge region. Assuming a model where polarisation is generated by multiple shocks, depolarisation can be due to Faraday dispersion in the foreground medium with respect to the shocks. We constrained its Faraday dispersion to be greater than 0.10 rad m−2 at 95% confidence level, which corresponds to an average magnetic field in the bridge region of greater than 0.46 nG (or 0.41 nG if we include regions of the Faraday spectrum that are contaminated by Galactic emission). This result is largely consistent with the predictions from numerical simulations for megaparsec regions where the gas density is about 300 times higher than the mean gas density.

A model for the infrared-radio correlation of main sequence galaxies at gigahertz frequencies and its variation with redshift and stellar mass

Context. The infrared-radio correlation (IRRC) of star-forming galaxies can be used to estimate their star formation rate (SFR) based on the radio continuum luminosity at MHz–GHz frequencies. For its practical application in future deep radio surveys, it is crucial to know whether the IRRC persists at high redshift z. Aims. Previous works have reported that the 1.4 GHz IRRC correlation of star-forming galaxies is nearly z-invariant up to z ≈ 4, but depends strongly on the stellar mass M⋆. This should be taken into account for SFR calibrations based on radio luminosity. Methods. To understand the physical cause behind the M⋆ dependence of the IRRC and its properties at higher z, we constructed a phenomenological model for galactic radio emission. Our model is based on a dynamo-generated magnetic field and a steady-state cosmic ray population. It includes a number of free parameters that determine the galaxy properties. To reduce the overall number of model parameters, we also employed observed scaling relations. Results. We find that the resulting spread of the infrared-to-radio luminosity ratio, q(z, M⋆), with respect to M⋆ is mostly determined by the scaling of the galactic radius with M⋆, while the absolute value of the q(z, M⋆) curves decreases with more efficient conversion of supernova energy to magnetic fields and cosmic rays. Additionally, decreasing the slope of the cosmic ray injection spectrum, αCR, results in higher radio luminosity, decreasing the absolute values of the q(z, M⋆) curves. Within the uncertainty range of our model, the observed dependence of the IRRC on M⋆ and z can be reproduced when the efficiency of supernova-driven turbulence is 5%, 10% of the kinetic energy is converted into magnetic energy, and αCR ≈ 3.0. Conclusions. For galaxies with intermediate to high (M⋆ ≈ 109.5 − 1011 M⊙) stellar masses, our model results in an IRRC that is nearly independent of z. For galaxies with lower masses (M⋆ ≈ 108.5 M⊙), we find that the IR-to-radio flux ratio increases with increasing redshift. This matches the observational data in that mass bin which, however, only extends to z ≈ 1.5. The increase in the IR-to-radio flux ratio for low-mass galaxies at z ≳ 1.5 that is predicted by our model could be tested with future deep radio observations.

Modeling the Galactic center gamma-ray emission with more realistic cosmic-ray dynamics

Context. Very-high-energy gamma-ray observations of the Galactic center (GC) show extended emission that is strongly correlated with the morphology of the central molecular zone (CMZ). The best explanation for that emission is a hadronic interaction between cosmic rays (CRs) and ambient gas, where a CR central and continuous source accelerates protons up to 1 PeV (“PeVatron”). However, current models assume very simplistic CR dynamics. Aims. Our goal is to verify if more realistic CR dynamics for the GC environment are consistent with current gamma-ray observations, and whether they could be constrained by upcoming observations with the Cherenkov Telescope Array (CTA). Methods. We generated synthetic gamma-ray maps using a CR transport model with spherical injection, different diffusion regimes (in and out of the CMZ), polar advection, and mono-energetic particles of 1 PeV, and including different CR populations injected from the Arches, Quintuplet, and nuclear clusters of young massive stars, plus supernova Sgr A East. We adopted two different 3D gas distributions consistent with the observed gas column density, either with or without an inner cavity. Results. In order to reproduce the existing observations detected by the High Energy Stereoscopic System (HESS), a ring-like gas distribution, with its mass set by the standard Galactic CO-to-H2 conversion factor, and CR acceleration from all relevant sources are required. For a conversion factor one order of magnitude lower, injection rates that are ten times higher are needed. We show that CTA will be able to differentiate between models with different CR dynamics, proton sources, and CMZ morphologies, owing to its unprecedented sensitivity and angular resolution. Conclusions. More realistic CR dynamics suggest that the CMZ has a large inner cavity and that the GC PeVatron is a composite CR population accelerated by the Arches, Quintuplet, and nuclear star clusters, and Sgr A East.

Galactic Diffuse γ-Ray Emission from GeV to PeV Energies in Light of Up-to-date Cosmic-Ray Measurements

Diffuse γ-ray emission between 10 and 1000 TeV from the Galactic plane was recently measured by the Large High Altitude Air Shower Observatory (LHAASO). These observations will help tremendously in constraining the propagation and interaction of cosmic rays in the Milky Way. Additionally, new measurements of cosmic-ray spectra reach a very high precision of up to 100 TeV energies, revealing multiple spectral structures of various species. In this work, based on up-to-date measurements of local cosmic-ray spectra and a simplified propagation setup, we confront a model prediction of diffuse γ-ray emission with measurements of diffuse γ-rays. To better constrain the low-energy part of the model, we analyze the 14.6 yr of Fermi Large Area Telescope (Fermi-LAT) data to extract the Galactic diffuse emission between 1 and 500 GeV from the same sky regions of LHAASO, after subtracting the contribution from known sources and the isotropic diffuse γ-ray background. The joint Fermi-LAT and LHAASO spectra thus cover a very wide energy range from 1 GeV to 1 PeV with small gaps from 0.5 to 10 TeV. Compared with the prediction, we find that clear excesses between several GeV and ∼60 TeV of the diffuse emission exist. Possible reasons to explain the excesses may include unresolved sources or more complicated propagation models. We illustrate that an exponential cutoff power-law component with an index of −2.40 and a cutoff energy of ∼30 TeV is able to account for such excesses.