Multireference perturbation theory (MRPT) with multiconfigurational self-consistent field (MCSCF) reference functions is applied to the calculations of core-electron binding energies (CEBEs) of atoms and molecules. Orbital relaxations in a core-ionized state and electron correlation are both taken into account in a conventional MCSCF-MRPT procedure. In the MCSCF calculation, the target core ionized state is directly optimized as an excited state and this treatment can completely prevent a variational collapse. Multireference Moller-Plesset perturbation theory and multiconfigurational self-consistent field reference quasidegenerated perturbation theory were used to treat electron correlation. The present method quite accurately reproduced the 1s CEBEs of CH4, NH3, H2O, and FH; the average deviation from the experimental data is 0.11 eV using Ahlrichs' VTZ basis set. The C 1s and O 1s CEBEs of formic acid and acetic acid were calculated and the results are consistent with the bonding characters of the atoms in these molecules. The present procedure can also be applied to CEBEs of higher angular momentum orbitals by including spin-orbit coupling. The calculated CEBEs of Ar 2p, HCl 2p, Kr 3d, and HBr 3d are in reasonable agreement with the available experimental values. In the calculation of the 3d CEBEs, a relativistic correction significantly improves the agreements. The effect of polarization functions is also discussed. (C) 2004 American Institute of Physics.