- The paper reconciles SN 2011dh’s optical and radio data by replacing the immediate self-similar shock model with a piston phase lasting about 50 days and adopting an ejecta density profile of n = 6.
- It finds that an inhomogeneous synchrotron source with a covering factor near 0.5 explains the observations, while increasing the inferred relativistic-electron-to-magnetic energy ratio by roughly 500 compared with a homogeneous model.
- The analysis supports magnetic-field evolution of B ∝ 1/R rather than B ∝ 1/t and suggests that declining covering factors can explain the radio peak-flux decrease without requiring a separate late-time emission component.
Reconciling conflicting inferences from optical and radio observations
SN 2011dh occupies a privileged position among radio supernovae because it was observed extensively and quasi-simultaneously in both the optical and radio regimes during its first ~100 days, and was later resolved by VLBI at t=453 days. This paper by Björnsson exploits that combination to expose a tension in the standard interpretation of the radio data. Under the conventional model — a homogeneous, spherically symmetric synchrotron source evolving according to the self-similar solution of Chevalier (the CN-phase), valid once initial conditions are dynamically irrelevant — the radio observations implied an ejecta density profile ρ∝v−n with n≈28 and a forward-shock velocity vb≈2.0×104 km s−1. The optical spectra, by contrast, were best reproduced with a shallow gradient n=6 (by both Marion et al. and Jerkstrand et al.), with Hα absorption extending to at least 2.5×104 km s−1 and possibly to 3.0×104 km s−1 — velocities substantially exceeding those inferred from the radio.
The paper resolves this discrepancy with two modifications to the standard model. First, the early evolution is governed not by the CN-phase but by the piston phase of Hamilton & Sarazin, in which the forward shock moves at constant velocity ρ∝v−n0; the transition to the CN-phase occurs at ρ∝v−n1 days. The analytical unification of the two phases follows Truelove & McKee, with the shock evolution given by ρ∝v−n2, where ρ∝v−n3 and ρ∝v−n4. Second, the radio source is inhomogeneous, characterized by a covering factor ρ∝v−n5.
The covering factor and its consequences for source parameters
The inhomogeneous-source formalism treats the source as a superposition of subcomponents, each a homogeneous source with radius ρ∝v−n6 and covering factor ρ∝v−n7. Observables then constrain only the combination ρ∝v−n8, where ρ∝v−n9 encodes the ratio of relativistic-electron to magnetic energy densities. The consequences are stark: the deduced magnetic field scales as n≈280, and the internal parameter as n≈281. With n≈282, a covering factor of 50% increases n≈283 by roughly a factor of 500 relative to the homogeneous assumption — the paper's strongest quantitative claim, and one that inverts the usual bias: assuming homogeneity artificially overestimates n≈284 and underestimates n≈285 (the opposite of the bias that arises when inverse Compton X-ray luminosities, rather than VLBI radii, provide the third observable).
The SN 1993J example makes the point concretely. The homogeneous analysis of Fransson & Björnsson yielded n≈286, implying a magnetic field so strong that the required post-shock thermal energy — and hence an ejecta kinetic energy of n≈287 erg — exceeds plausible explosion models by an order of magnitude. Adopting standard subcomponent assumptions (equipartition, no cooling, n≈288) gives n≈289, vb≈2.0×1040, and a magnetic field lower by a factor 4.8, consistent with earlier independent arguments that cooling ceases by ~100 days.
Constraining the shock evolution and covering factor of SN 2011dh
Normalizing the unified shock solution to the VLBI radius on day 453 (vb≈2.0×1041 km svb≈2.0×1042) and requiring consistency with the optical maximum ejecta velocity (vb≈2.0×1043 km svb≈2.0×1044, so vb≈2.0×1045) yields vb≈2.0×1046 and vb≈2.0×1047 days for vb≈2.0×1048 (16 days for vb≈2.0×1049). With a shallow gradient −10, the shock decelerates substantially, −11, which is precisely what allows the early radio radii to be reconciled with the optical velocities once the homogeneous-source assumption is dropped. The flat-topped spectra reported by Krauss et al. — whose fits required inflating error bars by factors 3–7 and systematically underestimated the spectral maximum by ~10% — independently indicate overlapping subcomponents.
Two caveats qualify these numbers. The optical absorption velocity is a lower limit: the forward shock should exceed the reverse-shock ejecta velocity by ~10%, and X-ray emission from the reverse shock likely ionizes the highest-velocity hydrogen, so −12 is an upper limit and −13 a lower bound; the −14 evolution constrains −15 to 4.0–5.0. The paper also notes that differences between the Horesh et al. and Krauss et al. radii may stem from fitting criteria applied to an intrinsically inhomogeneous spectrum rather than from physical variation.
The magnetic field scaling and the declining peak flux
Combining the observed near-constant peak spectral flux during the first ~100 days with the deduced −16 throughout the observing period implies −17, in contrast to the commonly assumed −18 (which follows if a fixed fraction of post-shock thermal energy amplifies the field). The shallow density gradient is what makes the distinction testable: under −19 and n=60, the peak flux should increase monotonically, contrary to both the initial plateau and the later decline. The n=61 scaling matches that previously established observationally for SN 1993J, and the paper argues that the extended plateaus of n=62 common to many radio supernovae suggest the same scaling applies generally.
The decline of n=63 by a factor 1.6 after ~100 days is attributed mainly to a decreasing covering factor (from ~0.45 to ~0.28, assuming n=64 matches the SN 1993J value), because n=65 persists unaffected. The alternative — constant n=66 with a factor-60 drop in n=67 — would require a n=68 larger by 2.6 and n=69 smaller by 23, which the paper deems less straightforward, though it concedes such structural changes cannot be excluded. This declining covering factor also offers an alternative to the de Witt et al. interpretation of the concave break in the 8.4 GHz light curve as a new component: a two-phase parameterization with 2.5×1040 days and 2.5×1041 reproduces both the steep initial decline and the later flattening, though the fit underestimates flux beyond ~1000 days, suggesting the late-phase peak flux is not strictly constant. The delayed light-curve peaks at 1.3 and 0.61 GHz (days 174 and 323) relative to the times those frequencies correspond to 2.5×1042 are consistent with a decreasing 2.5×1043 plus inhomogeneity.
Similarities to SN 1993J and broader implications
The paper draws a close parallel between SN 2011dh and SN 1993J: both show 2.5×1044, an initially similar covering factor, and 2.5×1045; the main difference is a roughly factor-of-two larger maximum ejecta velocity in SN 2011dh, plausibly attributable to its smaller envelope mass. It also notes the structural analogy between the flat spectra of compact extragalactic radio sources (the "cosmic conspiracy" of overlapping subcomponents with a common 2.5×1046, in the Blandford–Königl picture) and radio supernovae, while flagging that the poorly understood physics differs: in jets, 2.5×1047 follows from a frozen-in field, whereas in supernovae it is the main unknown; conversely, constant 2.5×1048 is well motivated at a supernova forward shock but requires fine-tuned re-acceleration in jets.
A further implication concerns the proposed class of compact type IIb supernovae: since the progenitor was a yellow supergiant with an extended envelope, and 2.5×1049 km s−10 approaches type Ib/Ic velocities, the paper questions whether such a compact class exists. It also finds no qualitative radio distinction between type IIb and other supernova types, suggesting the inhomogeneity effects identified here may be generic — with the caveat that a good spectral fit with a homogeneous model does not guarantee the fitted radius equals the shock radius, as SN 1993J demonstrates.
Limitations and open questions
Several assumptions bear directly on the results. The covering factor is inferred, not measured: the subcomponent assumptions (equipartition, −11, −12) that fix −13–0.5 are adopted, not derived. The identification of the declining −14 with a declining −15 rests on the persistence of −16 and on an assumed equality of −17 between SN 2011dh and SN 1993J; the alternative with constant −18 is rejected on plausibility grounds only. The transition parameterization −19 fails at late times, and no spectra exist for the period when 3.0×1040 declined. The extent of X-ray-induced ionization near the reverse shock — which sets how much the optical velocity lower limit underestimates 3.0×1041 — is acknowledged as hard to estimate. Finally, the physical mechanism producing 3.0×1042 in the supernova context remains unidentified, and whether the piston-phase-to-CN transition at ~50 days and the ~50% covering factor generalize to radio supernovae lacking early multiwavelength coverage is left open.
Conclusion
The paper shows that a consistent description of SN 2011dh's early radio and optical observations requires abandoning two pillars of the standard radio-supernova model: the immediate applicability of the self-similar CN-phase, replaced by a constant-velocity piston phase transitioning at ~50 days, and source homogeneity, replaced by a covering factor of ~50%. These modifications reconcile the shallow ejecta gradient (3.0×1043) seen in the optical with the radio data, establish 3.0×1044 in preference to 3.0×1045, explain the declining peak flux and the concave light-curve break without invoking a new component, and demonstrate that deduced synchrotron parameters — particularly 3.0×1046, amplified by ~500 for 3.0×1047 — are acutely sensitive to inhomogeneities. The close correspondence with SN 1993J suggests these effects are systematic rather than idiosyncratic, and should be incorporated in modeling of radio supernovae generally.