are the solvent effects critical in the modeling of polyoxoanions?

8
Are the Solvent Effects Critical in the Modeling of Polyoxoanions? XAVIER LO ´ PEZ, 1 JORGE A. FERNA ´ NDEZ, 1 SUSANNA ROMO, 1 JEAN FRANC ¸ OIS PAUL, 2 LEONID KAZANSKY, 3 JOSEP M. POBLET 1 1 Departament de Quı ´mica Fı ´sica i Inorga `nica, Universitat Rovira i Virgili, Imperial Ta `rraco 1, 43005 Tarragona, Spain 2 Laboratoire de Catalyse, Universite ´ de Lille1, 59655 Villeneuve d’Ascq, France 3 Institute of Physical Chemistry, Russian Academy of Sciences, 31 Leniski Prospekt, 119991 Moscow, Russia Received 7 January 2004; Accepted 5 May 2004 DOI 10.1002/jcc.20083 Published online in Wiley InterScience (www.interscience.wiley.com). Abstract: DFT calculations were driven for a set of differently charged polyoxoanions in the gas phase and in solution. We have calculated and analyzed their geometries and orbital energies to trace simple rules of behavior regarding the modeling of anions in isolated form. We discuss the quality of the results depending on the molecular charge, q, and the size of the cluster in terms of the number of metal centers, m. When the q/m ratio reaches a value of 0.8, DFT calculations for the isolated anion fail to describe the gap between the band of occupied oxo orbitals and the set of unoccupied orbitals delocalized among the metal atoms. In these cases the incorporation of the stabilizing external fields generated by the solvent through continuum models improves the geometries and orbital energies. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1542–1549, 2004 Key words: density-functional theory calculation; polyoxoanion; gas-phase calculation; solvent; COSMO; highly charged anions Introduction Polyoxometalates (POMs), or polyoxoanions, are inorganic spe- cies that often display high charges and only exist in condensed phases, where the external field stabilizes the anions. 1 However, due to technical limitations, most theoretical calculations driven on POMs assume molecules in the gas phase. 2–5 Since the first high- level theoretical studies, 6 greater accuracy in the description of POM’s properties, either by applying more exact functions or by introducing additional effects to the target systems, has been the main goal. The environment has been modeled in different ways to account for the external field generated by the solvent molecules and/or by counterions. In the last few years, some theoretical studies of several anions belonging to this family of compounds assumed that the external field generated on the anion is almost isotropic. 7 Consequently, the presence of the solvent 8 or the coun- terions 7 hardly modifies the properties derived from the sequence of molecular orbitals. It is commonly accepted that many POM properties, such as the relative acid-base properties of the external oxygen sites, the magnetic properties, the encapsulation, etc., can be well described by considering the isolated anion 2 —without including the secondary structure in the calculations. However, we have already noted that to study electron-transfer reactions be- tween differently charged polyoxoanions one needs to incorporate the stabilizing external field into the calculations. 8 This note is motivated from a routine calculation of the isolated Lindqvist hexaniobate, Nb 6 O 19 8 , carried out by our group using the BP86 functional. In these conditions, two results were unexpected: the large discrepancy between computed and experimental NbOO bond lengths (0.15 Å for NbOO t ) and the small gap obtained between HOMO and LUMO of 0.54 eV. The latter point is particularly important because niobates generally have similar electronic spectra (or are slightly shifted toward higher energies) than tungstates. 9 Taking into account that the HOMO–LUMO (H-L) gap for the homologous tungstate W 6 O 19 2 is larger than 3 eV, the electronic structure computed for the highly charged Nb 6 O 19 8 anion is totally unrealistic. This failure to reproduce the Correspondence to: J.M. Poblet; e-mail: [email protected] Contract/grant sponsor: the Spanish MCyT; contract grant number: BQU2002-04110-C02-01 Contract/grant sponsor: the CIRIT of the Generalitat de Catalunya (Autonomous Government of Catalonia); contract/grant number: SGR01- 00315 © 2004 Wiley Periodicals, Inc.

Upload: xavier-lopez

Post on 11-Jun-2016

212 views

Category:

Documents


0 download

TRANSCRIPT

Are the Solvent Effects Critical in the Modeling ofPolyoxoanions?

XAVIER LOPEZ,1 JORGE A. FERNANDEZ,1 SUSANNA ROMO,1 JEAN FRANCOIS PAUL,2

LEONID KAZANSKY,3 JOSEP M. POBLET1

1Departament de Quımica Fısica i Inorganica, Universitat Rovira i Virgili, Imperial Tarraco 1,43005 Tarragona, Spain

2Laboratoire de Catalyse, Universite de Lille1, 59655 Villeneuve d’Ascq, France3Institute of Physical Chemistry, Russian Academy of Sciences, 31 Leniski Prospekt,

119991 Moscow, Russia

Received 7 January 2004; Accepted 5 May 2004DOI 10.1002/jcc.20083

Published online in Wiley InterScience (www.interscience.wiley.com).

Abstract: DFT calculations were driven for a set of differently charged polyoxoanions in the gas phase and insolution. We have calculated and analyzed their geometries and orbital energies to trace simple rules of behaviorregarding the modeling of anions in isolated form. We discuss the quality of the results depending on the molecularcharge, q, and the size of the cluster in terms of the number of metal centers, m. When the q/m ratio reaches a valueof �0.8, DFT calculations for the isolated anion fail to describe the gap between the band of occupied oxo orbitals andthe set of unoccupied orbitals delocalized among the metal atoms. In these cases the incorporation of the stabilizingexternal fields generated by the solvent through continuum models improves the geometries and orbital energies.

© 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1542–1549, 2004

Key words: density-functional theory calculation; polyoxoanion; gas-phase calculation; solvent; COSMO; highlycharged anions

Introduction

Polyoxometalates (POMs), or polyoxoanions, are inorganic spe-cies that often display high charges and only exist in condensedphases, where the external field stabilizes the anions.1 However,due to technical limitations, most theoretical calculations driven onPOMs assume molecules in the gas phase.2–5 Since the first high-level theoretical studies,6 greater accuracy in the description ofPOM’s properties, either by applying more exact functions or byintroducing additional effects to the target systems, has been themain goal. The environment has been modeled in different ways toaccount for the external field generated by the solvent moleculesand/or by counterions. In the last few years, some theoreticalstudies of several anions belonging to this family of compoundsassumed that the external field generated on the anion is almostisotropic.7 Consequently, the presence of the solvent8 or the coun-terions7 hardly modifies the properties derived from the sequenceof molecular orbitals. It is commonly accepted that many POMproperties, such as the relative acid-base properties of the externaloxygen sites, the magnetic properties, the encapsulation, etc., canbe well described by considering the isolated anion2—withoutincluding the secondary structure in the calculations. However, we

have already noted that to study electron-transfer reactions be-tween differently charged polyoxoanions one needs to incorporatethe stabilizing external field into the calculations.8

This note is motivated from a routine calculation of the isolatedLindqvist hexaniobate, Nb6O19

8�, carried out by our group using theBP86 functional. In these conditions, two results were unexpected:the large discrepancy between computed and experimental NbOObond lengths (�0.15 Å for NbOOt) and the small gap obtainedbetween HOMO and LUMO of 0.54 eV. The latter point isparticularly important because niobates generally have similarelectronic spectra (or are slightly shifted toward higher energies)than tungstates.9 Taking into account that the HOMO–LUMO(H-L) gap for the homologous tungstate W6O19

2� is larger than 3eV, the electronic structure computed for the highly chargedNb6O19

8� anion is totally unrealistic. This failure to reproduce the

Correspondence to: J.M. Poblet; e-mail: [email protected]

Contract/grant sponsor: the Spanish MCyT; contract grant number:BQU2002-04110-C02-01

Contract/grant sponsor: the CIRIT of the Generalitat de Catalunya(Autonomous Government of Catalonia); contract/grant number: SGR01-00315

© 2004 Wiley Periodicals, Inc.

NbOO bond distances in this hexametalate, in comparison tometal–oxygen bond lengths in other POMs, was already observedby Bridgeman and Cavigliaso using a local functional.10

In this article, we analyze how the solvent affects the geometricand electronic structure for a set of differently charged well-knownPOMs to characterize in which conditions a gas-phase calculationproperly describes a given POM.

Results and Discussion

The importance of the negative charge of isolated clusters on theirelectronic properties may be crucial in some theoretical studies.We performed a series of calculations on a family of fully oxidizedPOMs to analyze a set of molecules with a variety of size andcharge: (1) the relatively small isopolyanions [W6�xNbxO19]q�

(x � 0–6 and q � 2–8) and W7O246�, some of which have high

negative charges; (2) the medium-sized Keggin clusters[XW12O40]q� (X � P, Si and q � 3, 4) and [SiM3W9O40]q� (M �Nb, Ti and q � 7, 10) and the large Preyssler anion,NaP5W30O110

14�, a cluster with a high negative charge but also witha relatively large number of metal atoms. The data discussed arisesfrom full geometry optimizations carried out for all the molecules.We used the ADF package to apply the DFT methodology.11 Thecalculations are based on the local density approximation charac-terized by the X� model12 for exchange with Becke’s gradient-corrected functional,13 and the VWN parameterization14 for cor-relation, corrected with Perdew’s functional.15 We used TZP-quality basis sets and the ZORA approximation16 for therelativistic effects of the most inert (core) electrons. The calcula-

tions were made in the vacuum and including the solvent effects bymeans of the conductor-like screening model (COSMO) (� �78).17,18 Additional calculations based on the B3LYP functional19

with and without the Polarizable Continuum Model (PCM)20,21 forthe solvent were also performed. For computing this series we usedthe Gaussian 98 package.22 The POM clusters presented here con-tain metal atoms free of valence electrons (d0 electronic configura-tion), so the molecular orbitals sequence consists of two well-sepa-rated sets of energy levels: the occupied oxo orbitals and the virtualmetal orbitals.23 Parameters related to the geometry and electronicstructure of the molecules are listed in Table 1. In this series, theresults for isolated anions are compared with those for anions insolution, where the solvent is modeled by a dielectric material simu-lating water. The central magnitude we refer to in this discussion isq/m. This represents the ratio between the total negative charge of theanion and the number of metal centers of the framework. The firstsection of the results concerns the relatively small Lindqvist structureand the heptatungstate anion, W7O24

6�. In the second section wediscuss the results for Keggin and Preyssler heteropolyanions.

Isopolyanions

For the set of W6�xNbx Lindqvist anions, it is known that formalsubstitution of W by Nb diminishes the electron affinity of theframework. This is a result of the larger H-L gaps in POMstructures containing Nb, which decrease the oxidizing powercompared to oxides of WVI. On the other hand, Ti-derivatives haveslightly smaller H-L gaps than analogous tungstates, but fairlysimilar ones to isostructural tungstates.24 The MOOt bond lengthscomputed either in the gas phase or in solution were systematically

Table 1. q/m Ratio, Calculated Bond Lengths (in Å) and HOMO-LUMO Gaps (in eV) for the Series of[W6�xNbxO19]q� Lindqvist Anions ( x � 0–6), Nb10O28

6�, W7O246�, XW12O40

q� (X � P, Si), SiM3W9O40q�

(M � Ti, Nb) and the Preyssler Anion NaP5W30O11014�, Computed in the Gas Phase and with Solvent.

q/ma sym

Gas phase Water

d(WOOt) d(MOOt)b

H-Lgap d(WOOt) d(MOOt)

b H-L gap

IsopolyanionsW6O19

2� 0.33 Oh 1.73 — 3.28 1.73 — 3.31W5NbO19

3� 0.50 C4v 1.75 1.77 3.31 1.74 1.77 3.39W4Nb2O19

4� 0.67 C2v 1.77 1.79 3.32 1.75 1.78 3.41W3Nb3O19

5� 0.83 C2v 1.78–1.79 1.82 3.23 1.76 1.79 3.42W2Nb4O19

6� 1.00 C2v 1.80 1.85 2.21 1.77 1.81 3.46WNb5O19

7� 1.17 C4v 1.82 1.87 1.27 1.78 1.82 3.49Nb6O19

8� 1.33 Oh — 1.90 0.54 — 1.83 3.62Nb10O28

6� 0.60 D2h — 1.78–1.80 3.42 — 1.77–1.78 3.65W7O24

6� 0.86 C2v 1.77–1.80 — 3.52 1.76–1.77 — 3.76Heteropolyanions

PW12O403� 0.25 Td 1.73 — 2.72 1.72 — 2.73

SiW12O404� 0.33 Td 1.74 — 2.73 1.73 — 2.76

NaP5W30O11014� 0.47 D5h 1.78–1.79 — 1.94 1.77–1.78 — 1.94

SiNb3W9O407� 0.58 C3v 1.77 1.80 2.84 1.75 1.79 2.95

SiTi3W9O4010� 0.83 C3v 1.80–1.81 1.76 2.43 1.76 1.74 2.60

aq/m indicates the (molecular charge)/(metal atoms) ratio.bM � Nb or Ti.

Are the Solvent Effects Critical in the Modeling of Polyoxoanions? 1543

longer at the DFT/BP86 level than the experimental values. Theoverestimation of the MAO distances obtained from calculationsperformed with the DFT methodology has been shown in theliterature for most polyoxoanions.2,25 On average, experimentalWOOt distances were about 1.69 Å and NbOOt distances were1.75–1.80 Å in Lindqvist anions, whereas in Keggin-type deriva-tives, the homologous bond lengths were 1.70–1.75 and 1.73 Å,respectively. TiOOt distances ranged from 1.65 to 1.70 Å inmixed-addenda Keggin clusters.26 These are structural parametersobtained from crystalline species. Even though it has been shownelsewhere that almost complete coincidence occurs between IRand Raman spectra of crystalline and dissolved POMs.27 Theinfluence of the counterions upon the geometries of such POMs isalso expected to be small.28

Table 1 shows that, within the gas phase approximation, thequality of the geometry described with DFT decreased steadily inparallel with the overall charge of the molecule. Figure 1 shows thecorrection introduced to the metal-terminal oxygen distances afteraccounting for the solvent effects. Two sets of values are plotted:&verbarl� d(MOO)� for M � W and M � Ti or Nb. For thehexatungstate anion, W6O19

2�, with q/m � 0.33 and a charge of �2,the description of the geometry was fairly accurate in the vacuum,with a deviation from the experimental value of 0.04 Å for d(WOOt).Incorporating the solvent effects via the continuum dielectric did notimprove the optimized distance. As the negative charge of the systemincreased, all the gas-phase MOOt istances increased fairly linearly.In the series of Lindqvist anions, the WOOt distance was 1.73 Å forW6. This converts into 1.82 Å in the WNb5 structure, which is verypoor. The trend was the same for d(NbOOt), but the deviationincreased faster because we obtained d(NbOOt) � 1.78 Å for q/m �0.5 and d(NbOOt) � 1.90 Å with q/m 1.33. The effects of the solvent

shortened the MAO bonds by 0.1–0.7 Å, depending on the value ofq/m. The other anions in Table 1 behaved similarly. We chosestructures with q/m parameters ranging from 0.25 to 0.86. Althoughthese values are not so large, there was also some improvement in thegeometries.

Analysis of the electronic structure in the series of Lindqvistanions produces some very interesting results. Figure 2 illustrates thegas-phase energies of the unoccupied (n�1)d and ns levels respectiveto the HOMO. For W6, the first unoccupied levels were, as expected,the d(W) orbitals, 3.28 eV above the HOMO, whereas the s(W) setwas located 7.4 eV from the oxo band. The W5Nb element of theseries behaved in a similar way, because the d orbitals (now a mixtureof W and Nb) were again 3.3 eV from the HOMO and the s(Nb)orbitals, although below s(W), were still 5.5 eV higher than the oxoband. Figure 2 shows the drop in the relative energy of the ns orbitalsas the q/m value increased. For larger values of q/m, this drop wasclear and in WNb5 and Nb6, the s(Nb) levels were 1 and 2 eV belowthe d(M) band, respectively. The energy of the empty d(M) orbitalsremained fairly constant in relation to the HOMO along the seriescomputed in the gas phase without solvent (see Fig. 1). Consideringthat series W6�xNbx is isoelectronic, the only change was the numberof protons in the nuclei. Therefore, the more external (diffuse) orbitalswere less destabilized when one proton was removed from the M6O19

framework. The more internal the orbital, the more sensitive to thenuclear charge. This explains the origin of the d/s crossing. Thisphenomenon is schematically shown in Figure 3. Introducing thesolvent effects modified this behavior because each set of orbitals wasaffected differently by the solvent. The occupied orbitals with a highparticipation of oxo ligands were the most stabilized because theywere the closest to the positively charged surface cavity.29 As Figure4 shows, the gap between d(W) and oxo bands grew as the q/m value

Figure 1. Absolute corrections introduced to the gas-phase bond lengths after accounting for the solvent,�� d(MOO)�, in Å, vs. the q/m parameter for the series of Lindqvist anions, [W6-xNbxO19]q�, listed inTable 1. Two sets of values for W (filled) and M � Ti or Nb (empty) are represented.

1544 Lopez et al. • Vol. 25, No. 12 • Journal of Computational Chemistry

increased. The largest difference occurred for Nb6 with a relativedestabilization of the d(Nb) orbitals by 0.40 eV. The energy of thes(M) orbitals was the least sensitive to the presence of the externalfield generated by the solvent (Fig. 3) and, in relative terms, theenergy of the d(M) orbitals varied more than that of the s(M) bandafter solvation. For example, the d and s(Nb) bands in Nb6, a moleculewith a negative charge of �8, appeared mixed �3.6 from the oxoband because d(M) had destabilized more than s(Nb) in relation to theHOMO. However, the d(Nb) levels were not completely below the sband, though the blue species contained a dxy-like metal electron. Thisis again because of the anion’s large negative charge. For the leastcharged W6 anion, the corresponding energies in solution of the d ands metal orbitals were 3.31 and 7.20 eV. The H-L gap agrees with theexperimental lowest charge transfer band.

The last two elements of this series, Nb10O266� and W7O24

6�, hadq/m ratios of 0.60 and 0.86, respectively. The decaniobate anion, withan H-L gap of 3.42 eV in the vacuum, was fairly well describedwithout the effect of the solvent. However, this orbital separationincreased considerably when the COSMO was introduced. This gapin solution was 3.65 eV, which was fairly similar to that of the Nb6

Lindqvist anion. Finally, for the heptatungstate, the H-L gap was 0.24eV smaller in the isolated form than in solution. In this case, theLUMO mainly comprised d(W) orbitals. The first s(W) levels werefound at higher energies, 4.2 eV above the HOMO. In solution thisvalue was larger, following the trend found for the Lindqvist series.Moreover, the relative energy of the s(W) orbitals in the heptatung-

state was similar to the 4.51 eV for this set in W3Nb3 [both had similarq/m ratios (�0.85)]. The packing in the W7O24 cluster was lesscompact than the other isopolyanions shown so far (Fig. 5), so thegeometry was somewhat more affected by the solvent effects (seeTable 2). X-ray data for W7O24

6� was taken from ref. 30 and forLindqvist anions from ref. 31. All the WOO distances improved afterthe solvent effects were accounted for. There was a decrease of 0.04Å in d(W2OO2), which was almost identical to the X-ray measure-ment. There was a generalized shortening (improvement) of theinteratomic distances caused by the external field. In the table, forcomparison purposes, we also list the homologous MOO distancesfor the hexatungstate and -niobate. As we discussed above, the com-putation of W6 in the gas phase produced very good results, whereasNb6, which carried a highly negative charge, was greatly affected bythe solvent, and therefore, with a clear contraction of the framework,approached the distances obtained experimentally. However, it wasthe metal-bridging oxygen distances that remained nearly unchangedin the majority of cases exposed. The terminal and central oxo sites inNb6 dramatically approached the metals (by up to 0.08 Å for NbOOc)to which they were coordinated. Comparatively, the hexatungstatewas almost unchanged by the effects of the solvent.

Heteropolyanions

A similar study was performed for several medium-sized and largeheteropolyanions. For the simplest molecules of this set, the XW12

Figure 2. Schematic view of the structure of unoccupied molecular orbitals in the gas phase. The energiesare relative to the HOMO. The first valence d-orbital, as well as the s(W) and s(Nb) levels of eachmolecule, are represented. The lowest level in each case corresponds to the LUMO.

Are the Solvent Effects Critical in the Modeling of Polyoxoanions? 1545

(X � P, Si) clusters (with charges �3 and �4), respectively, theexperimental distances were quite well reproduced in isolatedform. This may be attributed to a small q/m ratio (Table 1). Thecorrection introduced to the geometries with the COSMO was�0.01 Å for WOOt bonds. Both Keggin clusters roughly dis-played the same H-L gap in the gas phase and in the solution. Forlarger values of q/m, say �0.5, as in the Nb- and Ti-derivatives ofthe Keggin anion, the energy of the LUMO in the solution wassomewhat different from the energy in the gas phase. The relativeenergy of this d(W) orbital increased by 0.11 and 0.17 eV in SiNb3

and SiTi3, respectively. Despite the high charge of �10 carried bythe SiTi3W9 cluster, the framework contained 12 metal centers(q/m � 0.83) and the gas-phase description was still accurateenough. Neither the XW12 nor SiM3W9 anions suffered fromns/(n�1)d crossing described for the Lindqvist series.

A great computational effort was done to obtain the gas phase- andsolution-optimized structures of the Preyssler anion, NaP5W30O110

14�.This case is especially interesting because it enables us to checkanions with a high negative charge but a relatively small q/m value of0.47. The structure was optimized under the constraints of the D5h

symmetry group. The X-ray geometry was well reproduced32 with thehighest deviation in the MOOt distances by 0.05–0.06 Å in isolatedform.33 In agreement with this relatively small q/m value, the calcu-lations in aqueous solution hardly changed the geometry determinedfor the isolated anion. The energy of the d(W) band was also unalteredafter the effects of the solvent were included in the calculations. Theformation of big (or giant) POMs is generally accompanied by anincrease in the negative charge, but the values of q/m remain quitesmall. Take, for example, Muller’s wheel,34 whose formula is[Mo154(NO)14O420(OH)28(H2O)70](25�5)�, where q/m is less than

Figure 3. Molecular orbitals scheme for two differently charged polyanions studied in the gas phase andin solution. The gray boxes and the empty boxes denote sets of occupied and unoccupied orbitals,respectively. The gas-phase ns/(n�1)d orbital crossing is marked with dashed lines and discussed in thetext. Notice that this is a simplified scheme and that there are also unoccupied d orbitals above the s block.

1546 Lopez et al. • Vol. 25, No. 12 • Journal of Computational Chemistry

Figure 4. Energy of the first d(M)-like molecular orbital vs. the q/m parameter computed (1) in the gasphase (open diamonds) and (2) in solution (solid diamonds) for the series of W6�xNbx anions studied.Notice that for clarity the energy scale in the y-axis is much shorter than that in Figure 2.

Figure 5. Ball-and-stick views for (A) the M6O19q� and (B) W7O24

6� anions, with filled spheres represent-ing metal atoms. The point symmetry of the frameworks is Oh and C2v, respectively, and only thenecessary positions are labeled.

Are the Solvent Effects Critical in the Modeling of Polyoxoanions? 1547

0.20. This result is especially important because the COSMO treat-ment of the solvent did not involve a greater computational effort forsmall molecules. However, for a relatively large system, such as thePreyssler anion, the calculations in solution were �5 times moredemanding than for the isolated anion.

The results discussed up to this point were obtained by meansof the BP86 functional and the COSMO method. Let us put to thetest the conclusions drawn so far by changing the functional andthe solvation model. The series of Lindqvist derivatives wererecomputed, making use of the LDA and the B3LYP functionals.In the latter case, the PCM solvation method20,21 was used forcomparison to the COSMO model. Table 3 contains selectedinteratomic distances and H-L gaps, computed using LDA andB3LYP functionals for the isolated anions and for the anions in awater solution.35 The reader may find that there are great similar-ities between the sets of data contained in Tables 1 and 3. Inaddition to the well-known behavior of the LDA distances, ingeneral shorter than the BP86 or B3LYP ones, the only noticeabledifference is that the H-L gaps are slightly larger with the B3LYPfunctional. This behavior was already reported by Bridgeman,36

and does not modify the conclusions emerged from the systematicstudy using the BP86 functional.

To sum up, independently of the functional, we find that thegas-phase results show that higher charges of the anion lead to (1)longer MOO distances, and (2) smaller H-L gaps. In absoluteterms, this energy gap is characteristic of the functionals them-selves, but the tendency is to be smaller as the charge increases inall cases. The incorporation of the solvent in the calculations

Table 2. Selected Interatomic Distances (in Å) for W7O246�, W6O19

2�

and Nb6O198�, Computed in the Gas Phase and in Solution,

Together with X-ray Parameters.

Molecule

Computed

Gas phase Solution Experimental

HeptatungstateW1OO2 1.79 1.78 1.77W1OO5 1.92 1.92 1.89W1OO6 2.29 2.26 2.26W2OO2 2.51 2.47 2.48W2OO6 2.26 2.23 2.20W3OO3 2.01 1.98 1.98W3OO5 2.30 2.28 2.28W2OO1 1.80 1.77 1.74

W6O192�

WOOc 2.35 2.36 2.33WOOb 1.95 1.95 1.92WOOt 1.73 1.73 1.69

Nb6O198�

NbOOc 2.48 2.40 2.38NbOOb 2.04 2.03 2.01NbOOt 1.90 1.83 1.77

The labels in the left-hand column correspond to those shown in Figure 5.

Table 3. q/m Ratio, Calculated Bond Lengths (in Å) and HOMO-LUMO Gaps (in eV) for the Seriesof [W6�xNbxO19]q� Lindqvist Anions ( x � 0–6), Computed with the LDA and the B3LYPMethods in the Gas Phase and with Solvent.

q/ma sym

Gas phase Water

d(WOOt) d(NbOOt) H-L gap d(WOOt) d(NbOOt) H-L gap

LDA

W6O192� 0.33 Oh 1.71 — 3.44 1.71 — 3.47

W5NbO193� 0.50 C4v 1.73 1.76 3.48 1.72 1.75 3.52

W4Nb2O194� 0.67 C2v 1.75 1.78 3.48 1.73 1.76 3.56

W3Nb3O195� 0.83 C2v 1.77 1.80 3.14 1.74 1.78 3.56

W2Nb4O196� 1.00 C2v 1.79 1.83 2.13 1.75 1.79 3.60

WNb5O197� 1.17 C4v 1.81 1.85 1.18 1.76 1.80 3.63

Nb6O198� 1.33 Oh — 1.88 0.45 — 1.81 3.79

B3LYP

W6O192� 0.33 Oh 1.73 — 4.63 4.64

W5NbO193� 0.50 C4v 1.75 1.77 4.68 4.69

W4Nb2O194� 0.67 C2v 1.77 1.79 4.71 4.71

W3Nb3O195� 0.83 C2v 1.78 1.82 3.61 4.77

W2Nb4O196� 1.00 C2v 1.80 1.85 2.50 4.75

WNb5O197� 1.17 C4v 1.82 1.87 1.64 4.93

Nb6O198� 1.33 Oh — 1.90 1.03 4.84

LDA calculations were performed with the ADF set of programs and the solvent was modeled by the COSMOprocedure. Geometries were optimized in vacuum and in solution. For the B3LYP calculations the GAUSSIAN 98package was used, and single-point water solution calculations were carried out with the PCM solvation method. Thegas-phase geometries were not reoptimized in solution for the B3LYP/PCM calculations because the convergenceprocedure was extremely slow in GAUSSIAN 98 runs.

1548 Lopez et al. • Vol. 25, No. 12 • Journal of Computational Chemistry

always leads to an enlargement of the H-L gap. This broadening isproportional to the charge of the anion.

Conclusions

In this article, we have compared the electronic properties ofanions in isolated form and in solution. Bearing in mind thatpolyoxometalates (or polyoxoanions) only exist in condensedphases, their theoretical modeling must reproduce their behavior insolution. The gas-phase approximation may be satisfactory tostudy some chemical properties for low values of q/m; q � chargeand m � number of metal centers. However, it is strictly necessaryto incorporate the solvent effects in the calculations when q/m �0.8. For example, the HOMO–LUMO gap improves simply byaccounting for the external electrostatic effects generated by thesolvent molecules. These theoretical results also show that thestructures optimized in solution are always better than those for theisolated molecule. The iso- and heteropolyanions studied herefollow the same trend irrespective of their shape. In addition,different computational methods and solvation models have beentested, always leading to the same trend.

Acknowledgments

The authors thank the reviewers’ comments.

References

1. Hill, C. L. Chem Rev 1998, 98, 1.2. Borshch, S. A. Inorg Chem 1998, 37, 3116.3. Duclusaud, H.; Borshch, S. A. J Am Chem Soc 2001, 123, 2825.4. Bridgeman, A. J.; Cavligliasso, G. Faraday Discuss 2003, 124, 239.5. Poblet, J. M.; Lopez, X.; Bo, C. Chem Soc Rev 2003, 32, 297, and

references therein.6. Taketa, H.; Katsuki, S.; Eguchi, K.; Seiyama, T.; Yamazoe, N. J Phys

Chem 1986, 90, 2969; Kempf, J. Y.; Rohmer, M.-M.; Poblet, J. M.;Bo, C.; Benard, M. J Am Chem Soc 1992, 114, 1136.

7. (a) Rohmer, M.-M.; Benard, M.; Blaudeau, J.-P.; Maestre, J. M.;Poblet, J. M. Coord Chem Rev 1998, 179–180, 1019; (b) Maestre,M. M.; Poblet, M. M.; Bo, C.; Casan–Pastor, N.; Gomez–Romero, P.Inorg Chem 1998, 37, 3444.

8. Lopez, X.; Bo, C.; Poblet, J. M. J Am Chem Soc 2002, 124, 12574.9. So, H.; Pope, M. T. Inorg Chem 1972, 11, 1441.

10. Bridgeman, A. J.; Cavligliasso, G. Inorg Chem 2002, 41, 1761.11. ADF 2003.01; Department of Theoretical Chemistry, Vrije Univer-

siteit: Amsterdam; Baerends, E. J.; Ellis, D. E.; Ros, P. Chem Phys1973, 2, 41; Versluis, L.; Ziegler, T. J. Chem Phys 1988, 88, 322; TeVelde, G.; Baerends, E. J J Comput Phys 1992, 99, 84; FonsecaGuerra, C.; Snijders, J. G.; Te Velde, G.; Baerends, E. J Theor ChemAcc 1998, 99, 391.

12. Slater, J. C. Quantum Theory of Molecules and Solids; McGraw–Hill:New York, 1974, vol. 4.

13. Becke, A. D. J Chem Phys 1986, 84, 4524; Phys Rev 1988, A34, 3098.14. Vosko, S. H.; Wilk L.; Nusair, M. Can J Phys 1980, 58, 1200.15. Perdew, J. P. Phys Rev 1986, B33, 8822; 1986, B34, 7406.16. The reader is warned that the use of the ZORA in the present calcu-

lations may change some DFT results published before by our groups.

References: van Lenthe, E.; Ehlers, A. E.; Baerends, E. J. J Chem Phys1999, 110, 8943; van Lenthe, E.; Baerends, E. J.; Snijders, J. G.J Chem Phys 1993, 99, 4597; van Lenthe, E.; Baerends, E. J.; Snijders,J. G. J Chem Phys 1994, 101, 9783; van Lenthe, E.; Snijders, J. G.J Chem Phys 1996, 105, 6505.

17. Klamt, A.; Schuurmann, G. J Chem Soc, Perkin Trans 2, 1993, 799;Andzelm, J.; Kolmel, C.; Klamt, A. J Chem Phys 1995 103, 9312;Klamt, A. H Chem Phys 1995, 99, 2224.

18. Pye, C. C.; Ziegler, T. Theor Chem Acc 1999, 101, 396.19. Lee, C.; Yang, W.; Parr, R. G. Phys Rev B 1988, 37, 785; Miehlich,

B.; Savin, A.; Stoll, T.; Preuss, H. Chem Phys Lett 1989, 157, 200;Becke, A. D. Phys Rev A 1988, 38, 3098.

20. Miertus, S.; Scrocco, E.; Tomasi, J. Chem Phys 1981, 55, 117.21. Miertus, S.; Tomasi, J. Chem Phys 1982, 65, 239.22. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb,

M. A.; Cheeseman, J. R.; Zakrewski, J. A.; Montgenery, J. A.; Strat-man, R. E.; Burant, J. C.; Dapprich, S.; Milam, J. M.; Daniels, A. D.;Kudin, K. N.; Strain, M. C.; Farkas, O.; Tomasi, J.; Barone, V.; Cossi,M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Cifford, S.;Octerski, J.; Petersson, G. A.; Ayala, P. Y.; Cui, Q.; Morokuma, K.;Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.;Cioslowski, J.; Ortiz, J. V.; Stefano, B. B.; Liu, G.; Liashenko, A.;Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.; Fox, D. J.;Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Gonzalez,C.; Challacombe, M.; Gill, P. M. W.; Johnson, B. G.; Chen, W.; Wong,M. W.; Andres, J. L.; Head–Gordon, M.; Replogle, E. S.; Pople, J. A.Gaussian 98; Gaussian, Inc.: Pittsburgh, PA, 1998.

23. Maestre, J. M.; Lopez, X.; Bo, C.; Poblet, J. M.; Casan–Pastor, N.J Am Chem Soc 2001, 123, 3749.

24. Yamase, T.; Ozeki, T.; Motomura, S. Bull Chem Soc Jpn 1992, 65,1453.

25. Lopez, X.; Maestre, J. M.; Bo, C.; Poblet, J. M. J Am Chem Soc 2001,123, 9571; Lopez, X.; Bo, C.; Poblet, J. M. Sarasa, J. P. Inorg Chem2003, 42, 2634.

26. The structural information on polyoxometalates containing Ti is scarceand unclear. X-ray diffraction data normally show disorder, and TiAOdistances are, in general, not accurately measured.

27. Rocchiccioli–Deltcheff, C.; Thouvenot, R.; Rouassier, M. Inorg Chem1982, 21, 30; Rocchiccioli–Deltcheff, C.; Fournier, M.; Raymonde, F.;Thouvenot, R. Inorg Chem 1983, 22, 207.

28. Ichida, H.; A. Kobayashi, A.; Sasaki, Y. Acta Crystallogr B 1980, 36,1382.

29. In the COSMO method, the solvent is modeled by a set of pointcharges distributed over a surface cavity surrounding the solute mol-ecule.

30. Ikenoue, K.; Mikuriya, M.; Miyauchi, O.; Nukada, R.; Yagasaki, B.Bell Chem Soc Jpn 1994, 67, 2590.

31. Tytko, K. H.; Mehmke, J.; Fischer S. Struct Bonding 1999, 93, 129.32. The electronic structure of [XP5W30O110]q� will be discussed in detail

elsewhere.33. Alizadeh, M. H.; Harmalker, S. P.; Jennin, Y.; Martin–Frere, J.; Pope,

J. T. J Am Chem Soc 1985, 107, 2662.34. Muller, A.; Krickemeyer, E.; Meyer, J.; Bogge, H.; Peters, F.; Plass,

W.; Diemann, E.; Dillinger, S.; Nonnenbruch, F.; Randerath, M.;Menke, C. Angew Chem Int Ed Engl 1995, 34, 2122.

35. Calculations characterized by the LDA functional (both in the gasphase and with the COSMO method for the solvent) were carried outwith the ADF program. The B3LYP calculations were carried out withthe Gaussian 98 package. In this case, the PCM solvent method wasapplied in single-point calculations, so no structural parameters areincluded for this case in the table.

36. Bridgeman, A. J.; Cavigliasso, G. J Phys Chem A 2002, 106, 6114.

Are the Solvent Effects Critical in the Modeling of Polyoxoanions? 1549