First-principles calculations of zero-field splitting parameters

In this work, an implementation of an approach to calculate the zero-field splitting (ZFS) constants in the framework of ab initio methods such as complete active space self-consistent field, multireference configuration interaction, or spectroscopy oriented configuration interaction is reported. The spin-orbit coupling (SOC) contribution to ZFSs is computed using an accurate multicenter mean-field approximation for the Breit-Pauli Hamiltonian. The SOC parts of ZFS constants are obtained directly after diagonalization of the SOC operator in the basis of a preselected number of roots of the spin-free Hamiltonian. This corresponds to an infinite order treatment of the SOC in terms of perturbation theory. The spin-spin (SS) part is presently estimated in a mean-field fashion and appears to yield results close to the more complete treatments available in the literature. Test calculations for the first- and second-row atoms as well as first-row transition metal atoms and a set of diatomic molecules show accurate results for the SOC part of ZFSs. SS contributions have been found to be relatively small but not negligible (exceeding 1 cm(-1) for oxygen molecule). At least for the systems studied in this work, it is demonstrated that the presented method provides much more accurate estimations for the SOC part of ZFS constants than the emerging density functional theory approaches.     

Zero-field splittings and local tilting angles tauMn2+ for Mn2+ in ZnGeP2 and CdGeP2 crystals

The electron paramagnetic resonance (EPR) zero-field splittings (ZFSs) D of Mn2+ in ZnGeP2 and CdGeP2 crystals are calculated from both the microscopic spin-orbit coupling mechanism and the empirical superposition model. From the calculations, the ZFS D of ZnGeP2:Mn2+ is reasonably explained by using the local tilting angle tauMn2+ (rather than the corresponding angle tauZn2+ in the host crystal) and the local tilting angle tauMn2+ (which has not been reported) in CdGeP2:Mn2+ is estimated. The intrinsic ZFS parameter b2(R0) approximately -0.052(6)cm(-1) (with R0 approximately 2.43 angstroms) is suggested for Mn2+-P(3-) combination by using the local tilting angles tauMn2+. The value is quite unlike that (approximately 0.4(2)cm(-1)) obtained in the previous paper by combining the crystallographic data of host crystals with the ZFSs for Mn2+ ions in crystals. The reasonableness of these results is discussed.     

Ab initio and coupled-perturbed density functional theory estimation of zero-field splittings in MnII transition metal complexes

The paper presents a method comparison for the prediction of zero-field splitting (ZFS) parameters in a series of Mn (II) coordination complexes. The test set consists of Mn (II) complexes that are experimentally well-characterized by X-ray diffraction and high-field electron paramagnetic resonance. Their ZFS parameters have been calculated using density functional theory (DFT) as well as complete active space self-consistent field (CASSCF) methods. It is shown that the recently introduced coupled-perturbed spin-orbit coupling (CP-SOC) approach [ Neese, F. J. Chem. Phys. 2007, 127, 164112 ] together with hybrid-DFT functionals leads to a slope of the correlation line (plot of experimental vs calculated D values) that is essentially unity provided that the direct spin-spin interaction is properly included in the treatment. This is different from our previous DFT study on the same series of complexes where a severe overestimation of the D parameter has been found [ Zein, S. ; Duboc, C. ; Lubitz, W. ; Neese, F. Inorg. Chem. 2008, 47, 134 ]. CASSCF methods have been used to evaluate the ZFS in an "ab initio ligand-field" type treatment. The study demonstrates that a substantial part of the relevant physics is lost in such a treatment since only excitations within the manganese d-manifold are accounted for. Thus, a severe underestimation of the D parameter has been found. Because the CASSCF calculations in combination with quasidegenerate perturbation theory treats the SOC to all orders, we have nevertheless verified that second-order perturbation theory is an adequate approximation in the case of the high-spin d (5) configuration.     

Theoretical investigation of zero field splitting parameters for Mn2+ centres in L-asparagine monohydrate

Zero-field splitting (ZFS) parameters D and E are calculated using the point-charge model (PCM) and superposition model (SPM) for Mn(2+) centre in L-asparagine monohydrate (LAM) single crystal. The calculated ZFS parameters obtained using these two models are compared with the experimental values for interstitial site of Mn(2+). The SPM and PCM give ZFS parameters similar to those of experimental ones. This supports the notion that the impurity ion occupies interstitial site in LAM.     

Calculation of the zero-field splitting tensor on the basis of hybrid density functional and Hartree-Fock theory

The zero-field splitting (ZFS) (expressed in terms of the D tensor) is the leading spin-Hamiltonian parameter for systems with a ground state spin S>12. To first order in perturbation theory, the ZFS arises from the direct spin-spin dipole-dipole interaction. To second order, contributions arise from spin-orbit coupling (SOC). The latter contributions are difficult to treat since the SOC mixes states of different multiplicities. This is an aspect of dominant importance for the correct prediction of the D tensor. In this work, the theory of the D tensor is discussed from the point of view of analytic derivative theory. Starting from a general earlier perturbation treatment [F. Neese and E. I. Soloman, Inorg. Chem. 37, 6568 (1998)], straightforward response equations are derived that are readily transferred to the self-consistent field (SCF) Hartree-Fock (HF) or density functional theory (DFT) framework. The main additional effort in such calculations arises from the solution of nine sets of nonstandard coupled-perturbed SCF equations. These equations have been implemented together with the spin-orbit mean-field representation of the SOC operator and a mean-field treatment of the direct spin-spin interaction into the ORCA electronic structure program. A series of test calculations on diatomic molecules with accurately known zero-field splittings shows that the new approach corrects most of the shortcomings of previous DFT based methods and, on average, leads to predictions within 10% of the experimental values. The slope of the correlation line is essentially unity for the B3LYP and BLYP functionals compared to approximately 0.5 in previous treatments.     

Local structure of Mn4+ and Fe3+ spin probes in layered LiAlO2 oxide by modelling of zero-field splitting parameters

The zero-field splitting parameters (ZFS) of Mn(4+) and Fe(3+) ions in LiAlO(2) with a layered structure are analyzed experimentally and theoretically by using high-frequency electron paramagnetic resonance spectroscopy, Neuman superposition model (NSM), DFT and multiconfigurational calculations. The interpretation of ZFS is based on the comparison of the experimentally determined values with the calculated ones. This approach allows assessing the performance of different methods for computation of ZFS of Fe(3+) and Mn(4+) in layered oxide matrices. DFT and multiconfigurational calculations are used to analyze the effect of oxygen, aluminium, and lithium neighbours on ZFS of Fe(3+) and Mn(4+). These calculations are based on a cluster comprising Fe(3+) or Mn(4+) ions in a trigonally compressed octahedron with 6 metal ions (Al(3+) or Co(3+)) as first metal neighbours and 6 O(2-) and 2 Li(+) (above and below the layer) as second neighbours. A satisfactory agreement with the experimental data is achieved when the local structure of Mn(4+) and Fe(3+) deviates from the trigonal host-site geometry. The local structure of Fe(3+) comprises an axial distortion, while trigonal environment with reduced extent of distortion appears around Mn(4+).     

Investigations of zero-field splitting (ZFS) and local distortion parameters of ZnAl2S4:Mn2+ by theoretical analysis

Fourth-order perturbation formula on the basis of the dominant spin-orbit coupling mechanism is employed to investigate the local environment around Mn2+ centers in ZnAl2S4 single crystals. The zero-field splitting (ZFS) parameter D is calculated for the Mn2+ ions at the Al3+ site with local symmetry D3d using the different orbital reduction factors. Both the contributions of the lattice distortions to the crystal-field (CF) parameters and the D are examined by means of different cases. The comparison between the calculated results in this study and the previous experimental and theoretical values reveals a good agreement and reasonable distortion parameters for Mn2+ ions at Al3+ sites.     

Magnetic susceptibility and ground-state zero-field splitting in high-spin mononuclear manganese(III) of inverted N-methylated porphyrin complexes: Mn(2-NCH3NCTPP)Br

The crystal structures of diamagnetic dichloro(2-aza-2-methyl-5,10,15,20-tetraphenyl-21-carbaporphyrinato-N,N',N')-tin(IV) methanol solvate [Sn(2-NCH 3NCTPP)Cl 2.2(0.2MeOH); 6.2(0.2MeOH)] and paramagnetic bromo(2-aza-2-methyl-5,10,15,20-tetraphenyl-21-carbaporphyrinato-N,N',N')-manganese(III) [Mn(2-NCH 3NCTPP)Br; 5] were determined. The coordination sphere around Sn (4+) in 6.2(0.2MeOH) is described as six-coordinate octahedron ( OC-6) in which the apical site is occupied by two transoid Cl (-) ligands, whereas for the Mn (3+) ion in 5, it is a five-coordinate square pyramid ( SPY-5) in which the unidentate Br (-) ligand occupies the axial site. The g value of 9.19 (or 10.4) measured from the parallel polarization (or perpendicular polarization) of X-band EPR spectra at 4 K is consistent with a high spin mononuclear manganese(III) ( S = 2) in 5. The magnitude of axial ( D) and rhombic ( E) zero-field splitting (ZFS) for the mononuclear Mn(III) in 5 were determined approximately as -2.4 cm (-1) and -0.0013 cm (-1), respectively, by paramagnetic susceptibility measurements and conventional EPR spectroscopy. Owing to weak C(45)-H(45A)...Br(1) hydrogen bonds, the mononuclear Mn(III) neutral molecules of 5 are arranged in a one-dimensional network. A weak Mn(III)...Mn(III) ferromagnetic interaction ( J = 0.56 cm (-1)) operates via a [Mn(1)-C(2)-C(1)-N(4)-C(45)-H(45A)...Br(1)-Mn(1)] superexchange pathway in complex 5.     

Spin-spin contributions to the zero-field splitting tensor in organic triplets, carbenes and biradicals-a density functional and ab initio study

An evaluation study for the direct dipolar electron spin-spin (SS) contribution to the zero-field splitting (ZFS) tensor in electron paramagnetic resonance (EPR) spectroscopy is presented. Calculations were performed on a wide variety of organic systems where the SS contribution to the ZFS dominates over the second-order spin-orbit coupling (SOC) contribution. Calculations were performed using (hybrid) density functional theory (DFT), as well as complete active space self-consistent field (CASSCF) wave functions. In the former case, our implementation is an approximation, because we use the two-particle reduced spin-density matrix of the noninteracting reference system. In the latter case, the SS contribution is approximated by a mean-field method which, nevertheless, gives accurate results, compared to the approximation free computation of the SS part in a CASSCF framework. For the case of the triplet dioxygen molecule, it was shown that restricted open-shell density functional theory (RODFT), as well as CASSCF, can provide accurate spin-spin couplings while spin-unrestricted DFT leads to much larger errors. Furthermore, 15 organic radicals, including several 1,3 and 1,5 diradicals, dinitroxide biradicals, and even a chlorophyll a model system, were examined as test cases to demonstrate the accuracy and efficiency of our approach within a DFT framework. Accurate D values with root-mean-square deviations of 0.0035 cm(-1) were obtained. Furthermore, all trends, including those due to substituent effects, were correctly reproduced. In a different set of calculations, the polyacenes benzene, naphthalene, anthracene, and tetracene were studied. Applying DFT, the absolute D values were noticeably underestimated, but it was possible to correctly reproduce the trend to smaller D values with larger size of the systems. Finally, it was demonstrated that our approach is also well-suited for the study of carbenes. The smaller organic radicals of this work were also studied, through the use of CASSCF wave functions. This was a special advantage in the case of the triplet polyacenes, where the CASSCF approach gave better results than the DFT method. In comparing spin-restricted and spin-unrestricted results, it was shown through a natural orbital analysis and comparison to high-level ab initio calculations that even small amounts of spin polarization introduced by the unrestricted calculations lead to large deviations between the unrestricted Kohn-Sham (UKS) and restricted open-shell Kohn-Sham (ROKS) approaches. It is challenging to understand why the ROKS results show much better correlation with the experimental data.     

Ab initio and DFT studies of the spin-orbit and spin-spin contributions to the zero-field splitting tensors of triplet nitrenes with aryl scaffolds

Spin-orbit and spin-spin contributions to the zero-field splitting (ZFS) tensors (D tensors) of spin-triplet phenyl-, naphthyl-, and anthryl-nitrenes in their ground state are investigated by quantum chemical calculations, focusing on the effects of the ring size and substituted position of nitrene on the D tensor. A hybrid CASSCF/MRMP2 approach to the spin-orbit term of the D tensor (D(SO) tensor), which was recently proposed by us, has shown that the spin-orbit contribution to the entire D value, termed the ZFS parameter or fine-structure constant, is about 10% in all the arylnitrenes under study and less depends on the size and connectivity of the aryl groups. Order of the absolute values for D(SO) can be explained by the perturbation on the energy level and spatial distributions of π-SOMO through the orbital interaction between SOMO of the nitrene moiety and frontier orbitals of the aryl scaffolds. Spin-spin contribution to the D tensor (D(SS) tensor) has been calculated in terms of the McWeeny-Mizuno equation with the DFT/EPR-II spin densities. The D(SS) value calculated with the RO-B3LYP spin density agrees well with the D(Exptl) -D(SO) reference value in phenylnitrene, but agreement with the reference value gradually becomes worse as the D value decreases. Exchange-correlation functional dependence on the D(SS) tensor has been explored with standard 23 exchange-correlation functionals in both RO- and U-DFT methodologies, and the RO-HCTH/407 method gives the best agreement with the D(Exptl) -D(SO) reference value. Significant exchange-correlation functional dependence is observed in spin-delocalized systems such as 9-anthrylnitrene (6). By employing the hybrid CASSCF/MRMP2 approach and the McWeeny-Mizuno equation combined with the RO-HCTH/407/EPR-II//U-HCTH/407/6-31G* spin densities for D(SO) and D(SS), respectively, a quantitative agreement with the experiment is achieved with errors less than 10% in all the arylnitrenes under study. Guidelines to the putative approaches to D(SS) tensor calculations are given.     

Unusual structural types in manganese cluster chemistry from the use of N,N,N',N'-tetrakis(2-hydroxyethyl)ethylenediamine: Mn8, Mn12, and Mn20 Clusters

The syntheses, crystal structures, and magnetochemical characterization are reported for three new mixed-valent Mn clusters [Mn(8)O(3)(OH)(OMe)(O(2)CPh)7(edte)(edteH(2))](2)CPh) (1), [Mn(12)O(4)(OH)(2)(edte)(4)C(l6)(H(2)O)(2)] (2), and [Mn(20)O(8)(OH)(4)(O(2)CMe)(6)(edte)(6)](ClO(4))(2) (3) (edteH(4) = (HOCH(2)CH(2))(2)NCH(2)CH(2)N(CH(2)CH(2)OH)(2) = N,N,N',N'-tetrakis(2-hydroxyethyl)ethylenediamine). The reaction of edteH(4) with Mn(O(2)CPh)(2), MnCl(2), or Mn(O(2)CMe)(2) gives 1, 2, and 3, respectively, which all possess unprecedented core topologies. The core of 1 comprises two edge-sharing [Mn(4)O(4)] cubanes connected to an additional Mn ion by a micro(3)-OH- ion and two alkoxide arms of edteH(22-). The core of 2 consists of a [Mn(12)(micro(4-)O)(4)](24+) unit with S4 symmetry. The core of 3 consists of six fused [Mn(4)O(4)] cubanes in a 3 x 2 arrangement and linked to three additional Mn atoms at both ends. Variable-temperature, solid-state dc and ac magnetization (M) studies were carried out on complexes 1-3 in the 5.0-300 K range. Fitting of the obtained M/Nmicro(B) vs H/T data by matrix diagonalization and including only axial zero-field splitting (ZFS) gave ground-state spin (S) and axial ZFS parameter (D) of S = 8, D = -0.30 cm-1 for 1, S = 7, D = -0.16 cm-1 for 2, and S = 8, D = -0.16 cm-1 for 3. The combined work demonstrates that four hydroxyethyl arms on an ethylenediamine backbone can generate novel Mn structural types not accessible with other alcohol-based ligands.     

Origin of the zero-field splitting in mononuclear octahedral dihalide MnII complexes: an investigation by multifrequency high-field electron paramagnetic resonance and density functional theory

The synthesis, structural characterization, and electronic properties of a new series of high-spin six-coordinate dihalide mononuclear MnII complexes [Mn(tpa)X2] (tpa=tris-2-picolylamine; X=I (1), Br (2), and Cl (3)) are reported. The analysis of the crystallographic data shows that in all investigated complexes the manganese ion lies in the center of a distorted octahedron with a cis configuration of the halides imposed by the tpa ligand. By a multifrequency high-field electron paramagnetic resonance investigation (95-285 GHz), the electronic properties of 1-3 were determined (DI=-0.600, DBr=-0.360, DCl=+0.115 cm-1), revealing the important effect of (i) the nature of the halide and (ii) the configuration (cis/trans) of the two halides on the magnitude of D. The spin Hamiltonian parameters obtained by density functional theory calculations initiated from the crystal structure of 1-3 are in reasonable agreement with the experimental values. The absolute value of D is consistently overestimated, but the sign and the trend over the chemical series is well reproduced. Theoretical models (cis- and trans-[Mn(NH3)4X2], X=I, Br, Cl and F) have been used to investigate the different contributions to D and also to understand the origin of the experimentally observed changes in D within the series reported here. This study reveals that the spin-spin coupling contributions to the D tensor are non-negligible for the lighter halides (F, Cl) but become insignificant for the heavier halides (I, Br). The four different types of excitations involved in the spin-orbit coupling (SOC) part of the D tensor contribute with comparable magnitudes and opposing signs. The general trend observed for halide MnII complexes (DI>DBr>DCl) can be explained by the fact that the halide SOC dominates the D value in these systems with a major contribution arising from interference between metal- and halide-SOC contributions, which are proportional to the product of the SOC constants of Mn and X.     

Multifrequency electron paramagnetic resonance and theoretical studies of a Mn(II) (S = 5/2) complex: The role of geometrical elements on the Zero Field Splitting parameters

The Zero Field Splitting (ZFS) parameters of the mononuclear Mn(II) (S = 5/2) compound [Mn(Ac4Ph)₂], where HAc4Ph stands for 4-phenyl-2-acetylpyridine thiosemicarbazone, are determined by dual mode X-band and high field/high frequency electron paramagnetic resonance spectroscopy. Density functional theory (DFT) calculations reproduce both the sign and absolute value of the axial parameter D of the ZFS tensor, whose distribution is quantitatively correlated to distribution on geometrical elements of the complex.     

Studies of a linear single-molecule magnet

Reaction of the dinuclear complex [Mn2O2(bpy)4](ClO4)3 with H3cht (cis,cis-l,3,5-cyclohexanetriol) in MeCN produces the complex [Mn3(Hcht)2(bpy)4](ClO4)3.Et2O.2MeCN (1.Et2O.2MeCN). Dc magnetic susceptibility measurements reveal the existence of weak ferromagnetic exchange between the three Mn ions, leading to a spin ground state of S = 7, with D = -0.23 cm(-1). W-Band (94 GHz) EPR measurements on restrained powdered crystalline samples confirm the S = 7 ground state and determine the ground state zero-field splitting (ZFS) parameters of D = -0.14 cm(-1) and B4(0)= +1.5 x 10(-5) cm(-1). The apparent 4th order behaviour is due to a breakdown of the strong exchange limit approximation (J approximately d, the single-ion ZFS). Single crystal dc relaxation decay and hysteresis loop measurements reveal the molecule to have an appreciable energy barrier to magnetization relaxation, displaying low temperature sweep rate and temperature-dependent hysteresis loops. Density functional studies confirm the ferromagnetic exchange coupling between the Mn ions.     

Low-valent low-coordinated manganese(I) ion dimer: a temperature dependent W-band EPR study

W-Band EPR spectra of [[HC(CMeNAr)(2)]Mn](2) (Ar = 2,6-(i)Pr(2)C(6)H(3)) have been measured at different temperatures. The spectra show a behavior which is typical for an antiferromagnetically coupled dimer with excited states populating upon increasing temperature. By following the intensity variation of the different features of the spectra with temperature, we attributed different groups of resonances to the S = 1, 2, and 3 states of the dimer. Their corresponding spin Hamiltonian parameters were derived from simulations. The zero-field-splitting parameters measured in this way were D(S=1) = 1.57 cm(-1) and E(S=1) = 0.064 cm(-1), D(S=2) = 0.266 cm(-1) and E(S=2) = 0.0045 cm(-1), and D(S=3) = 0.075 cm(-1) and E(S=3) = 0. On the basis of the molecular structure of the system, we could estimate that zero-field splitting (ZFS) is the result of anisotropic exchange and single-ion anisotropic contributions of similar magnitude (|D| approximately 0.2 cm(-1)). These results allow a deeper insight into the electronic structure of the Mn(I) centers in low-coordination environments, further supporting the electronic structure of Mn(I) to be 4s(1)3d(5), as previously indicated by DFT calculations.     

Optical determination of the single-ion zero-field splitting in large spin clusters

The dodecametallic Cr(III) cluster has an S = 6 ground state with an axial zero-field splitting (ZFS) of DS=6 = +0.088 cm-1. Analysis of high-resolution optical data (MCD) allows us to determine the single-ion ZFS of the constituent Cr(III) ions directly (D = -1.035 cm-1). A vector coupling analysis demonstrates that the cluster ZFS is almost entirely single-ion in origin. Thus, the relative orientations of the local and cluster magnetic axes can lead to cluster ZFS of opposite sign to the single-ion even when this is the only significant contribution.     

Strain effects in electron spin resonance spectroscopy of quintet 2,6-bis(4'-nitrenophenyl)-4-phenylpyridine

Photolysis of 2,6-bis(4'-azidophenyl)-4-phenylpyridine in 2-methyltetrahydrofuran (2MTHF) glass at 7 K leads to quintet 2,6-bis(4'-nitrenophenyl)-4-phenylpyridine as a mixture of rotational isomers. The electron spin resonance (ESR) spectrum of this mixture of rotamers shows a considerable broadening of many transitions in the range of 0-5000 G and cannot be reproduced by computer simulations solely based on the tuning of the spin Hamiltonian parameters g, D(Q), and E(Q) alone or on predictions of DFT calculations. The best modeling of the experimental ESR spectrum is obtained only when the large line-broadening parameter of Γ(E(Q)) = 1200 MHz along with the spin Hamiltonian g = 2.003, D(Q) = 0.154 cm(-1), and E(Q) = 0.050 cm(-1) is used in the spectral simulations. The most accurate theoretical estimations of the magnetic parameters of the dinitrene in a 2MTHF glass are obtained from the B3LYP/6-311+G(d,p)+PBE/DZ/COSMO calculations of the spin-spin coupling parameters D(SS) and E(SS). Such calculations overestimate the E(Q) and D(Q) values of the dinitrene just by 1% and 10%, respectively, demonstrating that contributions of the spin-orbit coupling parameters D(SOC) and E(SOC) to the total D(Q) and E(Q) values are negligibly small. The research shows that ESR studies of polynuclear high-spin nitrenes, obtained by photolysis of rotational isomers of the starting azides, can only be successful if large E(Q) strain effects are taken into account in the spectral simulations.     

Spin ground state and magnetic properties of cobalt(II): relativistic DFT calculations guided by EPR measurements of bis(2,4-acetylacetonate)cobalt(II)-based complexes

The spin ground state of the core ion and structure of the bis(2,4-acetylacetonate)cobalt(II) model complex and its synthetic aqua and ethanol derivatives, Co(acac)(2)L(n), (L = EtOH, H(2)O), were examined by means of density functional theory (DFT) calculations supported by electron paramagnetic resonance (EPR) measurements. Geometry optimizations were carried out for low-spin (doublet) and high-spin (quartet) states. For the Co(acac)(2) complex two possible conformations, a square-planar and a tetrahedral one, were taken into account. For all structures relative energies were calculated with both "pure" and hybrid functionals. The calculated data were complemented with the results of the EPR investigations carried out at liquid helium temperature, allowing for definite assignment of the high-spin state for the Co(acac)(2)(EtOH)(2) complex. However, because of the unresolved spectral features, only effective g-values could be assessed, whereas the zero-field splitting parameters (ZFS) were calculated by means of the spin-orbit mean field (SOMF) relativistic DFT method for which direct spin-spin (SS) and spin-orbit coupling (SOC) contributions were quantified.     

The catalytic Mn2+ sites in the enolase-inhibitor complex: crystallography, single-crystal EPR, and DFT calculations

Crystals of Zn2+/Mn2+ yeast enolase with the inhibitor PhAH (phosphonoacetohydroxamate) were grown under conditions with a slight preference for binding of Zn2+ at the higher affinity site, site I. The structure of the Zn2+/Mn2+-PhAH complex was solved at a resolution of 1.54 A, and the two catalytic metal binding sites, I and II, show only subtle displacement compared to that of the corresponding complex with the native Mg2+ ions. Low-temperature echo-detected high-field (W-band, 95 GHz) EPR (electron paramagnetic resonance) and 1H ENDOR (electron-nuclear double resonance) were carried out on a single crystal, and rotation patterns were acquired in two perpendicular planes. Analysis of the rotation patterns resolved a total of six Mn2+ sites, four symmetry-related sites of one type and two out of the four of the other type. The observation of two chemically inequivalent Mn2+ sites shows that Mn2+ ions populate both sites I and II and the zero-field splitting (ZFS) tensors of the Mn2+ in the two sites were determined. The Mn2+ site with the larger D value was assigned to site I based on the 1H ENDOR spectra, which identified the relevant water ligands. This assignment is consistent with the seemingly larger deviation of site I from octahedral symmetry, compared to that of site II. The ENDOR results gave the coordinates of the protons of two water ligands, and adding them to the crystal structure revealed their involvement in a network of H bonds stabilizing the binding of the metal ions and PhAH. Although specific hyperfine interactions with the inhibitor were not determined, the spectroscopic properties of the Mn2+ in the two sites were consistent with the crystal structure. Density functional theory (DFT) calculations carried out on a cluster representing the catalytic site, with Mn2+ in site I and Zn2+ in site II, and vice versa, gave overestimated D values on the order of the experimental ones, although the larger D value was found for Mn2+ in site II rather than in site I. This discrepancy was attributed to the high sensitivity of the ZFS parameters to the Mn-O bond lengths and orientations, such that small, but significant, differences between the optimized and crystal structures alter the ZFS considerably, well above the difference between the two sites.