---
title: Cosmic-Ray Energy-Dependent Injection (CREDIT)
url: https://www.emergentmind.com/topics/cosmic-ray-energy-dependent-injection-time-credit
type: topic
---

# Cosmic-Ray Energy-Dependent Injection (CREDIT)

Cosmic-Ray Energy-Dependent Injection Time (CREDIT) denotes a class of non-stationary source models in which the time relevant to cosmic-ray release or entry into acceleration depends on particle energy or rigidity. In supernova-remnant (SNR) escape models, CREDIT means that higher-rigidity particles are released earlier and lower-rigidity particles later; in a non-stationary diffusive-shock-acceleration treatment, it means that the highest-energy particles observed at a later epoch must have entered acceleration at earlier times after the supernova explosion. The term should therefore be distinguished from work on the source spectral exponent alone, where “injection” refers to the power-law index of the emitted spectrum rather than to an energy-dependent time history [2409.11012, 1707.00136, 1703.02802].

## 1. Terminological scope and principal usages

The cited literature uses the term “injection” in distinct but related ways. In all cases, the central issue is temporal non-stationarity, but the time variable attached to particle production or release is not identical across models.

| Reference | Energy/time content | Relation to CREDIT |
|---|---|---|
| [1707.00136] | Higher-energy particles observed today must have been injected earlier into acceleration; gamma-ray energy maps to post-explosion epoch | Explicit CREDIT concept |
| [2409.11012], [2507.07088] | Higher-rigidity cosmic rays escape an SNR earlier; lower-rigidity particles escape later | Explicit CREDIT scenario |
| [2606.12390] | Continuous injection with energy-dependent effective escape from an evolving wind bubble | CREDIT-relevant, but not a full explicit CREDIT law |
| [1703.02802] | Reconstruction of the source spectral exponent $p$ from Galactic propagation | Not about energy-dependent injection time |

This terminological distinction is essential. In [2409.11012] and [2507.07088], CREDIT is a rigidity-dependent escape history from individual SNRs. In [1707.00136], CREDIT links gamma-ray energy to the epoch when parent particles entered the acceleration process. By contrast, [1703.02802] assumes an instantaneous point source with a power-law injection spectrum and studies how propagation through an inhomogeneous interstellar medium modifies the observed spectrum, not how injection time depends on energy.

## 2. CREDIT as a non-stationary acceleration history in supernova remnants

A formulation of CREDIT appears in the analysis of gamma-ray spectra from SNRs, where the highest-energy hadrons require the longest acceleration time and therefore preserve a memory of the earliest post-explosion injection history. In this picture, low-energy particles and low-energy gamma rays are mainly shaped by later injection, whereas TeV gamma rays originate from particles that began acceleration during the first months after the supernova explosion. The particle transport is modeled with the non-stationary acceleration equation
$$
\frac{\partial f}{\partial t}+u\frac{\partial f}{\partial x}
=
\frac{\partial}{\partial x}\left(D\frac{\partial f}{\partial x}\right)
+\frac{1}{3}\frac{du}{dx}\,p\frac{\partial f}{\partial p}+Q,
$$
and the injection term
$$
Q(t,x,p)=\frac{\eta n_1 u_1}{4\pi p_i^2}\,\delta(p-p_i)\,\delta(x)\,Q_t(t),
$$
where $Q_t(t)$ carries the time dependence of the injection efficiency [1707.00136].

In the test-particle regime and $t_1 \gg t_2$, the shock solution is written as
$$
f_o(t,p)=f_o(p)\int_0^\tau Q_t(\tau-\tau')\,\varphi_o(\tau')\,d\tau',
$$
with stationary spectrum
$$
f_o(p)=\frac{\eta n_1}{4\pi p_i^3}\frac{3\sigma}{\sigma-1}\left(\frac{p}{p_i}\right)^{-s_f},
\qquad
s_f=\frac{3\sigma}{\sigma-1}.
$$
The steady-state criterion is
$$
\mathcal{I}(\tau)\equiv \int_0^\tau \varphi_o(\tau')\,d\tau' \approx 1.
$$
For a power-law injection history, $Q_t\propto \tau^\beta$, the spectrum becomes approximately
$$
s=s_f+\alpha\beta,
$$
with $D(p)\propto p^\alpha$. Time-dependent injection therefore changes not only the normalization but also the slope of the accelerated particle spectrum.

The application to IC443 uses SN1987A as a proxy for the parent supernova. From ATCA radio data between 1517 and 8014 days after explosion, the inferred injection history rises smoothly with time and is approximately described during much of the observational interval by
$$
Q_t \sim \log(t/950\,\mathrm{days}),
$$
reaching the steady-state level $Q_t=1$ only around $10^4$ days. For the unobserved early phase, the adopted model assumes steady injection up to day 100, then $Q_t\propto t^{\beta_x}$ with $\beta_x=0.8$. With the same model, the fit to the IC443 gamma-ray spectrum reproduces the observed data from Fermi, MAGIC, and VERITAS without imposing an ad hoc proton break: the proton spectrum breaks around $p/m_pc\sim 50$, the gamma-ray spectrum steepens around $\sim 2$ GeV, the spectrum above $\sim 0.3$ TeV is controlled by the earliest injection phase, photons above $\sim 20$ GeV are influenced by injection before the radio observations of SN1987A began, and the gamma-ray spectrum above $1$ GeV reflects the first $\sim 20$ years of SNR evolution [1707.00136].

## 3. CREDIT as rigidity-dependent escape from supernova remnants

A second, and now more explicit, formulation of CREDIT treats the source term itself as rigidity dependent in time. In this SNR escape picture, cosmic rays do not all leave the remnant simultaneously; instead, the escape time depends on rigidity $\mathcal R$. The source term is
$$
Q(\mathcal R,t)\equiv \delta\!\left(t-t_{\rm esc}(\mathcal R)\right)\,Q_{\mathcal R}(\mathcal R),
$$
where $t_{\rm esc}(\mathcal R)$ is the rigidity-dependent escape time. The corresponding smooth-source null hypothesis averages over many remnants and yields
$$
\bar Q(\mathcal R)\equiv \nu \int_0^\infty dt\,\delta(t-t_{\rm esc}(\mathcal R))Q_{\mathcal R}(\mathcal R)=\nu Q_{\mathcal R}(\mathcal R).
$$
In that limit, the escape history is invisible because source discreteness is washed out [2409.11012, 2507.07088].

The physical picture is tied to the evolution of the maximum confining rigidity. Early in the remnant’s life, magnetic-field amplification is strongest and higher maximum rigidities can be confined. After the Sedov-Taylor transition, the shock slows and $\mathcal R_{\max}(t)$ drops. The escape delay is parameterized as
$$
t_{\rm esc}(\mathcal R)=t_0+\Delta t_{\rm esc}(\mathcal R),
$$
with
$$
\Delta t_{\rm esc}(\mathcal R)=t_{\rm Sed}\left(\frac{\mathcal R}{\mathcal R_{\rm knee}}\right)^a,
\qquad
a=\frac{\ln(t_{\rm life}/t_{\rm Sed})}{\ln(\mathcal R_b/\mathcal R_{\rm knee})},
$$
for $\mathcal R>\mathcal R_b$, and
$$
\Delta t_{\rm esc}(\mathcal R)=t_{\rm life}
$$
otherwise. The fiducial values are $t_{\rm Sed}=1\,\mathrm{kyr}$, $t_{\rm life}=100\,\mathrm{kyr}$, $\mathcal R_{\rm knee}=3\,\mathrm{PV}$, $\mathcal R_b\sim 10\,\mathrm{TV}$ as the physically motivated confinement scale, and source rate $\nu=0.03\,\mathrm{yr^{-1}}$. The injection spectrum is
$$
Q_{\mathcal R}(\mathcal R)\propto \mathcal R^{-2.2}.
$$

This construction interpolates between CREDIT and burst-like injection. For $\mathcal R>\mathcal R_b$, escape is energy dependent; in the limit $\mathcal R_b\to\infty$, all rigidities are released at the same time. The CREDIT prediction is that young and nearby SNRs produce rigidity-localized enhancements in the local proton spectrum, rather than a broad smooth tilt, because a source contributes sharply when its age and $t_{\rm esc}(\mathcal R)$ align [2409.11012, 2507.07088].

## 4. Galactic transport, Green’s functions, and stochastic source populations

The rigidity-dependent SNR CREDIT scenario is embedded in a simplified Galactic diffusion model. The transport equation for protons is
$$
\frac{\partial \psi_{\mathcal R}}{\partial t} -\kappa(\mathcal R)\nabla^2\psi_{\mathcal R} =Q(\mathcal R,t,\mathbf x),
$$
with diffusion coefficient
$$
\kappa(\mathcal R)=\kappa_0\,\beta(\mathcal R)\left(\frac{\mathcal R}{\mathcal R_0}\right)^\delta.
$$
The setup assumes that all sources lie in the Galactic disk $z=0$, diffusion occurs in a halo with free-escape boundaries at $z=\pm H$, no radial boundary is imposed, and only protons above a few GV or above about 10 GeV are considered, so inelastic losses, advection, reacceleration, and convection are neglected in the relevant implementations. Fiducial parameters are $\delta=0.6$, $\kappa_0=6.8\times10^{27}\,\mathrm{cm^2\,s^{-1}}$ at $\mathcal R_0=1\,\mathrm{GV}$, and halo height $H=4\,\mathrm{kpc}$ [2409.11012, 2507.07088].

For a single source at $(t_i,\mathbf x_i)$, the source term is
$$
Q_i(\mathcal R,t,\mathbf x)=Q_{\mathcal R}(\mathcal R)\,\delta(\mathbf x-\mathbf x_i)\,\delta\!\left(t-t_{{\rm esc},i}(\mathcal R)\right),
$$
and the Green’s-function solution is written with
$$
\sigma^2(\mathcal R,t;t_i)=2\,\kappa(\mathcal R)\left(t-t_{{\rm esc},i}(\mathcal R)\right).
$$
The factor $\vartheta$ enforces the free-escape boundary conditions by a mirror-charge construction, and the local flux is
$$
\Phi_{\mathcal R}=\frac{v}{4\pi}\psi_{\mathcal R}.
$$
A causality restriction is imposed by excluding source injections outside the observer’s past light cone.

The source population is stochastic. Ages are drawn from a uniform distribution consistent with the Galactic supernova rate, source positions are drawn from an axisymmetric disk distribution depending on Galactocentric radius, the Sun is placed at $R_\odot=8.3\,\mathrm{kpc}$, the Galaxy radius is $R=15\,\mathrm{kpc}$, and sources are treated as one SNR population injecting protons up to the knee. The break rigidity $\mathcal R_b$ is either fixed within realizations spanning $[1,25]\,\mathrm{TV}$ or drawn per source from
$$
p(\mathcal R_b)\propto \frac{1}{\mathcal R_b},
\qquad
\mathcal R_b\in[1,25]\,\mathrm{TV}.
$$
The spectra are generated by summing Green’s-function contributions from tens of millions of simulated sources, and one implementation uses GPU acceleration via **jax** [2409.11012, 2507.07088].

## 5. Predicted spectral structure and statistical discrimination

The principal observational prediction of the SNR CREDIT scenario is the appearance of narrow rigidity-dependent peaks in the local proton spectrum. Below $\sim 1\,\mathrm{TV}$, realizations are fairly smooth, with only percent-level deviations from the ensemble mean. Above $\sim 1\,\mathrm{TV}$, the CREDIT scenario generates pronounced, narrow peaks and rapid structure. For young and nearby sources, the flux can be enhanced by tens of percent in narrow rigidity intervals, and in some realizations even by a few times the average flux in narrow bins. These structures remain visible after rebinning to AMS-02 or DAMPE rigidity bins and exceed expected experimental uncertainties in many cases [2409.11012, 2507.07088].

The uncertainty model used for detectability estimates approximates statistical errors as
$$
\sigma_{\rm stat}[\%] = 0.42\,\big(\mathcal R[\mathrm{TV}]\big)^{0.6},
$$
with a $1\%$ uncorrelated systematic uncertainty added in quadrature. The key contrast is with the smooth-source null hypothesis, which assumes a continuous spatial source density, steady injection, and no observable imprint of individual source ages or distances. In that limit, the local spectrum is effectively the ensemble mean.

To test distinguishability, a decision tree classifier is trained on three hypotheses: smooth source distribution with uncorrelated experimental errors, burst-like discrete source injection, and CREDIT with $\mathcal R_b\in[1,25]\,\mathrm{TV}$. The classifier takes flux values in rigidity bins above $10\,\mathrm{GV}$ as a 58-dimensional input space. It is trained and validated with $4\times 10^6$ realizations per scenario and 10-fold cross-validation. The confusion matrices are almost diagonal, misclassifications are rare, and varying $\mathcal R_b$ either fixed per realization or randomly per source barely changes the performance. Halving the source rate $\nu$, halving the source lifetime $t_{\rm life}$, or smoothing the release time with a Gaussian kernel of width $1\,\mathrm{kyr}$ does not substantially alter the classification structure for $\mathcal R_b\gtrsim 1\,\mathrm{TV}$ [2409.11012, 2507.07088].

The interpretive use is explicit. A CREDIT-like classification could be used to infer or constrain the supernova time $t_i$ of the source associated with a feature, given assumptions about $\kappa(\mathcal R)$, $\mathcal R_b$, $t_{\rm Sed}$, and $t_{\rm life}$. A burst-like classification would imply a lower bound on $\mathcal R_b$. A smooth classification would challenge the SNR paradigm in the sense stated in the paper: non-detection of source signatures would be unlikely if SNRs are indeed the dominant Galactic cosmic-ray sources [2409.11012].

## 6. Related transport problems and non-CREDIT usages

A related but distinct time-dependent escape problem arises in wind bubbles. There, particles are continuously injected at the wind termination shock and propagate through advection and diffusion until escape at the time-dependent position of the forward shock, treated as a free escape boundary. The model is one-dimensional and spherically symmetric, with
$$
R_{\rm fs}\propto t^{3/5},
\qquad
R_{\rm ts}\propto t^{2/5},
$$
so the distance between termination shock and forward shock increases with time. The injected spectrum is $dN/dE\propto E^{-2}$ with a time-dependent high-energy cutoff $E_{\max}(t)$. There is no explicit energy-dependent injection-time formula, but the effective escape time becomes energy dependent because high-energy particles diffuse faster and are more likely to reach the moving boundary, while low-energy particles are delayed or trapped. The escaping spectra can be harder than $\sim E^{-2}$; for Kraichnan diffusion the average hardening is about $\Delta\gamma=\gamma_{\rm esc}-\gamma_{\rm injected}\simeq 0.1-0.2$, Bohm diffusion gives the strongest low-energy suppression, and Kolmogorov diffusion the weakest. The trapped low-energy population may contribute to multimessenger radiation and accumulated grammage within the bubble [2606.12390].

A separate source of confusion is the use of “injection” in studies of the source spectral exponent. In the analysis of the average cosmic-ray injection exponent at Galactic sources, the source is assumed to be instantaneous and point-like,
$$
S(r,t,E)=S_0 E^{-p}\delta(r)\delta(t),
$$
with diffusion coefficient
$$
D(E)=D_0E^\delta.
$$
In the homogeneous picture,
$$
N(r,E)\sim E^{-p-\delta},
$$
so if $N(E)\propto E^{-\eta}$ one recovers $p=\eta-\delta$, which for $\eta\approx 2.7$ and $\delta\sim 0.3$ gives $p\sim 2.4$. The paper’s main result is that, once interstellar-medium inhomogeneity is included, both anomalous diffusion and normal diffusion in a non-homogeneous medium yield a steeper average source index, $p\sim 2.8\text{--}3.0$. This is a result about the average energy spectral exponent at Galactic sources, not about a source model in which higher- or lower-energy cosmic rays are injected at different times [1703.02802].

Source: https://www.emergentmind.com/topics/cosmic-ray-energy-dependent-injection-time-credit