Pith. sign in

REVIEW 4 major objections 6 minor 15 references

SN 2014C: a metamorphic supernova exploded in the intricate and hydrogen-rich surroundings

T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read With daily-cadence observations from the first month, this paper shows that SN 2014C's light curve departs from pure radioactive decay about 20 days after explosion, requiring an early interaction with hydrogen-rich circumstellar material.

desk verdict New early-time data and a robust H-alpha identification, but the 20-day interaction timing is prior-driven and not uniquely established. read the letter →

arxiv 2411.17008 v3 pith:5SO6X4TQ submitted 2024-11-26 astro-ph.HE astro-ph.SR

classification astro-ph.HEastro-ph.SR
keywords supernovaTypeIbSN2014Ccircumstellarmediuminteractionbolometriclightcurvehigh-velocityH-alpha56Ninucleosynthesisstellarmassloss
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper reports nearly daily photometric and spectroscopic observations of SN 2014C during its first month and uses them to revise what powers this famous 'metamorphic' supernova. It argues that the explosion synthesized 0.14 solar masses of 56Ni, but that this nickel cannot sustain the observed brightness after about ten days past peak; the excess must come from the ejecta crashing into surrounding hydrogen-rich material. If right, the SN-circumstellar interaction starts around 20 days after explosion, much earlier than the roughly 100 days when hydrogen emission lines first appear in the spectra. That matters because it means the progenitor kept a small hydrogen envelope and lost mass violently shortly before exploding, and it makes SN 2014C a testbed for how Type Ib supernovae turn into Type IIn.

What carries the argument

The central tool is the bolometric light curve built from template-subtracted UBVRI photometry and compared with the Arnett radioactive-decay model. The residual excess over that model is interpreted with an ejecta-CSM interaction model, termed the CSIRD model, which adds shock-powered luminosity to radioactive decay and includes a gamma-ray leakage factor. Markov chain Monte Carlo fitting yields the ejecta mass, nickel mass, CSM density, and interaction start time. The velocity evolution of the 6200 Å absorption works as the diagnostic that separates high-velocity H-alpha from photospheric Si II, and the colors and light-curve shape identify the onset of the extra energy source.

What would settle it

A direct test would fit a magnetar spin-down model, with the same explosion date and reddening, to the bolometric light curve; if it reproduces the excess at least as well as the CSM-interaction model, the early-interaction timing claim would not be uniquely supported. Alternatively, a sensitive X-ray or radio observation designed to detect the shock at day 20 would confirm or rule out ongoing interaction.

Watch

Extended reading notes

Core claim

SN 2014C was a Type Ib supernova with an unusually high peak luminosity (about 4.3 × $10^{42}$ erg $s^{-1}$), a fast rise time near 11.6 days, and a slow decline compared with typical Type Ib events. The absorption feature near 6200 Å is high-velocity H-$\alpha$ in the outer ejecta, not Si II, as shown by its velocity evolution over the first 20 days; this means a small amount of hydrogen remained in the progenitor envelope. The bolometric light curve follows a radioactive-decay model with 0.14 solar masses of 56Ni for roughly the first 15 days, then shows an excess starting at about 20 days after explosion. A hybrid model combining radioactive decay and ejecta-CSM interaction fits the full light curve with an ejecta mass of about 2.2 solar masses and an interaction onset at about 19.9 days. The authors conclude that SN-CSM interaction began much earlier than the emergence of IIn-like nebular features, which are first seen near day 100.

Load-bearing premise

The load-bearing premise is that the light-curve excess beginning around 20 days is powered by ejecta hitting circumstellar material, rather than by an alternative energy source such as a magnetar; the authors note their model is only one possible explanation.

Editorial extensions

If this is right

  • The peak-luminosity nickel mass of 0.14 solar masses is an upper limit: if part of the peak brightness comes from interaction, less 56Ni is needed.
  • The inner edge of the dense hydrogen-rich shell lies at roughly 4 × 10^14 cm from the explosion, with denser material near 10^16 cm reached around day 100.
  • The progenitor was not fully stripped: high-velocity H-alpha in the outer ejecta shows a small hydrogen envelope remained at explosion.
  • The implied mass-loss rate of order 0.2 solar masses per year at 1000 km s^-1 points to an eruptive or binary mass-loss history rather than steady Wolf-Rayet or red-supergiant winds.
  • The interaction already supplies extra light by day 20, so classifying SN 2014C as a pure Ib during its first months misses the engine that shapes its luminosity.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • An untested alternative is that a magnetar, rather than CSM interaction, powers the post-peak excess; fitting a magnetar spin-down model to the same bolometric light curve would settle whether the 20-day interaction onset is unique.
  • Because the model assumes a spherical CSM, the inferred CSM mass and interaction time would shift if the shell is clumpy or torus-like, as some later observations suggest; testing with non-spherical geometry is a direct next step.
  • The same analysis applied to other fast-rising Type Ib/c supernovae with dense early photometry could reveal that hidden early interaction is common, turning 'metamorphic' events into an extreme of a continuous population.
  • A testable prediction is that sensitive X-ray and radio observations beginning near day 20 should catch the forward shock while it is still young; current non-detections only cover the first 20 days.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper presents new high-cadence photometric and spectroscopic observations of SN 2014C covering the first month after explosion, with template-subtracted light curves and a dense spectral series. The authors classify it as a luminous, fast-rising Type Ib SN, argue that the ~6200 Å absorption is high-velocity Hα based on its velocity evolution, and derive M_Ni = 0.14 Msun from the peak bolometric luminosity using the Arnett relation. They show that this radioactive-decay model cannot reproduce the observed luminosity after τ ≈ 20 days, and they fit a combined radioactive-decay plus CSM-interaction (CSIRD) model that yields an interaction start time ti ≈ 20 days. They conclude that the SN-CSM interaction in SN 2014C likely begins much earlier than the emergence of IIn-like nebular Hα features at ~100 days, and they discuss the progenitor's complex H-rich surroundings.

Significance. The central photometric result—the divergence of the bolometric light curve from an Arnett model with M_Ni = 0.14 Msun after τ ≈ 20 days—is a robust, reddening-independent comparison because both the observed peak and the model scale with the adopted extinction and distance. The early-time spectral series is valuable, and the identification of the 6200 Å feature as high-velocity Hα using the evolution of its absorption velocity is convincing. The paper also provides a useful comparison of SN 2014C with other metamorphic SNe. If the early interaction timing is correct, it would significantly revise the picture of the CSM geometry around this famous object. However, the interpretation of the excess as CSM interaction is not uniquely established, and the paper itself concedes in Section 6 that the CSIRD model is 'only one possible explanation.' The manuscript is honest about its limitations, but the load-bearing timing claim needs additional support from model comparison and a non-circular treatment of the interaction-start prior.

major comments (4)
  1. [Section 5, Eq. (2) and preceding paragraph] The prior on the interaction start time is set as ti = 20 ± 5 days because 'the detection of interaction signals starting 20 days post-explosion indicates that the SN ejecta caught up with the H-rich CSM at that moment.' Since the same light-curve excess is what motivates the prior, the recovered ti = 19.9 d is not independent evidence for the timing. Refit the CSIRD model with a wide, uninformative prior on ti (e.g., uniform from 5 to 50 days) and report whether the posterior remains peaked near 20 days; also discuss the degeneracy between ti and the CSM density normalization.
  2. [Section 5, Fig. 14 and Section 6] The excess over the radioactive-decay model is attributed solely to ejecta-CSM interaction, but a magnetar spin-down model is not tested, even though the paper cites SN 2005bf as a possible magnetar-powered case. A magnetar with P ≈ 10–30 ms and B ≈ 10^14 G can produce a slowly declining excess after the radioactive peak, and the flat-to-blue colors at τ ≳ 20 days do not uniquely select CSM interaction. The manuscript itself concedes in Section 6 that the CSIRD model is 'only one possible explanation' and that interaction 'could have commenced earlier and intensified rapidly.' To support the abstract's timing claim, either fit a magnetar (or another alternative) model to the same bolometric light curve and show it is disfavored, or explicitly temper the claim to state that the excess requires an additional energy source whose origin is not uniquely determined.
  3. [Section 3.1 and Section 5] The explosion epoch is derived from fireball fitting as MJD 56658.91, but the paper notes that Margutti et al. (2017) estimated MJD 56656, roughly 3 days earlier. All phases τ are measured from the adopted explosion date, so the onset of the excess (τ ≈ 20 d) and the inferred CSM radius shift by up to ~3 days. The paper should propagate this uncertainty into the CSIRD fit, for example by treating the explosion epoch as a nuisance parameter with a prior spanning both estimates, and should state how ti and the resulting interaction radius change.
  4. [Section 5, Eq. (1) and Section 6] The 56Ni mass is first derived from the peak luminosity using Eq. (1) under the assumption that the peak is entirely radioactive, and then it is floated in the MCMC fit, which returns 0.14 ± 0.01 Msun. This posterior is essentially the input value and does not test the assumption. In Section 6 the authors state that the nickel mass derived from the peak 'may represent an upper limit,' but the CSIRD fit does not allow for a substantially lower M_Ni to compensate with CSM power at peak. Run the fit with a broader prior on M_Ni (e.g., 0.05–0.20 Msun) and report the marginalized posterior and its correlation with the CSM parameters.
minor comments (6)
  1. [Section 5, Eq. (1)] The formula for the nickel mass is typeset ambiguously; please write it explicitly as M_Ni = (L_max / 10^43 erg s^-1) / [6.45 exp(-t_r/τ_Ni) + 1.45 exp(-t_r/τ_Co)] in solar masses.
  2. [Section 5, Eq. (2)] The gamma-ray leakage factor A is introduced without derivation or reference. Please clarify what is new relative to the CSMRD model of Chatzopoulos et al. (2012, 2013) and state the units and convention for the deposition term (1 - exp(-A t^{-2})).
  3. [Section 5] The CSM mass is fixed to 1 Msun based on a range of literature values, but the text later notes that Orlando et al. (2024) infer 2.5 Msun. Please quantify how the inferred ti and CSM density normalization vary with the assumed M_CSM.
  4. [Figure 12 caption] The spectrum at t ≈ 1016 d is from Shivvers et al. (2019), not from this work; the caption should attribute it explicitly to avoid implying that it is a new observation from this paper.
  5. [Figure 14] The 'BB_cor' curve is not defined in the text or caption; please specify whether it is the integral of the blackbody fit and how it differs from the direct UBVRI integrated flux.
  6. [Table 1 and Table A5] The V-band rise time is listed as 11.47 d in Table 1 and the R-band rise time as 12.13 d in Table A5, while the text quotes t_rise ≈ 11.6 d for the bolometric rise; please reconcile these values or state which quantity each refers to.

Circularity Check

1 steps flagged · score 6.0 of 10

Early-interaction start time is prior-fed: ti=19.9 d is recovered from a prior set to 20±5 d using the same τ≳20 d light-curve excess the model is meant to explain.

  1. fitted input called prediction [Section 5 (Bolometric Light Curve and Modeling), paragraph after Eq. (2)]
    "The detection of interaction signals starting 20 days post-explosion indicates that the SN ejecta caught up with the H-rich CSM at that moment. Hence, we set the prior on the initial time of interaction as ti = 20 ± 5 days."

    The 'interaction signals starting 20 days post-explosion' are the same bolometric excess at τ ≳ 20 d relative to the radioactive-decay model that the CSIRD model is invoked to explain. Setting the MCMC prior for ti to 20 ± 5 d from this same excess ensures the recovered value will cluster near 20 d; the later quoted result ti = 19.9+3.34−3.25 d is therefore a restatement of the input prior, not an independent measurement. The abstract's new claim that SN–CSM interaction began about 20 d after explosion, earlier than the nebular Hα emergence, is thus supported by construction rather than by independent evidence.

full rationale

The central derivation chain has one genuine circular step. Section 5 identifies a bolometric excess at τ ≳ 20 d, declares it an 'interaction signal,' sets the CSIRD model prior to ti = 20 ± 5 d on that basis, and then reports a fitted ti = 19.9 d as support for early interaction. Because the prior is derived from the same excess the model is used to interpret, the recovered start time is not independent evidence for the abstract's headline claim that SN–CSM interaction began about 20 d after explosion. The paper is candid: in the Summary it concedes the combined model is 'only one possible explanation' and that interaction 'could have commenced earlier and intensified rapidly,' and it flags the 56Ni mass as an upper limit. That candor mitigates the circularity but does not remove it, because the ti result is still presented as a fitted constraint. Other parts of the paper are self-contained: the high-cadence photometry, template-subtracted light curves, the Hα versus Si ii discrimination via velocity evolution, and the nebular spectroscopy are independent observations. The magnetar alternative (SN 2005bf) is not modeled, but that is a model-degeneracy and correctness concern, not a circularity; it reinforces that the early-timing claim is not uniquely determined. No load-bearing self-citation or imported uniqueness theorem is present. Overall, there is partial circularity in the timing claim, while the bulk of the observational analysis stands independently.

Assumptions & free parameters 7 free parameters · 6 assumptions · 0 invented entities

The central claims rest on standard supernova tools (Arnett, CSMRD, SYN++) and a handful of hand-chosen model parameters (n=7, x0=0.3, kappa_gamma, CSM mass = 1 Msun). No new physical entities are proposed. The load-bearing inputs are the reddening, the explosion epoch, and the model priors, especially t_i = 20 ± 5 d.

free parameters (7)
  • Host-galaxy reddening E(B-V)host = 0.62 ± 0.10 mag; adopted total E(B-V) = 0.70 ± 0.10 mag
    Adopted from Na I D EW and V-R color methods; directly scales all luminosities and absolute magnitudes.
  • Explosion epoch = MJD 56658.91
    Derived from fireball fitting to early photometry; the paper states this method has limitations for SNe Ib/c and notes Margutti et al. give MJD 56656, a ~3 day difference.
  • 56Ni mass = 0.14 ± 0.03 Msun
    Derived from peak bolometric luminosity via the Arnett relation; also returned as a fitted parameter in the CSIRD MCMC fit.
  • Ejecta mass Mej = 2.22 (+0.69, -0.77) Msun
    MCMC-fit parameter in the CSIRD model; affects diffusion timescales.
  • CSM density normalization rho_CSM,0 = 1.81 (+1.57, -0.78) x 10^-14 g cm^-3
    MCMC-fit parameter setting the CSM density at the inner radius.
  • Interaction start time t_i = 19.9 (+3.34, -3.25) d
    MCMC-fit parameter; prior set to 20 ± 5 d based on the same light-curve excess, so the posterior is not fully independent.
  • CSM interaction radiative efficiency = 0.2 (+0.3, -0.1)
    MCMC-fit parameter describing the efficiency of converting shock kinetic energy into radiation.
assumptions (6)
  • domain assumption The Arnett (1982) radioactive decay model is an appropriate baseline for the early bolometric light curve.
    Used to derive M_Ni from peak luminosity and to identify the excess after about 20 d; assumes spherical ejecta, constant opacity, and full gamma-ray trapping during the photospheric phase.
  • domain assumption The CSM follows a steady stellar-wind density profile (rho ~ r^-2).
    Standard assumption in the CSMRD model; the paper acknowledges a torus-shaped CSM is possible and that this simplification may bias the CSM mass.
  • ad hoc to paper The outer ejecta density is a power law with n=7 and transition radius x0=0.3.
    These values are chosen by hand, following common practice, rather than derived from observations or independent modeling.
  • standard math The gamma-ray opacity of the ejecta is kappa_gamma = 0.027 cm2/g.
    Taken from prior literature (Cappellaro et al. 1997 etc.) and used in the gamma-ray leakage factor; standard for Ni/Co decay gamma-rays.
  • domain assumption The bolometric luminosity can be recovered by black-body fitting to UBVRI photometry.
    Assumes the optical SED is representative and that UV/IR contributions are small; systematic errors in the bolometric correction are not propagated.
  • domain assumption The SYN++ spectral synthesis code provides a reliable description of line formation for this SN.
    Used to argue that H-alpha fits better than Si II at +6 d; the fits are visual, not quantitative.

how reviews work

0 comments
Cite this review

Pith. "Pith review of SN 2014C: a metamorphic supernova exploded in the intricate and hydrogen-rich surroundings." pith.science (2026). https://pith.science/paper/5SO6X4TQ

@misc{pith2026241117008,
  author       = {Pith},
  title        = {Pith review of: SN 2014C: a metamorphic supernova exploded in the intricate and hydrogen-rich surroundings},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5SO6X4TQ}},
  note         = {Machine review of arXiv:2411.17008}
}
abstract

We present photometric and spectroscopic observations of supernova (SN) 2014C, primarily emphasizing the initial month after the explosion at approximately daily intervals. During this time, it was classified as a Type Ib SN exhibiting a notably higher peak luminosity ($L_{\rm peak}\approx4.3\times10^{42}\rm erg\,s^{-1}$), a faster rise to brightness ($t_{\rm rise} \approx 11.6$ d), and a more gradual dimming ($\Delta m_{15}^{V} \approx 0.48$ mag) compared to typical SNe Ib. Analysis of the velocity evolution over the first $\sim$ 20 days after the explosion supports the view that the absorption near 6200\AA is due to high-velocity H$\alpha$ in the outer layers of the ejecta, indicating the presence of a small amount of hydrogen in the envelope of progenitor before the explosion. Assuming the peak luminosity is entirely attributed to radioactive decay, we estimate that 0.14 ${\rm M}_{\odot}$ of $^{56}$Ni was synthesized in the explosion. However, this amount of nickel could no longer maintain observed brightness approximately ten days after peak luminosity, suggesting additional energy sources beyond radioactive decay. This supplementary energy likely originates from interaction with the circumstellar medium (CSM). Consequently, the timing of the SN-CSM interaction in SN 2014C may occur much earlier than the emergence of IIn-like features during the nebular phase.

Figures

Figures reproduced from arXiv: 2411.17008 by the authors.

Figure 1
Figure 1. Finder chart of SN 2014C and its local reference stars, combined by BV R-bands images taken by the LJT and YFOSC in January 2014. The mean FWHM of the combined image is ∼ 1 ′′ .50 under the scale of ∼ 0 ′′ .28/pixel. The zoomed-in templates taken in November 2020 for the BV R-bands, which encompass the region containing SN 2014C, are presented on the right. about 100 days after the explosion (Margutti et al. 2017; M… view at source ↗
Figure 2
Figure 2. [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Spectra of SN 2014C obtained by LJT (black), Asiago 1.82 m telescope (blue), and NOT (orange-red). These spectra are shifted vertically for clarity and labeled according to the phase referring to the V -band maximum. A bin size = 6 pixels was applied to improve the signal-to-noise ratio of each spectrum. The region around the center wavelength of Na i D, lines at t ≈ 1 d is zoomed in. The red dotted line marks the H… view at source ↗
Figures from the paper (11 more)
Figure 4
Figure 4. Figure 4: UBV RI-band absolute light curves of SN 2014C, compared with the light curves of some well-studied SESNe. The two filled squares mark photometry from unfiltered KAIT images obtained on Jan. 5.09 (MJD 56659.10) and from earlier red magnitudes on a pre-discovery image ta…
Figure 5
Figure 5. Figure 5: Color curves of SN2014C comparing with other SNe Ib/c. All of these samples were corrected with reddening. brightness in t R 1/2 ≈ 7.3 d, half the typical rise time for SNe Ib. The rise time of SN depends on factors such as the ejecta mass and explosion energy, with a …
Figure 6
Figure 6. Figure 6: Spectra of SN 2014C at -5 d and +6 d reproduced by SYN++, including the element of H (panel a and c) or Si (panel b and d), respectively. -5 0 5 10 15 Days after the V band Maximum -18000 -16000 -14000 -12000 -10000 -8000 -6000 -4000 -2000 A b s o rp t i o n V e l o c …
Figure 7
Figure 7. Figure 7: Evolution of line velocities for Hα, Si ii λ 6355, He i λ5876, He i λ6678, and Ca iiIR triplet of SN 2014C, derived from the absorption minima of each spectral line. Note: velocities of Hα and Si ii λ 6355 are derived from the same absorption line around 6200˚A. better…
Figure 9
Figure 9. Figure 9: Velocity derived from the absorption minimum of He i λ 5876˚A, including SN 1998dt (Matheson et al. 2001); SN 2005bf (Tominaga et al. 2005); SN 2004gq, SN 2005hg and SN 2008D ( [PITH_FULL_IMAGE:figures/full_fig_p010_9.png]
Figure 8
Figure 8. Figure 8: Comparison of spectra of SN 2014C with similar epoch spectra of other well-studied SNe. Labels in paren￾theses to the right of each spectrum indicate days from the V band maximum time. References for each spectrum are listed in Table A6. These spectra have been correct…
Figure 11
Figure 11. Figure 11: The emission lines of [O i] 6300, 6364 dou￾ble peak of SNe 2014C, 2008D (Modjaz et al. 2009), 2008ax (Chornock et al. 2011; Taubenberger et al. 2011) in velocity space. The zero point is 6300˚A and is indicated by the gray dashed line. These spectra were corrected for…
Figure 12
Figure 12. Figure 12: Spectra of SN 2014C at t ∼ 303 d (this work) and t ∼ 1016 d (Shivvers et al. 2019) comparing to that of SN 2014C-like SNe 2019oys (Sollerman et al. 2020), 2019yvr (Ferrari et al. 2024), 2001em (Shivvers et al. 2019), 2004dk (Mauerhan et al. 2018) and SN II 1987A (from…
Figure 13
Figure 13. Figure 13: The bolometric light curve of SN 2014C compared with some well-studied SESNe, including: SNe IIb SN 2008ax (Pastorello et al. 2008), SN 2011fu (Kumar et al. 2013); SN IIb/Ib SN 2007Y (Stritzinger et al. 2009); SNe Ib SN 1999ex (Stritzinger et al. 2002), SN 2008D (Mazz…
Figure 14
Figure 14. Figure 14: Bolometric light curve of SN 2014C derived from the black body fit based on the UBV RI-band photometry. The integrated flux of UBV RI-band is also plotted. The dotted, dashed, and solid lines represent the radioactive decay (RD), CSM, and RD+CSM models, respectively. …
Figure 15
Figure 15. Figure 15: Joint confidence level contours of the CSIRD parameters inferred from the MCMC-based fitting. supports the hypothesis of delayed interaction between SN and surrounding CSM. Previous observational evi￾dence suggests that the interaction of SN 2014C began 45 d post-expl…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

15 extracted references · 14 canonical work pages

  1. [1]

    E., Horesh, A., Mooley, K

    Anderson, G. E., Horesh, A., Mooley, K. P., et al. 2017, MNRAS, 466, 3648, doi: 10.1093/mnras/stw3310 Arnett, W. D. 1982, ApJ, 253, 785, doi: 10.1086/159681 Arnett, W. D., Bahcall, J. N., Kirshner, R. P., & Woosley, S. E. 1989, ARA&A, 27, 629, doi: 10.1146/annurev.aa.27.090189.003213 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, ...

  2. [2]

    Bersten et al. (2014)

  3. [3]

    Fremling et al. (2014)

  4. [4]

    Folatelli et al. (2016)

  5. [6]

    Bianco et al. (2014)

  6. [7]

    Folatelli et al. (2006)

  7. [9]

    Stritzinger et al. (2018)

  8. [10]

    Stritzinger et al. (2009)

Show all 15 references
  1. [11]

    Mazzali et al. (2008)

  2. [12]

    Modjaz et al. (2009)

  3. [13]

    Valenti et al. (2011)

  4. [17]

    Modjaz et al. (2014)

  5. [18]

    Pastorello et al. (2008)

  6. [19]

    Tsvetkov et al. (2009)

  7. [20]

    (2024) T able A6.Reference of the SNe in Figure

    Ferrari et al. (2024) T able A6.Reference of the SNe in Figure. 8 Star (a) pre-maximum (b) maximum (c) 7 d post maximum (d) 14 days post maximum iPTF13bvn Cao et al. (2013) Cao et al. (2013) ... Srivastav et al. (2014) SN 2004gq Modjaz et al. (2014) Modjaz et al. (2014) ... Mo...

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.