Galactic Plane Research Articles

The FAST Galactic Plane Pulsar Snapshot Survey. II. Discovery of 76 Galactic Rotating Radio Transients and the Enigma of RRATs

We have carried out the Galactic Plane Pulsar Snapshot (GPPS) survey by using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), the most sensitive systematic pulsar survey in the Galactic plane. In addition to more than 500 pulsars already discovered through normal periodical search, we report here the discovery of 76 new transient radio sources with sporadic strong pulses, detected by using the newly developed module for a sensitive single-pulse search. Their small DM values suggest that they all are Galactic rotating radio transients (RRATs). They show different properties in the follow-up observations. More radio pulses have been detected from 26 transient radio sources but no periods can be found due to a limited small number of pulses from all FAST observations. The follow-up observations show that 16 transient sources are newly identified as being the prototypes of RRATs with a period already determined from more detected sporadic pulses, and 10 sources are extremely nulling pulsars, and 24 sources are weak pulsars with sparse strong pulses. On the other hand, 48 previously known RRATs have been detected by the FAST, either during verification observations for the GPPS survey or through targeted observations of applied normal FAST projects. Except for one RRAT with four pulses detected in a session of 5-minute observation and four RRATs with only one pulse detected in a session, sensitive FAST observations reveal that 43 RRATs are just generally weak pulsars with sporadic strong pulses or simply very nulling pulsars, so that the previously known RRATs always have an extreme emission state together with a normal hardly detectable weak emission state. This is echoed by the two normal pulsars J1938+2213 and J1946+1449 with occasional brightening pulses. Though strong pulses of RRATs are very outstanding in the energy distribution, their polarization angle variations follow the polarization angle curve of the averaged normal pulse profile, suggesting that the predominant sparse pulses of RRATs are emitted in the same region with the same geometry as normal weak pulsars.

Quasi-regular variations of subpulse drifting for PSR J1857+0057

ABSTRACT During observations of the Galactic Plane Pulsar Snapshot survey by the Five-hundred-meter Aperture Spherical radio Telescope, varying subpulse drifting of PSR J1857+0057 is detected. The follow-up observation confirms the quasi-regular changes of the drifting rate about every 50 periods. We determine the drift rate D through a linear fit to the pulse-central longitudes of subpulses in a drifting band, and determine P3 from the cross-points of two fitted lines at the zero longitude for two neighbouring drifting bands. The low frequency modulation of about every 50 periods is found on variations of not only pulse intensity but also drift parameters. In most of low frequency modulation cycles, the integrated pulse intensity I and the absolute drift rate |D| tend to increase first and then decrease, and the drifting periodicity P3 varies just in the opposite. In addition, the phase-forward intensity-enhancement is observed in many modulation cycles. Based on our polarization data, the average position angle curve for pulses with a smaller |D| and larger P3 is slightly steep in the leading edge of pulse profile compared with that of the fully averaged profile.

The FAST Galactic Plane Pulsar Snapshot Survey – III. Timing results of 30 newly discovered pulsars

ABSTRACT Timing observations are crucial for determining the basic parameters of newly discovered pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) with the L-band 19-beam receiver covering the frequency range of 1.0–1.5 GHz, the FAST Galactic Plane Pulsar Snapshot (GPPS) Survey has discovered more than 600 faint pulsars with flux densities of only a few or a few tens of μJy at 1.25 GHz. To obtain accurate position, spin parameters and dispersion measure of a pulsar, and to calculate derived parameters such as the characteristic age and surface magnetic field, we collect available FAST pulsar data obtained either through targeted follow-up observations or through coincidental survey observations with one of the 19 beams of the receiver. From these data we obtain time of arrival (TOA) measurements for 30 newly discovered pulsars as well as for 13 known pulsars. We demonstrate that the TOA measurements acquired by the FAST from any beams of the receiver in any observation mode (e.g. the tracking mode or the snapshot mode) can be combined to get timing solutions. We update the ephemerides of 13 previously known pulsars and obtain the first phase-coherent timing results for 30 isolated pulsars discovered in the FAST GPPS Survey. Notably, PSR J1904+0853 is an isolated millisecond pulsar, PSR J1906+0757 is a disrupted recycled pulsar, and PSR J1856+0211 has a long period of 9.89 s that can constrain pulsar death lines. Based on these timing solutions, all available FAST data have been added together to obtain the best pulse profiles for these pulsars.

Dissecting the Gaia HR diagram – II. The vertical structure of the star formation history across the solar cylinder

ABSTRACT Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic plane within a cylinder centred on the Sun with a 200 pc radius and spanning 1.3 kpc above and below the Galaxy’s midplane. We quantify both the concentration of the more recent star formation in the Galactic plane, and the age-related increase in the scale height of the Galactic disc stellar component, which is well-described by power laws with indices ranging from 1/2 to 2/3. The vertically-integrated star formation rate falls from $(1.147 \pm 0.039)\times 10^{-8}\, \text{M}_\odot \, \text{yr}^{-1} \, \text{pc}^{-2}$ at earlier times down to $(6.2 \pm 3.0) \times 10^{-9}\, \text{M}_\odot \, \text{yr}^{-1} \, \text{pc}^{-2}$ at present times, but we find a significant peak of star formation in the 2–3 Gyr age bin. The total mass of stars formed per unit area over time is $118.7 \pm 6.2\, \text{M}_\odot \, \text{pc}^{-2}$, which is nearly twice the present stellar mass derived from kinematics within 1 kpc from the Galactic plane, implying a high degree of matter recycling in successive generations of stars. The method is then modified by adopting an age-dependent correlation between the SFH across the different slices, which results in less noisy and more symmetrical results without significantly changing the previously mentioned quantities. This appears to be a promising way to improve SFH recovery in external galaxies.

Improvement of Pan-STARRS Photometric Calibration with LAMOST and Gaia

In this work, we perform the recalibration of Pan-STARRS1 (PS1) photometry by correcting for position-dependent systematic errors using the spectroscopy-based stellar color regression (SCR) method, the photometry-based SCR method, and the Gaia BP/RP synthetic photometry method. We confirm the significant large-scale and small-scale spatial variation of the magnitude offsets for all the grizy filters. We show that the PS1 photometric calibration precisions in the grizy filters are around 5–7 mmag when averaged over 14′ regions. We note a much larger calibration error up to 0.04 mag in the Galactic plane, which is probably caused by the systematic errors of the PS1 magnitudes in crowded fields. The results of the three methods are consistent with each other within 1–2 mmag or better for all the filters. We provide 2D maps and a Python package to correct for position-dependent magnitude offsets of PS1, which can be used for high-precision investigations and as a reference to calibrate other surveys.

Revisiting Galactic Disk and Spiral Arms Using Open Clusters

We use the largest catalog of open clusters in the post-Gaia era to provide an observational view of the Galactic disk. By compiling physical parameters such as age, distance, and kinematic information, we investigate the spatial distribution of open clusters and revisit the spiral arms and other asymmetries in the Galactic disk. Using young open clusters as a tracer of spiral arms, we map the spiral structure of the Galaxy and find that most of the clusters start migrating away from the spiral arms in about 10–20 Myr and fill the interarm regions as they age. Using the 3D kinematic information on 371 open star clusters, we derive different individual pattern speeds for spiral arms that closely follow the rotation curve of the Milky Way, hence favoring the transient nature of spiral arms in the Milky Way. The pattern rotation speeds of each spiral arm suggest that the spiral arms have not accelerated in the last 80 Myr. Based on the distribution of open clusters younger than 700 Myr above or below the Galactic plane, we found a solar offset of z ⊙ = 17.0 ± 0.9 pc north of the Galactic plane and estimated the scale height z h = 91.7 ± 1.9 pc from the Galactic plane.

FAST Pulsar Database. I. Polarization Profiles of 682 Pulsars

Pulsar polarization profiles form a very basic database for understanding the emission processes in a pulsar magnetosphere. After careful polarization calibration of the 19-beam L-band receiver and verification of beam-offset observation results, we obtain polarization profiles of 682 pulsars from observations by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during the Galactic Plane Pulsar Snapshot survey and other normal FAST projects. Among them, polarization profiles of about 460 pulsars are observed for the first time. The profiles exhibit diverse features. Some pulsars have a polarization position angle curve with a good S-shaped swing, some with orthogonal modes; some have components with highly linearly polarized components or strong circularly polarized components; some have a very wide profile, coming from an aligned rotator, and some have an interpulse from a perpendicular rotator; some wide profiles are caused by interstellar scattering. We derive geometric parameters for 190 pulsars from the S-shaped position angle curves or with orthogonal modes. We find that the linear and circular polarization or the widths of pulse profiles have various frequency dependencies. Pulsars with a large fraction of linear polarization are more likely to have a large Edot.

Probabilistic classification of infrared-selected targets for SPHEREx mission: in search of young stellar objects

ABSTRACT We apply machine learning algorithms to classify infrared (IR)-selected targets for NASA’s upcoming Spectro-Photometer for the History of the Universe, Epoch of Reionization and Ices Explorer (SPHEREx) mission. In particular, we are interested in classifying young stellar objects (YSOs), which are essential for understanding the star formation process. Our approach differs from previous works, which have relied heavily on broad-band colour criteria to classify IR-bright objects, and are typically implemented in colour–colour and colour–magnitude diagrams. However, these methods do not state the confidence associated with the classification and the results from these methods are quite ambiguous due to the overlap of different source types in these diagrams. Here, we utilize photometric colours and magnitudes from seven near- and mid-IR bands simultaneously and employ machine and deep learning algorithms to carry out probabilistic classification of YSOs, asymptotic giant branch (AGB) stars, active galactic nuclei (AGNs), and main-sequence (MS) stars. Our approach also subclassifies YSOs into Class I, II, III, and flat spectrum YSOs, and AGB stars into carbon-rich and oxygen-rich AGB stars. We apply our methods to IR-selected targets compiled in preparation for SPHEREx which are likely to include YSOs and other classes of objects. Our classification indicates that out of 8308 384 sources, 1966 340 have class prediction with probability exceeding 90 per cent, amongst which $\sim 1.7~{{\ \rm per\ cent}}$ are YSOs, $\sim 58.2~{{\ \rm per\ cent}}$ are AGB stars, $\sim 40~{{\ \rm per\ cent}}$ are (reddened) MS stars, and $\sim 0.1~{{\ \rm per\ cent}}$ are AGNs whose red broad-band colours mimic YSOs. We validate our classification using the spatial distributions of predicted YSOs towards the Cygnus-X star-forming complex, as well as AGB stars across the Galactic plane.

The kinematics, metallicities, and orbits of six recently discovered Galactic star clusters with Magellan/M2FS spectroscopy

ABSTRACT We present Magellan/M2FS spectroscopy of four recently discovered Milky Way star clusters (Gran 3/Patchick 125, Gran 4, Garro 01, and LP 866) and two newly discovered open clusters (Gaia 9 and Gaia 10) at low Galactic latitudes. We measure line-of-sight velocities and stellar parameters ([Fe/H], log g, Teff, and [Mg/Fe]) from high-resolution spectroscopy centred on the Mg triplet and identify 20–80 members per star cluster. We determine the kinematics and chemical properties of each cluster and measure the systemic proper motion and orbital properties by utilizing Gaia astrometry. We find Gran 3 to be an old, metal-poor (mean metallicity of [Fe/H] = −1.83) globular cluster located in the Galactic bulge on a retrograde orbit. Gran 4 is an old, metal-poor ([Fe/H] = −1.84) globular cluster with a halo-like orbit that happens to be passing through the Galactic plane. The orbital properties of Gran 4 are consistent with the proposed LMS-1/Wukong and/or Helmi streams merger events. Garro 01 is metal-rich ([Fe/H] = −0.30) and on a near-circular orbit in the outer disc but its classification as an open cluster or globular cluster is ambiguous. Gaia 9 and Gaia 10 are among the most distant known open clusters at $R_{\mathrm{GC}}\sim 18,~21.2~\mathrm{\, kpc}$ and most metal-poor with [Fe/H] ∼−0.50, −0.34 for Gaia 9 and Gaia 10, respectively. LP 866 is a nearby, metal-rich open cluster ([Fe/H] = +0.10). The discovery and confirmation of multiple star clusters in the Galactic plane shows the power of Gaia astrometry and the star cluster census remains incomplete.

The FAST Galactic Plane Pulsar Snapshot survey – IV. Discovery of five fast radio bursts

ABSTRACT We report five new fast radio bursts (FRBs) discovered from the Galactic Plane Pulsar Snapshot (GPPS) survey by the Five-hundred-meter Aperture Spherical radio Telescope (FAST): FRB 20210126, FRB 20210208, FRB 20210705, FRB 20211005, and FRB 20220306. To date, no repeating bursts from these FRB sources have been detected in the follow-up monitoring observations, leading to their classification as potential one-off events. We obtain the basic parameters for these bursts, including position, dispersion measure, pulse width, spectral index, scattering time-scale, etc. The fluences and flux densities are generally lower in comparison to the values observed in one-off bursts discovered by other telescopes. Among the observed bursts, polarization data for four bursts were recorded during observations. Consequently, we obtain polarization profiles and Faraday rotation measures for these bursts.

EMU Detection of a Large and Low Surface Brightness Galactic SNR G288.8–6.3

We present the serendipitous detection of a new Galactic supernova remnant (SNR), G288.8–6.3, using data from the Australian Square Kilometre Array Pathfinder (ASKAP) Evolutionary Map of the Universe (EMU) survey. Using multifrequency analysis, we confirm this object as an evolved Galactic SNR at high Galactic latitude with low radio surface brightness and typical SNR spectral index of α = − 0.41 ± 0.12. To determine the magnetic field strength in SNR G288.8–6.3, we present the first derivation of the equipartition formulae for SNRs with spectral indices α > − 0.5. The angular size is 1.°8 × 1.°6 (107.′6 × 98.′4), and we estimate that its intrinsic size is ∼40 pc, which implies a distance of ∼1.3 kpc and a position of ∼140 pc above the Galactic plane. This is one of the largest in angular size and closest Galactic SNRs. Given its low radio surface brightness, we suggest that it is about 13,000 yr old.

Resonant Effects of a Bar on the Galactic Disk Kinematics Perpendicular to Its Plane

Detailed analysis of kinematics of the Milky Way disk in the solar neighborhood based on the GAIA DR3 catalog reveals the existence of peculiarities in the stellar velocity distribution perpendicular to the galactic plane. We study the influence of resonances—the outer Lindblad resonance and the outer vertical Lindblad resonance—of a rotating bar with stellar oscillations perpendicular to the plane of the disk, and their role in shaping the spatial and the velocity distributions of stars. We find that the Z and VZ distributions of stars with respect to LZ are affected by the outer Lindblad resonance. The existence of bar resonance with stellar oscillations perpendicular to the plane of the disk is demonstrated for a long (large semi-axis 5 kpc) and fast rotating bar with Ωb=60.0kms−1kpc−1. We show also that, in the model with the long and fast rotating bar, some stars in the 2:1 OLR region deviate far from their original places, entering the bar region. A combination of resonance excitation of stellar motions at the 2:1 OLR region together with strong interaction of the stars with the bar potential leads to the formation of the group of ‘escapees’, i.e., stars that deviate in R and Z—directions at large distances from the resonance region. Simulations, however, do not demonstrate any noticeable effect on VZ-distribution of stars in the solar neighborhood.

Face-on map of the molecular disc and 3-kpc expanding ring of the galaxy based on a high-accuracy rotation curve

We analyze the longitude-velocity diagram (LVD) of 12CO-line emission from archival data and use the most accurate rotation curve (RC) of the Milky Way to transform radial velocity to face-on position in the galactic plane. We point out that the face-on transformation is highly sensitive to the adopted RC, especially in the inner Milky Way, in the sense that deviations of the RC from the true rotation velocity lead either to an artifact hole or overcrowded concentration along the tangent circle for over- or under-estimated RC, respectively. Even if the RC is sufficiently accurate, non-circular motion such as with the 3 kpc expanding ring introduces significant artifacts in the resulting face-on-map, as long as a circular rotation is assumed. On the other hand, if we properly take into account the non-circular motion, it can be used to solve the near-far degeneracy problem of determination of kinematic distance. We thus propose a new method to solve the degeneracy by incorporating the expanding motion of a ring or arms. We apply the method to the LVD of the 3-kpc expanding ring and present its face-on map projected onto the galactic plane for the first time.

Search for directional associations between baikal gigaton volume detector neutrino-induced cascades and high-energy astrophysical sources

ABSTRACT Baikal-GVD has recently published its first measurement of the diffuse astrophysical neutrino flux, performed using high-energy cascade-like events. We further explore the Baikal-GVD cascade data set collected in 2018–2022, with the aim to identify possible associations between the Baikal-GVD neutrinos and known astrophysical sources. We leverage the relatively high angular resolution of the Baikal-GVD neutrino telescope (2–3 deg.), made possible by the use of liquid water as the detection medium, enabling the study of astrophysical point sources even with cascade events. We estimate the telescope’s sensitivity in the cascade channel for high-energy astrophysical sources and refine our analysis prescriptions using Monte-Carlo simulations. We primarily focus on cascades with energies exceeding 100 TeV, which we employ to search for correlation with radio-bright blazars. Although the currently limited neutrino sample size provides no statistically significant effects, our analysis suggests a number of possible associations with both extragalactic and Galactic sources. Specifically, we present an analysis of an observed triplet of neutrino candidate events in the Galactic plane, focusing on its potential connection with certain Galactic sources, and discuss the coincidence of cascades with several bright and flaring blazars.

Swift Deep Galactic Plane Survey classification of Swift J170800−402551.8 as a candidate intermediate polar cataclysmic variable

ABSTRACT Here, we present the results of our multiwavelength campaign aimed at classifying Swift J170800−402551.8 as part of the Swift Deep Galactic Plane Survey (DGPS). We utilized Target of Opportunity (ToO) observations with the Neil Gehrels Swift Observatory, Chandra X-ray Observatory, Neutron star Interior Composition Explorer (NICER), X-ray Multi-Mirror Mission (XMM–Newton), Nuclear Spectroscopic Telescope Array (NuSTAR), and the Southern African Large Telescope (SALT), as well as multiwavelength archival observations from Gaia, VST Photometric Hα Survey, and VISTA Variables in the Via Lactea. The source displays a periodicity of 784 s in our XMM–Newton observation. The X-ray spectrum (XMM–Newton and NuSTAR) can be described by thermal bremsstrahlung radiation with a temperature of kT ≈ 30 keV. The phase-folded X-ray light curve displays a double-peaked, energy-dependent pulse profile. We used Chandra to precisely localize the source, allowing us to identify and study the multiwavelength counterpart. Spectroscopy with SALT identified a Balmer H α line, and potential He i lines, from the optical counterpart. The faintness of the counterpart (r ≈ 21 AB mag) favours a low-mass donor star. Based on these criteria, we classify Swift J170800−402551.8 as a candidate intermediate polar cataclysmic variable, where the spin period of the white dwarf is 784 s.

Geometry and Kinematics of a Dancing Milky Way: Unveiling the Precession and Inclination Variation across the Galactic Plane via Open Clusters

This Letter presents a study of the geometry and motion of the Galactic disk using open clusters in the Gaia era. The findings suggest that the inclination θ i of the Galactic disk increases gradually from the inner to the outer disk, with a shift in orientation at the Galactocentric radius of approximately 6 ± 1 kpc. Furthermore, this study brings forth the revelation that the mid-plane of the Milky Way may not possess a stationary or fixed position. A plausible explanation is that the inclined orbits of celestial bodies within our Galaxy exhibit a consistent pattern of elliptical shapes, deviating from perfect circularity; however, more observations are needed to confirm this. An analysis of the vertical motion along the Galactocentric radius reveals that the disk has warped with precession and that the line of node shifts at different radii, aligning with the results from the classical Cepheids. Although there is uncertainty for precession/peculiar motion in solar orbit, after considering the uncertainty, the study derives a median value of = 6.8 km s−1 kpc−1 in the Galaxy. This value for the derived precession in the outer disk is lower than those in the literature due to the systematic motion in solar orbit (θ i = 0.°6). The study also finds that the inclinational variation of the disk is significant and can cause systematic motion, with the variation rate decreasing along the Galactic radius with a slope of −8.9 μas yr−1 kpc−1. Moreover, the derived in solar orbit is 59.1 ± 11.2sample ± 7.7 μas yr−1, which makes it observable for high-precision astrometry.

SPYGLASS. IV. New Stellar Survey of Recent Star Formation within 1 kpc

Young stellar populations provide a powerful record that traces millions of years of star formation history in the solar neighborhood. Using a revised form of the SPYGLASS young star identification methodology, we produce an expanded census of nearby young stars (age < 50 Myr). We then use the HDBSCAN clustering algorithm to produce a new SPYGLASS Catalog of Young Associations, which reveals 116 young associations within 1 kpc. More than 25% of these groups are largely new discoveries, as 20 are substantively different from any previous definition, and 10 have no equivalent in the literature. The new associations reveal a yet undiscovered demographic of small associations with little connection to larger structures. Some of the groups we identify are especially unique for their high transverse velocities, which can differ from the solar velocity by 30–50 km s−1, and for their positions, which can reach up to 300 pc above the galactic plane. These features may suggest a unique origin, matching existing evidence of infalling gas parcels interacting with the disk interstellar medium. Our clustering also suggests links between often-separated populations, hinting to direct structural connections between Orion Complex and Perseus OB2, and between the subregions of Vela. The ∼30 Myr old Cepheus-Hercules association is another emerging large-scale structure, with a size and population comparable to Sco-Cen. Cep-Her and other similarly aged structures are also found clustered along extended structures perpendicular to known spiral arm structure, suggesting that arm-aligned star formation patterns have only recently become dominant in the solar neighborhood.

QUIJOTE scientific results – X. Spatial variations of Anomalous Microwave Emission along the Galactic plane

ABSTRACT Anomalous microwave emission (AME) is an important emission component between 10 and 60 GHz that is not yet fully understood. It seems to be ubiquitous in our Galaxy and is observed at a broad range of angular scales. Here we use the new QUIJOTE-MFI wide survey data at 11, 13, 17, and 19 GHz to constrain the AME in the Galactic plane (|b| &amp;lt; 10°) on degree scales. We built the spectral energy distribution between 0.408 and 3000 GHz for each of the 5309 0.9° pixels in the Galactic plane, and fitted a parametric model by considering five emission components: synchrotron, free–free, AME, thermal dust and CMB anisotropies. We show that not including QUIJOTE-MFI data points leads to the underestimation (up to 50 per cent) of the AME signal in favour of free–free emission. The parameters describing these components are then intercompared, looking for relations that help to understand AME physical processes. We find median values for the AME width, WAME, and for its peak frequency, νAME, respectively of $0.560^{+0.059}_{-0.050}$ and $20.7^{+2.0}_{-1.9}$ GHz, slightly in tension with current theoretical models. We find spatial variations throughout the Galactic plane for νAME, but only with reduced statistical significance. We report correlations of AME parameters with certain ISM properties, such as that between the AME emissivity (which shows variations with the Galactic longitude) and the interstellar radiation field, and that between the AME peak frequency and dust temperature. Finally, we discuss the implications of our results on the possible molecules responsible for AME.

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.

Neutrinos from the Galactic plane

Neutrinos from the Galactic plane