Density functional theory-based prediction of the formation constants of complexes of ammonia in aqueous solution: indications of the role of relativistic effects in the solution chemistry of gold(I)

Abstract

A prediction of the formation constants (log K1) for complexes of metal ions with a single NH3 ligand in aqueous solution, using quantum mechanical calculations, is reported. DeltaG values at 298 K in the gas phase for eq 1 (DeltaG(DFT)) were calculated for 34 metal ions using density functional theory (DFT), with the expectation that these would correlate with the free energy of complex formation in aqueous solution (DeltaG(aq)). [M(H2O)6]n+(g) + NH(3)(g) = [M(H2O)5NH3]n+(g) + H2O(g) (eq 1). The DeltaG(aq) values include the effects of complex changes in solvation on complex formation, which are not included in eq 1. It was anticipated that such changes in solvation would be constant or vary systematically with changes in the log K(1) value for different metal ions; therefore, simple correlations between DeltaG(DFT) and DeltaG(aq) were sought. The bulk of the log K1(NH3) values used to calculate DeltaG(aq) were not experimental, but estimated previously (Hancock 1978, 1980) from a variety of empirical correlations. Separate linear correlations between DeltaG(DFT) and DeltaG(aq) for metal ions of different charges (M2+, M3+, and M4+) were found. In plots of DeltaG(DFT) versus DeltaG(aq), the slopes ranged from 2.201 for M2+ ions down to 1.076 for M4+ ions, with intercepts increasing from M2+ to M4+ ions. Two separate correlations occurred for the M3+ ions, which appeared to correspond to small metal ions with a coordination number (CN) of 6 and to large metal ions with a higher CN in the vicinity of 7-9. The good correlation coefficients (R) in the range of 0.97-0.99 for all these separate correlations suggest that the approach used here may be the basis for future predictions of aqueous phase chemistry that would otherwise be experimentally inaccessible. Thus, the log K1(NH3) value for the transuranic Lr3+, which has a half-life of 3.6 h in its most stable isotope, is predicted to be 1.46. These calculations should also lead to a greater insight into the factors governing complex formation in aqueous solution. All of the above DFT calculations involved corrections for scalar relativistic effects (RE). Au has been described (Koltsoyannis 1997) as a "relativistic element". The chief effect of RE for group 11 ions is to favor linear coordination geometry and greatly increase covalence in the M-L bond. The correlation for M+ ions (H+, Cu+, Ag+, Au+) involved the preferred linear coordination of the [M(H2O)2]+ complexes, so that the DFT calculations of DeltaG for the gas-phase reaction in eq 2 were carried out for M = H+, Cu+, Ag+, and Au+. [M(H2O)2]+(g) + NH3(g) = [M(H2O)NH3]+(g) + H2O(g) (eq 2). Additional DFT calculations for eq 2 were carried out omitting corrections for RE. These indicated, in the absence of RE, virtually no change in the log K1(NH3) value for H+, a small decrease for Cu+, and a larger decrease for Ag+. There would, however, be a very large decrease in the log K1(NH3) value for Au(I) from 9.8 (RE included) to 1.6 (RE omitted). These results suggest that much of "soft" acid behavior in aqueous solution in the hard and soft acid-base classification of Pearson may be the result of RE in the elements close to Au in the periodic table.