We have studied the magnetic structure of the high symmetry vanadyl pyrophosphate ((VO) 2 P 2 O 7 , VOPO), focusing on the spin exchange couplings, using density functional theory (B3LYP) with the full three-dimensional periodicity. VOPO involves four distinct spin couplings: two larger couplings exist along the chain direction (a-axis), which we predict to be antiferromagnetic, J OPO = À156.8 K and J O = À68.6 K, and two weaker couplings appear along the c (between two layers) and b directions (between two chains in the same layer), which we calculate to be ferromagnetic, J layer = 19.2 K and J chain = 2.8 K. Based on the local density of states and the response of spin couplings to varying the cell parameter a, we found that J OPO originates from a super-exchange interaction through the bridging -O-P-O-unit. In contrast, J O results from a direct overlap of 3d x 2 À y 2 orbitals on two vanadium atoms in the same V 2 O 8 motif, making it very sensitive to structural fluctuations. Based on the variations in V-O bond length as a function of strain along a, we found that the V-O bonds of V-(OPO) 2 -V are covalent and rigid, whereas the bonds of V-(O) 2 -V are fragile and dative. These distinctions suggest that compression along the a-axis would have a dramatic impact on J O , changing the magnetic structure and spin gap of VOPO. This result also suggests that assuming J O to be a constant over the range of 2-300 K whilst fitting couplings to the experimental magnetic susceptibility is an invalid method. Regarding its role as a catalyst, the bonding pattern suggests that O 2 can penetrate beyond the top layers of the VOPO surface, converting multiple V atoms from the +4 to +5 oxidation state, which seems crucial to explain the deep oxidation of n-butane to maleic anhydride.