The recent observation of very high energy cosmic neutrinos by IceCube heralds the beginning of neutrino astronomy. At these energies, the dominant background to the astrophysical signal is the flux of `prompt' neutrinos, arising from the decay of charmed mesons produced by cosmic ray collisions in the atmosphere. In this work we provide predictions for the prompt atmospheric neutrino flux in the framework of perturbative QCD, using state-of-the-art Monte Carlo event generators. Our calculation includes the constraints set by charm production measurements from the LHCb experiment at 7 TeV, and has been recently validated with the corresponding 13 TeV data. Our results for the prompt flux are a factor of about 2 below the previous benchmark calculation, in general agreement with two other recent estimates, and with an improved estimate of the uncertainty. This alleviates the existing tension between the theoretical prediction and IceCube limits, and suggests that a direct direction of the prompt flux is imminent.

The recent observation by the IceCube experiment of cosmic neutrinos at energies up to a few PeV heralds the beginning of neutrino astronomy. At such high energies, the conventional neutrino flux is suppressed and the prompt component from charm meson decays is expected to become the dominant background to astrophysical neutrinos. Charm production at high energies is however theoretically uncertain, both since the charm mass is at the boundary of applicability of perturbative QCD, and also because the calculations are sensitive to the poorly-known gluon PDF at small-x. In this work we provide detailed perturbative QCD predictions for charm and bottom production in the forward region, and validate them by comparing with recent data from the LHCb experiment at 7 TeV. Finding good agreement between data and theory, we use the LHCb measurements to constrain the small-x gluon PDF, achieving a substantial reduction in its uncertainties. Using these improved PDFs, we provide predictions for charm and bottom production at LHCb at 13 TeV, as well as for the ratio of cross-sections between 13 and 7 TeV. The same calculations are used to compute the energy distribution of neutrinos from charm decays in pA collisions, a key ingredient towards achieving a theoretically robust estimate of charm-induced backgrounds at neutrino telescopes.