Nuclear spin-spin coupling over van der Waals bond has recently been observed via the frequency shift of solute protons in a solution containing optically hyperpolarized 129 Xe nuclei. We carry out a first-principles computational study of the prototypic van der Waals-bonded xenon dimer, where the spin-spin coupling between two magnetically non-equivalent isotopes, J( 129 Xe − 131 Xe), is observable. We use relativistic theory at the four-component Dirac-Hartree-Fock and Diracdensity-functional theory levels using novel completeness-optimized Gaussian basis sets and choosing the functional based on a comparison with correlated ab initio methods at the nonrelativistic level. J-coupling curves are provided at different levels of theory as functions of the internuclear distance in the xenon dimer, demonstrating cross-coupling effects between relativity and electron correlation for this property. Calculations on small Xe clusters are used to estimate the importance of many-atom effects on J( 129 Xe − 131 Xe). Possibilities of observing J( 129 Xe − 131 Xe) in liquid xenon are critically examined, based on molecular dynamics simulation. A simplistic spherical model is set up for the xenon dimer confined in a cavity, such as in microporous materials. It is shown that the on the average shorter internuclear distance enforced by the confinement increases the magnitude of the coupling as compared to the bulk liquid case, rendering J( 129 Xe − 131 Xe) in a cavity a feasible target for experimental investigation.