Dispersion of the Raman depolarization ratio of HDO in water and heavy water from 295 to 368 K, and from concentrated NaClO4D2OH2O

The dispersion of the Raman depolarization ratio rho(L) was measured for HDO in H(2)O and in D(2)O. rho(L) for the decoupled OD stretch displays a maximum at 2575 +/- 15 cm(-1) at 296 K and a minimum at 2675 +/- 15 cm(-1), in agreement with the isosbestic point 2570 +/- 10 cm(-1), and the enthalpy dispersion maximum, 2650-2675 cm(-1), respectively. However, three extrema were uncovered in rho(L) for the OH stretch of HDO in D(2)O, and their positions agree with the frequencies of a minimum and a maximum in the enthalpy dispersion and with the isosbestic frequency. The frequency of the rho(L) maximum (OH stretch) lies just above the frequency corresponding to the joint angle-frequency probability maximum. [Lawrence and Skinner, J. Chem. Phys. 118, 264 (2003)]. The low- and high-frequency minima in rho(L) (OH stretch), correspond, respectively, to very strong H-bonds, and extremely weak, long, bent H bonds. The frequencies of the maxima and minima in rho(L) for the decoupled OH and OD stretches are independent of temperature within experimental error between 295 and 368 K. rho(L) was also measured for the OD stretch from saturated NaClO(4) in D(2)OH(2)O; it displays a maximum at 2560 +/- 20 cm(-1) and a sharp minimum at 2650 +/- 5 cm(-1). The shape of the dispersion of (betaalpha)(2) approximately rho(L) for HDO in D(2)O was calculated with the aid of the molecular dynamics results of Lawrence and Skinner. beta(2) is the anisotropic polarizability and alpha is the isotropic polarizability. A maximum resulted in the calculated dispersion at 3400 +/- 10 cm(-1), in excellent agreement with the measured maximum of 3395 +/- 15 cm(-1). The H-bond angles decrease far below 180 degrees as the OH-stretching frequency increases to 3700 cm(-1) and above. Such small H-bond angles, and very large O-O distances, are tantamount to broken H-bonds and are thought to produce the minimum in rho(L) near 2650 cm(-1).     

Quantum vibrational analysis of hydrated ions using an ab initio potential

We present full-dimensional potential energy surfaces (PESs) for hydrated chloride based on the sum of ab initio (H(2)O)Cl(-), (H(2)O)(2), and (H(2)O)(3) potentials. The PESs are shown to predict minima and corresponding harmonic frequencies accurately on the basis of comparisons with previous and new ab initio calculations for (H(2)O)(2)Cl(-), (H(2)O)(3)Cl(-), and (H(2)O)(4)Cl(-). An estimate of the effect of the 3-body water interaction is made using a simple 3-body water potential that was recently fit to tens of thousands of ab initio 3-body energies. Anharmonic, coupled vibrational calculations are presented for these clusters, using the "local monomer model" for the high frequency intramolecular modes. This model is tested against previous "exact" calculations for (H(2)O)Cl(-). Radial distribution functions at 0 K obtained from quantum zero-point wave functions are also presented for the (H(2)O)(2)Cl(-) and (H(2)O)(3)Cl(-) clusters.     

First principles study on the solvation and structure of C2O4(2-)(H2O)n, n = 6-12

The structures and energies of hydrated oxalate clusters, C2O4(2-)(H2O)n, n = 6-12, are obtained by density functional theory (DFT) calculations and compared to SO4(2-)(H2O)n. Although the evolution of the cluster structure with size is similar to that of SO4(2-)(H2O)n, there are a number of important and distinctive futures in C2O4(2-)(H2O)n, including the separation of the two charges due to the C-C bond in C2O4(2-), the lower symmetry around C2O4(2-), and the torsion along the C-C bond, that affect both the structure and the solvation energy. The solvation dynamics for the isomers of C2O4(2-)(H2O)12 are also examined by DFT based ab initio molecular dynamics.     

Infrared spectra of M(OH)(1,2,4) (M = Pb, Sn) in solid argon

Infrared absorptions for the matrix-isolated lead and tin hydroxides M(OH), M(OH)2 and M(OH)4 (M = Pb, Sn) were observed in laser-ablated metal atom reactions with H2O2 during condensation in excess argon. The major M(OH)2 product was also observed with H2 and O2 mixtures, which allowed the substitution of 18O2. The band assignments were confirmed by appropriate D2O2, D2, 16O18O, and 18O2 isotopic shifts. MP2 and B3LYP calculations were performed to obtain molecular structures and to reproduce the infrared spectra. The minimum energy structure found for M(OH)2 has C(s) symmetry and a weak intramolecular hydrogen bond. In experiments with Sn, HD, and O2, the internal D bond is favored over the H bond for Sn(OH)(OD). The Pb(OH)4 and Sn(OH)4 molecules are calculated to have S4 symmetry and substantial covalent character.     

Overtone-induced dissociation and isomerization dynamics of the hydroxymethyl radical (CH2OH and CD2OH). I. A theoretical study

The dissociation of the hydroxymethyl radical, CH(2)OH, and its isotopolog, CD(2)OH, following the excitation of high OH stretch overtones is studied by quasi-classical molecular dynamics calculations using a global potential energy surface (PES) fitted to ab initio calculations. The PES includes CH(2)OH and CH(3)O minima, dissociation products, and all relevant barriers. Its analysis shows that the transition states for OH bond fission and isomerization are both very close in energy to the excited vibrational levels reached in recent experiments and involve significant geometry changes relative to the CH(2)OH equilibrium structure. The energies of key stationary points are refined using high-level electronic structure calculations. Vibrational energies and wavefunctions are computed by coupled anharmonic vibrational calculations. They show that high OH-stretch overtones are mixed with other modes. Consequently, trajectory calculations carried out at energies about ~3000 cm(-1) above the barriers reveal that despite initial excitation of the OH stretch, the direct OH bond fission is relatively slow (10 ps) and a considerable fraction of the radicals undergoes isomerization to the methoxy radical. The computed dissociation energies are: D(0)(CH(2)OH → CH(2)O + H) = 10,188 cm(-1), D(0)(CD(2)OH → CD(2)O + H) = 10,167 cm(-1), D(0)(CD(2)OH → CHDO + D) = 10,787 cm(-1). All are in excellent agreement with the experimental results. For CH(2)OH, the barriers for the direct OH bond fission and isomerization are: 14,205 and 13,839 cm(-1), respectively.     

Vibrational spectra of Na, K, Mn2+, Ni2+ and Zn2+ salts of 1,2,4,5-benzenetetracarboxylic (pyromellitic) acid--a short hydrogen bond evidence

Five salts of 1,2,4,5-benzenetetracarboxylic acid (pyromellitic acid), [C6H2(COO)4H4], have been synthesized and investigated by infrared and Raman spectroscopy and by single crystal X-ray diffraction methods: sodium salt [Na2(H2O)2][C6H2(COO)4H2], potassium salt [K(H2O)3][C6H2(COO)4H3] and transition metal salts [M(H2O)6][C6H2(COO)4H2], which M = Mn, Ni and Zn. Crystal structures of all five compounds show short intramolecular asymmetric hydrogen bonds (SHB) between adjacent carboxyl groups with O...O distance average of 2.40 A. The Raman and infrared spectra reported indicate the presence of short hydrogen bonds in all salts, in agreement with the X-ray data. The O-H stretching mode [nu(OH)] had been observed at about 2500 cm(-1). Deuterated analogues were synthesized and their Raman spectra show that nu(OH)/nu(OD) ratio average is about unit. The symmetric [nu(sym)(O..H..O)] and asymmetric [nu(asym)(O..H..O)] stretching modes have been attributed about 300 and 870 cm(-1), respectively, in all salts, and for deuterated analogues, the ratio nu(OH)/nu(OD) to nu(sym)(O..H..O, O..D..O) is close to unit like it occurs in nu(OH). The vibrational modes, mainly SHB modes, are tentatively assigned by molecular orbital ab initio calculations of pyromellitic acid and anions [C6H2(COO)4H3]- and [C6H2(COO)4H2]2-. Geometry optimizations showed a good agreement with experimental data. Frequency calculation confirms the assignment of specific vibrational modes. Ab initio calculations show that nu(C=O) and nu(sym)(COO) are strongly coupled with in plane OH bending [delta(OH)]. In Raman spectra of deuterated analogues is observed a frequency shift of these bands.     

Contribution of the asymmetric stretch, nu3B1, to the fundamental Raman spectrum of water

The OD-stretching overtone from liquid D2O, 2nu, and the fundamental OD stretch from dilute HDO, both display high-frequency depolarization ratio minima, but the fundamental OD stretch from neat D2O displays a maximum, at the equivalent position. The rhoL minima arises from the decreased depolarization ratio produced by the absence of B1 modes. The fundamentals of HDO are of A species, and the 2nu overtone of D2O only involves A1 species, e.g., 2nu3B1 has A1 species via B1 x B1 = A1. A and A1 modes display small rhoL values which produce minima in rhoL near 2665 cm(-1) for HDO, and near 5250 cm(-1) for the D2O overtone. These minima give way to a depolarization ratio maximum when the depolarized, rhoL = 34, nu3B1 fundamental, makes its appearance in D2O at 2650 cm(-1). Fundamental and overtone depolarization ratios were used to determine the nu3B1 contribution to the depolarization ratio of the fundamental OD stretch; a value of approximately 28% resulted at 2655 cm(-1). Liquid H2O displays completely analogous features; a value of approximately 20% resulted for it at 3660 cm(-1). Nonhydrogen-bonded nu3B1, and more strongly hydrogen-bonded nu3B1, modes are also indicated for D2O and H2O. A rigorous test of the current results can be accomplished by measuring the depolarization ratio of the extraordinarily weak second Raman overtone, 3nu, recently detected for D2O.     

A dodecanuclear Zn cluster sandwiched by polyoxometalate ligands

A dodecazinc silicotungstate K(20)Na(2)[Zn(6)(OH)(7)(H(2)O)(Si(2)W(18)O(66))](2)·34H(2)O (1) has been synthesized and characterized by X-ray crystallography, elemental analysis, infrared, UV-vis spectroscopy, cyclic voltammetry, acid-base titration, and DFT calculations. The twelve zinc atoms between the two [Si(2)W(18)O(66)](16-) frameworks make this complex more stable hydrolytically than the heteropolytungstate ligands, [Si(2)W(18)O(66)](16-), themselves. The structurally unique central Zn(12) core is formed by the fusion of two [Zn(6)(OH)(7)(H(2)O)](5+) units through two edge-sharing Zn6 atoms. DFT B3LYP calculations give HOMO-LUMO and (HOMO - 1)-LUMO energy gaps of ～3.65 and 3.91 eV, respectively, as compared to the band gap in ZnO of 3.35 eV.     

Ab initio potential and dipole moment surfaces for water. II. Local-monomer calculations of the infrared spectra of water clusters

We employ recent flexible ab initio potential energy and dipole surfaces [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] to the calculation of IR spectra of the intramolecular modes of water clusters. We use a quantum approach that begins with a partitioned normal-mode analysis of perturbed monomers, and then obtains solutions of the corresponding Schro?dinger equations for the fully coupled intramolecular modes of each perturbed monomer. For water clusters, these modes are the two stretches and the bend. This approach is tested against benchmark calculations for the water dimer and trimer and then applied to the water clusters (H(2)O)(n) for n = 6-10 and n = 20. Comparisons of the spectra are made with previous ab initio harmonic and empirical potential calculations and available experiments.     

Communication: quasiclassical trajectory calculations of correlated product-state distributions for the dissociation of (H2O)2 and (D2O)2

Stimulated by recent experiments [B. E. Rocher-Casterline, L. C. Ch'ng, A. K. Mollner, and H. Reisler, J. Chem. Phys. 134, 211101 (2011)], we report quasiclassical trajectory calculations of the dissociation dynamics of the water dimer, (H(2)O)(2) (and also (D(2)O)(2)) using a full-dimensional ab initio potential energy surface. The dissociation is initiated by exciting the H-bonded OH(OD)-stretch, as done experimentally for (H(2)O)(2). Normal mode analysis of the fragment pairs is done and the correlated vibrational populations are obtained by (a) standard histogram binning (HB), (b) harmonic normal-mode energy-based Gaussian binning (GB), and (c) a modified version of (b) using accurate vibrational energies obtained in the Cartesian space. We show that HB allows opening quantum mechanically closed states, whereas GB, especially via (c), gives physically correct results. Dissociation of both (H(2)O)(2) and (D(2)O)(2) mainly produces either fragment in the bending excited (010) state. The H(2)O(J) and D(2)O(J) rotational distributions are similar, peaking at J = 3-5. The computations do not show significant difference between the ro-vibrational distributions of the donor and acceptor fragments. Diffusion Monte Carlo computations are performed for (D(2)O)(2) providing an accurate zero-point energy of 7247 cm(-1), and thus, a benchmark D(0) of 1244 ± 5 cm(-1).     

Infrared spectroscopic observation of the group 13 metal hydroxides, M(OH)1,2,3 (M =Al, Ga, In, and Tl) and HAl(OH)2

Reactions of laser-ablated Al, Ga, In, and Tl atoms with H2O2 and with H2 + O2 mixtures diluted in argon give new absorptions in the O-H and M-O stretching and O-H bending regions, which are assigned to the metal mono-, di-, and trihydroxide molecules. Isotopic substitutions (D2O2, 18O2, 16,18O2, HD, and D2) confirm the assignments, and DFT calculations reproduce the experimental results. Infrared spectra for the Al(OH)(OD) molecule verify the calculated C2v structure. The trihydroxide molecules increase on annealing from the spontaneous reaction with a second H2O2 molecule. Aluminum atom reactions with the H2 + O2 mixtures favor the HAl(OH)2 product, suggesting that AlH3 generated by UV irradiation combines with O2 to form HAl(OH)2.     

Raman H-bond pair volume for water

The dispersion of the H-bond pair volume Delta V over the decoupled OD and coupled OH-stretching contours from HDO in H(2)O was determined from Raman intensities at pressures to 9700 bar at 301 K. The dispersion of Delta V was determined from -RT[partial differential ln(I(i)/I(REF))/ partial differential P](T) versus omega (in cm(-1)), where i refers to omega's over the stretching contours and I(REF) refers to the reference intensity at the isosbestic frequency. The maximum H-bond pair volume (defined for breakage) is 1.4+/-0.1 cm(3)/mol H-bond, which corresponds to the volume difference between a large dispersion maximum at 2,675 cm(-1) (near the OD stretch omega of HDO in dense supercritical water) and a large, broad minimum centered near 2,375 cm(-1) (just below the OD stretch omega of HDO in lda ice). The average DeltaV is 0.71+/-0.10 cm(3)/mol H-bond. Other minima near 2,625 cm(-1) (OD) and 3550 cm(-1) (OH) refers to bent H-bonds whose angles are approximately 150 deg. Isothermal pressurization of water lowers the molal volume by decreasing the concentration of long, weak H-bonds, and increasing the concentrations of bent H-bonds and short, strong, linear H-bonds. Such bending, shortening, and strengthening produces freezing to ice VI near 10 kbar at 301 K. The isobaric temperature derivative of the maximum H-bond volume is (partial differential Delta V/partial differential T)(P)< or =(2-5) x 10(-3) cm(3)/deg mol H-bond. The OH enthalpy dispersion curve for saturated NaBF(4) in water, yields a large maximum at 3,530-3,540 cm(-1) indicating that BF(4) (-) interacts preferentially with the dangling or "free" OH groups of water thus producing weak, strongly bent H-bonds having angles similar to those of the 3,550 cm(-1) high-pressure H-bond bending feature.     

H·(H2O)n clusters: microsolvation of the hydrogen atom via molecular ab initio gradient embedded genetic algorithm (GEGA)

A new version of the ab initio gradient embedded genetic algorithm (GEGA) program for finding the global minima on the potential energy surface (PES) of mixed clusters formed by molecules and atoms is reported. The performance of the algorithm is demonstrated on the neutral H·(H(2)O)(n) (n = 1-4) clusters, that is, a radical H atom solvated in 1-4 water molecules. These clusters are of a fundamental interest. The solvated hydrogen atom forms during photochemical events in water, or during scavenging of solvated electrons by acids, and transiently exists in biological systems and possibly in inclusion complexes in the deep ocean and in the ice shield of earth. The processes associated with its existence are intriguingly complex, however, and have been the subject of decades-long debates. Using GEGA, we explicate the apparently extreme structural diversity in the H·(H(2)O)(n) (n = 1-4) clusters. All considered clusters have four basic structural types: type I, where the H radical is weakly coordinated to the oxygen atom of one of the water molecules; type II, where H is weakly coordinated to a H atom of one of the water molecules; type III, consisting of H(2), the OH radical, and n - 1 H(2)O molecules; and type IV, consisting of H(3)O and n - 1 H(2)O. There are myriads of isomers of all four types. The lowest energy species of types I and II are the isoenergetic global minima. H·(H(2)O)(n) clusters appear to be a challenging case for GEGA because they have many shallow minima close in energy some of which are significantly less stable than the global minimum. Additionally, the global minima themselves have high structural degeneracy, they are only weakly bound, and they are prone to dissociation. GEGA performed exceptionally well in finding both the global and the low-energy local minima that were subsequently confirmed at higher levels of theory.     

Stoichiometry and structure of uranylVI hydroxo dimer and trimer complexes in aqueous solution

We studied the structure and stoichiometry of aqueous uranylVI hydroxo dimers and trimers by spectroscopic (EXAFS, FTIR, UV-vis) and quantum chemical (DFT) methods. FTIR and UV-vis spectroscopy were used for the speciation of uranyl complexes in aqueous solution. DFT calculations show that (UO2)2(OH)22+ has two bridging hydroxo groups with a U-U distance of 3.875 A. This result is in good agreement with EXAFS, where a U-U distance of 3.88 A was found. For the hydroxo trimer complex, DFT calculations show that the species (UO2)3(O)(OH)3+ with oxo bridging in the center is energetically favored in comparison to its stoichiometric equivalent (UO2)3(OH)5+. This is again in line with the EXAFS result, where a shorter U-U distance of 3.81-3.82 A and evidence for oxo bridging in the center were found. Several stable intermediates which lie several tens of kJ/mol above that of (UO2)3(O)(OH)3+ were identified, and their structures, energies, and intramolecular proton-transfer reaction are discussed.     

Computational study on the existence of organic peroxy radical-water complexes (RO2.H2O)

The existence of a series of organic peroxy radical-water complexes [CH3O2.H2O (methyl peroxy); CH3CH2O2.H2O (ethyl peroxy); CH3C(O)O2.H2O (acetyl peroxy); CH3C(O)CH2O2.H2O (acetonyl peroxy); CH2(OH)O2.H2O (hydroxyl methyl peroxy); CH2(OH)CH2O2.H2O (2-hydroxy ethyl peroxy); CH2(F)O2.H2O (fluoro methyl peroxy); CH2(F)CH2O2.H2O (2-fluoro ethyl peroxy)] is evaluated using high level ab initio calculations. A wide range of binding energies is predicted for these complexes, in which the difference in binding energies can be explained by examination of the composition of the R group attached to the peroxy moiety. The general trend in binding energies has been determined to be as follows: fluorine approximately alkyl < carbonyl < alcohol. The weakest bound complex, CH3O2.H2O, is calculated to be bound by 2.3 kcal mol-1, and the strongest, the CH2(OH)O2.H2O complex, is bound by 5.1 kcal mol-1. The binding energy of the peroxy radical-water complexes which contain carbonyl and alcohol groups indicates that these complexes may perturb the kinetics and product branching ratios of reactions involving these complexes.     

Communication: rigorous calculation of dissociation energies (D0) of the water trimer, (H2O)3 and (D2O)3

Using a recent, full-dimensional, ab initio potential energy surface [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] together with rigorous diffusion Monte Carlo calculations of the zero-point energy of the water trimer, we report dissociation energies, D(0), to form one monomer plus the water dimer and three monomers. The calculations make use of essentially exact zero-point energies for the water trimer, dimer, and monomer, and benchmark values of the electronic dissociation energies, D(e), of the water trimer [J. A. Anderson, K. Crager, L. Fedoroff, and G. S. Tschumper, J. Chem. Phys. 121, 11023 (2004)]. The D(0) results are 3855 and 2726 cm(-1) for the 3H(2)O and H(2)O + (H(2)O)(2) dissociation channels, respectively, and 4206 and 2947 cm(-1) for 3D(2)O and D(2)O + (D(2)O)(2) dissociation channels, respectively. The results have estimated uncertainties of 20 and 30 cm(-1) for the monomer plus dimer and three monomer of dissociation channels, respectively.     

Rationalization of the lanthanide-ion-driven magnetic properties in a series of 4f-5d cyano-bridged chains

Magnetic properties of new d-f cyanido-bridged 1D assemblies [RE(pzam)(3)(H(2)O)W(CN)(8)]·H(2)O (RE(III) = Gd, 1, Tb, 2, Dy, 3; pzam = pyrazine-2-carboxamide) were studied by temperature- and field-dependent magnetization measurements. No evidence for 3D interchain magnetic ordering is found above 2 K. Multiconfiguration ab initio calculations and subsequent modeling afforded simulation of the weak zero-field splitting effect in 1 and discussion of magnetic anisotropy in the f units of compounds 2 and 3. A semiquantitative corroboration with the experimental magnetic measurements is presented, performing the simulation of magnetic susceptibility vs temperature and magnetization vs field variation. The association into molecular and supramolecular architectures is analyzed by means of energy decomposition subsequent to the DFT calculations on idealized molecular models extracted from the experimental chain structure.     

The first photocrystallographic evidence for light-induced metastable linkage isomers of ruthenium sulfur dioxide complexes

Light-induced metastable linkage isomers of trans-[Ru(NH(3))(4)Cl(SO(2))]Cl and trans-[Ru(NH(3))(4)(H(2)O)(SO(2))](C(6)H(5)SO(3))(2) have been identified for the first time using photocrystallographic methods. In both linkage isomers the SO(2) ligand is side bound, but the Ru-O and Ru-S distances are considerably longer and almost equal in the trans-H(2)O isomer. DFT calculations confirm that both isomers correspond to minima on the ground-state potential energy surface and also predict the existence of a second oxygen-bound isomer for both compounds. The decay of the light-induced species has been studied by both DSC and IR. Activation energies for the thermal back-reaction, as derived from the temperature-dependent disappearance of light-induced IR bands, are 50.0 and 58.4 kJ/mol for the two isomers, which is larger than the corresponding numbers for photoinduced side-bound nitrosyl linkage isomers.     

Dimers of boroglycine and methylamine boronic acid: a computational comparison of the relative importance of dative versus hydrogen bonding

Boronic acids are widely used in materials science, pharmacology, and the synthesis of biologically active compounds. In this Article, geometrical structures and relative energies of dimers of boroglycine, H2N-CH2-B(OH)2, and its constitutional isomer H3C-NH-B(OH)2, were computed using second-order M?ller-Plesset perturbation theory and density functional theory; Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the MP2 calculations, and the Pople 6-311++G(d,p) basis set was employed for a majority of the DFT calculations. Effects of an aqueous environment were incorporated into the results using PCM and COSMO-RS methodology. The lowest-energy conformer of the H2N-CH2-B(OH)2 dimer was a six-membered ring structure (chair conformation; Ci symmetry) with two intermolecular B:N dative-bonds; it was 14.0 kcal/mol lower in energy at the MP2/aug-cc-pVDZ computational level than a conformer with the classic eight-centered ring structure (Ci symmetry) in which the boroglycine monomers are linked by a pair of H-O...H bonds. Compared to the results of MP2 calculations with correlation-consistent basis sets, DFT calculations using the PBE1PBE and TPSS functionals with the 6-311++G(d,p) basis set were significantly better at predicting relative conformational energies of the H2N-CH2-B(OH)2 and H3C-NH-B(OH)2 dimers than corresponding calculations using the BLYP, B3LYP, OLYP, and O3LYP functionals, particularly with respect to dative-bonded structures.     

DFT study on the mechanism for the substitution of F(-) into Al(III) complexes in aqueous solution

The mechanisms for the substitution of an aqua ligand with F(-) in monomeric Al complexes were studied with density functional theory (DFT). Typical mechanisms are modeled to determine the preferred substitution pathway according to the activation energy barriers. The present computational results are in favor of interchange associative (I(a)) mechanism for the substitution of F(-) into Al(H(2)O)(6)(3+), whereas interchange dissociative (I(d)) mechanism is preferred for the substitution into Al(H(2)O)(5)(OH)(2+), which is in agreement with the previous experimental findings. This implies the mechanistic changeover from I(a) to I(d) induced by the spectator hydroxyl ligand. Like the water-exchange reaction, the substitution rate is accelerated by OH(-) ligand. The difference of the computational and experimental activation enthalpy values is interpreted as the DFT errors in energy and the deviation of transmission coefficient from unity.