Forward neutrinos from charm at the Large Hadron Collider

The currently operating FASER experiment and the planned Forward Physics Facility (FPF) will detect a large number of neutrinos produced in proton-proton collisions at the LHC. In this work, we estimate neutrino fluxes at these detectors from charm meson decays, which will be particularly important for the nu e and nu tau channels. We make prediction using both the next-to-leading order collinear factorization and the kT-factorization approaches to model the production of charm quarks as well as different schemes to model their hadronization into charm hadrons. In particular, we emphasize that a sophisticated modeling of hadronization involving beam remnants is needed for predictions at FASER and FPF due to the sensitivity to the charm hadron production at low transverse momenta and very forward rapidity. As example, we use the string fragmentation approach implemented in PYTHIA 8. While both standard fragmentation functions and PYTHIA 8 are able to describe LHCb data, we find that PYTHIA 8 predicts significantly higher rate of high energy neutrinos, highlighting the importance of using the correct hadronization model when making predictions.