Highly efficient perturbative + variational strategy based on orthogonal valence bond theory for the evaluation of magnetic coupling constants. Application to the trinuclear Cu(II) site of multicopper oxidases