Theoretical discussion of inorganic reaction mechanisms includes the role of the solvent, requiring estimates of thermodynamics of solvation in reactant, product or transition states. These are often made without reference to a ''benchmark'' calculation for a relatively simple, well-characterized system, which could provide a stringent test of their reliability. A convenient test system for such a calculation is Pt(Cl) 2 (H 2 O) 2 , the cis and trans forms of which in solution exist in equilibrium and thus the relevant experimental thermodynamics is known very precisely [Inorg. Chem. 28 (1989) 3366, and earlier papers]. This paper reports a quantum chemical determination of the relative thermodynamic stabilities of cis- and trans-[Pt(Cl) 2 (H 2 O) 2 ] in aqueous solution by a range of methods and basis sets, concentrating largely on the role of solvation. The methods employed include density functional (DFT) and coupled cluster theories, the conductor-like screening model (COSMO) of solvation and the effective fragment potential (EFP) method. The basis sets range from double zeta to 6-31+G(2df,p) on the ligands and a multiply polarized basis on Pt, in combination with a relativistic effective core potential (ecp). In the gas phase at 0K, the trans - cis energy difference computed by coupled cluster methods is 0.1kcalmol - 1 . The analogous value at 298K in aqueous solution, obtained by DFT (B3LYP) and COSMO, is -0.8kcalmol - 1 , which is in good agreement with the experimental value of -0.1kcalmol - 1 . Modelling solvation by a combination of discrete (quantum chemical) and continuum (COSMO) methods, whereby four H 2 O molecules are explicitly accounted for, yielded a value of -1.5 +/- 1.0kcalmol - 1 . These calculations demonstrate that DFT and polarizable continuum methods (PCMs) with appropriate basis sets are capable of yielding relative thermodynamic stabilities in solution of near chemical accuracy.