Developing Re(V)-based therapeutic agents, where the ReO 3+ core is coordinated to a multidentate ligand, is of interest in the radiopharmaceutical sciences because of the desirable nuclear properties of 186 Re/ 188 Re. A reliable yet cost-effective computational method for evaluating the strength of each coordination bond while preserving the integrity of the metal-ligand complex would provide quantitative input for ligand design. A relaxed potential energy surface (PES) scan approach is assessed for trans-ReO(SH) 2 (NH 2)(NH 3), [ReO(SH) 3 (NH 2)] 1-, and [ReO(SH) 3 (N(H)CHO)] 1model complexes to calculate bond dissociation energies (BDEs) for Re-NH 3 , Re-NH 2 , Re-N(H)CHO and Re-SH bonds, common components of Re(V) coordination environments. The PES scans were performed using various combinations of DFT/coupled-cluster methods and basis sets, and the effect of bulk solvent was examined by using the integral equation formalism of the polarizable continuum model (IEF-PCM). BDEs obtained from the PES scans are compared to those obtained for infinite separation. In the gas phase, the BDE curves reach about 90% of the total BDE at 2 Å and plateau by 3.0 to 3.5 Å beyond the equilibrium bond length; in the presence of implicit solvent, the BDE water curves plateau at a shorter distance and more than 90% of the total BDE is recovered at 2 Å. The gas-phase PES scans follow the desired reaction coordinate for NH 3 , [N(H)CHO] 1and SH 1-, but not for the poor leaving group NH 2 1-. The desired heterolytic cleavage of the Re-NH 2 bond is achieved when the PES scans are performed in the presence of solvent. Elongating the Re-S/N bonds in a rigid, multidentate trans-ReO-N 2 S 2 complex yields BDE trends similar to those found for the model complexes (Re-NH 2 > Re-N(H)CHO > Re-SH > Re-NH 3).