Additional Computational Details PCM and NBO calculations Starting from the gas-phase BLYP/SDD/6-311+G** equilibrium structures, [UO 2 L 5 ] 2+ (L = H 2 O, NH 3) were re-optimised using the polarisable continuum model in its integral equation formalism (IEF-PCM), 1 as implemented in Gaussian 09. 2 For L = H 2 O, standard solvent parameters of water were used, for ammonia we employed the static dielectric constant of liquid ammonia at 20ºC ( = 16.6), 3 the same dynamic dielectric constant as water, 4 and a solvent radius of 2.5 Å. 5 Natural bond orbital analyses 6 employed the Gaussian NBO version 3.1 as implemented in Gaussian 09. CPMD simulations The same methods and basis sets as in our previous studies of uranyl complexes 7 were employed. Car-Parrinello molecular dynamics (CPMD) 8 simulations were performed using the BLYP functional 9 and norm-conserving pseudopotentials that had been generated according to the Troullier and Martins procedure 10 and transformed into the Kleinman-Bylander form. 11 For uranium, the semicore (or small-core) pseudopotential was employed that had been generated and validated in reference 12. Periodic boundary conditions were imposed using cubic supercells with a lattice contant of 13.22 Å. Kohn-Sham orbitals were expanded in plane waves at the -point up to a kinetic energy cutoff of 80 Ry. Simulations were performed in the NVT ensemble using a single Nosé-Hoover thermostat set to 300 K (frequency 1800 cm 1), a fictitious electronic mass of 600 a.u., and a time step of 0.121 fs. The boxes contained uranyl and and a total of 42 ammonia molecules, affording a density of 0.71. The system has two positive charges, neutralised by a uniform background charge. In order to maintain the time step, hydrogen was substituted with deuterium. Long-range electrostatic interactions were treated with the Ewald method. No electro-static decoupling between replicated cells was included. Constrained CPMD simulations were performed along a predefined reaction coordinate (a single U-N bond distance r), starting from the respective mean values for the five-coordinate minumum in gas and solution (as obtained from the unconstrained simulations. Changes in the Helmholtz free energes were evaluated by pointwise thermodynamic integration 13 of the mean constraint force f along this coordinate via A ab = a b f()d (1). Electronic Supplementary Material (ESI) for Dalton Transactions. This journal is