IR-IR double-resonance experiments were used to study the state-to-state rotational relaxation of CO with Ne as a collision partner. Rotational levels in the range J(i)=2-9 were excited and collisional energy transfer of population to the levels J(f)=2-8 was monitored. The resulting data set was analyzed by fitting to numerical solutions of the master equation. State-to-state rate constant matrices were generated using fitting law functions. Fitting laws based on the modified exponential gap (MEG) and statistical power exponential gap (SPEG) models were used; the MEG model performed better than the SPEG model. A rate constant matrix was also generated from scattering calculations that employed the ab initio potential energy surface of McBane and Cybulski [J. Chem. Phys. 110, 11 734 (1999)]. This theoretical rate constant matrix yielded kinetic simulations that agreed with the data nearly as well as the fitted MEG model and was unique in its ability to reproduce both the rotational energy transfer and pressure broadening data for Ne-CO. The theoretical rate coefficients varied more slowly with the energy gap than coefficients from either of the fitting laws. (C) 2004 American Institute of Physics.