Resolving Star Cluster Formation in Galaxy Simulations with Cosmic Ray Feedback

Star clusters host the massive stars responsible for feedback in star-forming galaxies. Stellar feedback shapes the interstellar medium (ISM), affecting the formation of future star clusters. To self-consistently capture the interplay between feedback and star formation, a model must resolve the parsec-scale star formation sites and the multiphase ISM. Additionally, the dynamical impact of cosmic rays (CRs) on star formation rates (SFRs) must also be considered. We present the first simulations of the formation of an ensemble of star clusters with dynamically-coupled CRs, near-individual star particles, and a feedback-regulated ISM. We analyze tallbox simulations performed using the CRISP model in the moving-mesh code AREPO. We apply varied implementations of CR transport under the theory of self-confinement. We find that CRs simultaneously reduce the SFR, the power law slope of the cluster mass function, and the cluster formation efficiency. Each simulation is compatible with observations, and CR feedback tends to move results along observed star cluster relations. We see only modest changes in cluster radius and velocity dispersions, but significant differences in the virial parameters. Ultimately, the primary impact of CRs is to reduce SFRs. Lower SFRs imply fewer supernovae, and consequently a lower turbulent energy budget for gas. Star clusters formed in a CR-regulated ISM have lower velocity dispersions, and are therefore more bound under self-gravity. The effective clustering of supernovae is unchanged by CRs. Despite the idealized setup, the CRISP feedback model recovers many key aspects of star cluster formation.