Bonding and Magnetic Response Properties of Several Toroid

Sep 1, 2011 - The Ni–ligand interaction has been studied by means projected density of states and energy decomposition analysis, which denotes its ...
0 downloads 0 Views 4MB Size
ARTICLE pubs.acs.org/JPCA

Bonding and Magnetic Response Properties of Several Toroid Structures. Insights of the Role of Ni2S2 as a Building Block from Relativistic Density Functional Theory Calculations Alvaro Mu~noz-Castro* Doctorado en Fisico-Quimica Molecular, Relativistic Molecular Physics (ReMoPh) group, Universidad Andres Bello, Republica 275, Santiago, Chile

bS Supporting Information ABSTRACT: Relativistic density functional calculations were carried out on several nickel toroid mercaptides of the general formula [Ni(μ-SR)2]n, with the aim to characterize and analyze their stability and magnetic response properties, in order to gain more insights into their stabilization and size-dependent behavior. The Niligand interaction has been studied by means projected density of states and energy decomposition analysis, which denotes its stabilizing character. The graphical representation of the response to an external magnetic field is applied for the very first time taking into account the spinorbit term. This map allows one to clearly characterize the magnetic behavior inside and in the closeness of the toroid structure showing the prescence of paratropic ring currents inside the Nin ring, and by contrast, diatropic currents confined in each Ni2S2 motif denoting an aromatic behavior (in terms of magnetic criteria). The calculated data suggests that the Ni2S2 moiety can be regarded as a stable constructing block, which can afford several toroid structures of different nuclearities in agreement with that reported in the experimental literature. In addition, the effects of the relativistic treatment over the magnetic response properties on these lighter compounds are denoted by comparing nonrelativistic, scalar relativistic, and scalar plus spinorbit relativistic treatments, showing their acting, although nonpronunced, role.

’ INTRODUCTION Crown ethers have received much attention since their discovery in 1967,1 prompted by their particular structures and properties giving rise to a rapid development of the hostguest chemistry. Their ability to selectively bind metal atoms or ions2 is of importance in several fields such as supramolecular chemistry, catalysis, material science,3 or as simple model systems for understanding noncovalent interactions.4 In this context, their inorganic counterparts have become an area of increasing interest driven by the versatility of a discrete number of redox centers into a molecular macrocyclic array with properties related to the organic crown ethers moieties, leading to the development of some inorganic cages5 and metallocrowns,6 toward the design of multimetallic catalysts7 and enzyme mimics,8 among other applications in materials science.9 Nickel mercaptides (Ni(SR)2) are a fundamental issue in the development of metallocrowns, since the pioneering work by Woodward10a in the middle of the 1960s. Initially considered as insoluble high polymers,11 their definitive structural characterization was achieved following efforts to include more soluble ligands, revealing their molecular toroid structure10 through the discrete hexanuclear [Ni(μ-SCH2CH3)2]6 complex,10a anticipating the characterization of related systems with different nuclearities and including other redox centers from the same group.1215 r 2011 American Chemical Society

The resulting robust structure allows one to include in a cyclic arrangement several redox centers in a molecular unit of four-,12 five-13 and six-10,14 nickel vertices, among others,15 where the NiNi intermetallic distance remains longer than the mean separation in nickel metal (2.49 Å),16 with values ranging from 2.65 to 3.05 Å,1014 suggesting the potential ability to selectively determinate a desired size of host domains by varying the number of vertices, which could give rise to materials with sizedependent properties leading to exceptional novelty and intriguing behavior.17 In this article we focus on the four-, five-, and six-membered nickel mercaptides, with the aim to characterize the energetic and electronic variations due to the increase in nuclearity, which in turn determine the diameter and bonding of the host, which ensure the stability of the toroid structure. In addition, we explore the magnetic behavior of these closed-shell systems by mapping the shielded and deshielded areas18 in an applied external field, via the nucleus-independent chemical shifts (NICS) procedure,1921 in order to clarify zones that exhibits diatropic (aromaticity) and paratropic (antiaromaticity) induced currents. Finally, we will compare nonrelativistic and relativistic calculations in order to Received: March 27, 2011 Revised: July 16, 2011 Published: September 01, 2011 10789

dx.doi.org/10.1021/jp2028438 | J. Phys. Chem. A 2011, 115, 10789–10794

The Journal of Physical Chemistry A

ARTICLE

characterize the dependence of the magnetic response to the scalar and spinorbit terms in the first-row transition metal systems as part of our current research.22,23

’ COMPUTATIONAL DETAILS Relativistic density functional theory calculations24 were done by using the ADF 2010.02 code,25 incorporating both scalar and spinorbit corrections via the two-component ZORA Hamiltonian.26,27 Triple-ξ Slater basis set plus polarization function (STO-TZP) were employed within the generalized gradient

Figure 1. Schematic representation of 1, 2, and 3.

’ RESULTS AND DISCUSSION

Table 1. Selected Distances (Å) and Angles (deg) of the Systems 1, 2, and 3 1 calc. (D2)

a

2 exp.a

calc. (Cs)

3 exp.a

calc. (D3)

exp.a

diam. Ni

3.883

4.876

5.848

5.848

diam. S

5.246

5.550

6.639

6.635

NiNi

2.737

2.924

2.928

NiNi BOb

0.22

NiS

2.253

NiS BOb

0.76

NiNiNi

90

89

110

108

120

120

NiSNi

75

74

80

80

83

84

SNiS

98

98

98

98

98

97

2.665

2.892

2.835

0.18 2.208

2.243

0.15 2.207

0.77

approximation (GGA) according to the PerdewBurkeErnzerhof (PBE)28a,b in order to account for the exchange-correlation (xc) effects due to its improved performance on long-range interactions and relative low computational cost for larger molecules.28c,d,e Geometry optimizations were done without any symmetry restrain, via the analytical energy gradient method implemented by Verluis and Ziegler.29 Bond order calculations were carried out according to Mayer as it is implemented in the ADF package, due to its valuable results in analyzing inorganic molecules and clusters.3032 The map of induced magnetic fields18 were calculated according to the NICS19,20 procedure of the Schleyer group, employing the GGA exchange expression proposed by Handy and Cohen28f and the correlation expression proposed by Perdew, Burke, and Ernzerhof28 (OPBE), incorporating the scalar relativistic (OPBE/ ZORA) and both scalar and spinorbit effects (OPBE/ZORA +SO) through the ZORA Hamiltonian and STO-TZP as the basis set. For comparison, nonrelativistic (NR) calculations were done at a similar level of theory.

2.205

2.210

0.77

holding anglec

110°

117°

122°

123°

twist angled

13°

10°





Averaged experimental values from refs 12a, 13b, and 14e for [Ni(μSC5H9N(CH3)2)2]4, [Ni(μ-SCH2SiCH3)2]5, and [Ni(μ-SCH2CH2OH)2]6, respectively. b Bond order calculations via Mayer analyses.3032 c Folding angle of the Ni2S2 moiety. d Twist angle of the NiS4 moiety.

Molecular Structures. Toroid structures are already obtained by using several bridging mercaptide ligands, such as NaSR and NaSeR,1015 and a nickel salt, resulting in cyclic arrays of different nuclearity. In order to obtain an overall picture of the four-, five-, and six- membered cyclic nickel mercaptides ([Ni(μSR)2]n), we consider methyl-mercaptide (R = CH3)-based compounds, namely, [Ni(μ-SCH3)2]4 (1), [Ni(μ-SCH3)2]5 (2), and [Ni(μ-SCH3)2]6 (3) (Figure 1). The optimized structures are close to the highest point symmetry possible, resulting in almost D2 structure for 1, Cs for 2, and D3 for 3. These structures exhibit only a few imaginary frequencies due to the free rotation of the methyl groups in the range of 120i cm1 to 140i cm1. Calculated geometrical parameters for these slightly simplified models are summarized in Table 1, denoting its reasonable agreement with the available experimental data for [Ni(μSC5H9N(CH3)2)2]4,12a [Ni(μ-SCH2SiCH3)2]5,13b and [Ni(μSCH2CH2OH)2]6.14e The closest NiNi distances vary from 2.737 Å for 1, to 2.924 Å for 3, leading to metallocrowns of diameter 3.883 Å, 4.876 Å, and 5.848 Å, depicted by the Nin ring, embeded between two Sn rings of 5.246 Å, 5.550 Å, and 6.639 Å diameter, respectively, for 1, 2 and 3. The NiS distance remain similar in the series. Moreover, the NiNiNi angle increases in relation to the circumference size, in contrast with the NiSNi angle of each Ni2S2 rhomboid varying to a lesser degree,15a suggesting that in the toroid structure, the Ni2S2 moiety is a

Figure 2. DOS, projected DOS for 3p of S and 3d of Ni atoms, and Mulliken orbital overlap population (MOOP) between both atomic shells of 1, 2, and 3. Red dashed line denotes the EF. 10790

dx.doi.org/10.1021/jp2028438 |J. Phys. Chem. A 2011, 115, 10789–10794

The Journal of Physical Chemistry A

ARTICLE

robust elementary building block to obtain toroid structures with different numbers of vertices, which vary slightly from 1 to 3. The nonplanarity of this rhomboid is denoted by the folding angle (Table 1),14a which indicates that, going from the four- to the sixmembered systems, the structure comes more relaxed from 110° to 122° toward a planar structure (180°). Accordingly, the twist angle of NiS4 denotes a more planar square for 3. It is important to note that the geometrical parameters of the Nin core of the metallocrowns 1, 2, and 3 are well reproduced by using methyl mercaptide ligands. This fact suggests that the robust toroid structure depicted by the NiS core remains similar despite the use of different ligands.1014 Electronic Structures. The resulting molecular electronic structure of the studied systems remains similar, as can be seen from the density of states (DOS) depicted in Figure 2, where the Fermi level (EF) is denoted by a dashed line. The energy gap between the highest occupied (HOMO) and lowest unoccupied Table 2. Energy Decomposition Analysis of the Interaction of a Ni Atom and the Rest of the Structure (kcal/mol), Charge Analysis (a.u.) via the Hirshfeld and Voronoi Deformation Density (VDD) Schemes, and HOMOLUMO gap (HL) 1

a

2

3

ΔEorb

435.64

42%

437.18

43%

464.29

45%

ΔEelstat

595.66

58%

582.72

57%

559.73

55%

ΔEpauli

+219.58

+211.27

+218.42

ΔEint Nia

811.72

808.63

805.6

Hirshfeld

+0.128

+0.131

+0.135

VDD

+0.092

+0.097

+0.110

HL

1.304

1.110

0.808

Charge analysis of each Ni nucleus.

Table 3. Values of δ(r) at Different Levels of Theory, Neglecting and Including Scalar and Spin-Orbit Terms, for 1, 2, and 3 at the Center of the Nin and Ni2S2 Rings, Respectively 1

a

2

3

Ni4b

Ni2S2c

Ni5b

Ni2S2c

Ni6b

Ni2S2c

NRa

+8.41

13.54

+7.89

13.71

+5.95

13.71

SRa

+8.70

13.24

+8.18

13.49

+6.19

13.59

SOa

+8.99

12.65

+8.29

12.85

+6.24

12.94

NR: Nonrelativistic; SR: scalar relativistic; SO: scalar + spinorbit calculations. b Center of Nin ring. c Center of Ni2S2 ring.

molecular orbital (LUMO) decreases with the increase of the toroid size (Table 2), where in all the studied systems the HOMO is mainly of Ni 3d antibonding character, and the LUMO is mainly of S 3p character (Supporting Information). To gain more insight into the bonding nature between the Ni center and the mercaptide ligand of the NiS4 moiety, we analyze the projected density of states (pDOS) of the 3d-Ni and 3p-S valence atomic levels for 1, 2, and 3. Figure 1 shows that the occupied molecular energy levels near the EF are mainly derived from the 3d-Ni shell, indicating that the Ni centers are largely involved in the chemical reactivity of these toroid systems.33 The bonding, nonbonding, and antibonding regions concerning the electronic structure are determined via the Mulliken orbital overlap population (MOOP) between the 3d-Ni shell and mercaptide ligands, where positive and negative values denote the bonding and antibonding regions, respectively, while the almost zero values denote the nonbonding regions (or almost nonbonding). The bonding block is located between 10.0 and 5.5 eV, in a region with mainly ligand character, while the antibonding region is located mainly at the LUMO. Thus, most of the antibonding combination of the NiS4 moiety remains unoccupied, ensuring the stability of the Niligand interaction, which in turn leads to a robust toroid structure. The mainly Ni 3d molecular levels (∼ 5.5 to 4.0 eV) are designated as nonbonding regions with a small antibonding character in the NiS4 fragment. Moreover, the increase of the amount of energy levels in the frontier band with orbitals of nonbonding character accounts for the decrease in the HOMOLUMO gap, making systems of higher nuclearities (n > 6) attractive targets for molecular electronics applications.34,35 The interaction between the nNi vertex was analyzed based on the overlap population between the 3d atomic shells of each redox center, and by using the Meyer bond-order procedure3032 at the scalar relativistic level. A decreasing NiNi interaction is concluded because of the small overlap between the basis functions going from 1 to 3, and by the Meyer bond orders values of 0.753 for 1, 0.172 for 2, and 0.143 for 3, suggesting that the stabilization of the toroid structure involves in part a NiNi bonding in 1, which decreases as a function of the NiNi bond length. Hence, the main stabilizing factor in 2 and 3 is the Ni2S2 building block, which is reinforced by the NiNi bond in 1. It is expected that these conclusions can be extended to other related systems with higher nuclearities. Binding Energy Analysis. Energy decomposition analysis (EDA) according to the MorokumaZiegler scheme36 was employed with the aim of obtaining a deeper description of the Niligand interaction, and how it varies with the number of Ni vertices. The interaction (binding energy) is calculated

Figure 3. Maps of δ(r) denoting shielded (negative values) and deshielded (positive values) areas in ppm, for 1, 2, and 3, in an xy contour-plane representation. 10791

dx.doi.org/10.1021/jp2028438 |J. Phys. Chem. A 2011, 115, 10789–10794

The Journal of Physical Chemistry A to be largely favorable by about 811.72 kcal/mol for 1, 808.63 kcal/mol for 2, and 805.60 kcal/mol for 3, as can be expected by the bonding metalligand interaction discussed in the Electronic Structures section (see above). Within this scheme, the interaction energy (ΔEint) can be further partitioned into three main components: ΔEint = ΔEpauli + ΔEelstat + ΔEorb. Here, the two first terms are computed by considering unperturbed fragments, and the latter is obtained when the densities are allowed to relax into the final molecular orbitals. The ΔEelstat term accounts for the stabilizing electrostatic interaction, and ΔEorb stands for the stabilizing covalent character of the fragmentfragment interaction, while the ΔEpauli term accounts for the Pauli (steric) repulsion. In this sense, the ΔEelstat represents more than the 55% of the total stabilizing energy (calculated as ΔEelstat /(ΔEelstat +ΔEorb) %) in the studied compounds. Hence, these results suggest that the electrostatic interaction plays the most important role, compared to the orbital interaction (∼45%), in the stabilization of the toroid complex. The variation of the stabilizing terms from 1 to 3 suggest that that a more strained structure leads to an interaction of more electrostatic character (58% for 1, with a NiNiNi angle of 90°, and 55% for 3, with a NiNiNi angle of 120°). In order to assign a total valence population to the Ni atoms (Table 2), we employed the Hirshfeld37a and Voronoi deformation density (VDD)37b partition schemes (Table 2). The Hirshfeld and VDD analyses show similar results, which suggests that the initially formal Ni2+ receives charge from the ligand, which decreases upon going to a more relaxed structure (3). Magnetic Response. The evaluation of the magnetic response of the studied toroid structures was carried out on the basis of NICS by Schleyer,19 mapping the induced magnetic field (Bind, in ppm units) over the space by application of an external magnetic field (Bext) of 1.0 T. These two quantities are related

Figure 4. Maps of δ(r) denoting shielded (negative values) and deshielded (positive values) areas in ppm, for 3, in an xz contour-plane.

ARTICLE

through the shielding tensor (σ), as is given elsewhere:18 ext Bind ðrÞ ¼  σ ðrÞ B

where the isotropic shielding value (δ(r)) can be determined by the negative of the trace of the second rank tensor σ(r)19 (δ(r) = NICS = (1/3)Trσ(r)) in units of ppm. It is noteworthy that in this type of system, deshielded areas (Table 3) alternate with a highly shielded region at the Ni2S2 ring (ca. 12 ppm at the center of Ni2S2), as can be seen from Figure 3. The diatropic currents are isolated in each Ni2S2 ring and decrease precipitously toward the center of the toroid (short-range induced field18b). In Figure 4, the xz cut-plane denotes the shape of the deshielded area through the center of the macrocyclic ring, which illustrates that the paratropic currents are retained inside the toroid structure until the Sn rings, where the values of δ(r) are close to 0.0 ppm. From 1 to 3, the value of δ(r) decreases in accordance with the increase of the diameter of the Nin ring, which decreases the interaction between the deshielded areas. From Figure 3, the isotropic magnetic response denotes the similar magnetic behavior of the studied systems. In order to gain more insight into the electronic delocalization inside the toroid structure, we focus on the magnetic criteria of induced electronic delocalization (aromaticity),38 where an applied external magnetic field through the z-axis in an aromatic ring system leads to diamagnetic ring currents characterized by shielded response in the z-axis (negative values of the zz component of δ(r)). On the contrary, deshielded areas denote paramagnetic rings currents usually characterized in antiaromatic ring systems,18,19 by the map representation of the zz component of δ(r) (δzz(r)) in Figure 5. At the center of each Nin ring, the positive values indicate deshielded areas, which in turn denote paratropic currents, as has been noted recently;21 however, by the map of δzz (r), it is possible to observe that the d aromaticity inside of the Nin ring leads to diatropic ring currents (aromaticity) at the closeness of the toroid structure, which in turn induce a paratropic response at the center of the toroid structure with values of δzz (r) of 1.56 ppm for 1, 3.87 ppm for 2, and 6.33 ppm for 3, at the spinorbit relatisvitic level. This behavior contrasts with the aromatic properties of organic molecules,18b which show pure shielded response inside the ring. This aromatic behavior suggests an enhanced stabilization,39 which is more pronounced in 1. In Table 3, the isotropic δ(r) value indicates that the Ni2S2 ring exhibits diatropic ring currents indepently of the external magnetic field (Bext) axis and the size of the system, supporting its size-independent stability and its consideration and use as a constructing block.

Figure 5. Maps of δzz (r) denoting shielding zones (negative values) close to the nNi ring and a deshielding zone (positive values) at the center. 10792

dx.doi.org/10.1021/jp2028438 |J. Phys. Chem. A 2011, 115, 10789–10794

The Journal of Physical Chemistry A In addition, the relativistic effects in these nonrelativistic molecules have been revisited by comparison of the magnetic response to an external magnetic field via nonrelativistic, scalar relativistic, and scalar + spinorbit calculations.24 It is quite important to note this variation, in the sense that the nickel is termed as a nonrelativistic atom; however, the titled effects vary in their magnetic behavior, as can be seen from Table 3. The different levels of theory vary the value of the induced magnetic field, increasing the unshielded areas (center of Nin) and decreasing the shielded areas (center of Ni2S2), which can be understood in terms of the relativistic valence atomic spinorbit shells 3d3/2 and 3d5/2 for Ni, and 3p1/2 and 3p3/2 for S,2224 which in fact take part in the magnetic and related properties.24

4. CONCLUSIONS Analysis of the electronic structure and induced magnetic fields have been achieved via the relativistic two-component ZORA Hamiltonian including scalar and spinorbit coupling. The calculated data suggest that the stability arises mainly from the bonding interaction between the mercaptide ligand and the metallic center in the Ni2S2 fragment, and the intermetallic interaction is almost negligible in the multimetallic toroid structures with n > 5. By calculating the induced magnetic fields by an applied external field, we achieve an understanding of the reason for the paratropic current behavior in terms of n paratropic sites (Ni vertex), which, in conjunction, leads to a paratropic current at the center. The zz-component of the shielding tensor suggests an aromatic behavior confined to the closeness of the Nin ring, which leads to deshielded zones at the center of structure. The Ni2S2 moiety exhibits aromatic character, a property that, in addition to the high stabilization of the NiS4 interaction, allows one to conclude that the rhomboid building block ensures the stability of the robust toroid structure. In addition, the relativistic effects vary the behavior of the induced magnetic fields, even in the systems termed as “nonrelativistic compounds” (Ni: Z = 28; S: Z = 16; C: Z = 12; H: Z = 1, where Z it is the atomic number). Which, denotes that the relativistic terms influence several molecular properties. ’ ASSOCIATED CONTENT

bS

Supporting Information. Isosurfaces of the HOMO and LUMO. This material is available free of charge via the Internet at http://pubs.acs.org.

’ AUTHOR INFORMATION Corresponding Author

*E-mail: [email protected].

’ ACKNOWLEDGMENT The author ackowledges the financial support of FONDECYT Grant 11100027 and PROJECT MILLENNIUM No. P07-006F. The author thanks the reviewers for their helpful comments. ’ REFERENCES (1) (a) Pedersen, C. J. J. Am. Chem. Soc. 1967, 89, 2495. (b) Pedersen, C. J. Angew. Chem., Int. Ed. Engl. 1988, 27, 1021. (2) (a) Slone, R. V.; Benkstein, K. D.; Belanger, S.; Hupp, J. T.; Guzei, I. A.; Rheingold, A. L. Coord. Chem. Rev. 1998, 171, 221.(b) Lehn, J.-M. Supramolecular Chemistry. Concepts and Perspectives; VCH: Weinhein,

ARTICLE

Germany, 1995. (c) Dietrich, B., Viout, P., Lehn, J.-M. Macrocyclic Chemistry: Aspects of Organic and Inorganic Supramolecular Chemistry; VCH: New York, 1993. (3) (a) Gokel, G. W. Crown Ethers and Cryptands: Monographs in Supramolecular Chemistry; Royal Society of Chemistry: Cambridge, U. K., 1991. (b) V€ogtle, F. Supramolecular Chemistry: An Introduction; Wiley: Chichester, U.K., 1991. (c) Beer, P. D. Adv. Inorg. Chem. 1992, 39, 79. (d) Van Veggel, F. C.; Verboom, W.; Reinhoudt, D. N. Chem. Rev. 1994, 94, 279. (4) (a) Basilio, N.; Garcia-Rio, L.; Mejuto, J. C.; Perez-Lorenzo, M. J. Org. Chem. 2006, 71, 4280. (b) Bourson, J.; Pouget, J.; Valeur, B. J. Phys. Chem. 1993, 97, 4552. (c) Wang, Z.; Chang, S. H.; Kang, T. J. Spectrochim. Acta A 2008, 70, 3135. (d) Horwitz, E. P.; Dietz, M. L.; Fisher, D. E. Solvent Extr. Ion Exch. 1991, 9, 1. (e) Gokel, G. W.; Leevy, W. M.; Weber, M. E. Chem. Rev. 2004, 104, 2723. (5) (a) Allan, R. E.; Bashall, A.; Palmer, J. S.; McPartlin, M.; Mosquera, M. E. G.; Rawson, J. M.; Wheatley, A. E. H.; Wright, D. S. Chem. Commun. 1997, 1975. (b) Lah, M. S.; Gibney, B. R.; Tierney, D. L.; Penner-Hahn, J. E.; Pecorato, V. L. J. Am. Chem. Soc. 1993, 115, 5857. (c) Huang, S.-P.; Kanatzidis, M. G. Angew. Chem., Int. Ed. Engl. 1992, 31, 787. (d) Muller, A.; Hovemeier, K.; Rohlfing, R. Angew. Chem., Int. Ed. Engl. 1992, 31, 1192. (e) Reuter, H. Angew. Chem., Int. Ed. Engl. 1992, 31, 1185. (6) (a) Pecoraro, V. L. Inorg. Chim. Acta 1989, 155, 171. (b) Pecoraro, V. L.; Stemmler, A. J.; Gibney, B. R.; Bodwin, J. J.; Wang, H.; Kampf, J. W.; Barwinski, A. Prog. Inorg. Chem. 1997, 45, 83. (c) Faust, T. B.; Heath, P. G.; Muryn, C. A.; Timco, G. A.; Winpenny, R. E. Chem. Commun. 2010, 46, 6258. (d) Song, L.-C.; Gao, J.; Wang, H.-T.; Hua, Y.-J.; Fan, H.-T.; Zhang, X.-G.; Hu, Q.-M. Organometallics 2006, 25, 5724. (e) Mimassi, L.; Guyard-Duhayon, C.; Noelle Rager, M.; Amouri, H. Inorg. Chem. 2004, 43, 6644. (7) Shibasaki, M., Yamamoto, Y., Eds. Multimetallic Catalysts in Organic Synthesis; Wiley-VCH: Weinheim, Germany, 2004. (8) (a) Ohki, Y.; Ikagawa, Y.; Tatsumi, K. J. Am. Chem. Soc. 2007, 129, 10457. (b) Laplaza, C. E.; Holm, R. H. J. Am. Chem. Soc. 2001, 123, 10255. (9) Mezei, G.; Zaleski, C. M.; Pecoraro, V. L. Chem. Rev. 2007, 107, 4933and references therein. (10) (a) Woodward, P.; Dahl, L. F.; Abel, E. W.; Crosse, B. C. J. Am. Chem. Soc. 1965, 87, 5251. (b) Abel, E. W.; Crosse, B. C. J. Chem. Soc. 1966, 1377. (11) Jensen, K. A. Z. Anorg. Chem. 1944, 252, 227. (12) (a) Gaete, W.; Ros, J.; Solans, X.; Font-Altaba, M.; Brianso, J. L. Inorg. Chem. 1984, 23, 39. (b) Kriege, M.; Henkel, G. Z. Naturforsch. B 1987, 42, 1121. (c) Kr€uger, T.; Krebs, B.; Henkel, G. Angew. Chem., Int. Ed. Engl. 1989, 28, 61. (d) M€uller, A.; Henkel, G. Naturforsch Z. B: Chem. Sci. 1997, 52, 1501. (13) (a) Mahmoudkhani, A. H.; Langer, V. Inorg. Chim. Acta 1999, 294, 83. (b) Koo, B.-K.; Block, E.; Kang, H.; Liu, S.; Zubieta, J. Polyhedron 1988, 7, 1397. (14) (a) Wark, T. A.; Stephan, D. W. Organometallics 1989, 8, 2836. (b) Schulbert, K.; Mattes, R. Z. Naturforsch. B 1994, 49, 770. (c) Gould, R. O.; Harding, M. M. J. Chem. Soc. A 1970, 875. (d) Mahmoudkhani, A. H.; Langer, V. Polyhedron 1999, 18, 3407. (e) Jian, F.-F.; Jiao, K.; Li, Y.; Zhao, P.-S.; Lu, L.-D. Angew. Chem., Int. Ed. 2003, 42, 5722. (15) (a) Dance, I. G.; Scudder, M. I.; Secomb, R. Inorg. Chem. 1985, 24, 1201. (b) Ivanov, S. A.; Kozee, M. A.; Merrill, W. A.; Agarwal, S.; Dahl, L. F. J. Chem. Soc., Dalton Trans. 2002, 4105. (c) Kunchur, N. R. Acta Crystallogr. 1968, B24, 1623. (d) Higgins, J. D., III; Suggs, J. W. Inorg. Chim. Acta 1988, 145, 247. (16) Baranov, A. I.; Kloo, L.; Olenev, A. V.; Popovkin, B. A.; Romanenko, A. I.; Shevelkov, A. V. J. Am. Chem. Soc. 2001, 123, 12375. (17) (a) Ozin, G. A., Arsenault, A. C. Nanochemistry; Royal Society of Chemistry: Cambridge, U.K., 2005. (b) Ozin, G. A. Adv. Mater. 1992, 4, 612. (c) Jin, P.; Li, F.; Chen, Z. J. Phys Chem. A 2011, 115, 2402. (18) (a) Sebastiani, D. ChemPhysChem 2006, 7, 164. (b) Merino, G.; Heine, T.; Seifert, G. Chem.—Eur. J. 2004, 10, 4367. (c) Sebastiani, D.; Kudin, K. N. ACS Nano 2008, 2, 661. (d) Klod, S.; Kleinpeter, E. J. Chem. 10793

dx.doi.org/10.1021/jp2028438 |J. Phys. Chem. A 2011, 115, 10789–10794

The Journal of Physical Chemistry A Soc., Perkin Trans. 2001, 2, 1893. (e) Aihara, J. I.; Oe, S. Bull. Chem. Soc. Jpn. 2003, 76, 1363. (19) (a) Schleyer, P.v. R; Maerker, C.; Dransfeld, A.; Jiao, H.; Hommes, H. J. R. J. Am. Chem. Soc. 1996, 118, 6317. (b) Chen, Z.; Wannere, C. S.; Corminboeuf, C.; Puchta, R.; Schleyer, P. v. R. Chem. Rev. 2005, 105, 3842. (c) Chen, Z.; Corminboeuf, C.; Heine, T.; Bohmann, J.; Schleyer, P. v. R. J. Am. Chem. Soc. 2003, 125, 13930. (20) (a) Zhang, C.-G.; Zheng, F.; Dixon, D. A. J. Am. Chem. Soc. 2002, 124, 14795. (b) Kuznetsov, A. E.; Boldyrev, A. I.; Li, X.; Wang, L.-S. J. Am. Chem. Soc. 2001, 123, 8825. (c) Kuznetsov, A. E.; Boldyrev, A. I. Inorg. Chem. 2002, 41, 532and references therein. (21) Datta, A.; John, N. S.; Kulkarni, G. U.; Pati, S. K. J. Phys. Chem. A 2005, 109, 11647. (22) (a) Mu~noz-Castro, A.; MacLeod Carey, D.; Arratia-Perez, R. Chem. Phys. Lett. 2009, 474, 290. (b) Alvarado-Soto, L.; Ramirez-Tagle, R.; Arratia-Perez, R. Chem. Phys. Lett. 2008, 467, 94. (c) Mu~ noz-Castro, A.; MacLeod Carey, D.; Arratia-Perez, R. J. Phys. Chem. A 2010, 114, 666. (d) Mu~noz-Castro, A.; MacLeod Carey, D.; Arratia-Perez, R. ChemPhysChem 2010, 11, 646. (e) Mu~noz-Castro, A.; Arratia-Perez, R. J. Phys. Chem. A 2010, 114, 5217. (23) Mu~noz-Castro, A.; MacLeod Carey, D.; Arratia-Perez, R. J. Chem. Phys. 2010, 132, 164308. (24) (a) Dyall, K. G.; Fægri, K. Introduction to Relativistic Quantum Chemistry; Oxford University Press: New York, 2007. (b) Malli, G. L.; Oreg, J. J. Chem. Phys. 1975, 63, 830. (c) Malli, G. L. Chem. Phys. Lett. 1980, 73, 510. (d) Malli, G. L.; Pyper, N. C. Proc. R. Soc. A 1986, 407, 377. (25) Amsterdam Density Functional (ADF) Code, release 2010; Vrije Universiteit Amsterdam: The Netherlands, 2011. (26) Van Lenthe, E.; Baerends, E. J.; Snijders, J. G. J. Chem. Phys. 1994, 101, 9783. (27) Te Velde, G.; Bickelhaupt, F. M.; Van Gisberger, S. J. A.; Guerra, C. F.; Baerends, E. J.; Snijders, J. G.; Ziegler, T. J. Comput. Chem. 2001, 22, 931. (28) (a) Perdew, J. P.; Burke, K.; Wang, Y. Phys. Rev. B 1996, 54, 16533. (b) Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1996, 77, 3865. (c) Swart, M.; Groenhof, A. R.; Ehlers, A. W.; Lammertsma, K. J. Phys. Chem. A 2004, 108, 5479. Swart, M. J. Chem. Theory Comput. 2008, 4, 2057. (d) van der Wijst, T.; Fonseca Guerra, C.; Swart, M.; Bickelhaupt, F. M. Chem. Phys. Lett. 2006, 426, 415. (e) Swart, M.; Sola, M.; Bickelhaupt, F. M. J. Chem. Phys. 2009, 131, 094103. (f) Handy, N. C.; Cohen, A. J. Mol. Phys. 2001, 99, 403. (29) Verluis, L.; Ziegler, T. J. Chem. Phys. 1988, 88, 322. (30) (a) Mayer, I. Chem. Phys. Lett. 1983, 97, 270. (b) Bridgeman, A. J.; Cavigliasso, G.; Ireland, L. R.; Rothery, J. J. Chem. Soc., Dalton Trans. 2001, 2095. (31) Bridgeman, A. J.; Cavigliasso, G. J. Phys. Chem. A 2003, 107, 4568. (32) Bridgeman, A. J.; Cavigliasso, G. Faraday Discuss. 2003, 124, 239. (33) (a) Parr, R. G.; Yang, W. Density Functional Theory for Atoms and Molecules; Oxford University Press: New York, 1982. (b) Pearson, R. G. Chemical Hardness; Wiley-VCH: Germany, 1997. (34) (a) Metzger, R. M. Acc. Chem. Res. 1999, 32, 950. (b) Perepichka, D. F.; Bryce, M. R. Angew. Chem., Int. Ed. 2005, 44, 5370. (c) Roncali, J. Chem. Rev. 1997, 97, 173. (d) Kobayashi, A.; Fujiwara, E.; Kobayashi, H. Chem. Rev. 2004, 104, 5243. (e) Mu~noz-Castro, A.; MacLeod Carey, D.; Morales-Verdejo, C.; Chavez, I.; Manríquez, J.-M.; Arratia-Perez, R. Inorg. Chem. 2010, 49, 4175. (35) (a) Sonmez, G.; Meng, H.; Wudl, F. Chem. Mater. 2003, 15, 4923. (b) Carroll, R. Ll.; Gorman, C. B. Angew. Chem., Int. Ed. 2002, 41, 4378. (36) (a) Bickelhaupt, F. M.; Baerends, E. J. In Reviews in Computational Chemistry; Lipkowitz, K. B., Boyd, D. B., Eds.; Wiley-VCH: New York, 2000; Vol. 15, p 1. (b) Morokuma, K. J. Chem. Phys. 1971, 55, 1236. (c) Ziegler, T.; Rauk, A. Theor. Chim. Acta 1977, 46, 1. (d) Vyboishchikov, S. F.; Krapp, A.; Frenking, G. J. Chem. Phys. 2008, 129, 144111.

ARTICLE

(37) (a) Hirshfeld, F. L. Theor. Chim. Acta 1977, 44, 129. (b) Fonseca Guerra, C.; Handgraaf, J.-W.; Baerends, E. J.; Bickelhaupt, F. M. J. Comput. Chem. 2004, 25, 189. (38) (a) Lazzaretti, P. Prog. Nucl. Magn. Reson. Spectrosc. 2000, 36, 1. (b) Carion, R.; Liegeois, V.; Champagne, B.; Bonifazi, D.; Pelloni, S.; Lazzeretti, P. J. Phys. Chem. Lett. 200, 1, 1563. (c) Havenith, R. W. A.; Fowler, P. W. Phys. Chem. Chem. Phys. 2006, 8, 3383. (d) Heine, T.; Corminboeuf, C.; Seifert, G. Chem. Rev. 2005, 105, 3889. (39) (a) Zhao, J.; Xu, J.; King, R. B. Inorg. Chem. 2008, 47, 9314. (b) Denk, M. K.; Gupta, S.; Lough, A. J. Eur. J. Inorg. Chem. 1999, 41. (c) King, R. B. Chem. Rev. 2001, 101, 1119. (d) Gudat, D.; Haghverdi, A.; Hupfer, H.; Nieger, M. Chem.—Eur. J. 2000, 6, 3414.

10794

dx.doi.org/10.1021/jp2028438 |J. Phys. Chem. A 2011, 115, 10789–10794