Time-Dependent Density Functional Theory Calculations of Ligand K

Dipartimento di Scienze Chimiche, UniVersita` di Trieste, Via L. Giorgieri 1, ... Nazionale per la Scienza e Tecnologia dei Materiali, INSTM, Unita' d...
0 downloads 0 Views 160KB Size
J. Phys. Chem. A 2004, 108, 8467-8477

8467

Time-Dependent Density Functional Theory Calculations of Ligand K Edge and Metal L Edge X-ray Absorption of a Series of Oxomolybdenum Complexes G. Fronzoni,* M. Stener, A. Reduce, and P. Decleva Dipartimento di Scienze Chimiche, UniVersita` di Trieste, Via L. Giorgieri 1, I-34127 Trieste, Italy, and Consorzio InteruniVersitario Nazionale per la Scienza e Tecnologia dei Materiali, INSTM, Unita’ di Trieste, and INFM DEMOCRITOS National Simulation Center, Trieste, Italy ReceiVed: May 13, 2004; In Final Form: July 22, 2004

X-ray absorption spectra (XAS) at the Cl and S K edge and Mo L edge have been calculated at the TDDFT level for a series of dioxomolybdenum complexes LMoO2X (L ) hydrotris(3,5-dimethyl-1-pyrazolyl)borate, X ) Cl, SCH2Ph, OPh), which play an important role in modeling the catalytic cycle of the sulfite oxidase enzyme. Also, the XAS spectra of model molecules of the Mo complexes have been simulated and interpreted in terms of the Mo 4d orbital splitting, in order to find possible correlations with the spectral pattern of the complexes. Comparison with the available experimental data allows us to assess the performances of the present computational scheme to describe the core excitations in large bioinorganic systems. The theoretical interpretation of the spectral features of both the metal and ligand core excitations in terms of the oscillator strength distribution provides important insight into the covalency of the metal-ligand bond.

1. Introduction Molybdenum is the only second-row transition metal known to be essential for living systems. Its active biological form, present in almost all Mo-containing enzymes, is known as molybdenum cofactor (Moco), and it appears that the coordination by this ligand is essential for the catalytic function of the Mo enzymes. The Mo active center is common to important human enzymes including xanthine oxidase, sulfite oxidase, and aldehyde oxidase. In particular, the sulfite oxidase enzyme catalyzes the conversion of sulfite to sulfate, the terminal step in the oxidative degradation of cysteine and methionine. Structural information on the coordination of Mo in the sulfite oxidase, provided by EXAFS1 and crystal studies,2 has revealed that in oxidized sulfite oxidase the Mo is pentacoordinated by three sulfur ligands and two oxo groups with approximate square pyramidal coordination geometry. Two of the sulfur ligands arise from the thiolene group of the molybdopterin unit, and the third is contributed by the side chain of a cysteine molecule:

The conversion of sulfite to sulfate, catalyzed by sulfite oxidase, occurs with the transfer of an oxygen atom from the dioxo-Mo(VI) center of the enzyme to sulfite to generate an oxo-Mo(IV) center and sulfate. Important insight into the mechanism of the oxygen-atom transfer can be greatly facilitated by studies of the chemistry of inorganic models for the Mo cofactor, designed to mimic the MoO2 center and simulate the cofactor ligand donor atom set. Among the several complexes investigated,3 an especially significant model is based on a tris(pyrazolyl)borate ligand (L) which is able to catalyze the transfer

of an oxygen atom from a suitable donor to an acceptor. It has been verified that a model of the form LMoO2X, with X as the coligand, reacts in the oxidized form with PPh3 (triphenyl phosphine), losing an oxo group which is then regenerated from water, therefore mimicking the key features of the catalytic cycle proposed for the sulfite oxidase.4 A series of dioxo-Mo(VI) complexes of the type LMoO2X (L ) hydrotris(3,5-dimethyl1-pyrazolyl)borate, X ) Cl, Br, OPh, SPh, SCH2Ph, SCH(CH3)2) have been recently synthesized and characterized.5 Electronic structure insight for these complexes has been obtained from the analysis of the pre-edge X-ray absorption spectra (XAS), which can provide useful information on the metal-ligand center covalency and its effects on the oxo transfer activity.4 It is well-known that XAS spectroscopy is one of the premier tools to elucidate the local structure environment of metal ions, even in low concentrations. XAS involves the photoexcitation of core electrons to unoccupied valence orbitals and to the continuum. Its great advantage is associated with the strongly localized nature of the core excitations, which makes the K and L edge spectra very sensitive to the electronic structure and to the local geometrical environment of the absorbing atom even when it is embedded in a very complex situation. XAS has been also applied to biological molecules,6-8 providing important insight into the electronic characterization of the metalloproteins’ active sites, which often exhibit unique spectroscopic features that reflect highly covalent sites. In particular, the metal L edge (2p f 3d + continuum) and the ligand K edge (1s f 3p + continuum) are particularly useful because of both the improved resolution of the experiments at lower energies and the electricdipole-allowed nature of the transitions involved, which can give important information on the covalency of ligand-metal bonds.9-11 In fact, the intensity of the metal L edge features provides a direct probe of the d metal contribution to the halfoccupied or unoccupied valence antibonding MOs, whereas the intensity of the ligand K edge structures allows us to estimate the 3(n)p ligand character in these MOs due to covalent bonding.

10.1021/jp047953u CCC: $27.50 © 2004 American Chemical Society Published on Web 09/09/2004

8468 J. Phys. Chem. A, Vol. 108, No. 40, 2004

Fronzoni et al.

Metal-ligand bonding in metallobiomolecules is dominated by the ligand’s σ and π donor interactions, with neutral ligands behaving mostly as σ donors. The metal centers with d orbitals that are relatively high in energy can be available for π bonding interaction with the unoccupied low-lying valence orbitals on the neutral ligands. In some instances, ligand-metal donor interactions can be highly covalent, with the covalency greatly influencing the reactivity of the metal center. Major advances have also been made in the field of theoretical approaches for the calculations of the electronic structure of large systems. In particular, the density functional theory (DFT) allows high-level calculations on large transition-metal systems, as required for metalloprotein active sites. Here, it should be emphasized that the detailed information obtainable from accurate calculations can be of significant support in the understanding of the virtual electronic states of the molecular system and in the assignment of the related measured spectral features. The time-dependent extension of the DFT (TDDFT) can be applied to electronic excited states, and from its origin, it has received a noteworthy impulse to describe optical absorption in atoms and molecules.12,13 TDDFT represents the state of the art for evaluation of the molecular electronic transitions in the DFT framework. The current TDDFT implementations14-17 allow the calculations of the outer valence excitation spectra in terms of excitation energies and oscillator strengths. Recently, an extension of the method in the ADF code18 to the treatment of core excitations has been implemented, and test calculations of the XAS spectra of the TiCl4 molecule have been performed at the K and L metal edges and ligand L edge, which gave excellent agreement between theory and experiment. It is important to underline that, at the TDDFT level also, the coupling between different excitation channels, which becomes allowed in the case of degenerate core holes (such as 2p), is included and proves to be very important in describing the 2p absorption edge. The results obtained have pointed out the reliability of the method to properly describe and interpret core excitation spectra, and its computational economy makes it a good candidate to treat larger and less symmetric systems with reasonable effort. In this work, we present an application of the TDDFT method to the calculations of the ligand K edge and metal L edge of the Mo(VI) complexes LMoO2X (X ) Cl, SCH2Ph, OPh), which mimic the [MoVIO2]2+ active site of the sulfite oxidase enzyme. The comparison of the TDDFT results with the XAS experimental spectra measured by Izumi et al.4 allows us to assess the performance of the present approach for simulation of the XAS spectra of large bioinorganic systems. The TDDFT calculations are also performed on the core excitations of model molecules for the three LMoO2X complexes in order to analyze in detail the electronic structure and the metal-ligand bond in simpler systems. The possibility to extend the estimate of metal-ligand covalency from the model analogous to the complex system through the comparison of their XAS spectra is pointed out. 2. Theoretical Method The TDDFT approach for electron excitations and its implementation in the ADF code has been described in detail in the literature;14,19 here, we just recall the salient steps and describe the new features. The general problem is cast in the following eigenvalue equation:

ΩFi ) ωi2Fi

(1)

where Ω is a four-index matrix with elements Ωiaσ,jbτ, and the indexes consist of products of occupied virtual (ia and jb) KS orbitals, while σ and τ refer to the spin variable. The eigenvalues ω2i correspond to squared excitation energies, while the oscillator strengths are extracted from the eigenvectors Fi.19 The Ω-matrix elements can be expressed in terms of KS eigenvalues () and the coupling matrix K:

Ωiaσ,jbτ ) δστδijδab(a - i)2 + 2x(a - i)Kiaσ,jbτx(b - j) (2) The elements of the coupling matrix K are given by the following expression:

Kijσ,klτ )

∫dr ∫dr′φiσ(r)φjσ(r)[|r -1 r′| +

]

fστ xc (r, r′, ω) φkτ(r′)φlτ(r′) (3) where the φ’s are the KS orbitals and fστ xc (r, r′, ω) is the exchange correlation kernel. In this work, the kernel is approximated according to the adiabatic local density approximation (ALDA).20 The space spanned by the solutions of eigenvalue eq 1 corresponds to the 1h-1p excited configurations, so it is possible to approximate this space operating a selection over the configurations, keeping only those necessary for an accurate description of the phenomenon, as is customary in ab initio CI calculations. So, in practice, according to the method in ref 18, the indexes which span the occupied orbital space (i and j) are limited to run only over the core shell. 3. Computational Details The calculations have been performed with the ADF program,21,22 modified according to the previous section. The allelectron triple-ζ plus polarization Slater-type orbitals (TZP STO) basis set has been used for all the atoms except hydrogen, for which the double-ζ (DZ) set without the polarization function has been chosen. All basis sets have been taken from the ADF database. The same tests previously performed with other basis sets have proven that the present choice is a good compromise between accuracy and computational effort. The LB9423 exchange correlation potential with the groundstate (GS) electron configuration has been employed. It has been chosen because of its correct asymptotic behavior, which is a necessary condition for a good description of virtual orbitals. The exchange correlation kernel is approximated within the usual ALDA scheme, and only excitations with singlet spin symmetry have been considered. In the model molecules, the nitrogen atoms of ligand L (reported in Chart 1) were replaced by three ammonia molecules, and the thiolate and alkoxide ligands were represented by SCH3 and OCH3, respectively (see Chart 1). For the geometry of the model molecules, the same strategy used in the work of Izumi et al.4 has been adopted without any geometry optimization or symmetry idealization. For the LMoO2X complexes, a DFT geometry optimization has been performed, employing the VWN potential24 using a double-ζ plus polarization (DZP) basis set for B, C, N, and O, a DZ set for the Mo atom, and a single-ζ (SZ-STO) basis set for H. All core levels (up to 1s for B, C, N, and O; 2p for S and Cl; and 3d for Mo) have been treated as frozen orbitals. The calculated Mo-X distances are 2.31, 2.38, and 1.94 Å for X ) Cl, SCH2Ph, and OPh, respectively, which are found in

X-ray Absorption of Oxomolybdenum Complexes

J. Phys. Chem. A, Vol. 108, No. 40, 2004 8469

CHART 1: Coordinate Systems Used in the DFT Calculations for the Model Molecules [(NH3)3MoO2X]+ (X ) Cl, SCH3, OCH3) (Upper Part) and the Dioxomolybdenum Complexes LMoO2X (X ) Cl, SCH2Ph, OPh) (Lower Part)

Figure 1. Energy diagram of the Mo 4d orbitals for the DFT ground state. Upper panel: [(NH3)3MoO2Cl]+, [(NH3)3MoO2SCH3]+, and [(NH3)3MoO2OCH3]+. Lower panel: LMoO2Cl, LMoO2SCH2Ph, and LMoO2OPh. The energy separations between the orbital energies are given in eV.

the range of the crystal distances reported in the literature5 for dioxo-Mo complexes with slightly different tridentate N-donor ligands. The calculated Mo-N distances are 2.19, 2.17, and 2.20 Å for X ) Cl, SCH2Ph, and OPh, respectively. All of the calculated spectra at the ligand 1s edge and Mo 2p edge have been convoluted with Gaussian functions with fwhm’s of 1.3 and 1.8 eV, respectively, to facilitate the comparison with the experimental spectra. 4. Results and Discussion The experimental spectra of the LMoO2X complexes4 are analyzed in detail in the lowest-energy pre-edge region; therefore, we report in figures and tables the calculated excitation energies and oscillator strengths relative to the same energy range, both for the model molecules [(NH3)3MoO2X]+ and the LMoO2X complexes. This pre-edge region involves the electronic transitions toward the lowest virtual valence orbitals, which are mostly characterized by the 4d metal contribution with only small np ligand participation. The intensity of these transitions maps the two different final orbital contributions through the excitations from the L metal and K ligand edges. Therefore, for a more complete discussion, it has appeared convenient to collect in a single table the results obtained for both edges of each compound. 4.1. Minimum Model Molecules. In the model molecules, the nitrogen atoms of the L ligand are replaced by three ammonia molecules, and the benzylthiolate and phenoxo ligands are substituted with SCH3 and OCH3 moieties, respectively. The coordinate system of the [(NH3)3MoO2X]+ (X ) Cl, SCH3, OCH3) model used for the DFT calculations is presented in

Chart 1. To correlate the calculated spectral features to the electronic structures of the model molecules, it can be useful to consider the MO scheme of the lowest virtual valence levels, which are characterized by strong 4d metal contribution, as already stressed (Figure 1). The splitting of the five 4d MOs can be interpreted as derived from a Mo esa-coordinated complex with a highly distorted Oh geometry. For each molecule, we observe a group of three MOs at lower energy (dxy, dyz, and dxz), formally derived from the t2g set in Oh symmetry, while the two MOs derived from the eg set (dz2 and dx2-y2), are destabilized. The dz2 MO has the highest energy in all molecules because of its strong antibonding interaction with the oxo groups, while the dxy and dyz MOs, which are diagonal with respect to the ligand’s position, are the most stabilized. Ligand K Edge. The TDDFT Cl 1s and S 1s pre-edge excitation spectra of [(NH3)3MoO2Cl]+ and [(NH3)3MoO2SCH3]+ are reported in Figure 2 and in Tables 1 and 2. The spectral features are associated with transitions from the ligand (Cl or S) 1s orbital, which has an essentially atomic character, and therefore, the dipole-allowed transitions gain intensity only if the final low-lying unoccupied MOs, generally dominated by the 4d metal character, also contain some 3p ligand contribution. This is likely to be small because the 3p AOs (of Cl or S atoms) also contribute to occupied orbitals; therefore, the relatively low values calculated for the pre-edge 1s features are not surprising. The low-energy Cl 1s spectrum of [(NH3)3MoO2Cl]+ (Figure 2 and Table 1) exhibits two main structures separated by about 2.7 eV and followed by a small tail at higher energy. The first peak derives from the Cl 1s excitations into the first three lower MOs (53a, 54a, and 55a), which have a predominant 4d character (dxy, dyz, and dxz, respectively). The highest intensity of the first transition toward the LUMO maps the larger

8470 J. Phys. Chem. A, Vol. 108, No. 40, 2004

Figure 2. TDDFT Cl 1s and S 1s pre-edge excitation spectra of [(NH3)3MoO2Cl]+and [(NH3)3MoO2SCH3]+. Calculated lines are convoluted with Gaussians of 1.3 eV fwhm.

contribution of the Cl 3p component in the 53a final state with respect to the next two MOs (54a and 55a), whose oscillator strengths decrease proportionally to their Cl 3p component (see Table 1). The second structure is contributed by two transitions toward the 56a and 57a MOs: The first one (at 2781.07 eV) is associated with the Cl 3p component of the dx2-y2 level (56a), and its intensity is comparable to that calculated for the 53a transition, reflecting the similar Cl 3p participation in these final states. The small Cl 3p component in the 57a final state is responsible for the intensity drop of the relative line at 2781.32 eV. It is interesting to note that this final state has only a low 4d metal contribution, while the highest MO with strong 4d character (dz2) corresponds to the 60a final state. Here, the Cl 3p component is negligible and the intensity of the relative transition very low, as for the following lines which contribute to the small tail in the higher-energy part of the spectrum. This analysis allows us to correlate the spectral features to the covalent bonding between the metal and the Cl ligand; the dominant contribution to the metal-ligand covalency is provided by the Cl 3pπ* and Cl 3pσ* (π and σ refer to the MoCl bond) orbitals which significantly overlap with the Mo 4dxy and 4dx2-y2 orbitals in the 53a and 56a levels, as probed by the higher intensity of the relative Cl 1s transitions in the spectrum. The interaction of the second Cl 3pπ* orbital with the 4dyz orbital (in the 54a level) is instead weaker because of the tilt of the Mo-Cl bond in the direction opposite to the O-Mo-O diagonal, with a consequent reduction of the orbital overlap in accord with the decrease of the spectral intensity of the relative Cl 1s transition.

Fronzoni et al. Consider now the results for [(NH3)3MoO2(SCH3)]+ (Figure 2 and Table 2). Note that the shape of the S 1s spectrum appears different from the Cl 1s spectrum, with three bands of comparable intensity: The first one is assigned as transitions to the lowest virtual orbitals with mainly 4d character, the second is contributed by a single transition toward the 61a MO (4dx2-y2), while the third and most intense band derives from many lines that are close in energy and of significant intensity. As concerns the first structure, the main difference with the X ) Cl ligand is the shift of the intensity from the LUMO transition (57a) to the (LUMO + 1) transition (58a, at 2434.32 eV), which maps the higher S 3p participation in the 58a final orbital with respect to the 57a and 59a orbitals (see Table 2). This indicates that the thiolate π* orbital (S 3pz) overlaps principally with the 4dyz orbital, in contrast to X ) Cl, being the second S π* orbital (3px) involved in the S-C bonding and therefore not available for significant interaction with the lowest 4dxy metal orbital. The third S 3p orbital (3py) can interact primarily with the 4dx2-y2 orbital, to give rise to a σ* orbital, which corresponds to the final state (61a) of the S 1s transition at 2436.84 eV and is responsible for the second strong peak of the spectrum. This shows that the dominant contribution to the metal-ligand covalency comes from the interaction between the S 3pσ* and 4dx2-y2 orbitals, as found also in the case of the Cl ligand. The first two peaks are separated by about 2.5 eV, a value very close to that found for the X ) Cl compound, therefore suggesting a similar splitting between the first three metal d orbitals and the dx2-y2 orbital in the two model molecules, as confirmed also by the energetic scheme of Figure 1. The higher-energy transitions which give rise to the third spectral feature sample the S 3p content of the higher virtual MOs. The oscillator strength distribution also reflects a certain mixing among 1h-1p configurations, in particular in the 67a and 68a final states. The S 3p content represents the virtual counterpart of the S 3px orbital involved in the S-C bond, which is pushed to higherlying virtual MOs and is responsible for the different spectral response to the bonding situation with respect to the X ) Cl ligand. It has to be noted that the highest 4d metal orbital, dz2, contributes to the 64a virtual orbital, whose S 3p component is quite small, as is apparent from the low intensity calculated for the corresponding transition. The results obtained for the ligand K edge show that there is a proportionality between the calculated intensity of the ligand features and the mixing of ligand orbitals into the metal d orbitals, confirming the capability of the XAS ligand spectra to probe the covalency of the metal-ligand bond. Molybdenum L Edge. Consider now the 2p metal edge: The relative excitations are expected to be dominated by the 2p4d transition moment, which directly probes the molybdenum 4d content of the different unoccupied levels and, therefore, the details of the d involvement into the bonding. The results are reported in Figure 3 and in Tables 1-3. The complexity of the spectra, with respect to the ligand K edge, derives from the large number of lines which contribute to the calculated features, as is apparent in Figure 3. Because of the lack of symmetry elements in these molecules, for each final orbital, there are three distinct transitions starting from the three Mo 2p orbitals, whose SCF energies are slightly different; furthermore, the effect of the 1h-1p configuration mixing included in the TDDFT calculations can further split the relative excitation energies. The theoretical spectra of [(NH3)3MoO2Cl]+ and [(NH3)3MoO2(SCH3)]+ exhibit three main features: a lower-energy and mostintense peak followed by a middle, less-resolved structure, and then a last, higher-energy peak. The second and third peaks

X-ray Absorption of Oxomolybdenum Complexes

J. Phys. Chem. A, Vol. 108, No. 40, 2004 8471

TABLE 1: Calculated TDDFT Excitation Energies (E(eV)) and Oscillator Strengths (f) for Chlorine 1s and Molybdenum 2p Excitation Spectra of [(NH3)3MoO2Cl]+ a Cl 1s

Mo 2p

E (eV)

f × 102

final state

Mo 4d and Cl 3p character

E (eV)

f × 102

2778.33

0.126

53a

Mo 4dxy (57%) + Cl 3p (12%)

2778.58

0.059

54a

Mo 4dyz (56%) + Cl 3p (6%)

2779.35

0.024

55a

Mo 4dxz (55%) + Cl 3p (2%)

2781.07

0.120

56a

Mo 4dx2-y2(41%) + Cl 3p (10%)

2781.32

0.033

57a

Mo 4dx2-y2 (10%) + Cl 3p (1.5%)

2782.46

0.004

58a

H 1s + N 2p

2782.70

0.005

59a

H 1s + N 2p

2783.28

0.006

60a

Mo 4dz2 (47%)

2783.85

0.0003

61a

Mo 4dxy (4%)

2784.29

0.001

62a

Mo 4dx2-y2 (3%)

2499.93 2500.06 2500.15 2500.22 2500.44 2500.56 2500.09 2501.18 2501.44 2502.68 2502.78 2502.85 2502.95 2503.04 2503.23 2504.09 2504.10 2504.11 2504.32 2504.33 2504.34 2504.95 2505.25 2505.31 2505.48 2505.49 2505.50 2505.92 2505.93 2505.93

0.001 0.095 0.364 0.002 0.379 0.880 0.316 0.487 1.225 0.043 0.308 0.119 0.010 0.491 1.025 0.004 0.026 0.017 0.004 0.013 0.043 ∼0 0.921 0.820 0.018 0.044 0.078 0.012 0.010 0.030

a

configuration mixing in Mo 2p transitions 53a 53a + 54a 53a 54a 54a 54a + 53a 55a 55a 55a 56a 56a + 57a 56a + 57a 57a 57a + 56a 57a + 56a 58a 59a 60a 61a 62a

Calculated Cl 1s and Mo 2p DFT-KS eigenvalues (ionization thresholds) are 2792.26 and 2513.89 eV, respectively.

coalesce into only one peak in the [(NH3)3MoO2OCH3]+ spectrum, reflecting the different energy splitting of the transition lines which distribute into two distinct groups, which is apparent as well in Figure 3. The spectral behavior can be rationalized by considering the diagram of the virtual valence levels with mainly Mo 4d character, reported in Figure 1. In the Cl complex, the five 4d levels are split into three groups by the ligand field: the lower-lying dxy, dyz, and dxz MOs which are close in energy, then the dx2-y2 orbital which lies at approximately 2 eV from both the lowest MO group, and the highest dz2 orbital. The splitting of the 4d orbitals is consistent with the shape of the calculated spectrum in Figure 3: The first peak (at about 2500.4 eV) is relative to the transitions toward the lowest 4d orbitals (53a, 54a, and 55a). The mixing among the 1h-1p configurations in the 53a and 54a final states is responsible for a nonequivalent intensity distribution between the two states, as expected on the grounds of their similar 4d content (see Table 1). The middle peak (around 2503 eV) is associated with the Mo 2p transitions to the 4dx2-y2 component present in the 56a and 57a final states. Also in this case, the intensity distribution between these two final states does not reflect the 4d participation to the MOs (41% in 56a and 10% in 57a) but, rather, reflects the mixing of the 2p f 56a and 2p f 57a excitations in the final states as reported in Table 1. The last calculated peak is attributed to the excitations to the higher-lying 4dz2 final orbital, which characterizes the 60a final state at approximately 2505 eV. The three-peaked shape is maintained also in the thiolate spectrum according to the similar energy splitting of the Mo 4d orbitals, as shown in Figure 1. The first peak (around 2500 eV) is assigned to the transitions toward the lower-lying 4d orbitals (57a, 58a, and 59a), and its intensity is comparable to that calculated for the Cl complex. The middle peak is essentially associated with the 2p f 61a (4dx2-y2) transitions around 2502.7

eV; it is characterized by smaller intensity with respect to that obtained for the corresponding peak of the previous molecule and a reversed shape for the presence of a higher-energy tail associated with the small 4dx2-y2 component present in the 63a final state. The most significant difference with the Cl complex spectrum is relative to the third peak, which appears to be contributed by several transitions of comparable intensity in SCH3. Although we can identify the 64a final orbital as mainly 4d (dz2) in character (about 37%), we do not calculate significant oscillator strength for the 2p transition toward this MO; instead, the intensity is spread over other excited states with less significant 4d metal participation. This is an effect of the strong mixing of 1h-1p configurations (2p f 64a, 2p f 65a, and 2p f 66a) in the final states lying in the 2504-2505 eV energy range. The spectrum of the alkoxide complex (Figure 3 and Table 3) shows only two features of similar intensity. The first one manifests the same nature as in the two previous compounds, derived from the transitions to the lowest 4d MOs of the alkoxide. Also, the energy range covered by these transitions is similar to that of the other two compounds (the peak is centered at 2500.8 eV), as well as the intensity. The strongest difference is the disappearance of the middle peak in the spectrum: This is an effect of different splitting of the metal d orbitals, as is apparent in Figure 1. We observe an increase of the energy separation between the lower-lying 4d orbitals and the most destabilized dx2-y2 and dz2 orbitals (in the alkoxide, this splitting is approximately 2.6 eV, while it is approximately 1.8 and 2 eV in X ) Cl and X ) SCH3 complexes, respectively) and a parallel decrease of the splitting between the last two d levels. This is responsible for the energy positions of the relative transitions and the appearance of only one higher-energy peak in the spectrum. This peak is centered at approximately 2505 eV and derives from several transitions toward final states with

8472 J. Phys. Chem. A, Vol. 108, No. 40, 2004

Fronzoni et al.

TABLE 2: Calculated TDDFT Excitation Energies (E(eV)) and Oscillator Strengths (f) for Sulfur 1s and Molybdenum 2p Excitation Spectra of [(NH3)3MoO2(SCH3)]+ a S 1s

Mo 2p

E (eV)

f × 102

final state

Mo 4d and S 3p character

E (eV)

f × 102

2433.61

0.034

57a

Mo 4dxy (58%) + S 3p (4%)

2434.32

0.137

58a

Mo 4dyz (50%) + S 3p (16%)

2434.81

0.033

59a

Mo 4dxz (54%) + S 3p (4%)

2436.22

0.003

60a

H 1s + Mo 5s

2436.84

0.168

61a

Mo 4dx2-y2 (43%) + S 3p (20%)

2437.47

0.008

62a

H 1s + N 2p

2437.63

0.008

63a

Mo 4dx2-y2 (3%) + S 3p (2%)

2438.56

0.027

64a

Mo 4dz2 (37%) + S 3p (7%)

2438.80

0.057

65a

Mo 4dz2 (11%) + S 3p (20%)

2439.00

0.055

66a

Mo 4dz2 (5%) + S 3p (8%)

2439.30

0.024

67a

Mo 4dz2 (3%) + S 3p (6%)

2439.42

0.072

68a

Mo 4dz2 (2%) + S 3p (7%)

2499.44 2499.63 2499.65 2500.13 2500.29 2500.38 2500.71 2500.85 2501.07 2502.02 2502.03 2502.05 2502.64 2502.78 2502.92 2503.28 2503.29 2503.30 2503.44 2503.46 2503.50 2504.38 2504.49 2504.52 2504.60 2504.71 2504.74 2504.78 2504.84 2504.87 2505.09 2505.11 2505.16 2505.18 2505.24 2505.25

0.003 0.383 0.282 0.0004 0.261 0.607 0.369 0.619 1.093 0.004 0.003 0.009 0.028 0.583 0.668 0.016 0.028 0.007 0.032 0.062 0.290 0.010 0.080 0.151 0.059 0.318 0.225 0.005 0.276 0.368 0.022 0.056 0.149 0.010 0.224 0.372

a

configuration mixing in Mo 2p transitions 57a 58a 58a 58a + 59a 59a 60a 61a 62a 63a 64a 64a + 65a 64a + 65a 65a 65a + 66a 65a + 66a 66a 66a 66a + 65a 67a 67a 67a + 68a 68a 68a 68a + 67a

Calculated S 1s and Mo 2p DFT-KS eigenvalues (ionization thresholds) are 2446.83 and 2512.64 eV, respectively.

different dx2-y2 and dz2 metal contribution and different intensities. As we can see in Table 3, the 2p f 4dx2-y2 transition involves the 57a and 59a final states (at 2503.8 and 2504.2 eV) which are mixed by the 1h-1p configuration interaction with a consequent redistribution of the intensity, which concentrates on the final state at 2504.3 eV. A quite large oscillator strength is calculated for the 2p f 60a transition which reflects the significant dz2 contribution to the 60a final orbital; the lowerintensity lines at higher energy derive from transitions toward final states, among which the dz2 content redistributes with lower weight. 4.2. Dioxomolybdenum Complexes. The energy diagram of the virtual MOs containing strong Mo 4d contribution (Figure 1) can be used as a guide for the interpretation of the calculated XAS spectra of the LMoO2X (X ) Cl, SCH2Ph, OPh) complexes. The coordinate systems used for the DFT calculations are shown in Chart 1. As we can see, all three molecules show the same ordering of the Mo 4d MOs: a group of three low-lying orbitals, which have the energy order dxy < dyz < dxz and are well separated by the two most destabilized dx2-y2 and dz2 orbitals, with the energy of dz2 greater than that of dx2-y2. We can note that the calculated 4d energy splitting for the complexes is analogous to that found for the model molecules. Ligand K Edge. The TDDFT-calculated Cl and S K edge spectra of LMoO2(SCH2Ph) and LMoO2(OPh) are reported in Tables 4 and 5 and in Figure 4 together with the experimental profiles.4 The theoretical Cl 1s spectrum shows two main features centered at 2778.9 and 2781.6 eV and reproduces correctly the experimental pattern (Figure 4). The main discrepancy is the

strong underestimate of the calculated excitation energies with respect to the experimental results (see caption to Figure 4), a typical problem ascribed to the self-interaction present in DFT, which is particularly important for core excitations. However, for the discussion of the spectra, it is more important to consider the energy splitting between the features than the absolute scale; therefore, to facilitate the comparison, we have shifted the experimental spectrum on the calculated energy scale. We can note the nice reproduction of the energy separation between the two first peaks, which amounts to 2.7 eV in the calculation and 2.3 eV in the experiment, while the intensity of the second structure appears to be underestimated in the calculation, producing a reversed intensity trend with respect to the experiment. The first peak is assigned to the Cl 1s transitions toward the three lower-lying virtual orbitals (118a, 119a, and 120a in Table 4) with predominant Mo 4d character (4dxy, 4dyz, and 4dxz, respectively) and only a small contribution from the Cl ligand. This is highest for the LUMO (Cl 3pπ component) and decreases for the next two MOs, as correctly mapped by the calculated oscillator strengths. The second peak corresponds to the transitions toward the three following final states (121a, 122a, and 123a) characterized by different 4dx2-y2 contributions and by smaller Cl 3p components. The similar intensity calculated for these lines does not reflect the different Cl 3p contributions to the final states (Table 4) and depends on the slight mixing of the Cl 1s f 121a and Cl 1s f 122a configurations, which redistributes the oscillator strengths over the final states. It has to be noted that the highest MO with 4d metal character (4dz2), corresponding to the 129a final state, has only a negligible Cl 3p participation and does not contribute to

X-ray Absorption of Oxomolybdenum Complexes

J. Phys. Chem. A, Vol. 108, No. 40, 2004 8473 TABLE 3: Calculated TDDFT Excitation Energies (E(eV)) and Oscillator Strengths (f) for Molybdenum 2p Excitation Spectrum of [(NH3)3MoO2(OCH3)]+ a final E (eV) f × 102 state 2499.98 2500.21 2500.24 2500.67 2500.85 2500.95 2501.19 2501.34 2501.53 2502.57 2502.59 2502.61 2503.68 2503.75 2503.90 2503.82 2503.84 2503.85 2504.09 2504.11 2504.35 2504.90 2505.18 2505.21 2505.28 2505.30 2505.31 2505.47 2505.49 2505.50 2505.67 2505.69 2505.71

0.002 0.278 0.388 0.0005 0.176 0.591 0.500 0.869 1.132 0.006 0.009 0.010 0.057 0.254 0.243 0.023 0.028 0.029 0.021 0.204 0.999 0.010 0.730 0.587 0.141 0.076 0.166 0.051 0.046 0.023 0.005 0.270 0.313

Mo 4d and O 2p characterb

configuration mixing

53a Mo 4dxy (62%)+ O 2p (3%)

53a

54a Mo 4dyz (57%)+ O 2p (9%)

54a 54a + 55a 54a + 55a 55a + 54a 55a + 54a 55a 56a

55a Mo 4dxz (54%) + O 2p (1%) 56a H 1s + Mo s 57a Mo 4dx2-y2 (41%) + O 2p (7%) 58a H 1s + N 2p 59a Mo 4dx2-y2 (11%)

57a 57a 57a + 59a 58a

60a Mo 4dz2 (47%)

59a 59a 59a + 57a 60a

61a Mo 4dxy (5%)

61a

62a Mo 4dx2-y2 (2%)

62a

63a Mo 4dz2 (4%)

63a

a Calculated Mo 2p DFT-KS eigenvalue (ionization threshold) is 2513.3 eV. b The O 2p contribution to the final state is relative to the oxygen atom of the alkoxide ligand.

Figure 3. TDDFT Mo 2p pre-edge excitation spectra of [(NH3)3MoO2Cl]+, [(NH3)3MoO2SCH3]+, and [(NH3)3MoO2OCH3]+. Calculated lines are convoluted with Gaussians of 1.8 eV fwhm.

any spectral feature. The higher-energy tail in the calculated spectrum, not considered in the experiment, derives from transitions to virtual MOs with ligand contributions. Apart from fine details, the interpretation of the two pre-edge features is analogous to that proposed for the corresponding model molecule as well as the considerations made for the covalent bonding between the Mo and Cl ligands: The most significant overlap is realized between the Cl 3pπ* and Mo 4dxy and between Cl 3pσ* and Mo 4dx2-y2 orbitals. The theoretical results for the S 1s spectrum of LMoO2(SCH2Ph) are reported in Table 5 and Figure 4, where the experimental data are also shown (shifted on the calculated energy scale). The reported experimental features have been obtained through a background subtraction from the normalized spectrum;4 this procedure has ruled out the intensity contribution from structures lying at higher energies. For this reason, we have chosen to compare our theoretical results with the difference spectrum of ref 4 and to limit the discussion only to the two lower-energy features. As we can see, the agreement between theory and experiment is very satisfactory, both concerning the energy separation of the first two peaks (2.3 eV in the calculation and 2.0 eV in the experiment) and the intensity distribution. The first peak is contributed by three transitions to the lowest 4d

virtual orbitals (142a, 143a, and 144a) with the highest intensity associated with the LUMO + 1 orbital (4dyz), which has the most significant overlap with the S 3pπ* component (see Table 5). The second band centered at 2436.79 eV is assigned to the 149a transition, whose high intensity reflects the significant S 3p contribution (S 3pσ*) to the final orbital dominated by the 4dx2-y2 component. The last and most intense calculated structure is mostly contributed by the transition to the 155a final state, which is characterized by a strong S 3p participation and only a minor 4dz2 contribution. The 4dz2 component characterizes the 156a MO whose S 3p content is, however, negligible. It is interesting to notice that our calculated S 1s ionization threshold (2442.77 eV, see Table 5) lies more than 8 eV above the first calculated structure. This suggests that the background subtraction performed on the experimental data might be revised with a shift to higher energies, to ascribe a wider energy range to the pre-edge region. In summary, the calculated S 1s spectral behavior for the LMoO2(SCH2Ph) complex closely resembles that obtained for its model analogues, also concerning the main differences with the X ) Cl complex. This behavior can be related to the similarity of the calculated 4d orbital energy splitting for the models and the LMoO2X molecules and of the extent of mixing between 4d metal and S 3p ligand components in the lower-lying virtual MOs. Molybdenum L Edge. The TDDFT results for the Mo L edge spectra of LMoO2Cl, LMoO2(SCH2Ph), and LMoO2(OPh) are reported in Tables 4, 5, and 6 and in Figure 5 together with the experimental profiles relative to the L3 edge spectra. These are

8474 J. Phys. Chem. A, Vol. 108, No. 40, 2004

Fronzoni et al.

TABLE 4: Calculated TDDFT Excitation Energies (E(eV)) and Oscillator Strengths (f) for Chlorine 1s and Molybdenum 2p Excitation Spectra of [LMoO2Cl]+ a Cl 1s

Mo 2p

E (eV)

f × 102

E (eV)

f × 102

2778.74

0.1079

2779.04

0.0603

2779.76

0.0180

2781.46

0.0363

2781.53

0.0417

2781.65

0.0488

2782.17

0.0030

2782.76

0.0002

2783.31

0.002

2783.35

0.0013

2499.99 2500.13 2500.21 2500.31 2500.53 2500.65 2501.13 2501.23 2501.47 2502.72 2502.75 2502.77 2502.79 2502.81 2502.86 2502.91 2503.06 2503.29 2503.43 2503.44 2503.45 2504.03 2504.03 2504.04 2504.58 2504.59 2504.66 2504.98 2505.03

0.001 0.115 0.357 ∼0 0.348 0.823 0.309 0.507 1.157 0.032 0.042 0.009 0.056 0.070 0.054 0.001 0.652 1.173 0.0004 0.005 0.053 0.002 0.011 0.004 0.004 ∼0 ∼0 1.065 0.961

a

final state

Mo 4d and Cl 3p character

configurations mixing in Mo 2p transition

118a

Mo 4dxy (56%) + Cl 3p (10%)

119a

Mo 4dyz (55%) + Cl 3p (6%)

120a

Mo 4dxz (55%) + Cl 3p (