---
title: 'G118.4+37.0 (Calvera’s SNR): Multiwavelength Study'
url: https://www.emergentmind.com/topics/g118-4-37-0-calvera-s-snr
type: topic
---

# G118.4+37.0 (Calvera’s SNR): Multiwavelength Study

Searching arXiv for recent papers on G118.4+37.0 / Calvera’s SNR to ground the article in the current literature.
G118.4+37.0, commonly referred to as Calvera’s SNR, is a candidate supernova remnant at high Galactic latitude, spatially associated with the X-ray pulsar Calvera, 1RXS J141256.0+792204. It was first reported as a faint, nearly circular radio ring in LOFAR Two-metre Sky Survey data, then identified as an extended GeV source in *Fermi*-LAT observations, and subsequently examined through deeper XMM-Newton and optical H$\alpha$ follow-up. Across these studies, the source is characterized by a $\sim 1^\circ$-scale shell-like morphology, non-thermal radio emission, extended $\gamma$-ray emission, and diffuse thermal X-rays interior to the ring, although its detailed emission mechanism and evolutionary state remain under active discussion [2207.14141; 2210.16340; 2211.10016; 2507.13210].

## 1. Discovery and basic identification

G118.4+37.0 was reported as a ring of low surface brightness radio emission around the Calvera pulsar in the LOFAR Two-metre Sky Survey (LoTSS) [2207.14141]. The ring is centered at $\alpha=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $\delta=+79^\circ23'15''$, corresponding to Galactic coordinates $l = 118.41^\circ$, $b = +37.03^\circ$, with inner and outer radii of $14.2'$ and $28.4'$ [2207.14141]. Its high latitude is central to its interpretation, since a confirmed association with Calvera would place the remnant far above the Galactic plane and make it one of the few known Galactic halo SNRs [2207.14141].

The designation G118.4+37.0 encodes the approximate Galactic coordinates and is used in the literature interchangeably with “Calvera’s SNR” [2210.16340; 2211.10016]. The radio ring is very round and of relatively uniform brightness, apart from a brighter, straight filament in the north-west [2207.14141]. Its angular diameter is approximately $0.95^\circ$, consistent with the later GeV extension measured with *Fermi*-LAT [2211.10016].

The discovery paper considered three interpretations for the ring: an H II region, a supernova remnant, or an Odd Radio Circle [2207.14141]. The authors preferred the SNR interpretation because of the positional coincidence of the ring, the pulsar, and an X-ray-emitting non-equilibrium ionisation plasma previously detected in the field [2207.14141]. Later GeV detection of extended emission consistent with the radio morphology strengthened that interpretation by confirming the presence of relativistic particles [2210.16340].

## 2. Radio properties and early multi-wavelength diagnostics

The radio ring was detected at 144 MHz in LoTSS with an integrated flux density of $1.08 \pm 0.15$ Jy, and also measured in WENSS at 325 MHz with $0.63 \pm 0.21$ Jy and in NVSS at 1.4 GHz with $0.22 \pm 0.04$ Jy [2207.14141]. A power-law fit to these three flux densities gave a radio spectral index $\alpha = -0.71 \pm 0.09$ for $S_\nu \propto \nu^\alpha$ [2207.14141]. In the later $\gamma$-ray modeling paper, the same radio behavior is described as a “flat” integrated radio spectrum with $\alpha_{\rm radio} = 0.71 \pm 0.09$ under the convention $S_\nu \propto \nu^{-\alpha_{\rm radio}}$ [2211.10016]. Both notations encode the same non-thermal synchrotron slope.

The mean surface brightness at 144 MHz was reported as $\Sigma_{144} \approx 6.7\times10^{-23}\ {\rm W\,m^{-2}\,Hz^{-1}\,sr^{-1}}$, and the corresponding 1 GHz surface brightness as $\Sigma_{1\,{\rm GHz}} \approx 1.74\times10^{-23}\ {\rm W\,m^{-2}\,Hz^{-1}\,sr^{-1}}$ [2207.14141]. These values place the object among low-surface-brightness remnants and were used for comparison with empirical $\Sigma$–$D$ relations [2207.14141]. The mean brightness temperature at 144 MHz was estimated as $\sim 11$ K [2207.14141].

No polarized emission associated with the ring was detected in LoTSS RM-synthesis analysis over Faraday depths of $\pm 100\ {\rm rad\ m^{-2}}$, with a polarized-intensity limit of $300\ \mu{\rm Jy\ beam^{-1}}$ at $3'$ resolution [2207.14141]. The resulting fractional polarization constraint is $p \lesssim 5\%$ on those scales, although strong low-frequency depolarization is expected [2207.14141]. This non-detection did not rule out an SNR interpretation.

The early optical and X-ray context was mixed but suggestive. No shell-like H$\alpha$ emission was detected around the radio ring in INT/WFC imaging, with a $5\sigma$ surface-brightness limit of $1.1\times10^{-16}\ {\rm erg\ cm^{-2}\ s^{-1}\ arcsec^{-2}}$ [2207.14141]. A small V-shaped H$\alpha$ “smudge,” about $23''$ long, was detected inside the ring at $\alpha = 14\mathrm{h}08\mathrm{m}28.3\mathrm{s}$, $\delta = +79^\circ22'54.0''$, but its relation to the radio shell was uncertain [2207.14141]. In X-rays, no diffuse shell emission was seen in the ROSAT All-Sky Survey, yet an interior extended source, 1RXS J140818.1+792113, had previously been modeled as a non-equilibrium ionization plasma with oxygen overabundance, a property described as a strong SNR diagnostic [2207.14141].

These combined diagnostics disfavored the H II region and Odd Radio Circle interpretations. The H II hypothesis was found inconsistent with the steep radio spectrum, the stringent H$\alpha$ non-detection, the lack of WISE emission, and the absence of a plausible ionizing star [2207.14141]. The ORC interpretation was disfavored by the $\sim 1^\circ$ angular size, which is much larger than the arcminute scales typical of known ORCs, and by the lack of an obvious host galaxy [2207.14141].

## 3. GeV detection and spatial morphology

Extended GeV emission from the region was reported independently in two 2022 *Fermi*-LAT studies, both based on long integrations of Pass 8 data [2210.16340; 2211.10016]. One study analyzed almost 14 years of observations and detected extended GeV emission consistent with the size and location of the radio source, arguing that these features and its similarities to other isolated SNRs establish the source as the remnant of a supernova [2210.16340]. The other used 13.6 years of Pass 8 data and provided detailed localization and extension tests [2211.10016].

In the 13.6-year analysis, the preferred morphology was a uniform disk [2211.10016]. The point-source localization in the 1–1000 GeV band gave $\mathrm{R.A.} = 212.454 \pm 0.021$ deg and $\mathrm{Dec.} = 79.344 \pm 0.028$ deg, while the best-fit disk had centroid $\mathrm{R.A.} = 212.649 \pm 0.038$ deg, $\mathrm{Dec.} = 79.389 \pm 0.038$ deg, with disk radius $R = 0.478$ deg and corresponding $R_{68} = 0.392 \pm 0.025$ deg [2211.10016]. A 2D Gaussian fit yielded $R_{68} = 0.394^{+0.068}_{-0.050}$ deg, and a uniform ring model fixed to the radio shell used inner radius $14'.2$ and outer radius $28'.4$ [2211.10016]. The uniform disk provided the best likelihood among the tested templates [2211.10016].

The same study reported a point-source TS of 23.8 and a uniform-disk TS of 54.1 in the 1–1000 GeV band, corresponding to an approximately $6.4\sigma$ detection [2211.10016]. The extension significance was quantified as
$$
TS_{\rm ext} = -2[\log L(\mathrm{point})-\log L(\mathrm{extended})] = 30.0,
$$
equivalent to an approximately $5.5\sigma$ preference for an extended source [2211.10016]. The GeV extension was described as closely matching the diffuse radio ring seen by LoTSS, WENSS, and NVSS [2211.10016].

A later reanalysis with 17 years of LAT data shifted the emphasis from single-template fitting to morphology on sub-degree scales [2507.13210]. Using a $7.5^\circ$ radius ROI centered at $\mathrm{RA}=212.8^\circ$, $\mathrm{Dec}=+79.4^\circ$, with 0.05° pixels and 10 energy bins per decade, that study found patchy extended residuals over $\approx 1^\circ$, with TS peaks at $3$–$4.5\sigma$ ($TS \approx 9$–20) after removing 4FGL J1409.8+7921 from the background model [2507.13210]. The highest TS peak was reported as colocated with the brightest soft X-rays and the isolated H$\alpha$ filament [2507.13210]. This suggests spatial non-uniformity in the high-energy emissivity rather than a perfectly uniform shell.

## 4. Spectral properties and emission scenarios

The *Fermi*-LAT spectrum between 100 MeV and 1 TeV has been modeled with both a single power law and a log-parabola [2211.10016]. In that analysis, the single power law,
$$
\frac{dN}{dE}\propto E^{-\alpha},
$$
gave $TS_{\rm PL} = 44.4$, photon index $\alpha = 1.67 \pm 0.10$, and integrated photon flux
$$
F_\gamma = (6.16 \pm 2.10)\times10^{-10}\ {\rm ph\ cm^{-2}\ s^{-1}}
$$
over 100 MeV–1 TeV [2211.10016]. The log-parabola,
$$
\frac{dN}{dE} = N_0 (E/E_0)^{-\alpha-\beta\log(E/E_0)},
$$
gave $TS_{\rm LogPb} = 55.8$, with $\alpha = 2.05 \pm 0.27$, $\beta = 0.31 \pm 0.14$, and integrated photon flux
$$
F_\gamma = (1.45 \pm 0.52)\times10^{-10}\ {\rm ph\ cm^{-2}\ s^{-1}}
$$
[2211.10016]. The improvement $\Delta TS = 11.4$ relative to the power law corresponds to $\sim 3.4\sigma$ evidence for spectral curvature concentrated at several tens of GeV [2211.10016].

The spectral energy distribution was constructed in six logarithmically spaced bins between 0.1 and 1000 GeV, with 95% upper limits assigned where $TS<5$ [2211.10016]. The highest-energy bins were reported as upper limits, reinforcing a turnover or softening at hundreds of GeV [2211.10016]. Assuming $d=3.3$ kpc, the 10–316 GeV luminosity was estimated as $L_\gamma(10$–$316\ {\rm GeV}) \approx 1.75\times10^{33}\ {\rm erg\ s^{-1}}$ [2211.10016].

Broadband modeling has not yielded a unique emission mechanism. In the leptonic scenario, synchrotron emission explains the radio flux while inverse Compton scattering on the cosmic microwave background, plus infrared and optical fields, accounts for the $\gamma$ rays; bremsstrahlung is negligible at the adopted halo density [2211.10016]. The electron distribution was parameterized as
$$
N_e(E) = K_e E^{-p}\exp(-E/E_{\rm cut,e}),
$$
with target photon fields including the CMB, an infrared component with $T_1 = 30$ K and $u_1 = 0.1\ {\rm eV\ cm^{-3}}$, and an optical component with $T_2 = 5000$ K and $u_2 = 0.1\ {\rm eV\ cm^{-3}}$ [2211.10016].

Two representative leptonic fits were highlighted [2211.10016]. A fit with $p=2.4$, $E_{\rm cut,e}=100$ TeV, $B\approx1.8\ \mu{\rm G}$, and $W_e(E>1\ {\rm GeV}) \approx 6.3\times10^{47}\ {\rm erg}$ reproduces the radio data but overpredicts low-energy $\gamma$ rays unless finely tuned [2211.10016]. A harder fit with $p=2.0$, $E_{\rm cut,e}\approx10$ TeV, $B\approx6\ \mu{\rm G}$, and $W_e(E>1\ {\rm GeV}) \approx 1.8\times10^{47}\ {\rm erg}$ better matches the hard GeV spectrum and respects the X-ray upper limits, but is in tension with the radio slope, since $p=2.0$ would imply $\alpha_{\rm radio}\approx0.5$ rather than the observed $\alpha_{\rm radio}\approx0.71$ [2211.10016]. The same study emphasized that the combination of a flat radio spectrum and a hard GeV spectrum is difficult to reconcile with a simple one-zone leptonic model [2211.10016].

The hadronic scenario parameterized the proton population as
$$
N_p(E)=K_p E^{-s}\exp(-E/E_{\rm cut,p}),
$$
with $\gamma$ rays produced by neutral-pion decay in $pp$ interactions [2211.10016]. To reproduce the hard GeV spectrum, the proton index had to be very hard, $s \lesssim 1.6$, with best fit $s\approx1.6$ and $E_{\rm cut,p}\lesssim20$ TeV [2211.10016]. For the halo density estimate $n_{\rm gas} \approx 4\times10^{-4}\ {\rm cm^{-3}}$, the required proton energy was
$$
W_p(E>1\ {\rm GeV}) \approx 4.6\times10^{52}\left(\frac{n_{\rm gas}}{4\times10^{-4}\ {\rm cm^{-3}}}\right)^{-1}\ {\rm erg},
$$
more than an order of magnitude above a canonical supernova kinetic energy of $10^{51}$ erg [2211.10016]. A broken-power-law proton spectrum still required $W_p \approx 6.4\times10^{52}(n_{\rm gas}/4\times10^{-4}\ {\rm cm^{-3}})^{-1}\ {\rm erg}$ [2211.10016]. On that basis, the hadronic model was judged energetically challenging unless the ambient density had been underestimated [2211.10016].

The 2025 XMM-Newton and LAT study altered this balance by inferring much larger local densities from thermal X-rays [2507.13210]. That work argued that the measured ambient density, together with the patchy morphology of the $\gamma$-ray emission and the detection of H$\alpha$ filaments, indicates that a hadronic origin is compatible with the $\gamma$-ray flux, though a mixed leptonic-hadronic scenario cannot be excluded [2507.13210]. This suggests that the viability of hadronic emission depends strongly on whether the remnant is expanding into a smooth halo medium or into a structured circumstellar environment.

## 5. X-ray plasma, optical filaments, and environmental structure

New XMM-Newton observations obtained in 2024 provided the first detailed characterization of the diffuse X-ray emission inside the ring [2507.13210]. Count-rate, O VII equivalent-width, and median photon energy maps isolated three diffuse soft regions—Central, North, and South—plus a harder “Clump” northeast of the Central region [2507.13210]. The diffuse emission was fitted in XSPEC with TBabs$\times$vcie and TBabs$\times$vnei models, grouped to 25 counts per bin, with two local background choices to bracket systematics [2507.13210].

For the fiducial Central region, the CIE fits gave $kT = 0.148\pm0.007$ keV for one background choice and $kT = 0.151\pm0.005$ keV for the other; across the bright region, $kT \approx 0.14$–0.16 keV was found to be robust [2507.13210]. Oxygen abundance was mildly sub-solar, with $O = 0.76^{+0.15}_{-0.11}$ or $0.61^{+0.08}_{-0.07}$ depending on background, while nitrogen was sub-solar but weakly constrained [2507.13210]. The best-fit thermal normalization was $(7.0\pm1.2)\times10^{-4}\ {\rm cm^{-5}}$ or $(8.3\pm1.1)\times10^{-4}\ {\rm cm^{-5}}$, and the absorbed 0.3–1 keV flux was $(4.8\pm0.3)\times10^{-13}\ {\rm erg\ s^{-1}\ cm^{-2}}$ [2507.13210]. NEI models did not improve the fits, although NEI could not be excluded because of a $kT$–$\tau$ degeneracy [2507.13210].

The authors derived the electron density from the XSPEC thermal normalization,
$$
K = \frac{10^{-14}}{4\pi d^2}\int n_e n_H\,dV,
$$
assuming $n_e \approx 1.2 n_H$ and an ellipsoidal geometry for the Central region [2507.13210]. For $d=4.5$ kpc, filling factor $f=1$, and $K=7.0\times10^{-4}\ {\rm cm^{-5}}$, they obtained $n_e \approx 0.28\ {\rm cm^{-3}}$ [2507.13210]. Allowing for the fitted normalization range, distances of 4–5 kpc, and $f\le1$ led to $n_e \approx 0.1$–$0.7\ {\rm cm^{-3}}$, corresponding to pre-shock density $n_0 \approx 0.02$–$0.15\ {\rm cm^{-3}}$ for compression ratio 4 [2507.13210]. These values are orders of magnitude larger than the smooth-halo estimate used in the earlier $\gamma$-ray modeling [2211.10016; 2507.13210].

The “Clump” region is spectrally harder than the diffuse soft component [2507.13210]. It can be fitted either by a thermal model with $kT \approx 1.6$–1.7 keV or by a power law with $\Gamma \approx 2.7 \pm 0.5$, with statistically comparable fit quality [2507.13210]. Its location close to the LAT TS maximum and to the isolated H$\alpha$ smudge was interpreted as suggestive of a site where the shock meets denser material or amplified turbulence [2507.13210].

Optical imaging with TNG/DOLORES reached a $5\sigma$ H$\alpha$ surface-brightness sensitivity of $2.5\times10^{-17}\ {\rm erg\ cm^{-2}\ s^{-1}\ arcsec^{-2}}$ [2507.13210]. No H$\alpha$ filaments were detected along the outer radio shell, but the previously noted smudge was recovered as an approximately $1'$-long filamentary complex with two bright heads and faint diffuse emission on the border of the X-ray bright patch [2507.13210]. Additional faint filaments were present nearby, with no continuum counterparts [2507.13210]. The morphology was described as clumpy rather than a smooth Balmer-dominated arc, which is consistent with radiative shocks in denser pockets, although spectroscopy is required for a definitive classification [2507.13210].

This X-ray and optical evidence supports a structured medium rather than a uniform halo. The 2025 study concluded that the remnant is expanding in a tenuous environment but has encountered a denser phase, likely the relic of wind activity from a massive progenitor star [2507.13210]. A plausible implication is that different sectors of the shell may sample different densities, magnetic fields, and radiative regimes, which would naturally complicate single-zone broadband modeling.

## 6. Distance, age, relation to Calvera, and broader significance

The relation between G118.4+37.0 and the Calvera pulsar is central to the source’s astrophysical significance. In the discovery paper, the pulsar lay $4.9'$ from the geometric center of the radio ring [2207.14141]. At the NICER-based distance estimate of $d \approx 3.3$ kpc for Calvera, this corresponds to a projected offset of about 4.7 pc, while the ring has outer radius about 27.3 pc and diameter about 54.5 pc [2207.14141]. The pulsar has period $P = 59.2$ ms, spin-down rate $\dot P = (3.19 \pm 0.08)\times10^{-15}\ {\rm s\ s^{-1}}$, characteristic age $\tau_c \approx 294$ kyr, spin-down power $\dot E = 6.1\times10^{35}\ {\rm erg\ s^{-1}}$, and surface dipole field $B_s \approx 4.4\times10^{11}$ G; it is radio-quiet and no confirmed LAT pulsations have been detected [2211.10016].

The distance and age of the SNR have evolved across the literature. The early work adopted $d \approx 3.3$ kpc, inherited from Calvera’s X-ray atmosphere modeling, and used a Sedov scaling
$$
R(t) \approx 12.9 (E_{51}/n_0)^{1/5} t_4^{2/5}\ {\rm pc}
$$
with $E_{51}=1$ and $n_0 = 4\times10^{-4}\ {\rm cm^{-3}}$ [2207.14141]. Using what that paper stated as an “SNR radius of 54 pc,” the inferred age was $t \approx 7700$ yr [2207.14141]. The same synthesis also noted that elsewhere the paper quoted 54 pc as a diameter, so the 7.7 kyr estimate depends on treating 54 pc as the shock radius [2207.14141]. In the 2022 $\gamma$-ray analysis, use of Sedov-Taylor arguments and low $n_{\rm gas}$ led instead to $t_{\rm age} \approx 1.4\times10^3$ yr, explicitly noting sensitivity to the density assumption [2211.10016].

The 2025 XMM-Newton study derived a different physical picture [2507.13210]. From the post-shock temperature relation
$$
kT_s = \frac{3}{16}\mu m_p v_s^2,
$$
the measured $kT \approx 0.15$ keV implies $v_s \approx 350\ {\rm km\ s^{-1}}$, with range $\approx 340$–370 km s$^{-1}$ for $kT = 0.14$–0.16 keV [2507.13210]. Using the observed angular radius $\theta \approx 28.4'$ and distance 4–5 kpc gives a physical radius $R \approx 33$–41 pc [2507.13210]. A naïve Sedov age from $v_s = (2/5)R/t$ then yields $\approx 35$–45 kyr, but the authors emphasized that the measured $kT$ pertains to a denser interior circumstellar patch rather than necessarily to the forward shock, so the true blast-wave speed is likely larger and the true age smaller [2507.13210]. Solving the Sedov relation with the X-ray-derived density range and canonical explosion energy gave a favored age of $t \approx 10$–20 kyr and a preferred distance of 4–5 kpc [2507.13210].

This later distance estimate remains consistent within uncertainties with that inferred for Calvera [2507.13210]. The pulsar’s proper motion was quoted as $\mu = 78.5 \pm 2.9\ {\rm mas\ yr^{-1}}$, directed away from the ring center [2507.13210]. A minimal kinematic estimate using angular offset $\Delta \theta \lesssim 5'$ and age 10–20 kyr gives
$$
v_t \approx 430 \left(\frac{d}{4.5\ {\rm kpc}}\right)\left(\frac{\Delta\theta}{5'}\right)\left(\frac{15\ {\rm kyr}}{t}\right)\ {\rm km\ s^{-1}},
$$
which was described as compatible with observed neutron-star kicks [2507.13210]. Using the directly measured proper motion gives
$$
v_t = 4.74\mu d \approx 1230(d/3.3\ {\rm kpc})\ {\rm km\ s^{-1}},
$$
so $v_t \approx 1500$–1800 km s$^{-1}$ for $d\approx4$–5 kpc, large but not unprecedented [2507.13210]. The association therefore remains plausible but not yet dynamically closed.

In the broader SNR context, G118.4+37.0 has been discussed as atypical. One 2022 study proposed that it and similar isolated remnants could belong to a radio-dim SNR population evolving in low-density environments and showing hard GeV emission of leptonic origin [2210.16340]. Another emphasized that its evident GeV curvature and absence of non-thermal X-ray emission make it an interesting source bridging young-aged SNRs with bright non-thermal X-ray emission and old-aged SNRs interacting with molecular clouds [2211.10016]. The 2025 multi-wavelength study instead classified it as a middle-aged remnant expanding in a tenuous medium but interacting with denser clumps, with hadronic $\gamma$-ray production becoming viable once the higher X-ray-derived density is taken into account [2507.13210].

The main unresolved issues follow directly from these differing inferences. The detailed extended-source $\gamma$-ray spectrum remains uncertain in the more recent LAT reanalysis, CIE and NEI X-ray descriptions are not decisively distinguished, and H$\alpha$ imaging alone cannot determine whether the filaments are Balmer-dominated or radiative [2507.13210]. Even so, the combination of a $\sim 1^\circ$ radio shell, extended GeV emission, soft thermal X-rays interior to the ring, and close association with Calvera establishes G118.4+37.0 as a significant laboratory for studying SNR evolution, particle acceleration, and neutron-star birth environments at high Galactic latitude [2207.14141; 2210.16340; 2211.10016; 2507.13210].

Source: https://www.emergentmind.com/topics/g118-4-37-0-calvera-s-snr