PubMed Health⌕ Search

Biomedical subjects

Frank Neese

Publications and source records attributed to Frank Neese.

At least 19 recordsLinked to original sources
PubMed2006

Multireference ab initio calculations on reaction intermediates of the multicopper oxidases.

The multicopper oxidases (MCOs) couple the four-electron reduction of dioxygen to water with four one-electron oxidations of various substrates. Extensive spectroscopic studies have identified several intermediates in the MCO catalytic cycle, but they have not been able to settle the structures of three of the intermediates, viz. the native intermediate (NI), the peroxy intermediate (PI), and the peroxy adduct (PA). The suggested structures have been further refined and characterized by quantum mechanical/molecular mechanical (QM/MM) calculations. In this paper, we try to establish a direct link between theory and experiment, by calculating spectroscopic parameters for these intermediates using multireference wave functions from the multistate CASPT2 and MRDDCI2 methods. Thereby, we have been able to reproduce low-spin ground states (S = 0 or S = 1/2) for all the MCO intermediates, as well as a low-lying (approximately 150 cm-1) doublet state and a doublet-quartet energy gap of approximately 780 cm-1 for the NI. Moreover, we reproduce the zero-field splitting (approximately 70 cm-1) of the ground 2E state in a D3 symmetric hydroxy-bridged trinuclear Cu(II) model of the NI and obtain a quantitatively correct quartet-doublet splitting (164 cm-1) for a mu3-oxo-bridged trinuclear Cu(II) cluster. All results support the suggestion that the NI has an O2- atom in the center of the trinuclear cluster, whereas both the PI and PA have an O22- ion in the center of the cluster, in agreement with the QM/MM results and spectroscopic measurements.

Copper↗
PubMed2006

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.

Algorithms↗
PubMed2006

Electronic structure of bis(imino)pyridine iron dichloride, monochloride, and neutral ligand complexes: a combined structural, spectroscopic, and computational study.

The electronic structure of a family of bis(imino)pyridine iron dihalide, monohalide, and neutral ligand compounds has been investigated by spectroscopic and computational methods. The metrical parameters combined with Mössbauer spectroscopic and magnetic data for ((i)PrPDI)FeCl(2) ((i)PrPDI = 2,6-(2,6-(i)Pr(2)C(6)H(3)N=CMe)(2)C(5)H(3)N) established a high-spin ferrous center ligated by a neutral bis(imino)pyridine ligand. Comparing these data to those for the single electron reduction product, ((i)PrPDI)FeCl, again demonstrated a high-spin ferrous ion, but in this case the S(Fe) = 2 metal center is antiferromagnetically coupled to a ligand-centered radical (S(L) = (1)/(2)), accounting for the experimentally observed S = (3)/(2) ground state. Continued reduction to ((i)PrPDI)FeL(n) (L = N(2), n = 1,2; CO, n = 2; 4-(N,N-dimethylamino)pyridine, n = 1) resulted in a doubly reduced bis(imino)pyridine diradical, preserving the ferrous ion. Both the computational and the experimental data for the N,N-(dimethylamino)pyridine compound demonstrate nearly isoenergetic singlet (S(L) = 0) and triplet (S(L) = 1) forms of the bis(imino)pyridine dianion. In both spin states, the iron is intermediate spin (S(Fe) = 1) ferrous. Experimentally, the compound has a spin singlet ground state (S = 0) due to antiferromagnetic coupling of iron and the ligand triplet state. Mixing of the singlet diradical excited state with the triplet ground state of the ligand via spin-orbit coupling results in temperature-independent paramagnetism and accounts for the large dispersion in (1)H NMR chemical shifts observed for the in-plane protons on the chelate. Overall, these studies establish that reduction of ((i)PrPDI)FeCl(2) with alkali metal or borohydride reagents results in sequential electron transfers to the conjugated pi-system of the ligand rather than to the metal center.

Algorithms↗
PubMed2006

Octahedral non-heme oxo and non-oxo Fe(IV) complexes: an experimental/theoretical comparison.

Electron-transfer series are described for three ferric complexes of the pentadentate ligand 4,8,11-trimethyl-1,4,8,11-tetraazacyclotetradecane-1-acetate (Me(3)cyclam-acetate) with axial chloride, fluoride, and azide ligands. These complexes can all be reduced coulometrically to their Fe(II) analogs and oxidized reversibly to the corresponding Fe(IV) species. The Fe(II), Fe(III), and Fe(IV) species have been studied spectroscopically and their UV-vis, Mössbauer, EPR, and IR spectra are presented. The fluoro species [(Me(3)cyclam-acetate)FeF](n+) (n = 0, 1, 2) have been studied computationally using density functional theory (DFT), and the electronic structure of the Fe(IV) dication [(Me(3)cyclam-acetate)FeF](2+) is compared with that of the isoelectronic Fe(IV) oxo cation [(Me(3)cyclam-acetate)FeO](+); the different properties of the two species are mainly due to the significantly covalent Fe=O pi bonds in the latter.

Acetates↗
PubMed2006

Electronic structure of mononuclear bis(1,2-diaryl-1,2-ethylenedithiolato)iron complexes containing a fifth cyanide or phosphite ligand: a combined experimental and computational study.

A series of mononuclear square-based pyramidal complexes of iron containing two 1,2-diaryl-ethylene-1,2-dithiolate ligands in various oxidation levels has been synthesized. The reaction of the dinuclear species [Fe(III)2(1L*)2(1L)2]0, where (1L)2- is the closed shell di-(4-tert-butylphenyl)-1,2-ethylenedithiolate dianion and (1L*)1- is its one-electron-oxidized pi-radical monoanion, with [N(n-Bu)4]CN in toluene yields dark green crystals of mononuclear [N(n-Bu)4][Fe(II)(1L*)2(CN)] (1). The oxidation of 1 with ferrocenium hexafluorophosphate yields blue [Fe(III)(1L*)2(CN)] (1ox), and analogously, a reduction with [Cp2Co] yields [Cp2Co][N(n-Bu)4][Fe(II)(1L*)(1L)(CN)] (1red); oxidation of the neutral dimer with iodine gives [Fe(III)(1L*)2I] (2). The dimer reacts with the phosphite P(OCH3)3 to yield [Fe(II)(1L*)2{P(OCH3)3}] (3), and [Fe(III)2(3L*)2(3L)2] reacts with P(OC6H5)3 to give [Fe(II)(3L*)2{P(OC6H5)3}] (4), where (3L)2- represents 1,2-diphenyl-1,2-ethylenedithiolate(2-). Both 3 and 4 were electrochemically one-electron oxidized to the monocations 3ox and 4ox and reduced to the monoanions 3red and 4red. The structures of 1 and 4 have been determined by X-ray crystallography. All compounds have been studied by magnetic susceptibility measurements, X-band EPR, UV-vis, IR, and Mössbauer spectroscopies. The following five-coordinate chromophores have been identified: (a) [Fe(III)(L*)2X]n, X = CN-, I- (n = 0) (1ox, 2); X = P(OR)3 (n = 1+) )3ox, 4ox) with St = 1/2, SFe = 3/2; (b) [Fe(II)(L*)2X]n, X = CN-, (n = 1-) (1); X = P(OR)3 (n = 0) (3, 4) with St = SFe = 0; (c) [Fe(II)(L*)(L)X]n <--> [Fe(II)(L)(L*)X]n, X = CN- (n = 2-) (1red); X = P(OR)3 (n = 1-) (3red, 4red) with St = 1/2, SFe = 0 (or 1). Complex 1ox displays spin crossover behavior: St = 1/2 <--> St = 3/2 with intrinsic spin-state change SFe = 3/2 <--> SFe = 5/2. The electronic structures of 1 and 1(ox) have been established by density functional theoretical calculations: [Fe(II)(1L*)2(CN)]1- (SFe = 0, St = 0) and [Fe(III)(1L*)2(CN)]0 (SFe = 3/2, St = 1/2).

Computer Simulation↗
PubMed2006

Importance of direct spin-spin coupling and spin-flip excitations for the zero-field splittings of transition metal complexes: a case study.

This work reports the evaluation of several theoretical approaches to the zero-field splitting (ZFS) in transition metal complexes. The experimentally well-known complex [Mn(acac)3] is taken as an example. The direct spin-spin contributions to the ZFS have been calculated on the basis of density functional theory (DFT) or complete active space self-consistent field (CASSCF) wave functions and have been found to be much more important than previously assumed. The contributions of the direct term may exceed approximately 1 cm(-1) in magnitude and therefore cannot be neglected in any treatment that aims at a realistic quantitative modeling of the ZFS. In the DFT framework, two different variants to treat the spin-orbit coupling (SOC) term have been evaluated. The first approach is based on previous work by Pederson, Khanna, and Kortus, and the second is based on a "quasi-restricted" DFT treatment which is rooted in our previous work on ZFS. Both approaches provide very similar results and underestimate the SOC contribution to the ZFS by a factor of 2 or more. The SOC is represented by an accurate multicenter spin-orbit mean-field (SOMF) approximation which is compared to the popular effective DFT potential-derived SOC operator. In addition to the DFT results, direct "infinite order" ab initio calculations of the SOC contribution to the ZFS based on CASSCF wave functions, the spectroscopy-oriented configuration interaction (SORCI), and the difference-dedicated CI (DDCI) approach are reported. In general, the multireference ab initio results provide a more realistic description of the ZFS in [Mn(acac)3]. The conclusions likely carry over to many other systems. This is attributed to the explicit treatment of the multiplet effects which are of dominant importance, since the calculations demonstrate that, even in the high-spin d4 system MnIII, the spin-flip excitations make the largest contribution to the SOC. It is demonstrated that the ab initio methods can be used even for somewhat larger molecules (the present calculations were done with more than 500 basis functions) in a reasonable time frame. Much more economical but still fairly reasonable results have been achieved with the INDO/S treatment based on CASSCF and SOC-CI wave functions.

Journal Article↗
PubMed2006

Molecular and electronic structures of tetrahedral complexes of nickel and cobalt containing N,N'-disubstituted, bulky o-diiminobenzosemiquinonate(1-) pi-radical ligands.

The reaction of 2 equiv of the bulky ligand N,N'-bis(3,5-di-tert-butylphenyl)-1,2-phenylenediamine, H2[3L(PDI)], excess triethylamine, and 1 equiv of M(CH3CO2)2.4H2O (M = Ni, Co) in the presence of air in CH3CN/CH2Cl2 solution yields violet-black crystals of [Ni(II)(3L(ISQ))2] CH3CN (1) or violet crystals of [Co(3L)2] (3). By using Pd(CH3CO2)2 as starting material, green-blue crystals of [Pd(II)(3L(ISQ))2].CH3CN (2) were obtained. Single-crystal X-ray crystallography revealed that 1 and 3 contain (pseudo)tetrahedral neutral molecules [M(3L)2] (M = Ni, Co) whereas in 2 nearly square planar, neutral molecules [Pd(II)(3L(ISQ))2] are present. Temperature-dependent susceptibility measurements established that 1 and 2 are diamagnetic (S = 0) whereas 3 is paramagnetic with an S = 3/2 ground state. It is shown that 1 contains two pi radical benzosemiquinonate(1-)-type monoanions, ((3L(ISQ))(1-*), S(rad) = 1/2), and a central Ni(II) ion (d8; S = 1) which are antiferromagnetically coupled yielding the observed S(t) = 0 ground state. This result has been confirmed by broken symmetry DFT calculations of 1. In contrast, the S(t) = 3/2 ground state of 3 is more difficult to understand: the two resonance structures [Co(III)(3L(ISQ))(3L(PDI))] <--> [Co(II)(3L(PDI))(3L(IBQ))] might be invoked (for tetrahedral [Co(II)(3L(ISQ))2] containing an S(Co) = 3/2 with two antiferromagnetically coupled pi-radical ligands an S(t) = 1/2 is anticipated). Complex 2 is diamagnetic (S = 0) containing a Pd(II) ion (d8, S(Pd) = 0 in an almost square planar ligand field) and two antiferromagnetically coupled ligand radicals (S(rad) = 1/2). The electrochemistry and spectroelectrochemistry of 1, 2, and 3 have been studied, and electron-transfer series comprising the species [M(L)2]z (z = 2+, 1+, 0, 1-, 2-) have been established. All oxidations and reductions are ligand centered.

Benzoquinones↗
PubMed2006

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.

Journal Article↗
PubMed2006

An octahedral coordination complex of iron(VI).

The hexavalent state, considered to be the highest oxidation level accessible for iron, has previously been found only in the tetrahedral ferrate dianion, FeO4(2-). We report the photochemical synthesis of another Fe(VI) compound, an octahedrally coordinated dication bearing a terminal nitrido ligand. Mössbauer and x-ray absorption spectra, supported by density functional theory, are consistent with the octahedral structure having an FeN triple bond of 1.57 angstroms and a singlet d2(xy) ground electronic configuration. The compound is stable at 77 kelvin and yields a high-spin Fe(III) species upon warming.

Chemical Phenomena↗
PubMed2006

Molecular and electronic structures of oxo-bis(benzene-1,2-dithiolato)chromate(V) monoanions. A combined experimental and density functional study.

Two oxo-bis(benzene-1,2-dithiolato)chromate(V) complexes, namely, [CrO(L(Bu))2]1- and [CrO(L(Me))2]1-, have been synthesized and studied by UV-vis, EPR, magnetic circular dichroism (MCD), and X-ray absorption spectroscopy and by X-ray crystallography; their electro- and magnetochemistries are reported. H2L(Bu) represents the pro-ligand 3,5-di-tert-butylbenzene-1,2-dithiol, and H2L(Me) is the corresponding 4-methyl-benzene-1,2-dithiol. A structural feature of interest for both the complexes is the folding of the dithiolate ligands about the S-S vector providing Cs symmetry to the complexes. Geometry optimizations using all-electron density functional theory with scalar relativistic corrections at the second-order Douglas-Kroll-Hess (DKH2) and zeroth-order regular approximation (ZORA) levels result in excellent agreement with the experimentally determined structures and electronic and S K-edge X-ray absorption spectra. From DFT calculations, the Cs instead of C2v symmetry for the complexes is attributed to the strong S(3p) --> Cr(3d(x2-y2)) pi-donation in Cs geometry providing additional stability to the complexes.

Anions↗
PubMed2006

Vibrational markers for the open-shell character of transition metal bis-dithiolenes: an infrared, resonance raman, and quantum chemical study.

Transition metal complexes involving the benzene-1,2-dithiol (L(2-)) and Sellmann's 3,5-di-tert-butylbenzene-1,2-dithiol(L(Bu 2-)) ligands have been studied by UV-vis, infrared (IR), and resonance Raman (rR) spectroscopies. Raman spectra were obtained in resonance with the intervalence charge transfer (IVCT) bands in the near-infrared region and ligand-to-metal charge transfer (LMCT) bands in the near-UV region. Geometry optimization and frequency calculations using density functional theory (DFT) have been performed for [M(L)(2)](z) and [M(L(Bu))(2)](z) species (M = Ni, Pd, Pt, Co, Cu, Au, z = -1; M = Au, z = 0). On the basis of frequency calculations and normal-mode analysis, we have assigned the most important totally symmetric vibrations as well as corresponding overtone and combination bands that appear in rR spectra of compounds [Ni(L)(2)](1-), [M(L(Bu))(2)](1-) (M = Ni, Pt, Co, Cu). Experimental values of dimensionless normal coordinate displacements in excited states have been determined by fitting of rR spectra together with the absorption band shape, based on the time-dependent theory of Heller. Time-dependent density functional theory (TD-DFT) and multireference post-Hartree-Fock ab initio calculations, using the difference dedicated configuration interaction (MR-DDCI) method, were carried out to evaluate dimensionless normal coordinate displacements quantum chemically. The calculations show encouraging agreement with the experimental values. The large distortions along several normal modes led to significant vibronic broadening of IVCT and LMCT bands, and the broadening was accounted for in the deconvolution of the absorption spectra. The presence of an intense rR band around approximately 1100 cm(-1) was found to be a reliable marker for the presence of sulfur-based radicals.

Journal Article↗
PubMed2006

Effect of N-methylation of macrocyclic amine ligands on the spin state of iron(III): a tale of two fluoro complexes.

Syntheses and characterization of [(cyclamacetate)FeF]PF6 (1) and the corresponding N-methylated complex [(trimethylcyclamacetate)FeF]PF6 (3) are presented. Compound 1 is prepared in good yields from the analogous chloro complex, whereas 3 is prepared by hydrolysis of the oxo-bridged diiron compound (Me3cyclam-acetate)Fe-O-FeCl3 (2) in the presence of PF(6) anions. Magnetic susceptibility and spectroscopic data including electron paramagnetic resonance and Mössbauer spectra indicate that 1 contains low-spin Fe(III) (S = 1/2), while 3 is high spin (S = 5/2). Both octahedral fluoro complexes were investigated theoretically by density functional theory in order to determine why the spin states of the two molecules are different. Energies calculated using the B3LYP functional correctly predict 1 to have a low-spin S = 1/2 ground state and 3 to be high spin, regardless of whether a solvation model is included. The difference between 1 and 3 is most likely a combination of steric effects caused by the N-methyl groups, which compel the Fe-N bond distances to be longer in 3 than they ordinarily would be, and also electronic effects, which cause the N-methylated ligand to be a weaker sigma donor than its nonmethylated counterpart.

Amines↗
PubMed2006

Theoretical spectroscopy of model-nonheme [Fe(IV)OL5]2+ complexes in their lowest triplet and quintet states using multireference ab initio and density functional theory methods.

The structure, energies and spectroscopic properties of a simple [FeO(NH(3))(5)](2+) model with ground states (3)A(2g) and (5)A(1g) (in approximate C(4v) symmetry) have been studied in some detail using density functional (DFT) and simplified correlated multireference ab initio methods. The results reveal similarities as well as some pronounced differences in the properties of the molecule in the two alternative spin states.

Iron↗
PubMed2006

Calculation of solvent shifts on electronic g-tensors with the conductor-like screening model (COSMO) and its self-consistent generalization to real solvents (direct COSMO-RS).

The conductor-like screening model (COSMO) was used to investigate the solvent influence on electronic g-values of organic radicals. The previously studied diphenyl nitric oxide and di-tert-butyl nitric oxide radicals were taken as test cases. The calculations employed spin-unrestricted density functional theory and the BP and B3LYP density functionals. The g-tensors were calculated as mixed second derivative properties with respect to the external magnetic field and the electron magnetic moment. The first-order response of the Kohn-Sham orbitals with respect to the external magnetic field was determined through the coupled-perturbed DFT approach. The spin-orbit coupling operator was treated using an accurate multicenter spin-orbit mean-field (SOMF) approach. Provided that important hydrogen bonds are explicitly modeled by a supermolecule approach and that the basis set is sufficiently saturated, the COSMO calculations lead to accurate predictions of isotropic g-shifts with deviations of not more than 100 ppm relative to experiment. Very accurate results were obtained by employing a recently developed self-consistent modification of the COSMO method to real solvents (COSMO-RS), which we briefly introduce in this paper as direct COSMO-RS (D-COSMO-RS). This model gives isotropic g-shifts of similar high accuracy for water without using the supermolecule approach. This is an important result because it solves many of the problems associated with the supermolecule approach such as local minima and the choice of a suitable model system. Thus, the self-consistent D-COSMO-RS incorporates some specific solvation effects into continuum models, in particular it appears to successfully model the effects of hydrogen bonding. Although not yet widely validated, this opens a novel approach for the calculation of properties which so far only could be calculated by the inclusion of explicit solvent molecules in continuum solvation methods.

Algorithms↗
PubMed2006

Electronic structure of binuclear mixed valence copper azacryptates derived from integrated advanced EPR and DFT calculations.

Binuclear, mixed valence copper complexes with a [Cu(+1)(.5), Cu(+1)(.5)] redox state and S = (1)/(2) can be stabilized with rigid azacryptand ligands. In this system the unpaired electron is delocalized equally over the two copper ions, and it is one of the very few synthetic models for the electron mediating Cu(A) site of nitrous oxide reductase and cytochrome c oxidase. The spatial and electronic structures of the copper complex in frozen solution were obtained from the magnetic interactions, namely the g-tensor and the (63,65)Cu, (14)N, (2)H, and (1)H hyperfine couplings, in combination with density functional theory (DFT) calculations. The magnetic interactions were determined from continuous wave (CW) electron paramagnetic resonance (EPR), pulsed electron nuclear double resonance (ENDOR), two-dimensional TRIPLE, and hyperfine sublevel correlation spectroscopy (HYSCORE) carried out at W-band or/and X-band frequencies. The DFT calculated g and Cu hyperfine values were in good agreement with the experimental values showing that the structure in solution is indeed close to that of the optimized structure. Then, the DFT calculated hyperfine parameters were used as guidelines and starting points in the simulations of the various experimental ENDOR spectra. A satisfactory agreement with the experimental results was obtained for the (14)N hyperfine and quadrupole interactions. For (1)H the DFT calculations gave good predictions for the hyperfine tensor orientations and signs, and they were also successful in reproducing trends in the magnitude of the various proton hyperfine couplings. These, in turn, were very useful for ENDOR signals assignments and served as constraints on the simulation parameters.

Journal Article↗
PubMed2006

QM/MM calculations with DFT for taking into account protein effects on the EPR and optical spectra of metalloproteins. Plastocyanin as a case study.

A detailed study of the influence of the surrounding protein on magnetic and optical spectra of metalloproteins is presented using the quantum-mechanical/molecular mechanical (QM/MM) approach. The well-studied type I copper site in plastocyanin in the cupric oxidation state is taken as a test case because its spectroscopic properties have been extensively studied and are well understood. The calculations have been performed using nonrelativistic and scalar relativistic (at the level of the zeroth order regular approximation, ZORA) calculations (B3LYP functional). Linear response theory has been used to calculate first- and second-order properties, namely the EPR g-tensor, the central metal hyperfine couplings (HFCs), the HFCs of the directly coordinating ligands, as well as superhyperfine couplings (1H, 14N) from remote nuclei, transition energies, and oscillator strengths. Two different model systems have been defined that do not and do include important amino acids from the second coordination sphere, respectively. For comparison, calculations have been carried out in the gas phase and in a dielectric continuum (conductor like screening model, COSMO) with a dielectric constant of four. The best results were obtained at the scalar relativistic ZORA level for the largest model in conjunction with explicit modeling of the protein environment through the QM/MM procedure, which is also considered to be the highest level of theory used in this work. The protein effects beyond the second coordination sphere were found to be quite substantial (up to 30% changes on some properties), and were found to require an explicit treatment of the protein beyond the second coordination sphere. In addition, the embedding water cage was found to have a nonnegligible influence on the calculated spectroscopic data, which is of the same order as the influence of the protein backbone charges. However, while qualitatively satisfactory, the errors in the calculated spectroscopic parameters are still substantial, and can all be traced back to the fact that the linear-response of the presently available functionals is "too stiff" with respect to the external perturbations at least for the model systems studied here. Ligand field-based approaches are used to correct for systematic errors in the DFT procedures. As a consequence, we propose a new breakdown of the copper hyperfine interaction into Fermi-contact, spin-dipolar and spin-orbit contributions.

Binding Sites↗