Introduction of a New Theory for the Calculation of Magnetic Coupling Based on Spin-Flip Constricted Variational Density Functional Theory. Application to Trinuclear Copper Complexes which Model the Native Intermediate in Multicopper Oxidases.