ABSTRACTSeveral coupled-cluster approaches with an approximate treatment of connected triply or quadruply excited clusters, including the CCSD(T), , and CR-CC(2,3) corrections to CCSD, the active-space CCSDt and CCSDTq methods, and the CC(t;3) and CC(q;4) corrections to CCSDt and CCSDTq derived using the CC(P;Q) framework, and their CCSDT and CCSDTQ parents are applied to the ground-state potential curve and vibrational term values of the magnesium dimer. The correlation-consistent aug-cc-pV(n+d)Z and aug-cc-pwCVnZ ( and Q) basis sets are used. Among the noniterative triples corrections to CCSD, the CR-CC(2,3) approach performs the best, but, in analogy to the previously studied beryllium dimer [I. Magoulas et al., J. Phys. Chem. A 122, 1350 (2018)], the CC(t;3) and CC(q;4) methods outperform other employed approaches in reproducing the CCSDT and CCSDTQ data. Composite calculations combining the nearly all-electron CCSDT and valence CC(q;4) and CCSDTQ computations reproduce the experimentally resolved part of the vibrational spectrum of to within ∼1 when the vibrational spacings are examined, while predicting five extra vibrational states near the dissociation threshold, which have not been experimentally resolved. The corresponding binding energies and equilibrium bond lengths are within a few and 0.01 Å from their experimentally derived values.