Diffusion of PeV Cosmic Rays in the Turbulent and Multiphase Interstellar Medium

Galactic cosmic rays (CRs) are a fundamental non-thermal component of the interstellar medium (ISM). Understanding the transport of super-high-energy particles is essential for interpreting observations of Galactic PeVatrons. Classical diffusion models assuming a homogeneous and isothermal medium oversimplify the multiphase ISM. We utilize high-resolution 3D MHD simulations to self-consistently generate a multiphase ISM, comprising the warm (WNM), unstable (UNM), and cold neutral medium (CNM), and investigate 1.5-15 PeV particle transport using a test-particle approach. We find that thermal phase transitions induce steep magnetic field strength gradients at phase boundaries, creating localized magnetic fluctuations that act as efficient sites for adiabatic mirror reflections and non-adiabatic pitch-angle scattering, strongly enhancing cross-field transport at these interfaces. However, because phase boundaries occupy only a small volume fraction and particles spend most of their trajectory in the weakly scattering WNM and UNM, the global pitch-angle scattering coefficient in the multiphase ISM is smaller than in an equivalent isothermal medium. This locally strong scattering nevertheless drives both parallel and perpendicular spatial diffusion coefficients to $\sim10^{30} {\rm cm^2 s^{-1}$ at 1.5~PeV, with the perpendicular component exceeding its isothermal counterpart ($\sim 10^{28}{\rm cm^2 s^{-1}$) by two orders of magnitude. Using a phase--phase diffusion matrix decomposition, we show that global CR transport is governed by the volume-filling, trans-Alfvénic WNM and UNM, where particles stream along stochastically wandering field lines. Cross-phase displacement correlations are universally positive, indicating cooperative transport between thermal phases. In contrast, the super-Alfvénic CNM acts as an efficient confinement that substantially suppresses local diffusion.