We compute next-to-next-to-leading order (NNLO) QCD corrections to neutral vector boson production in association with a charm jet at the LHC. This process is studied in the forward kinematics at s = 13  TeV, which may provide valuable constraints on the intrinsic charm component of the proton. A comparison is performed between fixed order and NLO predictions matched to a parton shower showing mutual compatibility within the respective uncertainties. NNLO corrections typically lead to a reduction of theoretical uncertainties by a factor of two and the perturbative convergence is further improved through the introduction of a theory-inspired constraint on the transverse momentum of the vector boson plus jet system. A comparison between these predictions with data will require an alignment of a flavour-tagging procedure in theory and experiment that is infrared and collinear safe.