Upper left panel: ζ_CR as a function of the radius, r. The two dotted black lines at 270 and 500 au identify the radius ranges used to compute the slopes (s₁ and s₂, respectively). The gray shaded region shows the expected X-ray ionization in the outflow cavity for a typical and a flaring T Tauri star (Rab et al. 2017). Lower left panel: hydrogen column density, , from the source center outwards as a function of radius. Upper right panel: ζ_CR as a function of the H₂ column density. The upper x-axis shows the minimum energy that CR protons must have in order not to be thermalized (Padovani et al. 2018). The orange, magenta, and green filled regions in the three upper panels represent the envelope of the ζ_CR and profiles (see Sect. 5.2), while the red, purple, and dark green solid lines show their median value. Lower right panel: ζ_CR map of the blueshifted component (zoom of the upper left panel of Fig. 5) in logarithmic scale (colored map). The purple hatched region shows a circle of radius 270 au, where CRs propagate according to the diffusive regime. For illustration purposes, the black arrow at the position angle # shows a direction used to extract the ζ_CR profile from the map. The black dashed lines identify the radius as well as the corresponding values of and ζ_CR, where the slope of ζ_CR changes from s₁ to s₂.