The aim of this work was to determine and understand the origin of the electronic properties of Mn(IV) complexes, especially the zero-field splitting (ZFS), through a combined experimental and theoretical investigation on five well-characterized mononuclear octahedral Mn(IV) compounds, with various coordination spheres (N6, N3O3, N2O4 in both trans (trans-N2O4) and cis configurations (cis-N2O4) and O4S2). High-frequency and -field EPR (HFEPR) spectroscopy has been applied to determine the ZFS parameters of two of these compounds, MnL(trans-N2O4) and MnL(O4S2). While at X-band EPR, the axial-component of the ZFS tensor, D, was estimated to be +0.47 cm(-1) for MnL(O4S2), and a D-value of +2.289(5) cm(-1) was determined by HFEPR, which is the largest D-magnitude ever measured for a Mn(IV) complex. A moderate D value of -0.997(6) cm(-1) has been found for MnL(trans-N2O4). Quantum chemical calculations based on two theoretical frameworks (the Density Functional Theory based on a coupled perturbed approach (CP-DFT) and the hybrid Ligand-Field DFT (LF-DFT)) have been performed to define appropriate methodologies to calculate the ZFS tensor for Mn(IV) centers, to predict the orientation of the magnetic axes with respect to the molecular ones, and to define and quantify the physical origin of the different contributions to the ZFS. Except in the case of MnL(trans-N2O4), the experimental and calculated D values are in good agreement, and the sign of D is well predicted, LF-DFT being more satisfactory than CP-DFT. The calculations performed on MnL(cis-N2O4) are consistent with the orientation of the principal anisotropic axis determined by single-crystal EPR, validating the calculated ZFS tensor orientation. The different contributions to D were analyzed demonstrating that the d-d transitions mainly govern D in Mn(IV) ion. However, a deep analysis evidences that many factors enter into the game, explaining why no obvious magnetostructural correlations can be drawn in this series of Mn(IV) complexes.