The global ground and first three excited electronic state adiabatic as well as the corresponding quasidiabatic potential energy surfaces is reported as a function of nuclear geometries in the Jacobi coordinates (R →, r →, γ) using Dunning's cc-pVTZ basis set at the internally contracted multi-reference (single and double) configuration interaction level of accuracy. Nonadiabatic couplings, arising out of relative motion of proton and the vibrational motion of CO, are also reported in terms of coupling potentials. The quasidiabatic potential energy surfaces and the coupling potentials have been obtained using the ab initio procedure [Simah et al., J. Chem. Phys. 111, 4523 (1999)] for the purpose of dynamics studies. © 2016 AIP Publishing LLC.