Modified equipartition calculation for supernova remnants

Determination of the magnetic field strength in the interstellar medium is one of the most complex tasks of contemporary astrophysics. We can only estimate the order of magnitude of the magnetic field strength by using a few very limited methods. Besides Zeeman effect and Faraday rotation, the equipartition or the minimum-energy calculation is a widespread method for estimating magnetic field strength and energy contained in the magnetic field and cosmic ray particles by using only the radio synchrotron emission. Despite of its approximate character, it remains a useful tool, especially when there is no other data about the magnetic field in a source. In this paper we give a modified calculation which we think is more appropriate for estimating magnetic field strengths and energetics in supernova remnants (SNRs). Finally, we present calculated estimates of the magnetic field strengths for all Galactic SNRs for which the necessary observational data are available. The web application for calculation of the magnetic field strength of SNRs is available at http://poincare.matf.bg.ac.rs/~arbo/eqp/.

💡 Research Summary

The paper presents a revised equipartition (or minimum‑energy) method specifically tailored for supernova remnants (SNRs), addressing the shortcomings of the classical Pacholczyk (1970) formulation and the later Beck & Krause (2005) revision. The authors begin by reviewing the traditional approach, which estimates the total energy of relativistic particles by integrating a power‑law radio spectrum over frequency and assumes a fixed energy ratio between protons and electrons (k‑3) together with a homogeneous magnetic field. While widely used, this method neglects the detailed particle acceleration physics that governs SNRs and ignores the contribution of heavier ions.

To overcome these limitations, the authors adopt Bell’s (1978) theory of diffusive shock acceleration (DSA). They assume that all accelerated particles are injected at a common energy (E_{\rm inj}\approx \frac12 m_p v_s^2), where (v_s) is the shock velocity. The particle momentum distribution is taken as a pure power law (f(p)=k,p^{-\gamma}) with (2<\gamma<3). This single‑power‑law description naturally yields a smooth transition from the non‑relativistic regime (spectral index (\gamma)) to the ultra‑relativistic regime (spectral index (\gamma+1)) without the artificial break imposed in the BK05 scheme.

The total cosmic‑ray energy density (\varepsilon_{\rm CR}) is derived by integrating over momentum rather than frequency. The integration leads to an expression involving the Gamma function and a dimensionless integral (I(x)) (equations 1–3). Crucially, the authors introduce a new parameter (\kappa) that quantifies the energy ratio of ions to electrons, explicitly incorporating elemental abundances, mass‑to‑charge ratios, and the injection energy. In the limit (E_{\rm inj}\ll m_ec^2) the parameter reduces to the familiar constant used in earlier works, confirming consistency.

The synchrotron emissivity (\varepsilon_\nu) is expressed in terms of the electron normalization constant (K_e), magnetic field strength (B), and the averaged pitch‑angle factor (\langle\sin\Theta\rangle). By relating the observed flux density (S_\nu) to (\varepsilon_\nu) through the source volume (characterized by radius (R), filling factor (f), distance (d), and angular radius (\theta)), the authors obtain two coupled equations: one for the synchrotron emission and another for the total energy (magnetic plus particle). Minimizing the total energy with respect to (B) yields an explicit formula for the magnetic field (equation 11). This formula has the same functional dependence on observable quantities as the BK05 result but contains the modified (K_e) and (\kappa) that reflect the underlying DSA physics.

Because the integral (I(x)) can be written in terms of the hypergeometric function (2F_1), the authors develop a practical approximation (equations 18–20) that reproduces the exact result to within 3 % across the entire range of (x=p{\rm inj}/(mc)). This enables a fast, accurate computation of (\kappa) as a function of the injection energy, which is itself derived from the measured shock velocity.

The theoretical development is complemented by a web‑based calculator implemented in PHP (available at http://poincare.matf.bg.ac.rs/~arbo/eqp/). Users input the radio spectral index (or (\gamma)), observing frequency, flux density, distance, angular size, filling factor, shock speed, and elemental abundances. The tool then returns the magnetic field strength and the minimum total energy. Default parameters assume typical SNR values (spectral index 0.5–1.0, filling factor (f=1) or derived from a shell thickness (\delta), and ISM abundances H:He = 10:1). The Gamma function is evaluated using Nemes’ (2010) approximation to keep the computation lightweight.

In the discussion, the authors emphasize that strict equipartition ((\varepsilon_{\rm CR}=\varepsilon_B)) need not be imposed; any constant ratio (\beta=\varepsilon_B/\varepsilon_{\rm CR}) can be used if independent constraints (e.g., X‑ray inverse‑Compton, gamma‑ray pion decay) are available. This flexibility allows the method to be applied to both mature SNRs (low shock speeds, (\alpha) between 0.5 and 1) and younger, faster remnants, because the injection‑energy dependence is explicitly accounted for. Figures in the paper demonstrate that the new (\kappa(E_{\rm inj})) matches Bell’s original calculations and that, for low (E_{\rm inj}), the simplified constant‑(\kappa) approximation (equation 3) is justified.

Overall, the paper delivers a physically motivated, mathematically robust, and user‑friendly equipartition framework for SNRs. By integrating DSA theory, heavy‑ion contributions, and a smooth momentum‑space treatment, it provides more reliable magnetic‑field estimates than previous methods, while the accompanying web tool makes the approach readily accessible to the broader astrophysical community. This advancement is expected to improve our understanding of particle acceleration, magnetic‑field amplification, and energy budgets in supernova remnants across the Galaxy.