Reliable density functional theory (DFT) calculations can be performed in conjuction with spectroscopic measurements to elucidate the structural properties of physiologically important bis(amino acidato)copper(II) compounds in solutions. They can provide insight into the influence of intermolecular interactions on the molecular geometry in the crystal lattice or solution when compared with a DFT gas-phase minimum. Our previous paper [Marković et al. (2014) Eur J Inorg Chem 198] reported the DFT-determined geometries and Raman spectra for different conformers of physiological bis(L-histidinato)copper(II) with 20 explicit water molecules, as calculated using the B3LYP functional. The present study examined the reliability of those B3LYP results by applying the M06 functional instead, as it should better account for noncovalent interactions. The water molecules were positioned more compactly around the complex by M06 than by B3LYP. The accuracies of the two functionals when compared to relevant experimental data showed that M06 was better at reproducing in-plane Cu−N bond lengths but B3LYP gave more accurate axial Cu−O distances. Both functionals reproduced the experimental Raman spectrum at pH 8 to similar levels of accuracy and provided precise information on the Cu(II) coordination mode and conformation in aqueous solution. Additionally, we assessed several DFT and DFT-D functionals (BP86, B3LYP, B3LYP-D, M06, M06 L, wB97XD, mPW2PLYPD) by using them to model the geometries of experimental bis(L-histidinato)copper(II) crystalline conformations as isolated systems, and then benchmarking the results against those from high-level second-order pertubation Møller–Plesset (MP2) calculations. Although this assessment resulted in an equivocal conclusion because the MP2 results for the isolated complex were inconsistent with the corresponding DFT outcomes, it does provide new information on future benchmark options.