Abstract
We apply periodic domain-based local pair natural orbital second-order Møller-Plesset perturbation theory (DLPNO-MP2) to probe the adsorption energy of CO on MgO(001), the consensus model system for surface adsorption. A number of robust correlated wavefunction methods now achieve excellent agreement with experiment for the adsorption of a single CO molecule onto the MgO surface. However, studies probing denser coverage ratios are scarce because of the increased computational expense and the larger configuration space needed for optimization. We leveraged the computational efficiency of periodic DLPNO-MP2 to perform simulations beyond a single unit cell. By using large supercells, we highlight the importance of accurately representing the thermodynamic limit of the surface and demonstrate in turn that different coverage ratios can be consistently probed. In the dilute regime, we show that the adsorption energies obtained from periodic DLPNO-MP2 agree with existing benchmarks. We then obtain adsorption energies at increasing densities, approaching full monolayer coverage. Our results show a reduction in binding strength at full coverage, agreeing with experimental observations, which is explained by the increasing lateral repulsions between the COs. This study demonstrates the efficacy of periodic DLPNO-MP2 for probing increasingly sophisticated adsorption systems at the thermodynamic limit.