Mid-plane CR ionisation rate per H₂, ζ, plotted against radius ϖ in the TT, LMP, FU Ori, and HMP models from Shu et al. (2007) together with the standard (unmagnetised) MMSN model from Hayashi (1981). The upper and lower borders of the shaded areas correspond to the unmodulated Galactic CR spectra and , respectively. The kinks seen at the level of ζ ≈ 2 × 10⁻¹⁸ s⁻¹ occur at the transition surface density Σ_tr = 130 g cm⁻² in Eq. (46). The horizontal dotted line at 1.4 × 10⁻²² s⁻¹ shows the value of the ionisation rate set by LLR.