Dependence of ion hydration on the sign of the ion's charge

The solvation of simple ions in water is studied using molecular dynamics simulations with a polarizable force field. Previous simulations using this potential demonstrated that anions are more favorably solvated in water than cations. The present work is an attempt to explain this result by examining the effects of ions on the surrounding water structure, with particular focus on the first solvation shell and its interactions with the surrounding water. We conclude that while the first solvation shell surrounding cations is frustrated by competition between ion-water and water-water interactions, solvation of anions is compatible with good water-water interactions.     

Li+ cation environment, transport, and mechanical properties of the LiTFSI doped N-methyl-N-alkylpyrrolidinium+TFSI- ionic liquids

Molecular dynamics (MD) simulations have been performed on N-methyl-N-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide (mppy(+)TFSI(-)) and N,N-dimethyl- pyrrolidinium bis(trifluoromethanesulfonyl)imide (mmpy(+)TFSI(+)) ionic liquids (ILs) doped with 0.25 mol fraction LiTFSI salt at 303-500 K. The liquid density, ion self-diffusion coefficients, and conductivity predicted by MD simulations were found to be in good agreement with experimental data, where available. MD simulations reveal that the Li(+) environment is similar in mppy(+)TFSI(-) and mmpy(+)TFSI(+) ILs doped with LiTFSI. The Li(+) cations were found to be coordinated on average by slightly less than four oxygen atoms with each oxygen atom being contributed by a different TFSI(-) anion. Significant lithium aggregation by sharing up to three TFSI(-) anions bridging two lithiums was observed, particularly at lower temperatures where the lithium aggregates were found to be stable for tens of nanoseconds. Polarization of TFSI(-) anions is largely responsible for the formation of such lithium aggregates. Li(+) transport was found to occur primarily by exchange of TFSI(-) anions in the first coordination shell with a smaller (approximately 30%) contribution also due to Li(+) cations diffusing together with their first coordination shell. In both ILs, ion self-diffusion coefficients followed the order Li(+) < TFSI(-) < mmpy(+) or mppy(+) with all ion diffusion in mmpy(+)TFSI(-) being systematically slower than that in mppy(+)TFSI(-). Conductivity due to the Li(+) cation in LiTFSI doped mppy(+)TFSI(-) IL was found to be greater than that for a model poly(ethylene oxide)(PEO)/LiTFSI polymer electrolyte but significantly lower than that for an ethylene carbonate/LiTFSI liquid electrolyte. Finally, the time-dependent shear modulus for the LiTFSI doped ILs was found to be similar to that for a model poly(ethylene oxide)(PEO)/LiTFSI polymer electrolyte on the subnanosecond time scale.     

Structural investigations of polymer electrolyte poly(propylene oxide)-LiClO4 using diffraction experiments and reverse Monte Carlo simulation

The structure of an amorphous polymer electrolyte, poly(propylene oxide) (PPO) complexed with LiClO4, has been studied using reverse Monte Carlo (RMC) simulations. The simulations require no force field but are based on experimental data only, in this case from x-ray and neutron diffraction experiments. Excellent agreement between the experimental data and the structures resulting from the RMC simulation is obtained. Samples with ether-oxygen to lithium concentrations (molar ratios) O:Li=16:1 and 5:1 were studied and compared to results of pure PPO from a previous study. We focus on the effects of the solvated salt on the structure of the polymer matrix, the spatial distribution of ions, and the correlations between the anions and the polymer chains. Analyzing the structures produced in the simulations, we find that for a concentration 16:1, the interchain distance is approximately the same as in pure PPO but more well defined. For a concentration 5:1, we find a larger and less well-defined interchain distance compared to the 16:1 concentration. This signifies that at the 16:1 salt concentration, there is enough free volume in the polymer host to accommodate the ions, and that the solvation of salt induces ordering of the polymer matrix. At the higher salt concentration 5:1, the polymer network must expand and become less ordered to host the ions. We also note, in accordance with previous studies, that the solvation of salt changes the conformation of the polymer chain towards more gauche states. The simulations furthermore reveal marked correlations between the polymer chains and the anions, which we suggest arise predominantly from an interaction mediated via cations, which can simultaneously coordinate both ether oxygens in the polymer chains and anions. Interanionic distances at 5 A, which are consistent with two or more anions being coordinated around the same cation, are also observed. On a larger scale, the RMC structure of PPO-LiClO4 16:1 clearly indicates the presence of salt-rich and salt-depleted domains having a length scale of <20 A. In view of such a heterogeneous structure of PPO-LiClO4 16:1, it is plausible that the increased ordering of the polymer matrix is due to rather well-defined structural arrangements within the salt-rich domains, and that the characteristic interchain distance in the salt-rich domains is similar to that of the pure polymer.     

Solvation of cations in the LiNO3-Ca(NO3)2-H2O system at 25°C. A molecular dynamics simulation study

A series of solutions of the ternary LiNO₃-Ca(NO₃)₂-H₂O system is simulated under standard conditions by means of classical molecular dynamics (MD). Several variants of the potentials of interparticle interactions predicting different structures and compositions of the lithium cation hydration shell are used in the calculations. The coordination numbers of ions, the mean residence times of water molecules and nitrate anions in the solvation shells of cations, and the self-diffusion coefficients of solution components are estimated for all of the investigated systems. The structural features of the solvation shells of cations in multicomponent aqueous solutions are described on the basis of the obtained data.     

盐对甲醇微观结构的影响

利用拉曼光谱研究盐对甲醇微观结构的影响.比较了不同盐/甲醇体系的O—H伸缩谱段和C—O伸缩谱段的超额拉曼光谱,对比给出了阴、阳离子与甲醇的相互作用.O—H伸缩谱段的超额拉曼光谱明显地显示了阴离子与甲醇形成弱氢键,氢键强度排序为CH3OH-CH3OHCl--CH3OHNO3--CH3OHClO4--CH3OH,在这个波段内,基本观察不到阳离子与甲醇的相互作用.在C—O伸缩谱段内,阴阳离子均有显著的体现,且与它们作用的甲醇C—O伸缩振动频率有如下的关系:CH3—OH(阴离子)CH3—OH(体相)CH3—OH(阳离子).根据C—O伸缩谱段的超额拉曼光谱,拟合了该谱段的拉曼光谱,由分解的谱峰强度得到阴、阳离子第一溶剂化层中甲醇分子的数目,结果显示在该浓度(～0.005)下离子对第一溶剂化层以外的甲醇氢键网络结构没有明显影响.     

Solvation of uranium hexachloro complexes in room-temperature ionic liquids. A molecular dynamics investigation in two liquids

We report a molecular dynamics study of the solvation of UCl(6)(-), UCl(6)(2-), and UCl(6)(3-) complexes in the [BMI][Tf(2)N] and [MeBu(3)N][Tf(2)N] ionic liquid cations based on the same anion (bis(trifluoromethylsulfonyl)imide (Tf(2)N-)) and the butyl-3-methyl-imidazolium+ (BMI+) or methyl-tri-n-butyl-ammonium (MeBu(3)N+) cation, respectively. The comparison of two electrostatic models of the complexes (ionic model with -1 charged halides versus quantum mechanically derived charges) yields similar solvation features of a given solute. In the two liquids, the first solvation shell of the complexes is positively charged and evolves from purely cationic in the case of UCl(6)(3-) to a mixture of cations and anions in the case of UCl(6)(-). UCl(6)(3-) is exclusively "coordinated" to BMI+ or MeBu(3)N+ solvent cations that mainly interact via their CH aromatic protons or their N-Me group, respectively. Around the less charged UCl(6)(-) complex, the cations interact via the less polar moieties (butyl chains of BMI+ or MeBu(3)N+) and the anions display nonspecific interactions. In no case does the uranium atom further coordinate solvent ions. According to an energy components analysis, UCl(6)(3-) interacts more attractively with the [BMI][Tf(2)N] liquid than with [MeBu(3)N][Tf(2)N], while UCl(6)(-) does not show any preference, suggesting a significant solvation effect of the redox properties of uranium, also supported by free energy perturbation simulations. The effect of ionic liquid (IL) humidity is investigated by simulating the three complexes in 1:8 water/IL mixtures. In contrast to the case of "naked" ions (e.g., lanthanide(3+), UO2(2+), alkali, or halides), water has little influence on the solvation of the UCl(6)(n-) complexes in the two simulated ILs, as indicated by structural and energy analysis. This is in full agreement with the experimental observations (Nikitenko, S. I.; et al. Inorg. Chem. 2005, 44, 9497).     

LiTFSI structure and transport in ethylene carbonate from molecular dynamics simulations

Molecular dynamics (MD) simulations using a many-body polarizable force field were performed on ethylene carbonate (EC) doped with lithium bistrifluoromethanesulfonamide (LiTFSI) salt as a function of temperature and salt concentration. At 313 K Li+ was coordinated by 2.7-3.2 EC carbonyl oxygen atoms and 0.67-1.05 TFSI- oxygen atoms at EC:Li = 10 and 20 salt concentrations. In completely dissociated electrolytes, however, Li+ was solvated by approximately 3.8 carbonyl oxygen atoms from EC on average. The probability of ions to participate ion aggregates decreased exponentially with an increase in the size of the aggregate. Ion and solvent self-diffusion coefficients and conductivity predicted by MD simulations were in good agreement with experiments. Approximately half of the charge was transported by charged ion aggregates with the other half carried by free (uncomplexed by counterion) ions. Investigation of the Li+ transport mechanism revealed that contribution from the Li+ diffusion together with its coordination shell to the total Li+ transport is similar to the contribution arising from Li+ exchanging solvent molecules in its first coordination shell with solvents from the outer shells.     

Investigation of fundamental transport properties and thermodynamics in diglyme-salt solutions

Ionic mobility, the thermodynamics of ionic association, and the structure of associated species are studied in solutions of diglyme containing either lithium triflate or tetrabutylammonium triflate. Infrared spectroscopic, PFG NMR, thermodynamic, and crystallographic data suggest that the solute species existing in diglyme-lithium triflate are "free" ions, contact ion pairs, and dimers. Equilibrium constants, S(o), deltaH(o), and deltaG(o) are calculated for processes occurring between these species. In particular, the equilibrium constant, corrected for nonideality using a modified Debye-Hückel expression, is calculated for the dissociation of contact ion pairs into "free" cations and anions. A second equilibrium constant for the formation of dimers from contact ion pairs is also calculated; these constants do not significantly vary with salt concentration up to about 1.3 x 10(-3) mol cm(-3). The measured temperature dependence of equilibrium constants was used to calculate deltaH(o) and deltaS(o) for the two processes. The value of deltaS(o) = -102 J mol(-1) K(-1) for the dissociation of contact ion pairs shows that the large entropy decrease due to cation solvation outweighs the entropy increase due to dissociation of a contact ion pair. Ionic mobilities are calculated in lithium triflate-diglyme solutions using conductivity data in conjunction with information about the nature and concentrations of solute species obtained from IR spectroscopy. Mobilities in tetrabutlyammonium triflate-diglyme solutions are calculated directly from conductivity data. It was concluded that the concentration dependence of the molar conductivity is due in large part to the variation of the ion mobilities with concentration.     

Europium(III) and its halides in anhydrous room-temperature imidazolium-based ionic liquids: a combined TRES, EXAFS, and molecular dynamics study

Combining spectroscopic techniques (TRES and EXAFS) and molecular dynamics simulations, we have investigated the state of trivalent europium dissolved in room-temperature ionic liquids (RTILs), as a function of the RTIL anion and in the presence of added chloride anions. The studied RTILs are based on the 1-butyl-3-methyl-imidazolium (Bumim+) cation and differ by their anionic counterparts: BF4-, PF6-, Tf- (triflate, CF3SO3-), and Tf2N- [(CF3SO2)2N-]. The results show the strong influence of the RTIL nature on the first solvation shell of europium and on its complexation with chloride. Depending on the RTIL, europium(III), which was introduced in solution as a triflate salt, is found to be solvated either by RTIL anions only or as neutral undissociated EuTf3 moieties completed by solvent anions. Kinetic effects, related to the viscosity of the RTIL and the nature of the europium salt, also markedly influence the coordination of added Cl- or F- anions to the metal.     

Properties of Apolar Solutes in Alkyl Imidazolium‐Based Ionic Liquids: The Importance of Local Interactions

The solvation and the dynamic properties of apolar model solutes in alkyl imidazolium‐based ionic liquids (IL) are studied by using all‐atom molecular dynamics simulations. In regards to specific IL effects, we focused on the often used 1‐ethyl‐3‐methyl imidazolium cation in combination with the anions tetrafluoroborate, acetate, and bis(trifluoromethanesulfonyl)imide. Our findings reveal that the size of the anion crucially influences the accumulation behavior of the cations, which results in modified IL solvation properties. Deviations between the different alkyl imidazolium‐based IL combinations can be also observed with regard to the results for the radial distribution functions, the number of surrounding molecules, and the molecular orientation. The analysis of the van Hove function further shows pronounced differences in the dynamic behavior of the solutes. The simulations verify that the solute mobilities are mainly influenced by the composition of the local solvent shell and the properties of the underlying Lennard–Jones interactions. Additional simulations with regard to modified short‐range dispersion energies for alkyl imidazolium‐based ILs validate our conclusions.     

Solvation of uranyl(II), europium(III) and europium(II) cations in "basic" room-temperature ionic liquids: a theoretical study

We report a molecular dynamics study of the solvation of UO2(2+), Eu3+ and Eu2+ ions in two "basic" (Lewis acidity) room-temperature ionic liquids (IL) composed of the 1-ethyl-3-methylimidazolium cation (EMI+) and a mixture of AlCl4- and Cl- anions, in which the Cl-/AlCl4- ratio is about 1 and 3, respectively. The study reveals the importance of the [UO2Cl4]2- species, which spontaneously form during most simulations, and that the first solvation shell of europium is filled with Cl- and AlCl4- ions embedded in a cationic EMI+ shell. The stability of the [UO2Cl4]2- and [Eu(III)Cl6]3- complexes is supported by quantum mechanical calculations, according to which the uranyl and europium cations intrinsically prefer Cl- to the AlCl4- ion. In the gas phase, however, [Eu(III)Cl6]3- and [Eu(II)Cl6]4- complexes are predicted to be metastable and to lose two to three Cl- ions. This contrasts with the results of simulations of complexes in ILs, in which the "solvation" of the europium complexes increases with the number of coordinated chlorides, leading to an equilibrium between different chloro species. The behavior of the hydrated [Eu(OH2)8]3+ complex is considered in the basic liquids; the complex exchanges H2O molecules with Cl- ions to form mixed [EuCl3(OH2)4] and [EuCl4(OH2)3]- complexes. The results of the simulations allow us to better understand the microscopic nature and solvation of lanthanide and actinide complexes in "basic" ionic liquids.     

Anion-promoted cation motion and conduction in zeolites

The motion of sodium cations in sodalite and cancrinite has been investigated by force field calculations, solid-state NMR, and impedance spectroscopy. Special emphasis is dedicated to the influence of anions on sodium mobilities. Local cation motion is promoted when they interact with anions. However, not all systems with high local mobilities exhibit good ion conductivities, as cooperativity of the motion appears to be an important factor, as well. The activation barrier for local sodium motion (calculations) and long-range transport (dc conductivities) is lowered in sodalite when halogenide anions, Cl(-), Br(-), or I(-), are present. The activation barriers increase with increasing size of the anion and decreasing coordination in the transition state. On the basis of (23)Na solid-state NMR data, all the sodium ions in the dense sodalite structure are rather rigid up to 470 K. All the cations in chromate sodalite, and Na(+) in the small cancrinite epsilon-cages without anion interactions, show a restricted local motion at higher temperatures. There is a selective high local motion of Na(+) in the neighborhood of chromate anions in the more open channel system of cancrinite. These results suggest that sodium migration can be enhanced, at least locally, in open channel systems by anion interactions. A dynamics coupling between anion reorientation and cation mobility was not observed.     

Effects of ion solvation on phase equilibrium and interfacial tension of liquid mixtures

We study the bulk thermodynamics and interfacial properties of electrolyte solution mixtures by accounting for electrostatic interaction, ion solvation, and inhomogeneity in the dielectric medium in the mean-field framework. Difference in the solvation energy between the cations and anions is shown to give rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The ion solvation affects the phase equilibrium of the solvent mixture, depending on the dielectric constants of the solvents, reflecting the competition between the solvation energy and translation entropy of the ions. Miscibility is decreased if both solvents have low dielectric constants and is enhanced if both solvents have high dielectric constant. At the mean-field level, the ion distribution near the interface is determined by two competing effects: accumulation in the electrostatic double layer and depletion in a diffuse interface. The interfacial tension shows a nonmonotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations and decreases approximately as the square root of the salt concentration for dilute solutions, reaching a minimum near 1 mM. We also find that, for a fixed cation type, the interfacial tension decreases as the size of anion increases. These results offer qualitative explanations within one unified framework for the long-known concentration and ion size effects on the interfacial tension of electrolyte solutions.     

Molecular dynamics simulation of LiTFSI-acetamide electrolytes: structural properties

The liquid structures of nonaqueous electrolytes composed of lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) and acetamide, with LiTFSI/acetamide molar ratios of 1:2, 1:4, and 1:6, were studied by molecular dynamics simulations. The simulations indicate that the Li+ cations prefer to be six-coordinate by the sulfonyl oxygen atoms of the TFSI- anions and the carbonyl oxygen atoms of the acetamide molecules, rather than by the most electronegative nitrogen atom of the TFSI- anion. Therefore, close Li+-TFSI- contact pairs exist in the system. The TFSI- anion prefers to provide only one of four possible oxygen atoms to coordinate to the same Li+ cation. Three conformations (cis, trans, and gauche) of the TFSI- anions were found to coexist in the liquid electrolyte. At high salt concentrations, the TFSI- anions mainly adopt the gauche conformation in order to provide more oxygen atoms to coordinate to different Li+ cations, while simultaneously reducing the repulsion among the Li+ cations. On the other hand, the fraction of TFSI- anions adopting the cis conformation is largest for the system with the molar ratio of 1:6, in which many clusters, mainly composed of the Li+ cations and the TFSI- anions, are immersed in the acetamide molecules. The size and charge distribution of clusters were also investigated. In the system with the molar ratio of 1:2, nearly all of the ions in the PBC (periodic boundary conditions) box aggregate into a bulky cluster that gradually disassembles into small clusters with decreasing salt concentration. The addition of acetamide molecules was found to effectively relax the liquid electrolyte structure, and the system with the molar ratio of 1:4 was found to exhibit a more homogeneous liquid structure than the other two electrolyte systems with molar ratios of 1:2 and 1:6.     

On the structure of aqueous LiCl solutions

The structure of highly concentrated aqueous lithium chloride solutions was investigated by the Reverse Monte Carlo method. Two total structure factors, obtained from neutron and x-ray diffraction experiments, were applied as input information. From the resulting particle configurations, partial pair correlation functions, coordination numbers and cosine distributions of bond angles have been determined. It was found that, in accordance with common-sense expectations, the hydrogen bonded network of water molecules is breaking up continuously as the concentration of the electrolyte increases. The hydration shell of the cations becomes more and more distorted as concentration grows whereas the hydration structure of the anions appears to be nearly invariant. Ion-pairing was not detected even at the highest salt concentration.     

室温熔盐二(三氟甲基磺酸酰)亚胺锂-尿素体系的分子动力学模拟

应用分子动力学模拟了25 ℃和50 ℃时新型室温熔盐二(三氟甲基磺酸酰)亚胺锂[LiN(SO2CF3)2, LiTFSI]与尿素(摩尔比为1:3.6)体系的结构与动力学性质. 在两个温度下体系的微观结构基本相同, Li+的配位数约为5, 且都是与溶剂和阴离子中的氧原子发生配位. 对TFSI-的研究表明, 每个TFSI-只提供四个氧中的一个与Li+配位; 而且在Li+的配位层中, TFSI-具有顺、反和gauche 等不同的构象, 并且不同构象出现的几率会随着温度的改变而改变.     

Molecular simulation study of temperature effect on ionic hydration in carbon nanotubes

Molecular dynamics simulations have been performed to investigate the hydration of Li(+), Na(+), K(+), F(-), and Cl(-) inside the carbon nanotubes at temperatures ranging from 298 to 683 K. The structural characteristics of the coordination shells of ions are studied, including the ion-oxygen radial distribution functions, the coordination numbers, and the orientation distributions of the water molecules. Simulation results show that the first coordination shells of the five ions still exist in the nanoscale confinement. Nevertheless, the first coordination shell structures of cations change more significantly than those of anions because of the preferential orientation of the water molecules induced by the carbon nanotube. The first coordination shells of cations are considerably less ordered in the nanotube than in the bulk solution, whereas the change of the first coordination shell structures of the anions is minor. Furthermore, the confinement induces the anomalous behavior of the coordination shells of the ions with temperature. The first coordination shell of K(+) are found to be more ordered as the temperature increases only in the carbon nanotube with the effective diameter of 1.0 nm, implying the enhancement of the ionic hydration with temperature. This is contrary to that in the bulk solution. The coordination shells of the other four ions do not have such behavior in the carbon nanotube with the effective diameter ranging from 0.73 to 1.00 nm. The easier distortion of the coordination shell of K(+) and the match of the shell size and the nanotube size may play roles in this phenomenon. The exchange of water molecules in the first coordination shells of the ions with the solution and the ion diffusion along the axial direction of the nanotube are also investigated. The mobility of the ions and the stability of the coordination shells are greatly affected by the temperature in the nanotube as in the bulk solutions. These results help to understand the biological and chemical processes at the high temperature.     

Molecular dynamics study of the coordination sphere of trivalent lanthanum in a highly concentrated LiCl aqueous solution: a combined classical and ab initio approach

The first coordination sphere of trivalent lanthanum in a highly concentrated (14 M) lithium chloride solution is studied with a combination of classical molecular dynamics and density functional theory based first principle molecular dynamics. This method enables us to obtain a solvation shell of La3+ containing 2 chloride ions and 6 water molecules. After refinement using first principle molecular dynamics, the resulting cation-water and cation-anion distances are in very good agreement with experiment. The 2Cl- and the 6 water molecules arrange in a square antiprism around La3+. Exchange of water molecules was also observed in the first-principle simulation, with an intermediate structure comprising 7 water molecules stable for 2.5 ps. Finally, evaluation of dipole moments using maximally localized Wannier functions shows a substantial polarization of the choride anions and the water molecules in the first solvation shell of trivalent lanthanum.     

Structure and dynamics of dodecaborate clusters in water

We have studied using molecular dynamics simulations the interaction of the dodecaborate anion, B(12)H(12)(2-), and its amino, trimethyl, and triethyl derivatives with water molecules. We found peculiar organization of the water molecules in the first solvation shell with the formation of a dihydrogen bond between the hydrogen atoms of the anions and the hydrogen atoms of the water molecules. The simulations also show that the organization of the hydration shell is strongly influenced by the substituents in the anions. These differences are likely to play an important role in understanding the interaction of the anions with biological systems like membranes and proteins in aqueous environments.     

Chloride complexation by uranyl in a room temperature ionic liquid. A computational study

The stepwise addition of 1 to 4 Cl(-) anions to the uranyl cation has been studied via potential of mean force (PMF) calculations in the [BMI][Tf 2N] ionic liquid based on the 1-butyl-3-methylimidazolium cation (BMI(+)) and the bis(trifluoromethylsulfonyl)imide anion (Tf2N(-)). According to these calculations, the four Cl(-) complexation reactions are favored and UO2Cl4(2-) is the most stable chloride complex in [BMI][Tf2N]. The solvation of the different chloro-complexes is found to evolve from purely anionic (ca. 5 Tf2N(-) ions around UO2(2+)) to purely cationic (ca. 8.5 BMI (+) cations around UO2Cl4(2-)), with onion-type alternation of solvent shells. We next compare the solvation of the UO2Cl4(2-) complex to its reduced analogue UO2Cl4(3-) in the [BMI][Tf2N] and [MeBu3N][Tf2N] liquids that possess the same anion, but differ by their cation (imidazolium BMI(+) versus ammonium MeBu3N(+)). The overall solvation structure of both complexes is found to be similar in both liquids with a first solvation shell formed exclusively of solvent cations (about 9 BMI(+) cations or 7 MeBu3N(+) cations). However, a given complex is better solvated by the [BMI][Tf2N] liquid, due to hydrogen bonding interactions between Cl(-) ligands and imidazolium-ring C-H protons. According to free energy calculations, the gain in solvation energy upon reduction of UO2Cl4(2-) to UO2Cl4(3-) is found to be larger in [BMI][Tf2N] than in [MeBu3N][Tf2N], which is fully consistent with recent experimental results (Inorg. Chem. 2006, 45, 10419).