Molecular properties govern how molecules interact with one another or external fields. Accurate molecular properties are essential for constructing intermolecular interaction models, and to achieve high accuracy we need to go beyond leading-order (dipolar) terms. This work presents CCSD(T) references for 73 small non-spin-polarized molecules for molecular dipole and quadrupole moments, and dipole-dipole and quadrupole-quadrupole polarizabilities. Using these results, we provide a holistic performance analysis of molecular properties computed with HF, MP2, CCSD, and a wide range of density-functional methods. Additionally, we investigate in detail what levels of basis sets are required to simultaneously describe all properties in the data set. The best-performing density functionals for all four properties are the asymptotically cor- rected hybrid GGAs B97-3-AC, closely followed by PBE0-AC. Surprisingly, modern meta- GGAs and range-separated functionals perform inconsistently, with large errors in the quadrupo- lar properties. These results have direct implications for methods for intermolecular interac- tions, such as symmetry-adapted perturbation theory based on density functional theory, or for the quality of properties and interactions in machine-learning datasets. Finally, for DFT methods, the Jensen aug-pcseg-2 basis set is a computationally efficient alternative to more established basis sets, but for generating correlated references, third-row elements benefit from core polarization, and double augmentation is strictly necessary for quadrupole-quadrupole polarizabilities.