The structure of [Cu(aq)]2+ has been investigated by using full multiple-scattering theoretical (MXAN) analysis of the copper K-edge X-ray absorption (XAS) spectrum and density functional theory (DFT) to test both ideal Td and square-planar four-coordinate, five-coordinate square-pyramidal, and six-coordinate octahedral [Cu(aq)]2+ models. The best fit was an elongated five-coordinate square pyramid with four Cu-O(eq) bonds (2 x 1.98 +/- 0.03 A and 2 x 1.95 +/- 0.03 A) and a long Cu-O(ax) bond (2.35 +/- 0.05 A). The four equatorial ligands were D2d-distorted from the mean equatorial plane by +/-(17 +/- 4) degrees, so that the overall symmetry of [Cu(H2O)5]2+ is C2v. The four-coordinate MXAN fit was nearly as good, but the water ligands (4 x 1.96 +/- 0.02 A) migrated +/-(13 +/- 4) degrees from the mean equatorial plane, making the [Cu(H2O)4]2+ model again D2d-distorted. Spectroscopically calibrated DFT calculations were carried out on the C2v elongate square-pyramidal and D2d-distorted four-coordinate MXAN copper models, providing comparative electronic structures of the experimentally observed geometries. These calculations showed 0.85e spin on Cu(II) and 0.03e electron spin on each of the four equatorial water oxygens. All covalent bonding was restricted to the equatorial plane. In the square-pyramidal model, the electrostatic Cu-O(ax) bond was worth only 96.8 kJ mol(-1), compared to 304.6 kJ mol(-1) for each Cu-O(eq) bond. Both MXAN and DFT showed the potential well of the axial bond to be broad and flat, allowing large low-energy excursions. The irregular geometry and D2d-distorted equatorial ligand set sustained by unconstrained [Cu(H2O)5]2+ warrants caution in drawing conclusions regarding structural preferences from small molecule crystal structures and raises questions about the site-structural basis of the rack-induced bonding hypothesis of blue copper proteins. Further, previously neglected protein folding thermodynamic consequences of the rack-bonding hypothesis indicate an experimental disconfirmation.

In this paper a new software procedure is presented, named MXAN, able to fit the XANES part (from the edge to about 200 eV) of experimental X-ray absorption data. The method is based on the fitting between the experimental spectrum and several theoretical calculations generated by changing the relevant geometrical parameters of the site around the absorbing atom. The X-ray photoabsorption cross section is calculated using the full multiplescattering scheme; different choices of the exchange correlation part of the potential can be utilised. To show the potentialities of the method the analysis of the nickel K-edge of the nickel ion in aqueous solution is presented. The procedure is able to recover the correct information on the symmetry and atomic distances from the experimental Ni K-edge XANES spectrum. The recovered structure is found to be independent of the starting conditions, showing the theoretical uniqueness of the structural solution.